信道编码Huffman编码信息论与编码

信道编码Huffman编码信息论与编码

MATLAB做的信息论与编码的实验一小部分 结果编码是数字数组，还可计算平均码长，编码效率，编码后信息传输速率 话不多说，直接代码

function [W,L,q,r]=huffman(P) %哈夫曼编码 n=length(P); p=P; mark=zeros(n-1,n); for i=1:n-1 [p,num]=sort(p); mark(i,:)=[num(1:n-i+1),zeros(1,i-1)]; p=[p(1)+p(2),p(3:n),1]; end for i=1:n-1 table(i,:)=blanks(n*n); end table(n-1,n)='1'; table(n-1,2*n)='0'; for i=2:n-1 table(n-i,1:n-1)=table(n-i+1,n*(find(mark(n-i+1,:)==1))-(n-2):n*(find(mark(n-i+1,:)==1))); table(n-i,n)='1'; table(n-i,n+1:2*n-1)=table(n-i,1:n-1); table(n-i,2*n)='0'; for j=1:i-1 table(n-i,(j+1)*n+1:(j+2)*n)=table(n-i+1,n*(find(mark(n-i+1,:)==j+1)-1)+1:n*find(mark(n-i+1,:)==j+1)); end end for i=1:n W(i,1:n)=table(1,n*(find(mark(1,:)==i)-1)+1:find(mark(1,:)==i)*n); l(i)=length(find(abs(W(i,:))~=32)); end L=sum(P.*l); sum1=0; for i=1:n %使用for循环进行信息熵、平均码长求解 a(i)=-log2(P(i)); %a(i)表示单个信源的自信息量 K(i)=ceil(a(i)); %K(i)表示对自信息量向上取整 c(i)=a(i)*P(i); sum1=sum1+c(i); %信息熵 end H=sum1; r=H/L; q=H/L; [m,n]=size(W); W=reshape(W',1,m*n); W=double(W)-double('0') W(find(W==-16))=[]; s1 = 'Huffman 编码所得码字 W:'; s2=' Huffman 编码的平均码字长度L:'; s3=' Huffman 编码的编码效率 Q:'; s4='Huffman 编码后信息传输速率 R:'; disp(s1),disp(W); disp(s2),disp(L); disp(s3),disp(q); disp(s4),disp(r);

上面是哈夫曼的编码函数，在主程序里面，直接调用就行

[w,l,q,r]=huffman(S);%w为编码，L为平均长度，q为编码效率,r编码后信息传输速率

最后显示一下结果 用到的离散无记忆信源S=[0.4,0.2,0.2,0.1,0.1];%离散无记忆信源

Huffman 编码所得码字 W: 1 0 0 0 0 1 0 0 1 1 0 0 1 0 Huffman 编码的平均码字长度L: 2.2000 Huffman 编码的编码效率 Q: 0.9645 Huffman 编码后信息传输速率 R: 0.9645

我没写单位，要加的话在函数中每一行disp后面再加一个disp就行 算了，我加上吧

disp(s1),disp(W); disp(s2),disp(L,'(bit/sign)'); disp(s3),disp(q); disp(s4),disp(r,'(bit/CodeSign)');