步骤：
（1）计算原图像f(x,y)的DFT，
（2） 讲频谱的零频点移动到频谱图的中心位置；
（3）计算滤波器函数H(U,V)与F(U,V)的乘积G（U,V）；
（4）讲频谱G（U,V）的零频点移回到频谱图的坐上角。
（5）计算（4）的结果的傅立叶反变换g(x,y)；
（6）取g(x,y)的实部作为最终的滤波后的结果图像。
代码：
大家别激动的啦 代码来啦 哈哈哈 ~~~
clc;
clear all;
close all;
if(ndims==3 && size(I,3)==3)
I=rgb2gray(I);
end
f=fft2(I);
s=fftshift(f);
out=s.*ff
out=ifftshift(out);
out=ifft2(out);
out=abs(out);
out=out/max(out(:));
原理很简单的，当然代码的构思也so easy!!!!!!