在前年暑假的时候,用C实现了哈夫曼编译码的功能,见文章《哈夫曼树及编译码》。不过在通信仿真中,经常要使用到Matlab编程,所以为了方便起见,这里用Matlab实现的哈夫曼编码的功能。至于哈夫曼编译码的基本原理,我们可以参考之前的文章《哈夫曼树及编译码》,里面有详细的说明及图解过程。下面直接给出具体的Matlab实现的哈夫曼编码函数,由于程序中注释还算比较详细,在此就不予与说明:
function [ h,e ] = Huffman_code( p )
%p为概率分布,此函数功能是进行哈夫曼编码
% 此处显示详细说明
% h为各个元素的麻子
% e为输出的平均码长
if length(find(p<0))~=0
error('概率不应该小于0!')
end
if abs(sum(p)-1)>10e-10
error('概率之和大于1,请检查输入!')
end
n=length(p);
p=sort(p)
q=p;
m=zeros(n-1,n);
for i=1:n-1
[q,e]=sort(q);
m(i,:)=[e(1:n-i+1),zeros(1,i-1)]; %由数组l 构建一个矩阵,该矩阵表明概率合并时的顺序,用于后面的编码
q=[q(1)+q(2),q(3:n),1];
end
for i=1:n-1
c(i,1:n*n)=b