cl;
img=imread('15.bmp');
%img=double(img);
f=fft2(img);        %傅立葉變換
f=fftshift(f);      %使影象對稱
r=real(f);          %影象頻域實部
i=imag(f);          %影象頻域虛部
margin=log(abs(f));      %影象幅度譜，加log便於顯示
phase=log(angle(f)*180/pi);     %影象相位譜
l=log(f);           
subplot(2,2,1),imshow(img),title('源影象');
%subplot(2,2,2),imshow(l,[]),title('
 
影象頻譜');
subplot(2,2,3),imshow(margin,[]),title('影象幅度譜');
subplot(2,2,4),imshow(phase,[]),title('影象相位譜');