应用Matlab对图像信号进行频谱分析及滤波
选取一张彩色图片，建议把像素设置成200200，提取图像的
灰度值，并显示出灰度图像；在图像中增加正弦噪声信号自己设置
几个频率的正弦信号，画出加入噪声信号后的灰度图像；给定滤波
器的性能指标，采用窗函数法或者双线性变换设计数字低通滤波器，
并画出滤波器的频率响应；然后用自己设计的滤波器对含噪声图像信
号进行滤波，画出滤波后图像信号的灰度图像。利用频谱分析对原始
灰度图像、加入噪声信号的灰度图像、和滤波后的灰度图像进行频谱
分析和对比，分析信号的变化。
ximreadD1jpgx1rgb2grayxMNsizex1x2im2doublex1figure1subplot131imshowx2x3zeros1MNfori1M
forj1Nx3Mi1jx2ij
e
de
d
LMNfs40000dt1fs
0L1x7si
2pi15000
dt10si
2pi18000
dt6x4x3x7
fori1Mforj1Nx5ijx4Mi1je
d
fe
dfigure1subplot132imshowx5
wp05piws075pirp3rs50
w
buttordwppiwspirprssBzAzbutter
wppizpkbuttap
BapAapzp2tfzpkBAlp2lpBapAapw
BzAzbili
earBAFhwfreqzBzAzLfsfigure3subplot121plotwabsh
x6zeros1Lfori1M
forj1Nx6Mi1jx5ij
e
de
dx8zeros1Lx8filterBzAzx6
x3kfftx3Lfigure2subplot131plotw2absx3k
x4kfftx4Lfigure2subplot132plotw2absx4k
fx7kfftx7Lfigure3subplot122plotw2absx7k
x8kfftx8Lfigure2subplot133plotw2absx8k
fori1Mforj1Nx9ijx8Mi1je
d
e
dfigure1subplot133imshowx9
fr