I?=?imread（‘Lena.bmp‘）;
J?=?dct2（I）;?%2維DCT變換
%figure（1）;imshow（log（abs（J））[]）colormap（jet（64））colorbar
[ab]?=?size（J）;???%ab保存矩陣大小
J1?=?J.*（abs（J）>10）;?%將小于10的值令為0
k1?=?idct2（J1）/255;??%DCT反變換重建圖像
Yasuolv1?=?（sum（sum（abs（J）>10）））./（a*b）?%計算壓縮率
J2?=?J.*（abs（J）>100）;?%將小于100的值令為0
k2?=?idct2（J2）/255;