1. 先载入一幅灰度图像，如下：
（非灰度图）
2. 利用函数fft2,对其进行快速傅立叶变换, 并利用函数fftshift 将变换后的图像原点移动到频率矩形的中心。
3. 利用abs()函数来得到傅立叶频谱;angle()函数得到相位图；
4. 利用imshow 来可视化图像，观察图像的特点；
一． 结果（图像）：
二． 分析说明：
1. 在载入图片的时候要注意图片要存放到该文件的文件夹中。
2. 下载的图片看似是灰度图，实际上并非是灰度图，所以加多了一句代码“A=rgb2gray(A)”
三． 附件程序
A=imread(‘1.jpg‘); %载入图片
A=rgb2gray(A)
B=fftshift(fft2(A)); % 进行傅立叶变换
subplot(231)
imshow(A);
title(‘原始图像‘);
subplot(232)
imshow(abs(B),[ ]);
title(‘原始频谱图‘);
subplot(233)
imshow(log(abs(B)),[ ]);
title(‘取对数后的频谱图‘);
subplot(234)
imshow(angle(B),[ ]);
title(‘相位图‘);
subplot(235)
imshow(real(B),[ ]);
title(‘实部图‘);
subplot(236)
imshow(imag(B),[ ]);
title(‘虚部图‘);
colormap(jet(64)) %给图片上色
转自：http://blog.163.com/luoyun_dreamer/blog/static/215529070201322612858437/
原文：http://www.cnblogs.com/sweeeper/p/4687550.html