%图像高斯平滑滤波处理 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'); 处理前后比照: