%DCT
x=imread（‘frame11.png‘）;?
x=rgb2gray（x）;?
x=double（x）;%變成雙精度
t=dctmtx（8）;%生成DCT變換矩陣
%亮度量化表，乘以一個系數，用來控制壓縮率。系數越大，壓縮率越高，得到的圖像質量越差。
a1=1*[1611101624405161;
???1212141926586055;
???1413162440576956;
???1417222951878062;
???182237566810910377;
???243555648110411392;
???49647887103121120101;
???7292959811210010399];
y?=?blkproc（x?[8?8]?‘P1?*?x?*?P2‘?t?t‘）;?
%量化
y10=?blkproc（y?[8?8]?‘round（x./P1）‘a1）;?
figure（1）;imshow（uint8（y10）[]）;title（‘DCT系數‘）;?
%?figure（1）;imshow（log（abs（y10））[]）;title（‘DCT系數‘）;?
%掩膜矩陣用來壓縮DCT系數，保留的AC系數越多，得到的圖像質量越高。
mask1=[1?1?0?0?0?0?0?0?
1?0?0?0?0?0?0?0?
0?0?0?0?0?0?0?0?
0?0?0?0?0?0?0?0?
0?0?0?0?0?0?0?0?
0?0?0?0?0?0?0?0?
0?0?0?0?0?0?0?0?
0?0?0?0?0?0?0?0];%保留了左上3個?
mask2=[1?1?1?1?0?0?0?0?
1?1?1?0?0?0?0?0?
1?1?0?0?0?0?0?0?
1?0?0?0?0?0?0?0?
0?0?0?0?0?0?0?0?
0?0?0?0?0?0?0?0?
0?0?0?0?0?0?0?0?
0?0?0?0?0?0?0?0];%保留了左上10?
mask3=[1?1?1?1?1?1?1?1?
???????1?1?1?1?1?1?1?0?
???????1?1?1?1?1?1?0?0?
???????1?1?1?1?1?0?0?0?
???????1?1?1?1?0?0?0?0?
???????1?1?1?0?0?0?0?0?
???????1?1?0?0?0?0?0?0?
???????1?0?0?0?0?0?0?0];%保留了左上36?
mask4=[?1?1?1?1?1?1?1?1?
1?1?1?1?1?1?1?1?
1?1?1?1?1?1?1?1?
1?1?1?1?1?1?1?1?
1?1?1?1?1?1?1?1?
1?1?1?1?1?1?1?1?
1?1?1?1?1?1?1?1?
1?1?1?1?1?1?1?1];%保留了全部?
%%?反量化與加掩膜
y1=blkproc（y10[8?8]‘P1.*P2.*x‘mask1a1）;?
y2=blkproc（y10[8?8]‘P1.*P2.*x‘mask2a1）;?
y3=blkproc（y10[8?8]‘P1.*P2.*x‘mask3a1）;?
y4=blkproc（y10[8?8]‘P1.*P2.*x‘mask4a1）;?

y11=uint8（blkproc（y1[8?8]‘P1*x*P2‘t‘t））;?
y22=uint8（blkproc（y2[8?8]‘P1*x*P2‘t‘t））;?
y33=uint8（blkproc（y3[8?8]‘P1*x*P2‘t‘t））;?
y44=uint8（blkproc（y4[8?8]‘P1*x*P2‘t‘t））;?

figure（2）;imshow（uint8（x））title（‘原圖‘）;?
figure（3）;imshow（y11）title（‘還原保留左上3‘）;?
figure（4）;imshow（y22）title（‘還原保留左上10‘）;?
figure（5）;imshow（y33）title（‘還原保留左上36‘）;?
figure（6）;imshow（y44）title（‘還原全部‘）;?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2019-04-22?16:39??DCT\
?????文件????????2055??2019-04-22?16:39??DCT\DCT.m
?????文件??????402536??2007-07-28?03:16??DCT\frame11.png