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("影象相位譜");