1 cl;
 2 img=imread(\'15.bmp\');
 3 %img=double(img);
 4 f=fft2(img);        %傅里叶变换
 5 f=fftshift(f);      %使图像对称
 6 r=real(f);          %图像频域实部
 7 i=imag(f);          %图像频域虚部
 8 margin=log(abs(f));      %图像幅度谱，加log便于显示
 9 phase=log(angle(f)*180/pi);     %图像相位谱
10 l=log(f);           
11 subplot(2,2,1),imshow(img),title(\'源图像\');
12 %subplot(2,2,2),imshow(l,[]),title(\'图像频谱\');
13 subplot(2,2,3),imshow(margin,[]),title(\'图像幅度谱\');
14 subplot(2,2,4),imshow(phase,[]),title(\'图像相位谱\');