最新信息論與編碼課程設(shè)計報告統(tǒng)計信源熵與香農(nóng)編碼-資料下載頁
2025-06-28 05:30本頁面
【正文】 next。 printf(num\tgailv\tsum\tlb(p(ai))\tlenth\tcode\n)。 printf(\n)。 for(i=0。imax_PN。i++) { printf(a%d\t%\t%\t%f\t%d\t,i,ppb,pp_sum,*log10(ppb),pkl)。 j=0。 for(j=0。jpkl。j++) printf(%d,pcode[j])。 printf(\n)。 hx=hxppb**log10(ppb)。 /*計算消息序列的熵*/ KL=KL+pkl*ppb。 /*計算平均碼字*/ p=pnext。 } printf(H(x)=%f\tKL=%f\nR=%fbit/code,hx,KL,hx/KL)。 /*計算編碼效率*/ } shnolist *setnull() { shnolist *head。 head=(shnolist *)malloc(sizeof(shnolist))。 headnext=NULL。 return(head)。 } shnolist *my_creat(datatype a[],int n) { shnolist *head,*p,*r。 int i。 head=setnull()。 r=head。//87 for(i=0。in。i++) { p=(shnolist *)malloc(sizeof(shnolist))。 ppb=a[i]。 pnext=NULL。 rnext=p。 r=p。 } return(head)。 } void valuelist(shnolist *L) { shnolist *head,*p。 int j=0。 int i。 datatype temp,s。 head=L。 p=headnext。 temp=0。 while(jmax_PN) { pp_sum=temp。 temp=temp+ppb。 pkl=*log10(ppb)+1。 /*編碼*/ { s=pp_sum。 for(i=0。ipkl。i++) pcode[i]=0。 for(i=0。ipkl。i++) { pcode[i]=2*s。 if(2*s=1) s=2*s1。 else if(2*s==0) break。 else s=2*s。 } } j++。 p=pnext。 } } int main(void) { shnolist *head。 pb_scan()。 pb_sort()。 head=my_creat(sym_arry,max_PN)。 valuelist(head)。 codedisp(head)。 }
研究報告相關(guān)推薦
鄂ICP備17016276號-1