2021年整理huffman的MATLAB程序.pptxVIP

%%%%%%%%数字图像Huffman 压缩编码实现%%%%%%%%%%%%%%%%%%% clc clear I=[0 1 3 2 1 3 2 1;0 5 7 6 2 5 6 7;1 6 0 6 1 6 3 4;2 6 7 5 3 5 6 5;3 2 2 7 2 6 1 6; 2 6 5 0 2 7 5 0;1 2 3 2 1 2 1 2;3 1 2 3 1 2 2 1]; %读入一幅图像的灰度值 % load I; [m,n]=size(I); %%%%%%将矩阵的不同数统计在数组c 的第一列中%%%%%%%%%%%% p1=1;s=m*n; for k=1:m for L=1:n f=0; for b=1:p1-1 if(c(b,1) == I(k,L)) f=1;break;end end if(f==0)c(p1,1)=I(k,L);p1=p1+1;end end end %%%%%%%%%%%将相同的数占站整个数组总数的比例统计在数组p 中%%%%%%% for g=1:p1-1 p(g)=0;c(g,2)=0; for k=1:m for L=1:n if(c(g,1) == I(k,L)) p(g)=p(g)+1;end end end p(g)=p(g)/s; end p11=p; %%%%%%%%找到最小的概率,相加直到等于1,把最小概率的序号存在tree 第一列中,次小放在第二列,和放 在p 像素比例之后%%%%%%% pn=0;po=1; while(1) if(pn = 1.0)break; else [pm,p2]=min(p(1:p1-1));p(p2)=1.1; [pm2,p3]=min(p(1:p1-1));p(p3)=1.1; pn=pm+pm2;p(p1)=pn; tree(po,1)=p2;tree(po,2)=p3; po=po+1;p1=p1+1; end end %%%%%%%%%C 数组第一维表示值,第二维表示代码数值大小,第三维表示代码的位数po=8 g=8 for k=1:po-1;tt=k;m1=1; if(or(tree(k,1) = g,tree(k,2) = g)) if(tree(k,1) = g) c(tree(k,1),2)=c(tree(k,1),2)+m1; m2=1; while(tt po-1) m1=m1*2; for L=tt:po-1 if(tree(L,1) == tt+g) c(tree(k,1),2)=c(tree(k,1),2)+m1; m2=m2+1;tt=L;break; elseif(tree(L,2) == tt+g) m2=m2+1;tt=L;break; end end end c(tree(k,1),3)=m2; end tt=k;m1=1; if(tree(k,2) g) m2=1; while(tt po-1) m1=m1*2; for L=tt:po-1 if(tree(L,1) == tt+g) c(tree(k,2),2)=c(tree(k,2),2)+m1; m2=m2+1;tt=L;break; elseif(tree(L,2) == tt+g) m2=m2+1;tt=L;break; end end end c(tree(k,2),3)=m2; end end end %%%%%%%%%把概率小的值为 1 标识，概率大的值为 0 标识%%%%%% [M,N]=size(c); disp(编码) A1=dec2bin(c(1,2),c(1,3)) %%%说明：这里可以把编码存在高维数组或构架数组、元胞数组同时显示。 A2=dec2bin(c(2,2),c(2,3)) A3=dec2bin(c(3,2),c(3,3)) A4=dec2bin(c(4,2),c(4,3)) A5=dec2bin(c(5,2),c(5,3));A6=dec2bin(c(6,2),c(6,3)) A7=dec2bin(c(7,2),c(7,3)) A8=dec2bin(c(8,2),c(8,3)) for m=1:M if (p11(m)~=0) H(m)=-p11(m)*log2(p11(m)); end end disp( 信 源 的 熵 ) H1=sum(H) %信源的熵 NN=0; for i=1:M NN=NN+p11(1,i)*c(i,3); %平均码长 end disp(平均码长) NN disp( 编 码 效 率 ) yita=H1/(NN*log2(2)) %效率 disp(冗余度) Rd=1-yita %冗余度