代码及运行结果

clear ;
close all;
gr= imread('cameraman.tif');

subplot(1,3,1);
imshow(gr);
subplot(1,3,2);
D = dct2(gr);
F = log(abs(D));
imshow(F,[]);
subplot(1,3,3);
D(abs(D)<200) =0;
%会丢掉很多点 ，但是原图还可识别，适合做信息提取
M = idct2(D)/255;
imshow(M);