Huffman编码与译码代码.docVIP

% 1 完成对输入的序列进行各个码元的概率统计； %完成对字符串中的各字符的统计，并列出其概率分布矩阵，返回pro矩阵 %s:待编码序列，S：所含的码元序列 function pro=getpro(s) pro=[]; a=length(s); S=unique(s); b=length(S); c=zeros(1,b); %用以存放个序列中各个码元的个数； %///////////////////////////////////////////// %进行概率计算； for i=1:b for j=1:a if S(i)==s(j) c(i)=c(i)+1; else continue; end; end; end; pro=c./a; disp(S); disp(pro); %完成对已知编码序列的译码，以及在改变码表中的某一位值得情况下，再一次译码，计算其误码率； %Codenumber:已编码序列； huffmantable:码表；Code2:各码元的码长；pro2：各码元的概率分布矩阵； %s:原始序列；decodenumber： 译码序列； function decodenumber=huffmandecode(Codenumber,huffmantable,Code2,pro2,bit) mm=unique(bit); %码元序列 mm=mm(pro2(2,:)); [lx,ly]=size(huffmantable); LL=size(Codenumber,2); decodenumber=[]; ZF=[Codenumber,-ones(1,max(Code2))]; for j=1:length(bit) for i=1:lx k=Code2(i); while(ZF(1:k)==huffmantable(i,1:k)) decodenumber=[decodenumber,mm(i)]; ZF=ZF(k+1:end); end end end disp(译码序列如下：); disp(decodenumber); disp(原始序列如下：); disp(bit); end %对于译码部分所用到的部分主要是编码时生成的码表以及huffmantree，在进行编码的时候， %通过筛选后的源字符串的字符序列的下标，与码表中的每行相对应的原则，遍历编码序列； %在遍历的时候，通过码表中各行的码长，控制遍历的长度，与每行中的码表进行比较，输出相对应的字符，即完成了译码。 % 3 编码，完成构成完整的基本树，形成huffmantree，并对其遍历，得到码表，通过码表对序列进行编码，并求其平均码长，信息熵，编码效率； % HuffmanTree 哈夫曼树，pro 码元概率分布矩阵，S，待编码序列； function [Codenumber,huffmantable,Code2]=huffmanencode(HuffmanTree,pro,bit) p=pro; a1=unique(bit); a2=bit; len1=length(a1); %码元序列的长度; len2=length(a2); %所要编码序列的长度; a3=zeros(1,len2); %生成与编码序列长度一样的零矩阵，用以存放编码对应的下标; Code=[]; %存放其遍历一个码元所对应的编码; Code2=[]; %记录各个码元的码长； Lastnumber=1; Sumnumber=0; %累积计算编码总长度； huffmantable=-ones(len1,len1); %建立len1*len1的单位负矩阵，用以存放其码表； %遍历二叉树，生成码表huffmantable; for i=1:len1; %循环完成len1个符号的编码； k=pro(2,i); key=1;