mat lab实现hafman三元编码_专业指导代码类资源

% 哈弗曼编码 clc clear % A=input('Enter matrix A=>'); % A=[0.2,0.19,0.18,0.17,0.15,0.1,0.01]; A=[0.4,0.2,0.1,0.1,0.05,0.05,0.05,0.05,0]; l=length(A); sum0=0; for x=1:l if A(x)>1 disp('输入矩阵概率'); break; else sum0 = sum0 +A(x); end end if sum0 ~= 1 disp('输入矩阵不能归一化'); else A=fliplr(sort(A)); T=A; [m,n]=size(A); q=ceil((n-3)/2); B=zeros(n,q+1); for i=1:n B(i,1)=T(i);%生成编码表的第一列 end r=B(i,1)+B(i-1,1)+B(i-2,1);%最后三个元素相加 T(n-2)=r; T(n-1)=0; T(n)=0; T=fliplr(sort(T)); t=n-2; for j=2:q+1%生成编码表的其他各列 for i=1:t B(i,j)=T(i); end K=find(T==r); B(n,j)=K(end);%从第二列开始，每列的最后一个元素记录特征元素在该列的位置 r=(B(t-2,j)+B(t-1,j)+B(t,j));%最后三个元素相加 T(t-2)=r; T(t-1)=0; T(t)=0; T=fliplr(sort(T)); t=t-1; end END1=sym('[0,1,2]');%给最后一列的元素编码 END=END1; t=5; d=1;%树的广度 for j=q:-1:1%从倒数第二列开始依次对各列元素编码 for i=1:t-3 if i>1 & B(i,j)==B(i-1,j) d=2; else d=1; end B(B(n,j+1),j+1)=-1;%将组合码赋值为-1从而在编码时将其省略 temp=B(:,j+1); x=find(temp==B(i,j));%确定待编码的B(i,j)在左支还是右支 END(i)=END1(x(d));%x是记录前一列元素在后一列元素中左右位置的矩阵 % p=[d,x(d),t] % a=B(i,j) end %衍生处再生成新码字 y=B(n,j+1);%为初始化后面的元素做准备 END(t-2)=[char(END1(y)),'0']; END(t-1)=[char(END1(y)),'1']; END(t)=[char(END1(y)),'2']; t=t+2; END1=END; end end A%排序后的原概率序列 % B%输出编码表 % END%编码结果 for i=1:n [a,b]=size(char(END(i))); L(i)=b; end avlen=sum(L.*A)%平均码长 der=0;%初始化方差 for i=1:n der=der+A(i)*power((L(i)-avlen),2); end H1=log2(A); H=-A*(H1')%熵 P=H/avlen/log2(3)%编码效率 der