huffman編碼與譯碼 代碼

《huffman編碼與譯碼 代碼》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(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之