function [h,l]=huffman(p)
if (length(find(p<0))~=0) %判断输入P矩阵各元素是否全为大于零的有效概率值
error('Not a prob,negative component');
end
if (abs(sum(p)-1)>10e-10) %判断输入矩阵的概率和是否为合理值1
error('Not a prob.vector,component do not add to 1')
end
n=length(p); %N取输入行向量P的长度,即需要编码元素个数。
q=p;
m=zeros(n-1,n);%M为N-1行、N列矩阵,用来记录每行最小两概率叠加后概率排列次序。
%确定概率大小值的排列,得到M矩阵。
for i=1:n-1
[q,l]=sort(q);
m(i,:)=[l(1:n-i+1),zeros(1,i-1)];
q=[q(1)+q(2),q(3:n),1];
end
%生成一个N-1行、N2(N×N)列矩阵C,每行可看作N个段,每段长为N,记录一个码字
%(每个码字的长度不会超过N)。给C矩阵的N-1行的第一个段赋值0,第二个段赋值1,
%这两个码字对应编码中最后相加为一的两个概率。
for i=1:n-1
c(i,:)=blanks(n*n);
end
c(n-1,n)='0';
c(n-1,2*n)='1';
% 循环N-2次,决定矩阵C从倒数第二行开始到第一行的每段的码字值。每一行值都从下一行值得到,
% 找到在下一行码字中相加本行最小两个概率得到的概率的对应码字,本行两个最小概率对应码字分
%别为此码字最后加“0”,加“1”。
for i=2:n-1
c(n-i,1:n-1)=c(n-i+1,n*(find(m(n-i+1,:)==1))...
-(n-2):n*(find(m(n-i+1,:)==1)));
c(n-i,n)='0';
c(n-i,n+1:2*n-1)=c(n-i,1:n-1);
c(n-i,2*n)='1';
%找到其余的本行在下一行对应的码字,该码字保持不变。循环结束后,C矩阵第一行的N段对应
%输入N个概率所对应符号的码字。该码字按码字长短排列。
for j=1:i-1
c(n-i,(j+1)*n+1:(j+2)*n)=c(n-i+1,...
n*(find(m(n-i+1,:)==j+1)-1)+1:n*find(m(n-i+1,:)==j+1));
end
end
%根据M矩阵第一行记录的概率排序位置分配给每个概率对应符号的码字。
for i=1:n
h(i,1:n)=c(1,n*(find(m(1,:)==i)-1)+1:find(m(1,:)==i)*n);
ll(i)=length(find(abs(h(i,:))~=32));
end
l=sum(p.*ll);