实验四 用双线性变换法设计IIR数字滤波器.pdfVIP

实验四 用双线性变换法设计 IIR 数字滤波器 一、 实验目的 1、熟悉用双线性变换法设计IIR设计数字滤波器的原理与方法。 2 、掌握数字滤波器的计算机仿真方法。 3、通过观察对实际心电图信号的滤波作用,观察数字滤波的感性知识。 二、 实验内容及步骤 (1) 用双线性变换法设计一个巴特沃斯低通 IIR 数字滤波器。设计指标参数为: 在通带内频率低于0.2π 时,最大衰减小于 1dB；在阻带内[0.3π , π ] 频率区间上, 最小衰减大于15dB。 (2) 打印出数字滤波器在频率区间[ 0, 0.5π ] 上的幅频衰减曲线，和[−4π ,4π ] 上的 幅频响应曲线。 (3) 用所设计的滤波器对实际心电图信号采样序列(在本实验后面给出)进行仿 真滤波处理,并分别打印出滤波前后的心电图波形图,观察总结滤波作用与效果。 3．实验步骤 （1）复习有关巴特沃斯模拟滤波器设计和用双线性变换法设计 IIR 数字滤波器 的内容，用双线性变换法设计满足设计指标的数字滤波器系统函数 H (z ) 。 0.0007378(1 + z − 1 )6 H (Z ) − 1 − 2 − 1 − 2 − 1 −2 (1 − 1.2686z + 0.705z )(1 − 1.0106z + 0.3583z )(1− 0.904z + 0.2155z ) 3 ( ) =∏ H K z (1.1) k 1 A(1 + 2z − 1 + z − 2 ) H k (z ) − 1 − 2 , κ =1,2,3 (1.2) 1− B z − C z k k A 0.09036 B1 1.2686, C 1 0.7051 = =− B 2 1.0106,= C 3 =−0.3583 Β3 =0.9044, C1 =− 0.2155 由（1.1）式和（1.2）式可见，滤波器H (z ) 由三个二阶滤波器H (z ) 、H (z ) 和H (z ) 1 2 3 级联组成，如图1,1 所示。 x(n) y (n) y (n) y3(n)=y(n) 1 2 H (z) H (z) H (z) 1 2 3 图 1.1 （2）打印出数字滤波器在频率区间[ 0, 0.5π ] 上的幅频衰减曲线，和[−4π ,4π ] 上 的幅频响应曲线。主要使用的 MATLAB 函数：freqz() ( 参考电脑桌面上数字信 号处理文件夹内的文件名为 matlab 的 pdf 文档第 4、5 页的内容) 。 下面的 M 文件举例说明如何绘制第一个二阶滤波器的幅频衰减曲线和幅频响应 曲线。 部分函数功能简介： max(x): 向量x 的元素的最大值 eps