matlab滤波器设计.doc

电子科技大学生命科学与技术学院标 准 实 验 报 告（实验）课程名称 数字信号处理 2012-2013-第2学期电子科技大学教务处制表电 子 科 技 大 学实 验 报 告学生姓名： 学 号： 指导教师：李永杰实验地点：清水河校区 实验时间：2013年 5月 22日一、实验室名称：科二 504 机房 二、实验名称：滤波器传输函数的零点和极点对滤波特性的影响三、实验学时：2学时四、实验原理：五、实验目的：（详细填写） 1. 直观地了解滤波器传输函数的零点和极点（的个数和位置）对滤波特性的影响。2. 利用设计的滤波器进行滤波。 六、实验内容：（详细填写） 给定某个滤波器的传输函数H(z):增加其零点和极点的个数，分析新的滤波器的滤波特性（幅度谱）。1）根据前面得到的增加零点/极点后的滤波器幅度谱特性，设计含噪声和有用信号的复杂信号（比如不同频率的cos或sin信号）；2）利用Matlab函数filtfilt( )对信号进行滤波，分析滤波结果。 七、实验器材（设备、元器件）： 八、实验步骤： 九、实验数据及结果分析：（详细填写）（包括程序、图、结果等） b=conv(0, 0.0798, 1, 0.6566); %a=conv(1, -0.5335, 1, -(0.625+0.3773i); %a=conv(a, 1, -(0.625-0.3773i); % ra = 0.5; %rwwa = 0.1*pi; %(02*pi)Rrea = ra*cos(wa);Rima = ra*sin(wa);a1 = Rrea + i*Rima; a2 = Rrea - i*Rima; rb = 0.2; %rwwb = 1.2*pi; %(02*pi)Rreb = rb*cos(wb);Rimb = rb*sin(wb);b1 = Rreb + i*Rimb; b2 = Rreb - i*Rimb; a=conv(a,1,a1);a=conv(a,1,a2);b=conv(b,1,b1);b=conv(b,1,b2);coef = sum(b)/sum(a); %b = b/coef; figure(1);zplane(b,a); %title();h,w = freqz(b,a);figure(2);plot(w/pi, abs(h);title(); x=0:100;y=sin(0.1*pi*x); %noise=1.2*rand(1,length(x)-0.8; %z=y + noise; %figure(3);title() ();plot(x,y,r,x,z,b);figure(4);Y=fft(y);plot(abs(Y);title();figure(5);Z=fft(z);plot(abs(Z);title(); q=filtfilt(b,a,z);figure(6);title();plot(x,q);十、实验结论：（详细填写） 通过观察原信号的频谱图来确定零点极点的位置。通过改造零点极点的位置来设计滤波器。极点位置越接近原信号的频谱峰值效果越好，零点的位置越接近原信号的谷值效果越好。十一、总结及心得体会：（详细填写） 了解了如何设计滤波器，懂得如何自己设计含噪声的有用信号，利用上面得到的滤波器对信号进行滤波。对于给定的滤波器传输函数H