霍夫曼码编码方式matlab,【Matlab编程】哈夫曼编码的Matlab实现

在前年暑假的时候，用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