Question

MATLAB高手帮我看看程序

高手帮我看看下面的程序！请逐行解释一下每句的含义。。。不用很专业！能看懂就行，在线等答案，速度快的追分！
clear all
N=100;
f=zeros(50,50);
f(15:35,23:28)=1;
figure(1)
imshow(f,'notruesize')
F=fft2(f,N,N);
F2=fftshift(abs(F));
figure(2)
x=1:N;y=1:N;
mesh(x,y,F2(x,y));colormap(gray);colorbar

Answer 1

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两个维度的频率上都是采样函数