《哈夫曼编码的MATLAB实现.pdf》由会员分享,可在线阅读,更多相关《哈夫曼编码的MATLAB实现.pdf(3页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、%哈夫曼编码的 MATLAB 实现(基于 0、1 编码):clc;clear;A=,;信源消息的概率序列A=fliplr(sort(A);%按降序排列T=A;m,n=size(A);B=zeros(n,n-1);%空的编码表(矩阵)for i=1:nB(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:tB(i,j)=T(i);endK=find(T=r);B(n,j)=K(end);%从第二列开始,每列的最
2、后一个元素记录特征元素在%该列的位置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-2if i1&B(i,j)=B(i-1,j)d=d+1;elsed=1;endB(B(n,j+1),j+1)=-1;temp=B(:,j+1);x=find(temp=B(i,j);END(i)=END1(x(d);endy=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:na,b=size(char(END(i);L(i)=b;endavlen=sum(L.*A)%平均码长H1=log2(A);H=-A*(H1)%熵P=H/avlen%编码效率