噪后的图像
sfftshiftfft2I4将灰度图像的二维不连续Fourier变换的零频率成分
移到频谱的中心
MNsizes
分别返回s的行数到M中，列数到N中

1floorM2
对M2进行取整

2floorN2
对N2进行取整
d040
fori1M
forj1N
dsqrti
12j
22点（ij）到傅立叶变换中心的距离
hexplog1sqrt2dd02
sijhsij
ILPF滤波后的频域表示
e
d
e
d
sifftshifts
对s进行反FFT移动
sim2ui
t8realifft2s对s进行二维反离散的Fourier变换后，取
复数的实部转化为无符号8位整数
subplot133
创建图形图像对象
imshows
显示ILPF滤波后的图像
titleELPF滤波后的图像（d40）
运行结果：
f14梯形低通滤波器TLPF
I1imreadeighttif
读取图像
I2im2doubleI1
I3im
oiseI2gaussia
001
I4im
oiseI3saltpepper001
figuresubplot131
imshowI2
显示灰度图像
title原始图像
为图像添加标题
subplot132
imshowI4
加入混合躁声后显示图像
title加噪后的图像
sfftshiftfft2I4将灰度图像的二维不连续Fourier变换的零频率成分
移到频谱的中心
MNsizes
分别返回s的行数到M中，列数到N中

1floorM2
对M2进行取整

2floorN2
对N2进行取整
d010d1160
fori1M
forj1N
dsqrti
12j
22点（ij）到傅立叶变换中心的距离
ifdd0
h1
elseifd0d1
hdd1d0d1
elseh0
e
d
e
d
sijhsije
d
ILPF滤波后的频域表示
fe
d
sifftshifts
对s进行反FFT移动
sim2ui
t8realifft2s对s进行二维反离散的Fourier变换后，取
复数的实部转化为无符号8位整数
subplot133
创建图形图像对象
imshows
显示ILPF滤波后的图像
titleTLPF滤波后的图像
为图像添加标题
运行结果：
15高斯低通滤波器GLPF
I1imreadeighttif
读取图像
I2im2doubleI1
I3im
oiseI2gaussia
001
I4im
oiseI3saltpepper001
figuresubplot131
imshowI2
显示灰度图像
title原始图像
为图像添加标题
subplot132
imshowI4
加入混合躁声后显示图像
title加噪后的图像
sfftshiftfft2I4将灰度图像的二维不连续Fourier变换的零频率成分
移到频谱的中心
MNsizes
分别返回s的行数到M中，列数到N中

1floorM2
对M2进行取整

2floorN2
对N2进行取整
fd040
fori1M
forj1N
dsqrti
12j
22点（ij）到傅立叶变换中心的距离
h1exp12d2d02
GLPF滤波函数
sijhsij
ILPF滤波后的频域表示
e
d
e
d
sifftshifts
对s进行反FFT移动
sim2ui
t8realifft2s对s进行二维反离散的Fourier变换后，取
复数的r