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('圖像相位譜');