费诺和香农编码源程序及运行结果

香农编码源程序#include<iostream.h>#include<math.h>#include<iomanip.h>#include<stdlib.h>classT{public:T(){}~T();voidCreate();voidCoutpxj();voidCoutk();voidCoutz();voidPrint();protected:intn;double*p;double*pxj;int*k;double*mz;};voidT::Create(){cout<<"请输入信源符号个数:"；cin>>n;p=newdouble[n]；cout<<"请分别输入这"<<n<<"个概率:\n";for(inti=0；i<n；i++)cin>>p[i];pxj=newdouble[n];k=newint[n];mz=newdouble[n];doublesum=0.0;for(i=0;i<n;i++)sum+=p[i];if(sum!=1.0)throw1;else{for(i=0;i<n;i++){intk=i;for(intj=i+1;j<n;j++)if(p[k]<p[j])k=j;doublem=p[i];p[i]=p[k];p[k]=m;}}}T::~T(){deletep;deletepxj;deletek;deletemz;}voidT::Coutpxj(){pxj[0]=0;for(inti=1;i<n;i++){pxj[i]=0;for(intj=0;j<i;j++)pxj[i]+=p[j];}}voidT::Coutk(){for(inti=0;i<n;i++){doubled=(-1)*(log(p[i])/log(2));if(d-(int)d>0)k[i]=(int)d+1;elsek[i]=(int)d;}}voidT::Print(){cout<<"Xi"<<setw(8)<<"P(xi)"<<setw(8)<<"Pa(xj)"<<setw(8)<<"Ki"<<setw(8)<<"码字"<<endl;for(inti=0;i<n;i++){cout<<"X"<<i+1<<setw(8)<<setprecision(2)<<p[i]<<setw(8)<<setprecision(2)<<pxj[i]<<setw(8)<<k[i]<<"";mz[i]=pxj[i];for(intj=0;j<k[i];j++){if(2*mz[i]-1>=0){cout<<1;mz[i]=2*mz[i]-1;}else{cout<<0;mz[i]=2*mz[i];}}cout<<endl;}doubleK=0.0,H=0.0,Y;for(i=0;i<n;i++){K+=(double)p[i]*k[i];H+=(-1)*p[i]*(log10(p[i])/log10(2.0));}Y=H/K;cout<〈"平均码长:"<<K<<endl;cout<<"信源熵："<<H<<endl;cout<<"编码效率:"<<Y<<endl;}voidmain(){Tt;inte;try{t.Create();t.Coutpxj();t.Coutk();t.Print();}catch(inte){if(e==1)cout<<"输入错误，请重新运行";}}费诺编码源程序#include<iostream.h>#include<math.h>#defineN15intpa[N][N];voidfano(floatp[],inta[N][N],intn,intm,intk){floatg=0.0,h=0.0,d,b,c;inti,j,flase;if(n<m){for(i=n;i<=m;i++){g=p[i]+g;}g=g/2;for(i=n;i<=m;i++){h=h+p[i];if(h>g){d=h-p[i];b=h-g;c=g-d;if(c>b){for(j=n;j<=i;j++)a[j][k]=0;fano(p,a,n,i,k+1);for(j=i+1;j<=m;j++)a[j][k]=1;fano(p,a,i+1,m,k+1);}else{for(j=n;j<=i-1;j++)a[j][k]=0;fano(p,a,n,i-1,k+1);for(j=i;j<=m;j++)a[j][k]=1;fano(p,a,i,m,k+1);}break;}}}}voidmain(){inti,j,k[N],n,flase=0;floatp[N],m,H=0.0,K=0.0,sum=0.0;cout<<"输入信源符号个数"<<endl;cin>>n;cout<<"输入各信源符号概率"<<endl;for(i=1;i<=n;i++){cin>>p[i];}for(i=1;i<=n;i++){sum=sum+p[i];}for(i=1;i<=n;i++){if(p[i]<0.0||p[i]>1.0||sum!=1.0){cout<<"inputgailverror!";flase=1;break;}if(flase==0){for(i=0;i<=n;i++)for(j=0;j<=n;j++){pa[i][j]=10;}fano(p,pa,1,n,1);cout<<"信源费诺编码如下：\n";for(i=1;i<=n;i++){k[i]=0;cout<<"x"<<i<<"="<<p[i]<<"\t码字为\t";for(j=1;j<=n;j++){if(pa[i][j]!=10){cout<<pa[i][j];k[i]++;}}cout<<"\t码长为\t"<<k[i]<<endl;}for(i=1;i<=n;i++){H=-(p[i]*log(p[i])/log(2))+H;}cout<<"信源熵 H(X)="<<H<<"(比特/符号)"<<endl;for(i=1;i<=n;i++){K=p[i]*k[i]+K;}cout<<"平均码长K="<<K<<"(比特/符号)"<<endl;cout<<"编码效率为"<<(H/K)*100<<"%"<<endl;}//if(flase==0)cout<<endl;}//main()借□IX攵s.号r\55/75211辙4辑B,B,B,0点1233>号为为为为W