数字信号处-MATLAB设计IIR滤波器.docVIP

浙江万里学院实验报告 课程名称： 数字信号处理 实验名称： 用双线性变换设计IIR数字滤波器 专业班级： 电子104 姓名： 黄文龙 学号： 2010014139 实验日期： 2012/12/5 实验目的： 熟悉用双线性变换法设计IIR数字滤波器的原理和方法。 掌握数字滤波器的计算机仿真方法。 通过观察对心电信号的滤波作用，获得数字滤波器的感性知识。 实验内容： 人体心电图信号在测量过程中往往受到工业高频干扰所以必须经过低通滤波处理后才能作为判断心脏功能的有用信息下面给出一实际心电图信号采样序列样式本x(n)其中存在高频干扰在实验中以x(n)作为输入序列滤除其中的干扰成分{ x(n) } = { -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}为滤除x(n)的噪声，设计一IIR滤波器，其系统函数H(z)： H(z)=（0.0007378(1+z^(-1))^6）/((1-1.268z^(-1)+0.705z^(-2))*(1-1.0106z^(-1)+0.3583z^(-2))*(1-0.904z^(-1)-0.215z^(-2))) 根据H(z)的结构可知，该滤波器由三级IIR滤波器级联而成， Hk(z)=(A(1+2z^-1+z^-2))/(1=Bkz^-1-Ckz^-2),k=1,2,3 其中A=0.09036 B1=1.2686,C1=-0.7051 B2=1.0106,C2=-0.3583 B3=0.9044,C3=-0.2155 输入序列x(n)依次通过三级滤波器即可得到最终输出。 实验结果： 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]; t=0:55; figure(1) plot(t,x,g); title(心电图信号); x1=[-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]; t=0:55; B1=[.09036,.18072,.09036]; A1=[1,-1.268,0.7051]; y1=filter(B1,A1,x1); subplot(2,2,1); plot(t,x1,g,t,y1,k-); xlabel(n); ylabel(x(n)); title(x(n)的心电脉冲函数); B2=[0.09036,.18072,.09036]; A2=[1,-1.0106,0.3583]; y2=filter(B2,A2,y1); subplot(2,2,2); plot(t,y1,g,t,y2,k-); xlabel(n); ylabel(y(n)); title(通过滤波器H1(Z)后的y1(n)函数); B3=[0.09036,.18072,.09036]; A3=[1,-0.904,0.215]; y3=filter(B3,A3,y2); subplot(2,2,3); plot(t,y2,g,t,y3,k-); xlabel(n); ylabel(y(n)); title(通过滤波器H1(Z),H2(Z)后的y2(n)函数); subplot(2,2,4); plot(t,x1,g,t,y3,k-); xlabel(n); ylabel(y(n)); title(通过滤波器H1(Z),H2(Z),H3(Z)后的y3(n)函数); A=fft(x1,200); B=fft(y3,200); magA=abs(A); magB=abs(B); figure(2) subplot(1,2,1) plot(magA); title(未滤