标准的低通滤波

I=imread('picture.tif');%将要滤波的图像读到I中，当然要换成你的图片名字了
imshow(I);
f=double(I);
g=fft2(f);
g=fftshift(g);
[N1,N2]=size(g);
n=2;
d0=50; %此处d0为截止频率，你可以依次换成60、70、80进行滤波比较
n1=fix(N1/2);
n2=fix(N2/2);
for i=1:N1
for j=1:N2
d=sqrt((i-n1)^2+(j-n2)^2);
h=1/(1+0.414*(d/d0)^(2*n));
result(i,j)=h*g(i,j);
end
end
result=ifftshift(result);
X2=ifft2(result);
J=uint8(real(X2));
figure,imshow(J);