clear all ％把命令行清空，准备开始
N=100; ％设定FFT的点数
f=zeros(50,50); ％设置一个50×50的全零矩阵
f(15:35,23:28)=1; ％把上面那个矩阵的15～35行的23～28列全设为1
figure(1) ％画第一个图
imshow(f,'notruesize') ％画f的图，这个图就是采样的窗口
F=fft2(f,N,N); ％对f矩阵执行FFT变换
F2=fftshift(abs(F)); ％将变换后的矩阵的模值中心转到原先矩阵的中心
figure(2) ％画第二个图
x=1:N;y=1:N; ％设定横纵坐标
mesh(x,y,F2(x,y));colormap(gray);colorbar ％画三维图，灰色，显示颜色柱
这个程序应该是反映图像处理中常用的采样。与一维的采样类似，二维矩形采样的频谱图在xy两个维度的频率上都是采样函数
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024