huffman的MATLAB程序參考模板

1、1 / 4%數(shù)字圖像Huffman壓縮編碼實(shí)現(xiàn)% 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); %將矩陣的不同數(shù)統(tǒng)計(jì)在數(shù)組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 %將相同的數(shù)占站整個(gè)數(shù)組總數(shù)的比例統(tǒng)計(jì)在數(shù)組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,把最小概率的序號(hào)存

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數(shù)組第一維表示值,第二維表示代碼數(shù)值大小,第三維表示代碼的位數(shù)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標(biāo)識(shí)，概率大的值為0標(biāo)識(shí)% M,N=size(c); disp(''編碼'') A1=dec2bin(c(1,2),c(1,3)    %說(shuō)明：這里可以把編碼存在高維數(shù)組或構(gòu)架數(shù)組、元胞數(shù)組同時(shí)顯示。  A2=dec2bin(c(2,2),c(2,3) A3=dec2bin(c(3,2)