資源描述:
《huffman編碼與譯碼 代碼》由會(huì)員上傳分享,免費(fèi)在線(xiàn)閱讀,更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫(kù)。
1、%1完成對(duì)輸入的序列進(jìn)行各個(gè)碼元的概率統(tǒng)計(jì);%完成對(duì)字符串中的各字符的統(tǒng)計(jì),并列出其概率分布矩陣,返回pro矩陣%s:待編碼序列,S:所含的碼元序列functionpro=getpro(s)pro=[];a=length(s);S=unique(s);b=length(S);c=zeros(1,b);%用以存放個(gè)序列中各個(gè)碼元的個(gè)數(shù);%/////////////////////////////////////////////%進(jìn)行概率計(jì)算;fori=1:bforj=1:aifS(i)==s(j)c(i)=c(i)+1;elsecontinue;end;end
2、;end;pro=c./a;disp(S);disp(pro);%完成對(duì)已知編碼序列的譯碼,以及在改變碼表中的某一位值得情況下,再一次譯碼,計(jì)算其誤碼率;%Codenumber:已編碼序列;huffmantable:碼表;Code2:各碼元的碼長(zhǎng);pro2:各碼元的概率分布矩陣;%s:原始序列;decodenumber:譯碼序列;functiondecodenumber=huffmandecode(Codenumber,huffmantable,Code2,pro2,bit)mm=unique(bit);%碼元序列mm=mm(pro2(2,:));[lx,l
3、y]=size(huffmantable);LL=size(Codenumber,2);decodenumber=[];ZF=[Codenumber,-ones(1,max(Code2))];forj=1:length(bit)fori=1:lxk=Code2(i);while(ZF(1:k)==huffmantable(i,1:k))decodenumber=[decodenumber,mm(i)];ZF=ZF(k+1:end);endendenddisp('譯碼序列如下:');disp(decodenumber);disp('原始序列如下:');disp
4、(bit);end%對(duì)于譯碼部分所用到的部分主要是編碼時(shí)生成的碼表以及huffmantree,在進(jìn)行編碼的時(shí)候,%通過(guò)篩選后的源字符串的字符序列的下標(biāo),與碼表中的每行相對(duì)應(yīng)的原則,遍歷編碼序列;%在遍歷的時(shí)候,通過(guò)碼表中各行的碼長(zhǎng),控制遍歷的長(zhǎng)度,與每行中的碼表進(jìn)行比較,輸出相對(duì)應(yīng)的字符,即完成了譯碼。%3編碼,完成構(gòu)成完整的基本樹(shù),形成huffmantree,并對(duì)其遍歷,得到碼表,通過(guò)碼表對(duì)序列進(jìn)行編碼,并求其平均碼長(zhǎng),信息熵,編碼效率;%HuffmanTree哈夫曼樹(shù),pro碼元概率分布矩陣,S,待編碼序列;function[Codenumber,huf
5、fmantable,Code2]=huffmanencode(HuffmanTree,pro,bit)p=pro;a1=unique(bit);a2=bit;len1=length(a1);%碼元序列的長(zhǎng)度;len2=length(a2);%所要編碼序列的長(zhǎng)度;a3=zeros(1,len2);%生成與編碼序列長(zhǎng)度一樣的零矩陣,用以存放編碼對(duì)應(yīng)的下標(biāo);Code=[];%存放其遍歷一個(gè)碼元所對(duì)應(yīng)的編碼;Code2=[];%記錄各個(gè)碼元的碼長(zhǎng);Lastnumber=1;Sumnumber=0;%累積計(jì)算編碼總長(zhǎng)度;huffmantable=-ones(len1,
6、len1);%建立len1*len1的單位負(fù)矩陣,用以存放其碼表;%遍歷二叉樹(shù),生成碼表huffmantable;fori=1:len1;%循環(huán)完成len1個(gè)符號(hào)的編碼;k=pro(2,i);key=1;m=find(HuffmanTree(6,1:len1)==k);while(HuffmanTree(5,m)~=1)%判斷是否遍歷到根結(jié)點(diǎn);Code(key)=HuffmanTree(4,m);key=key+1;m=HuffmanTree(3,m);%指向父節(jié)點(diǎn);endlc=length(Code);huffmantable(i,1:lc)=fliplr
7、(Code);%將Code矩陣中的編碼左右翻轉(zhuǎn),完成倒序排列;Code2=[Code2,lc];%將各個(gè)碼元的編碼長(zhǎng)度賦予矩陣Code2;end%顯示碼表,即只輸出矩陣huffmantable中非-1的部分;disp('碼表如下:');Code4=[];fori=1:len1%顯示其碼表;disp(a1(pro(2,i)));flag=1;while(huffmantable(i,flag)~=-1)Code4(flag)=huffmantable(i,flag);flag=flag+1;end;Code4,disp('碼長(zhǎng)='),Code2(i)end;%
8、通過(guò)以上生成的碼表,然后遍歷源字符串與unique之