哈夫曼编码的MATLAB实现（基于0、1编码）：

%哈夫曼编码的MATLAB实现（基于0、1编码）：clc;clear;A=[0.5,0.19,0.19,0.12];%信源消息的概率序列A=fliplr(sort(A));%按降序排列T=A;[m,n]=size(A);B=zeros(n,n-1);%空的编码表（矩阵）for i=1:n    B(i,1)=T(i);%生成编码表的第一列endr=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;endB;%输出编码表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;endA%排序后的原概率序列END%编码结果for i=1:n    [a,b]=size(char(END(i)));    L(i)=b;endavlen=sum(L.*A)%平均码长H1=log2(A);H=-A*(H1')%熵P=H/avlen%编码效率