基于一维滤波器的图像处理_solution.pdfVIP

基于一维滤波器的图像处理 一个二维（2D ）离散时间信号x(m,n)有两个自变量，一幅图像就是一个2D 离散信号。在过 滤2D 信号中普遍采用的方法是用两个 1D 滤波器的级联来处理信号，每个用于一个方向， 本设计所需图像在plus.mat 中。 1）用butter 和remez 产生3 个离散时间滤波器的系数： Wc=0.4; n1=10; n2=4; n3=12; [b1,a1]=butter(n1,Wc); a2=1; b2=remez(n2, [0 Wc-0.04 Wc+0.04 1], [1 1 0 0]); a3=1; b3=remez(n3, [0 Wc-0.04 Wc+0.04 1], [1 1 0 0]); 2 ）考察上述三个滤波器频率响应的幅值和相位，哪个滤波器具有线性相位？ 3 ）画出每个滤波器的阶跃响应，哪个滤波器具有最大的超调（稳态值与最大值的差）？ 4 ）用image （64*x ）画出图像x ，用这三个滤波器分别对x 进行行、列滤波，给出滤波后的 结果 （用filter 函数），比较哪个输出图像具有更多的振荡？ 5 ）哪个滤波器在原图像的形状上引入更多的失真？要密切注意在输入图像的对称性上所 受到的任何破坏，这种失真是由于滤波器的非线性相位带来的。 6 ）画出带有噪声的图像xn ，重复4 ，对于消除噪声，哪个滤波器的效果最好？讨论图像的 行列为什么相对原图发生了变化。 （2 ） [b1,a1]=butter(n1,Wc); figure(1); freqz(b1,a1); a2=1; b2=remez(n2, [0 Wc-0.04 Wc+0.04 1], [1 1 0 0]); figure(2); freqz(b2,a2); a3=1; b3=remez(n3, [0 Wc-0.04 Wc+0.04 1], [1 1 0 0]); figure(3); freqz(b3,a3); 0 ) B -100 d ( e d u -200 t i n g a M -300 -400 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Normalized Frequency ( rad/sample) 0 ) s e e r g e d -500 ( e s a h P -1000 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Normalized Frequency ( rad/sample) 50 ) B d 0 ( e d u t i n g -50 a M -100 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Normalized Frequency ( rad/sample) 0 ) s e e -100 r g e d ( e s -200 a h P -300 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Normalized Frequency ( rad/sample) 50 ) B d