频域的图像滤波1.实验目的(1) 学习MATLAB的相关使用方法。(2) 了解如何给图像添加噪声。(3) 学习如何在频域对图形进行滤波处理。2. 实验原理(1) 陷波滤波器指的是一种可以在某一个频率点迅速衰减输入信号，以达到阻 碍此频率信号通过的滤波效果。从通过信号的频率范围的角度讲，陷波滤波器属 于带阻滤波器的一种，只是他的阻带非常狭窄。既然陷波滤波器属于带阻滤波器， 它的阶数必须是二阶(含二阶)以上.(2) 使用快速傅立叶变换(fft)和快速傅立叶逆变换(IFFT)。而fftshift的作用是对fft的输出进行重新排列，将零频分量移到频谱的中心。调用格式：y=fftshift(x)当x为向量时，fftshift(x)直接将x中左右两半交换而产生y。当x为矩阵时，fftshift(x)直接将x中左右、上下进行交换而产生y。3. 实验内容(1) 对原图像添加已知频率的噪声。(2) 对含噪声的图像进行频域陷波滤波。(3) 原图像有用成分一般位于低频部分，可对图像添加纹理噪声，形成噪声图 像；分析纹理的频率，设计陷波滤波器滤除噪声。(4) 实验方法：对沿x轴方向的波纹加性噪声进行陷波滤波。参考程序如下：%1生成波纹噪声图像img 二double(rgb2gray(imread(peppers.bmp)；figure; imshow(img,);sizec =size(img);w =0.4*2*pi;%噪声的数字频率N=2*pi/w;%噪声每一周期的采样点数img_noice =img+20*ones(sizec(1),1)*sin(w*1:sizec(2);figure; imshow(img_noice,);%图像频谱F0=fft2(img);F0=fftshift(F0);figure; imshow(log(abs(F0),);F=fft2(img_noice);F=fftshift(F);figure; imshow(log(abs(F),);%2.设计理想陷波滤波器H=ones(sizec(1),sizec(2);%图像中心点x0 =sizec(1)/2+1;y0 =sizec(2)/2+1;%噪声所处频率点(x,y)x =x0;y =y0-round(sizec(2)/N);H(x,y-3:y+3)=0;H(x,(y0-y)+y0-3:(y0-y)+y0+3)=0;%3.滤波结果I=ifftshift(F.*H);img1 =ifft2(I);figure; imshow(img1,);实验结果：（a）原图（b）波纹噪声图像（c）原图频谱（d）含噪声的图像频谱(e)陷波滤波结果4.思考题(1) 对参考程序给出功能注释。(2) 对实验结果进行分析。(3) 该实验可以应用到哪些实际问题中？