matlab课后实验六.docxVIP

实 验 六实验项目名称：IIR滤波器的设计与信号滤波实验日期：2015.6.8实验学时：2实验地点：S1-326姓名：王文赛学号： 201201100126实验成绩：指导教师签字：１、实验目的（１）熟悉用双线性变换法设计IIR数字滤波器的原理与方法。（２）掌握数字滤波器的计算机仿真方法。（３）通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。２、实验原理利用双线性变换设计IIR滤波器（只介绍巴特沃斯数字低通滤波器的设计），首先要设计出满足指标要求的模拟滤波器的传递函数，然后由通过双线性变换可得所要设计的IIR滤波器的系统函数。如果给定的指标为数字滤波器的指标，则首先要转换成模拟滤波器的技术指标，这里主要是边界频率的转换，对指标不作变化。边界频率的转换关系为。接着，按照模拟低通滤波器的技术指标根据相应设计公式求出滤波器的阶数和截止频率；根据阶数查巴特沃斯归一化低通滤波器参数表，得到归一化传输函数；最后，将代入去归一，得到实际的模拟滤波器传输函数。之后，通过双线性变换法转换公式，得到所要设计的IIR滤波器的系统函数。利用所设计的数字滤波器对实际的心电图采样信号进行数字滤波器。3、心电图信号采样序列人体心电图信号在测量过程中往往受到工业高频干扰，所以必须经过低通滤波处理后，才能作为判断心脏功能的有用信息。下面给出的数据是一实际心电图信号采样序列样本，其中存在高频干扰。本实验中，以作为输入序列，滤除其中的干扰成分。4、实验步骤及内容（１）复习有关巴特沃斯模拟滤波器的设计和用双线性变换法设计IIR数字滤波器的内容，用双线性变换法设计一个巴特沃斯IIR低通数字滤波器。设计指标参数为：在通带内频率低于时，最大衰减小于；在阻带内频率区间上，最小衰减大于。（２）绘制出数字滤波器在频率区间上的幅频响应特性曲线。（３）用所设计的滤波器对实际心电图信号采样序列（实验数据在后面给出）进行仿真滤波处理，并分别绘制出滤波前后的心电图信号波形图，观察总结滤波作用与效果。（４）编写程序完成各部分实验内容。5、实验程序及运行结果1.双线性变换法设计一个巴特沃斯IIR低通数字滤波器wp=0.2*pi;ws=0.3*pi;ap=1;as=15;T=1;Wp=2/T*tan(wp/2); Ws=2/T*tan(ws/2);[N,Wc]=buttord(Wp,Ws,ap,as,s); [BS,AS]=butter(N,Wc,s);[BZ,AZ]=bilinear(BS,AS,1/T);[H,W]=freqs(BS,AS);plot(W/pi,20*log10(abs(H)));hold on;grid;[H,W]=freqz(BZ,AZ);plot(W/pi,20*log10(abs(H)),r); 2.wp=0.2*pi;ws=0.3*pi;ap=1;as=15;T=1;Wp=2/T*tan(wp/2); %求模拟原型滤波器参数Ws=2/T*tan(ws/2);[N,Wc]=buttord(Wp,Ws,ap,as,s); %求滤波器的阶数[BS,AS]=butter(N,Wc,s);%得到模拟原型滤波器的系统函数[BZ,AZ]=bilinear(BS,AS,1/T);%双线性变换法转换到数字滤波器的系统函数[H,W]=freqs(BS,AS);hold on;grid;[H,W]=freqz(BZ,AZ);axis([0,0.5,-120,50])plot(W/pi,20*log10(abs(H)),r); %画数字滤波器的频响3.用所设计的滤波器对实际心电图信号采样序列（实验数据在后面给出）进行仿真滤波处理T=1;Fs=1/T;wpz=0.2;wsz=0.3;wp=2*tan(wpz*pi/2);ws=2*tan(wsz*pi/2);rp=1;rs=15;[N,wc]=buttord(wp,ws,rp,rs,s);[B,A]=butter(N,wc,s);[Bz,Az]=bilinear(B,A,1/T);[Nd,wdc]=buttord(wp,ws,rp,rs,s);[Bdz,Adz]=butter(Nd,wdc,s);[H,w]=freqz(B,A);n=0:55;x=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66 ,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0];figure(1);stem(n,x,.);y=filter(Bz,Az,x);n1=0:55;figure(2);stem(n1,y,.);（３）由绘制的特性曲线及设计过程简