《实验三用双线性变换法设计IIR数字滤波器》由会员分享，可在线阅读，更多相关《实验三用双线性变换法设计IIR数字滤波器（9页珍藏版）》请在金锄头文库上搜索。
1、实验三题目用双线性变换法设计IIR数字滤波器第16周星期_3第6,7节一.实验目的(1) 熟悉用双线性变换法设计IIR数字滤波器的原理与方法。(2) 掌握数字滤波器的计算机仿真方法。(3) 通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。二、实验内容、方法、设计程序及实验结果(1)复习有关巴特沃斯模拟滤波器设计和用双线性变换法设计IIR数字滤波器的内容，用双线性变换法设计数字滤波器系统函数H0。其中满足本实验要求的数字滤波器系统函数为：H(z)0.00073781+zt力1.2686z，i+0.7051z，2丿(1.0106z，i+0.3583z，2丿10.9044z，i+0.21
2、55z，23.1)=1Hz)kk=1式中：+2z-1+z-2H(z)=A1+2+z,k=1,2,3k1Bz_1Cz-2kk3.2)A0.09036B=1.2686，1B=1.0106，2B=0.9044，3C=0.70511C=-0.35832C=-0.21553根据设计指标，调用MATLAB信号处理工具箱buttord和butter，也可以得到H(z)。由公式(3.1)和(3.2)可见，滤波器H(z)由三个二阶滤波器H0、1H(z)和H(z)级联而成，如图3-1所示。23图3-1滤波器H0的组成此参数下的程序如下：%实验三，用双线性变换法设计IIR数字滤波器x=-4,-2,0,-4,-6,-
3、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,.%控制滤波循环变量%关闭全部绘图窗口%创建绘图窗口%定位子图1%横坐标%画出枝干图%调整坐标%标注横坐标%标注纵坐标0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0;closeall;k=1;figure(1);subplot(2,2,1);n=0:55;stem(n,x,.);axis(0,56,-100,50);xlabel(n);ylabel(x(n);ti
4、tle(心电图信号采集序列x(n);%命名该子图%H1滤波器的分母系数矩阵%H2滤波器的分母系数矩阵%H3滤波器的分母系数矩阵%进行滤波%重新赋值X进行下一次滤波%控制循环变量B=0.09036,2*0.09036,0.09036;%H1滤波器的分子系数矩阵A=1.2686,-0.7051;A1=1.0106,-0.3583;A2=0.9044,-0.2155;while(k=3)y=filter(B,A,x);x=y;k=k+1;ifk=2A=A1;elseA=A2;endendsubplot(2,2,3);%定位子图3stem(n,y,.);axis(0,56,-100,50);xlabe
5、l(n);ylabel(y(n);title(三级滤波后的心电图信号(原坐标)；subplot(2,2,2)stem(n,y,.);axis(0,56,-15,5);xlabel(n);ylabel(y(n);title(调整坐标后的心电图信号)；%求数字滤波器的幅频特性A=0.09036,2*0.09036,0.09036;%滤波器的分子系数矩阵B1=1,-1.2686,0.7051;%H1滤波器的分母系数矩阵B2=1,-1.0106,0.3583;%H2滤波器的分母系数矩阵B3=l,-0.9044,0.2155;Hl,w=freqz(A,Bl,100);H2,w=freqz(A,B2,10
6、0);H3,w=freqz(A,B3,100);H4=H1.*(H2);H=H4.*(H3);db=20*log10(abs(H)+eps);subplot(2,2,4)plot(w/pi,db);axis(0,0.5,-50,10);%H3滤波器的分母系数矩阵%进行滤波器幅频特性分析%点积%显示方格xlabel(w);ylabel(IH(eAjAw)l);gridon;title(滤波器的幅频响应曲线);心电图信号采集序列剧町调整坐标后的心电图信号1-.三级滤液后的心电圈信号i原坐标1滤波器的幅频响应曲线（2）用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。设计指标参数为:在通带内频率低
7、于0.2n时，最大衰减小于ldB；在阻带内0.3n,n频率区间上，最小衰减大与15dB。参数如下截图:N,Wn=BUTrOED(0.2,0.3,1,15)Wn=0.2329比2二BUT7ER(U,n)Warning:Funct.oncallBUTTERinvokesinexactnatchD:VJUAILABrXtoolboKXsignalXsignalVbutter.m.0.00070.00440.01110.OHS0.0111D.0044O.OOOfA=1.0000-3.LS364.6222-3.7r951.3136-D.4S0D0.0544程序设计如下：x二-4,-2,0,-4,-6,-
8、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;k=1;%控制滤波循环变量%关闭全部绘图窗口figure(2);subplot(2,2,1);n=0:55;stem(n,x,.);axis(0,56,-100,50);xlabel(n);ylabel(x(n);%创建绘图窗口%定位子图1%横坐标%画出枝干图%调整坐标%标注横坐标%标注纵坐标title(心电图信号
9、采集序列x(n);%命名该子图B=0.0007,0.0044,0.0111,0.0148,0.0111,0.0044,0.0007;%H1滤波器的分子系数矩阵A=1.0000,-3.183&4.6228,-3.7800,1.8139,-0.4801,0.0545;%H1滤波器的分母系数矩阵y二filter(B,A,x);%定位子图3subplot(2,2,3);stem(n,y,.);axis(0,56,-100,50);xlabel(n);ylabel(y(n);title(滤波后的心电图信号(原坐标)；subplot(2,2,2)stem(n,y,.)；axis(0,56,-100,30)
10、；xlabel(n)；ylabel(y(n)；title(调整坐标后的心电图信号)；%求数字滤波器的幅频特性A=0.0007,0.0044,0.0111,0.0148,0.0111,0.0044,0.0007;%滤波器的分子系数矩阵B=1.0000,-3.183&4.6228,-3.7800,1.8139,-0.4801,0.0545;%H1滤波器的分母系数矩阵H,w=freqz(A,B,100);%进行滤波器幅频特性分析db=20*log10(abs(H)+eps);subplot(2,2,4)plot(w/pi,db);axis(0,0.5,-50,10);xlabel(w);ylabel
11、(|H(ejw)|);gridon;%显示方格title(滤波器的幅频响应曲线)；调整坐标后的心电图信号滤波后的心电图信号原坐标)n(1)、双线性变换法的特点：1. 模拟滤波器经过双线性变换后，不存在频率特性的混叠失真，因而对模拟滤波器的频率响应函数Ha(s)无限带要求，而且能够直接用于设计低通、高通、带通、带阻等各种类型的数字滤波器。2. 与冲激响应不变法中模拟频率与数字频率之间的线性关系s=QTs不同的是，双线性变换法中模拟滤波器的频率与所转换成的数字滤波器的频率之间是非线性关系，但是，如果事先进行频率预畸变，这种非线性关系不会使所设计的数字滤波器的幅频特性受到影响。由于频率的非线性关系会
12、产生相频特性失真，所以若对数字滤波器的相位特性要求较严，则不宜采用这种变换方法。3. 双线性变换方法比较容易，不需要将模拟系统函数进行部分分式分解。、数字滤波的过程：顺序迭代.设yk(n)为第k级二阶滤波器Hk(z)的输出序列，yk(n)为输入序列。由式31kkk-1可得到差分方程：y(n)=Ay(n)+2Ay(n-1)+Ay(n-2)+B|y(n-1)+Cy(n-2)当k=1时，yk-1(n)=x(n)o所以H(z)对x(n)的总响应序列y(n)可以用顺序迭代算法得到。即依次对k=1，2,3,求解差分方程(3.3)，最后得到y3(n)=y(n)。数字滤波的作用：通过对有噪声的心电图采集信号波形图和经过三级二阶滤波器滤波后的心电图信号波形图的对比分析，可以看出低通滤波器滤除信号中高频噪声的滤波效果，信号变得更平缓。四、思考题答：T取值对结果没有影响。设数字滤波器W处衰减为51，w2处衰减为62双线性变换公式：Q=ctan由于是低通滤波器:得:模拟巴特沃思低通滤波器的幅度平方函数为:令:a,N是一个与T无关的量。查表得归一化模拟低通滤波器为:则去归一化模拟低通滤波器为:数字系统函数为:由以上各式得:a是一个与T无关的量。H(z)是一个与T无关的函数，所以T取值对结果没有影响。
