霍夫曼编码代码matlab,matlab 实现霍夫曼编码

编码内容包括：用matlab实现霍夫曼编码，并且求出相应的信源熵，平均码长，和编码效率。

以下是代码：

%霍夫曼编码

%huffman_code

%编码思路：根据位置矩阵的变化过程反推生成霍夫曼编码

n=input('请输入编码个数:\n');

c=[];

% for i=1:n

%     fprintf('请输入第%d个编码:\n',i);

%     c(i)=input('');

% end

c=rand(1,n);%********随机生成数据********************

c=c/sum(c);

p=c;

disp(p);

w=[;];%***************位置矩阵，记录数据的变化过程*****

for i=1:n-1

[p,l]=sort(p);

w(i,:)=[l(1:n-i+1),zeros(1,i-1)];

p=[p(1)+p(2),p(3:n),1];

end

a={;};%**********字符数组，生成霍夫曼码****************

for i=1;2*(n-1)

for j=1:n

a{i,j}=' ';

end

end%*****************初始化*********************

a{2*n-2,1}='0';

a{2*n-2,2}='1';

for i=1:n-1

for j=1:i+1

a{2*(n-i)-1,j}=a{2*(n-i),find(w(n-i,:)==j)};

end

if((n-i-1)<=0)

break;

end

a{2*(n-i-1),1}=strcat(a{2*(n-i)-1,1},'0');

a{2*(n-i-1),2}=strcat(a{2*(n-i)-1,1},'1');

for j=3:i+2

a{2*(n-i-1),j}=a{2*(n-i)-1,j-1};

end

end

h=-c.*log2(c);

hx=sum(h);

for i=1:n

len(i)=length(a{1,i});

end

plen=len.*c;

averagelen=sum(plen);

fprintf('霍夫曼编码为:\n');

% disp(a);

for i=1:n

disp(a(1,i));

end

fprintf('平均码长为:\n');

disp(averagelen);

fprintf('信源熵为:\n');

disp(hx);

fprintf('编码效率为:\n');

disp(hx/averagelen);