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

clear 
a=imread('E:\1.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)); 
但是程序实现不了，请高手改正指点，谢谢～～！