能使霍夫曼编码在matlab上得到实现,从而进一步了解霍夫曼编码的原理

				%哈夫曼编码的MATLAB实现（基于0、1编码）：
				clc;
				clear;
				A=[5,3,1,6,2];%原概率序列
				A=A/sum(A);
				A=fliplr(sort(A));%按降序排列
				T=A;
				[m,n]=size(A);
				B=zeros(n,n-1);%空的编码表（矩阵）
				for i=1:n
				    B(i,1)=T(i);%生成编码表的第一列
				end
				r=B(i,1)+B(i-1,1);%最后两个元素相加
				T(n-1)=r;
				T(n)=0;
				T=fliplr(sort(T));
				t=n-1;
				for j=2:n-1%生成编码表的其他各列
				    for i=1:t
				        B(i,j)=T(i);
				    end
				        K=find(T==r);
				        B(n,j)=K(end);%从第二列开始，每列的最后一个元素记录特征元素在
				%该列的位置
				        r=(B(t-1,j)+B(t,j));%最后两个元素相加
				        T(t-1)=r;
				        T(t)=0;
				        T=fliplr(sort(T)); 
				        t=t-1;
				end
				B;%输出编码表
				END1=sym('[0,1]');%给最后一列的元素编码
				END=END1;
				t=3;
				d=1;
				for j=n-2:-1:1%从倒数第二列开始依次对各列元素编码
				    for i=1:t-2
				        if i>1 & B(i,j)==B(i-1,j)
				            d=d+1;
				        else
				            d=1;
				        end
				        B(B(n,j+1),j+1)=-1;
				        temp=B(:,j+1);
				
				        x=find(temp==B(i,j));
				        END(i)=END1(x(d));
				    end
				    y=B(n,j+1);
				    END(t-1)=[char(END1(y)),'0'];
				    END(t)=[char(END1(y)),'1'];
				    t=t+1;
				    END1=END;
				end
				    A%排序后的原概率序列
				    END%编码结果
				for i=1:n
				    [a,b]=size(char(END(i)));
				    L(i)=b;
				end
				avlen=sum(L.*A)%平均码长  
				H1=log2(A);
				H=-A*(H1')%熵
				P=H/avlen%编码效率