信息論仙農(nóng)編碼實(shí)驗(yàn)報(bào)告new

仙農(nóng)編碼實(shí)驗(yàn)報(bào)告姓名:資雁兵學(xué)號(hào):20091616310029班級(jí):信息安全一班題目：仙農(nóng)編碼實(shí)驗(yàn)要求用仙農(nóng)編碼實(shí)現(xiàn)對(duì)信源的編碼實(shí)驗(yàn)步驟將信源符號(hào)按概率從大到小的排列，為方便起見，令P(s1)≥p(s2)≥…≥p(sn)令p(s0)=0,用pa(sj)表示第i個(gè)碼字的累加概率，其中j=i+1，則：pa(sj)=p(s0)+…+p(sj-1)確定滿足下列不等式的整數(shù)li，并令li為第i個(gè)碼字的長度為-logp(si)≤li＜1-logp(si)將pa(si)和j=i+1用二進(jìn)制表示，并取小數(shù)點(diǎn)后li位作為符號(hào)si的編碼設(shè)計(jì)函數(shù)voidGetOrder(shannongsn[50],intn)//從大到小排序voidGetCode(shannongsn[50],intn)//獲得編碼程序使用方法運(yùn)行程序，先輸入所要編碼的個(gè)數(shù)，然后一次輸入一個(gè)字符作為信源名，以及一個(gè)信源的概率，即可得編碼源程序#include<stdio.h>#include<math.h>#defineMAX100typedefstruct{ floata;//概率 floatb;//碼長intcode[MAX];//編碼數(shù)組}xiannong;//結(jié)構(gòu)體voidScanfNode(xiannongxiannode[MAX],intn){inti;//輸入各個(gè)概率 for(i=0;i<n;i++){ printf("請(qǐng)輸入編碼的概率\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[]數(shù)組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;//把小數(shù)轉(zhuǎn)換到二進(jìn)制 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++){//生成二進(jìn)制 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;} } }}/***************主函數(shù)****************/voidmain(){xiannongxiannode[MAX];floatpa[MAX];intn,i,j;printf("請(qǐng)輸入編碼的個(gè)數(shù)\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]);}測(cè)試專用*/printf("\n");for(i=0;i<n;i++){ printf("第%d個(gè)仙農(nóng)編碼",i+1); for(j=0;j<xiannode[i