huffman的MATLAB程序参考模板

1、1 / 4%数字图像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、#160;  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  &

3、#160;  for L=1:n         f=0;         for b=1:p1-1             if(c(b,1) = I(k,L) f=1;break;end      

4、;   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       

5、60; 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,把最小概率的序号存

6、在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;  

7、;       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-

8、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,

9、1),2)+m1;             m2=1;             while(tt < po-1)                 m1=m1*2;  

10、;               for L=tt:po-1                     if(tree(L,1) = tt+g)       

11、60;                 c(tree(k,1),2)=c(tree(k,1),2)+m1;                         m2=m2+1;tt=L;brea

12、k;                     elseif(tree(L,2) = tt+g)                        

13、0;m2=m2+1;tt=L;break;                     end                 end          

14、   end             c(tree(k,1),3)=m2;         end         tt=k;m1=1;         if(tree(k,2) < 

15、g)             m2=1;             while(tt < po-1)                 m1=m1*2;   

16、              for L=tt:po-1                     if(tree(L,1) = tt+g)        

17、0;                c(tree(k,2),2)=c(tree(k,2),2)+m1;                         m2=m2+1;tt=L;break; 

18、0;                   elseif(tree(L,2) = tt+g)                         m2=m2

19、+1;tt=L;break;                     end                 end           &

20、#160; 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)