多元霍夫曼编码的matlab实现

%哈夫曼编码：
clcclear
P=[0.36 0.14 0.13 0.12 0.10 0.09 0.04 0.02]; 
r=2; 
p=P;
n=length（p）;
s=ceil（（n-r）/（r-1））;
a=（r-1）*s+r-n;
p=[pzeros（1a）];
n=length（p）;
p=fliplr（sort（p））;
T=p;
B=zeros（ns+1）;
B（:1）=T;
 temp=sum（B（n-r+1:end1））;
T（n-r+1）=temp;
 for i=n-r+2:n
     T（i）=0;
 end
 T=fliplr（sort（T））;
 t=n-r+1;
for j=2:s+1 
    for i=1:t
        B（ij）=T（i）;
    end
        K=find（T==temp）;
        B（nj）=K（end）; 
        temp=sum（B（t-r+1:tj））;
        T（t-r+1）=temp;
        T=fliplr（sort（T））; 
        t=t-r+1;
end