费诺与香农编码源程序及运行结果.doc

香农编码源程序 #includeiostream.h #includemath.h #includeiomanip.h #includestdlib.h class T { public: T() {} ~T(); void Create(); void Coutpxj(); void Coutk(); void Coutz(); void Print(); protected: int n; double *p; double *pxj; int *k; double *mz; }; void T::Create() { cout请输入信源符号个数:; cinn; p=new double[n]; cout请分别输入这n个概率:\n; for(int i=0;in;i++) cinp[i]; pxj=new double[n]; k=new int[n]; mz=new double[n]; double sum=0.0; for(i=0;in;i++) sum+=p[i]; if(sum!=1.0) throw 1; else { for(i=0;in;i++) { int k=i; for(int j=i+1;jn;j++) if(p[k]p[j]) k=j; double m=p[i]; p[i]=p[k]; p[k]=m; } } } T::~T() { delete p; delete pxj; delete k; delete mz; } void T::Coutpxj() { pxj[0]=0; for(int i=1;in;i++) { pxj[i]=0; for(int j=0;ji;j++) pxj[i]+=p[j]; } } void T::Coutk() { for(int i=0;in;i++) { double d=(-1)*(log(p[i])/log(2)); if(d-(int)d0) k[i]=(int)d+1; else k[i]=(int)d; } } void T::Print() { coutXisetw(8)P(xi) setw(8)Pa(xj) setw(8)Ki setw(8)码字 endl; for(int i=0;in;i++) { coutXi+1 setw(8)setprecision(2)p[i] setw(8)setprecision(2)pxj[i] setw(8)k[i] ; mz[i]=pxj[i]; for(int j=0;jk[i];j++) { if(2*mz[i]-1=0) { cout1; mz[i]=2*mz[i]-1; } else { cout0; mz[i]=2*mz[i]; } } coutendl; } double K=0.0,H=0.0,Y; for(i=0;in;i++) { K+=(double)p[i]*k[i]; H+=(-1)*p[i]*(log10(p[i])/log10(2.0)); } Y=H/K; cout平均码长:Kendl; cout信源熵:Hendl; cout编码效率:Yendl; } void main() { T t;int e; try { t.Create(); t.Coutpxj(); t.Coutk(); t.Print(); } catch(int e) {if(e==1) cout输入错误，请重新运行;} } 费诺编码源程序 #includeiostream.h #includemath.h #define N 15 int pa[N][N]; void fano(float p[],int a[N][N],int n,int m,int k) { float g=0.0,h=0.0,d,b,c; int i,j