MATLAB 图像傅里叶变换，幅度谱，相位谱

 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('图像相位谱');