一、编写程序完成不同滤波器的图像频域降噪和边缘增强的算法
并进行比较，得出结论。
1、不同滤波器的频域降噪
11理想低通滤波器ILPF
I1imreadeighttif
读取图像
I2im2doubleI1
I3im
oiseI2gaussia
001
I4im
oiseI3saltpepper001
figuresubplot131
imshowI2
显示灰度图像
title原始图像
为图像添加标题
subplot132
imshowI4
加入混合躁声后显示图像
title加噪后的图像
sfftshiftfft2I4将灰度图像的二维不连续Fourier变换的零频率成分
移到频谱的中心
MNsizes
分别返回s的行数到M中，列数到N中

1floorM2
对M2进行取整

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

2
对
赋初值
f
1floorM2
对M2进行取整

2floorN2
对N2进行取整
d040
初始化d0
fori1M
forj1N
dsqrti
12j
22点（ij）到傅立叶变换中心的距离
h11dd02

BLPF滤波函数
sijhsij
ILPF滤波后的频域表示
e
d
e
d
sifftshifts
对s进行反FFT移动
sim2ui
t8realifft2s对s进行二维反离散的Fourier变换后，取
复数的实部转化为无符号8位整数
subplot133
创建图形图像对象
imshows
显示ILPF滤波后的图像
titleBLPF滤波后的图像（d40）
实验结果：
13指数型低通滤波器ELPF
I1imreadeighttif
读取图像
I2im2doubleI1
I3im
oiseI2gaussia
001
I4im
oiseI3saltpepper001
figuresubplot131
fimshowI2
显示灰度图像
title原始图像
为图像添加标题
subplot132
imshowI4
加入混合躁声后显示图像
title加r