实验六-IIR数字滤波器的设计与实现

一、实验目的1、熟悉用双线性变换法设计IIR设计数字滤波器的原理与方法。2、掌握数字滤波器的计算机仿真方法。3、通过观察对实际心电图信号的滤波作用,观察数字滤波的感性知识。二、实验内容及步骤(1)用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。设计指标参数为:在通带内频率低于0.2π时,最大衰减小于1dB；在阻带内[0.3,]ππ频率区间上,最小衰减大于15dB。(2)打印出数字滤波器在频率区间[0,0.5]π上的幅频衰减曲线，和[4,4]ππ−上的幅频响应曲线。(3)用所设计的滤波器对实际心电图信号采样序列(在本实验后面给出)进行仿真滤波处理,并分别打印出滤波前后的心电图波形图,观察总结滤波作用与效果。3．实验步骤（1）复习有关巴特沃斯模拟滤波器设计和用双线性变换法设计IIR数字滤波器的内容，用双线性变换法设计满足设计指标的数字滤波器系统函数。()Hz()1612121310.0007378(1)(11.26860.705)(11.01060.3583)(10.9040.2155)()KkzHZzzzzzHz−−−−−−=+=−+−+−+=∏(1.1)2z−()1212k(12),1BkkAzzHzzCz−−−−++=κ=1,2,3(1.2)−−0.0903611.2686,10.705121.0106,30.358310.2155ABCBCC===−==−Β3=0.9044,=−由（1.1）式和（1.2）式可见，滤波器由三个二阶滤波器、和级联组成，如图1,1所示。()Hz1()Hz2()Hz3()Hzy3(n)=y(n)y1(n)y2(n)x(n)H1(z)H2(z)H3(z)图1.1（2）打印出数字滤波器在频率区间[0,0.5]π上的幅频衰减曲线，和[4,4]ππ−上的幅频响应曲线。主要使用的MATLAB函数：freqz()(参考电脑桌面上数字信号处理文件夹内的文件名为matlab的pdf文档第4、5页的内容)。下面的M文件举例说明如何绘制第一个二阶滤波器的幅频衰减曲线和幅频响应曲线。部分函数功能简介：max(x):向量x的元素的最大值eps：系统的浮点（Floating-point）精确度axis:坐标轴的控制函数，调用格式如下axis([xmin,xmax,ymin,ymax])gridon在图形中绘制坐标网格(3）编写滤波器仿真程序，计算对心电图信号采样序列()Hz()xn的响应序列()yn。设()kyn为第k级二阶滤波器的输出序列，()kHz1()kyn−为输入序列，如图4.1所示。由（1.2）式可得到差分方程：111()()2(1)(2)(1)(2)kkkkkkkkynAynAynAynBynCyn−−−=+−+−+−+−(1.3)当时，。所以对κ=1k-1y(n)=x(n)()Hz()xn的总响应序列()yn可以用顺序选代算法得到。即依次对，求解差分方程(1.3)，最后得到1,2,3k=3()()ynyn=。仿真程序就是实现上述求解差分方程和顺序迭代算法的通用程序。也可以直接调用MATLABfilter函数实现仿真。调用MATLABfilter()函数对实际心电图信号滤波(参考电脑桌面上数字信号处理文件夹内的文件名为matlab的pdf文档第3、4页的内容)。并调用DFT()函数绘制滤波前、后信号的波形图，和幅频特性曲线。下面的M文件举例说明如何绘制调用MATLABfilter()函数，完成第一个二阶滤波器对实际心电图信号滤波。注意：…符号将一条语句拆分在不同行,如果语句在同一行则不要。附录：心电图信号采样序列x(n)人体心电图信号在测量过程中往往受到工业高频干扰,所以必须经过低通滤波处理后,才能作为判断心脏功能的有用信息.下面给出一实际心电图信号采样序列样式本x(n),其中存在高频干扰.在实验中,以x(n)作为输入序列,滤除其中的干扰成分。xn=[-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];请大家注意：本学期数字信号处理实验只写前三次的实验报告。请大家在本学期15周周末在开放实验中心（202.202.43.114）网站上查询、核对自己的实验成绩。