matlab 图像离散余弦变换（源代码）

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;