实践2：一个自定义二维空域函数的FFT变换
f=zeros(30,30);
f(5:24,13:17)=1;   %创建矩形函数
F=fft2(f,256,256); %进行补零填充
figure;
imshow(log(abs(F)),[-1,5]);
colormap(jet);
colorbar

实践3：计算和显示线性过滤器的频率响应
freqz2
freqz2其实是傅立叶变换在图像处理中的一个简单应用，它可以同时计算和显示线性过滤器的频率响应。在过滤器理论中，线性过滤器脉冲响应的傅立叶变换给出了该过滤器的频率响应。
例如：
h=fspecial('laplacian');   %拉普拉斯过滤器
mesh(freqz2(h));          %低通频率响应
xlabel('Frequency');
ylabel('Frequency');
zlabel('Magnitude');