信息論與編碼課程設(shè)計(jì)-資料下載頁
2025-07-24 07:20本頁面
【正文】 0 j=j+1。 str2(j)=sort_str1(i)。 end str2(j+1)=sort_str1(i+1)。endfor i = 1:length(str2)%length函數(shù)獲取字符串長度 str_num =strfind(sort_str1,str2(i))。%strfind(S1,S2):尋找S2是否匹配S1，并返回S2的位置 count1(i) = length(str_num)。endstr2count=count1(3:end)p=count./sum(count)sum(p.*log2(p))%計(jì)算信源熵香農(nóng)編碼程序function c=shannon(p)% p=[ ]。% shannon(p)。[p,index]=sort(p)。p=fliplr(p)。%從大到小n=length(p)。pa=0。%累加概率for i=2:n pa(i)=pa(i1)+p(i1)。 endk=ceil(log2(p))。%碼長計(jì)算c=cell(1,n)。%生成元胞數(shù)組，存碼字，是cell，跟上一行不一樣for i=1:n c{i}= 39。39。 tmp=pa(i)。 for j=1:k(i) tmp=tmp * 2。 if tmp=1 tmp=tmp 1。 c{i}(j)= 39。139。 else c{i}(j)= 39。039。 end endend%p%pa%交換回原來的順序c=fliplr(c)。c(index)=c。fprintf(39。信源信息熵:\n39。)。H=sum(p.*log2(p))%計(jì)算信源熵fprintf(39。平均碼長:\n39。)。K=sum(p.*k)%計(jì)算平均碼長fprintf(39。編碼效率:\n39。)。w=H./K%計(jì)算編碼效率fprintf(39。碼字:\n39。)。shannon(p)。信源信息熵:H = 費(fèi)諾編碼程序主程序function c=fano1(p)% p=[ ]% c=fano1(p)n=size(p,2)。if n==1 c=cell(1,1)。 c{1}=39。39。 returnend[p,index]=sort(p)。%按概率排序p=fliplr(p)。total=sum(p)。%總概率acc=0。%累積概率flag=0。%是否到達(dá)尾部的標(biāo)志for i=1:n1 newacc=acc+p(i)。 if abs(total2 * newacc)=abs(total 2*acc) flag=1。 break。 end acc=newacc。endif ~flag i=n。endsplit=i。%從分界點(diǎn)對(duì)兩邊的碼遞歸做fanoc1=fano1(p(1:split1))。c2=fano1(p(split:n))。c=cell(1,n)。%添加前綴0，1for i=1:split1 c{i}=strcat(39。039。,c1{i})。endfor i = split:n c{i}=strcat(39。139。,c2{isplit+1} )。end%將順序調(diào)整回去c=fliplr(c)。c(index)=c。子程序function []=fano2(c,p)for i=1:length(c)%求平均碼長 count(i)=length(cell2mat(c(i)))。 endfprintf(39。信源信息熵:\n39。)。H=sum(p.*log2(p))%計(jì)算信源熵fprintf(39。平均碼長:\n39。)K=sum(count.*p)%計(jì)算平均碼長fprintf(39。編碼效率:\n39。)w=H./K%計(jì)算編碼效率fprintf(39。碼字:\n39。)c16
語文相關(guān)推薦
鄂ICP備17016276號(hào)-1