信息论仙农编码实验报告new

仙农编码实验报告姓名:资雁兵学号:20091616310029班级:信息安全一班题目：仙农编码实验要求用仙农编码实现对信源的编码实验步骤将信源符号按概率从大到小的排列，为方便起见，令P(s1)≥p(s2)≥…≥p(sn)令p(s0)=0,用pa(sj)表示第i个码字的累加概率，其中j=i+1，则：pa(sj)=p(s0)+…+p(sj-1)确定满足下列不等式的整数li，并令li为第i个码字的长度为-logp(si)≤li＜1-logp(si)将pa(si)和j=i+1用二进制表示，并取小数点后li位作为符号si的编码设计函数voidGetOrder(shannongsn[50],intn)//从大到小排序voidGetCode(shannongsn[50],intn)//获得编码程序使用方法运行程序，先输入所要编码的个数，然后一次输入一个字符作为信源名，以及一个信源的概率，即可得编码源程序#include<stdio.h>#include<math.h>#defineMAX100typedefstruct{ floata;//概率 floatb;//码长intcode[MAX];//编码数组}xiannong;//结构体voidScanfNode(xiannongxiannode[MAX],intn){inti;//输入各个概率 for(i=0;i<n;i++){ printf("请输入编码的概率\n"); scanf("%f",&xiannode[i].a);printf("\n"); }}voidGetNode(xiannongxiannode[MAX],intn){ inti;//从大到小排序 intj; floatk=1.0; for(i=0;i<n;i++){ for(j=i;j<n;j++) if(xiannode[j].a>xiannode[i].a){ k=xiannode[j].a; xiannode[j].a=xiannode[i].a; xiannode[i].a=k; } }}voidAddOrder(xiannongxiannode[MAX],intn,floatpa[MAX]){//累加到pa[]数组intj;pa[1]=0;for(j=1;j<=n;j++){pa[j+1]=pa[j]+xiannode[j-1].a; printf("%f",pa[j]);}}voidGetOrder(xiannongxiannode[MAX],intn,floatpa[MAX]){inti,k,m;//把小数转换到二进制 floatj; for(i=0;i<n;i++){//得到编码的码长 j=-log(xiannode[i].a)/log(2.0); xiannode[i].b=j; //printf("%f\n",xiannode[i].b);m=j; if(m==j) xiannode[i].b=m; else xiannode[i].b=m+1; }// for(i=1;i<=n;i++){//生成二进制 for(k=0;k<xiannode[i-1].b;k++){ pa[i]=pa[i]*2; if(pa[i]<1) xiannode[i-1].code[k]=0; else{xiannode[i-1].code[k]=1; pa[i]=pa[i]-1;} } }}/***************主函数****************/voidmain(){xiannongxiannode[MAX];floatpa[MAX];intn,i,j;printf("请输入编码的个数\n");scanf("%d",&n);ScanfNode(xiannode,n);GetNode(xiannode,n);AddOrder(xiannode,n,pa);GetOrder(xiannode,n,pa);/*for(i=0;i<n;i++){printf("%3f",xiannode[i].a);}printf("\n");for(i=0;i<n;i++){printf("%f",xiannode[i].b);}for(i=1;i<=n;i++){printf("%f",pa[i]);}测试专用*/printf("\n");for(i=0;i<n;i++){ printf("第%d个仙农编码",i+1); for(j=0;j<xiannode[i