傅立叶变换－MATLAB

理论

1⃣️物体图像的生成可以直接由矩阵运算生成，也可利用windows下的画图工具，生成一幅黑白图像，并调用命令函数imread()输入图像。
2⃣️输入的图像是一个巨大的二维矩阵，利用matlab函数库中的fft2()命令对该矩阵进行二维离散傅里叶变换，得到图像的频谱。
3⃣️该频谱是一个复数矩阵，然后用取模函数abs()对该复数矩阵取模，得到振幅谱矩阵。
4⃣️再利用函数fftshift()对取模后的矩阵进行频谱位移，这是因为变换后的二维矩阵的直流分量位于图像的周边角，该函数交换矩阵的第1，3象限和第2，4象限，使直流分量移动到频谱中心，从而使FFT频谱可视效果与实际图像相吻合。
5⃣️最后利用imshow()函数将图像显示出来。

clear;
a=imread('D:\*.bmp');
grid on
figure(1)
imshow(a,[])
afft=fft2(a);
aabs=abs(afft);
aabss=fftshift(aabs);
figure(2)
imshow(aabss,[])
colormap(gray);
colorbar
figure(3)
plot(aabss)
colormap(gray);
figure(4)
meshc(aabss)
maxx1=max(max(aabss));

说明

这是MATLAB下的傅立叶变换。

参考书

光学实验与仿真 孙绪宝 北京理工大学出版社