%图像高斯平滑滤波处理

img=imread('1.jpg');

f=rgb2gray(img);

subplot(1,2,1);

imshow(f);

f=double(f);

f=fft2(f);

f=fftshift(f);

[m,n]=size(f);  %

d0=80;

m1=fix(m/2);

n1=fix(n/2);

for i=1:m

    for j=1:n

        d=sqrt((i-m1)^2+(j-n1)^2);

        h(i,j)=exp(-d^2/2/d0^2);

    end

end

g=f.*h;

g=ifftshift(g);

g=ifft2(g);

g=mat2gray(real(g));

subplot(1,2,2);

imshow(g);

imwrite(img,'2.jpg');

处理前后比照：

http://home.eeworld.com.cn/attachment/2009/02/99525_200902041924011VZVF.jpg

***************************************************************************************************