图像处理 高斯低通滤波 MATLAB

%高斯低通滤波器
RGB = imread（‘C:\Users\Administrator\Desktop\lengleng.jpg‘）;  
I0 = rgb2gray（RGB）; 
subplot（231）imshow（I0）;title（‘原图‘）;
I1 = imnoise（I0‘gaussian‘）;         %对原图像加噪声
subplot（232）imshow（I1）;title（‘加入噪声后‘）
%将灰度图像的二维不连续Fourier 变换的零频率成分移到频谱的中心
s=fftshift（fft2（I1））;
subplot（233）imshow（log（1+abs（s））[]）;title（‘fftshift‘）;
[MN]=size（s）;                       %分别返回s的行数到M中，列数到N中
%GLPF滤波 
d0=10;                               %初始化d0
n1=floor（M/2）;                       %对M/2进行取整
n2=floor（N/2