基于DCT图像有损压缩MATLAB仿真

基于DCT图像有损压缩MATLAB仿真，按照JEPG标准，使用了霍夫曼编码，适合初学者学习，适合课程作业。

function bit_seq=ACEncoding（xy） 
Z=x;
v0=y;
R=abs（y）;
if R==1;amplen=1;
elseif（R >=   2 & R <=    3）;amplen = 2;
elseif（R >=   4 & R <=    7）;amplen = 3;
elseif（R >=   8 & R <=   15）;amplen = 4;
elseif（R >=  16 & R <=   31）;amplen = 5;
elseif（R >=  32 & R <=   63）;amplen = 6;
elseif（R >=  64 & R <=  127）;amplen = 7;
elseif（R >= 128 & R <=  255）;amplen = 8;
elseif（R >= 256 & R <=  511）;amplen = 9;
elseif（R >= 512 & R <= 1023）;amplen = 10;
end 
if （Z==0 & amplen==1）; codelen=2;accode=0;
elseif  （Z==0 & amplen==2）;codelen=2;accode=1;
elseif  （Z==0 & amplen==3）;codelen=3;accode=4;
elseif  （Z==0 & amplen==4）;codelen=4;accode=11;
elseif  （Z==0 & amplen==5）;codelen=5;accode=26;
elseif  （Z==0 & amplen==6）;codelen=7;accode=120;
elseif  （Z==0 & amplen==7）;codelen=8;accode=248;
elseif  （Z==0 & amplen==8）;codelen=10;accode=1014;
elseif  （Z==0 & amplen==9）;codelen=16;accode=65410;
elseif （Z==0 & amplen==10）;codelen=16;accode=65411;
    
elseif  （Z==1 & amplen==1）; codelen=4;accode=12; 
elseif  （Z==1 & amplen==2）; codelen=5;accode=27;  
elseif  （Z==1 & amplen==3）; codelen=7;accode=121; 
elseif  （Z==1 & amplen==4）; codelen=9;accode=502;  
elseif  （Z==1 & amplen==5）;codelen=11;accode=2038; 
elseif  （Z==1 & amplen==6）;codelen=16;accode=65412;
elseif  （Z==1 & amplen==7）;codelen=16;accode=65413;
elseif  （Z==1 & amplen==8）;codelen=16;accode=65414;
elseif  （Z==1 & amplen==9）;codelen=16;accode=65415;
elseif （Z==1 & amplen==10）;codelen=16;accode=65416;
    
elseif （Z==2 & amplen==1）;codelen=5;accode=28;
elseif （Z==2 & amplen==2）;codelen=8;accode=249;
elseif （Z==2 & amplen==3）;codelen=10;accode=1015;
elseif （Z==2 & amplen==4）;codelen=12;accode=4084;
elseif （Z==2 & amplen==5）;codelen=16;accode=65417;
elseif （Z==2 & amplen==6）;codelen=16;accode=65418;
elseif （Z==2 & amplen==7）;codelen=16;accode=65419;
elseif （Z==2 & amplen==8）;codelen=16;accode=65420;
elseif （Z==2 & amplen==9）;codelen=16;accode=65421;
elseif （Z==2 & amplen==10）;codelen=16;accode=65422;
    
    
elseif （Z==3 & amplen==1）;codelen=6;accode=58;
elseif （Z==3 & amplen==2）;codelen=9;accode=503;
elseif （Z==3 & amplen==3）;codelen=12;accode=4085;
elseif （Z==3 & amplen==4）;codelen=16;accode=65423;
elseif （Z==3 & amplen==5）;codelen=16;accode=65424;
elseif （Z==3 & amplen==6）;codelen=16;accode=65425;
elseif （Z==3 & amplen==7）;codelen=16;accode=65426;
elseif （Z==3 & amplen==8）;codelen=16;accode=65427;
elseif （Z==3 & amplen==9）;codelen=16;accode=65428;
elseif （Z==3 & amplen==10）;codelen=16;accode=65429;
    
    
elseif （Z==4 & amplen==1）;codelen=6;accode=59;
elseif （Z==4 & amplen==2）;codelen=10;accode=1016;
elseif （Z==4 & amplen==3）;codelen=16;accode=65430;
elseif （Z==4 & amplen==4）;codelen=16;accode=65431;
elseif （Z==4 & amplen==5）;codelen=16;accode=65432;
elseif （Z==4 & amplen==6）;codelen=16;accode=65433;
elseif （Z==4 & amplen==7）;codelen=16;accode=65434;
elseif （Z==4 & amplen==8）;codelen=16;accode=65435; 
elseif （Z==4 & amplen==9）;codelen=16;accod

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     173687  2018-11-01 15:47  my\1.jpg

     文件      25951  2018-11-01 15:59  my\2.jpg

     文件    3011863  2018-11-12 20:01  my\3.jpg

     文件       9088  2018-11-12 19:48  my\ACEncoding.m

     文件         58  2018-11-01 16:02  my\compression.m

     文件       4828  2018-11-12 19:21  my\compression1.m

     文件       1023  2018-11-12 19:48  my\DCEncoding.m

     文件     196662  2011-05-14 10:30  my\lena.bmp

     文件        309  2013-03-21 08:49  my\zigzag.m

     文件       4544  2018-11-12 21:00  my\Zprogram.m

     目录          0  2018-11-12 21:32  my

----------- ---------  ---------- -----  ----

              3428013                    11