快速傅里叶变换

2.3图像的傅里叶变换

%定义N=100

N = 100

%生成一个50x50的全零矩阵

f=zeros(50,50);

%生成一个长方形

f(15:35,23:28)=1;

%将含长方形的全零矩阵显示在从左到右的第一个位置上，题目为空间与图像

subplot(121),imshow(f),title('空间域图像');

%调用fft2函数实现二维图像的快速傅里叶变换，输入为含f（长方形）的50x50全零矩阵，输出为傅里叶变换结果F

F=fft2(f,N,N);

%是变换后的零频率分量位于中心

F2=fftshift(abs(F));

%将变换后的图像放在从左到右的第二个位置上

subplot(122);

%取x前N个元素，取y前N个元素

x=1:N;y=1:N;

%绘制由x,y,F2构成的曲面,设置一个测控条，题目为傅里叶变换幅值

mesh(x,y,F2(x,y));colorbar,title('傅里叶变换幅值')';