|
#include<stdio.h> #include<math.h> //使用各种数学函数 #define N 7 main() { int i,j; double p[N]={0},max,sum=0,R,H=0,K1=0,n; double q[N]={0};//概率数组 int length[N];//码长数组 double K[N]={0},I[N]={0};//求码长的中间数组 int a[N][N]={0};//二进制存储数组 for (i=0;i<N;i++) scanf("%lf",&p[ i]); printf("\n"); //判断输入概率是否正确 for(i=0;i<N;i++) { sum+=p[ i]; } printf("概率之和=%lf\n",sum); printf("\n"); if(sum==1.000000) { printf("输入概率p(i)：\n"); for(i=0;i<N;i++) { printf("p[%d]=%lf\n",i,p[ i]); } printf("\n"); } else {printf("error\n");return 0;} //将概率进行排序操作 for(i=0;i<N;i++) { for(j=i+1;j<N;j++) { if(p[ i]<p[j]) { max=p[ i]; p[ i]=p[j]; p[j]=max; } } } //输出排序后的数组 printf("排序：\n"); for(i=0;i<N;i++) { printf("p[%d]=%lf\n",i,p[ i]); } printf("\n"); //求累加概率 for(i=0;i<N;i++) { q[i+1]=q+p[ i];//从第二个开始加 } //输出累加概率数组 printf("累加概率：\n"); for(i=0;i<N;i++) { printf("q[%d]=%lf\n",i,q); } printf("\n"); //求-logp（i） for(i=0,j=0;i<N;i++,j++) { K[j]=-log(p[ i])/log(2); } printf("-logp(i): \n"); for(j=0;j<N;j++) printf("-logp(%d)=%.2lf\n",j,K[j]); printf("\n"); //求码字长度 for(i=0,j=0;i<N;i++,j++) { K[j]=-log(p[ i])/log(2); I[j]=K[j]+1; if(I[j]>(int)(I[j]))//香农编码的第三个步骤 length[j]=(int)I[j]; else length[j]=K[j]; K1+=p[ i]*length[j];//求平均码长 } //输出码字长度 printf("码长：\n"); for(i=0;i<N;i++) { printf("length[%d]=%d\n",i,length); } printf("\n"); //求累加概率的二进制形式，即码字 for(i=0;i<N;i++) { for(j=0;j<length;j++) { q=q*2; if(q>=1) { q=q-1; a[j]=1; } else a[j]=0; } } //输出码字 printf("码字：\n"); for(i=0;i<N;i++) { for(j=0;j<length;j++) printf("%d",a[j]); printf("\n"); } //求信源熵 编码效率 for(i=0,j=0;i<N;i++,j++) { K[j]=-log(p[ i])/log(2); H+=p[ i]*K[j]; n=H/K1; } printf("H(x)=%.2lf,n=%lf",H,n); }
| 
QQ好友和群
QQ空间
腾讯微博
腾讯朋友
收藏
淘帖
顶
踩
