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

实验三用双线性变换法设计IIR数字滤波器实验项目名称：用双线性变换法设计IIR数字滤波器实验项目性质：验证性实验所属课程名称：数字信号处理实验计划学时：2一.实验目的(1)熟悉用双线性变换法设计IIR数字滤波器的原理与方法。(2)掌握数字滤波器的计算机仿真方法。(3)通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。二.实验内容和要求(1)用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。设计指标参数为：在通带内频率低于0.2π时，最大衰减小于1dB；在阻带内[0.3π，π]频率区间上，最小衰减大与15dB。(2)以0.02π为采样间隔，打印出数字滤波器在频率区间[0，π/2]的幅频响应特性曲线。(3)用所设计的滤波器对实际心电图信号采样序列（在本实验后面给出）进行仿真滤波处理，并打印出滤波前后的心电图信号波形图，观察总结滤波作用与效果。三.实验主要仪器设备和材料计算机，MATLAB6.5或以上版本四.实验方法、步骤及结果测试(1)复习有关巴特沃斯模拟滤波器设计和用双线性变换法设计IIR数字滤波器的内容，用双线性变换法设计数字滤波器系统函数H(z)。其中满足本实验要求的数字滤波器系统函数为：根据设计指标，调用MATLAB信号处理工具箱buttord和butter，也可以得到H(z)。由公式（3.1）和（3.2）可见，滤波器H(z)由三个二阶滤波器H1(z)、H2(z)和H3(z)级联而成，如图3-1所示。(2)编写滤波器仿真程序，计算H(z)对心电图采样序列x(n)的响应序列y(n)。设yk(n)为第k级二阶滤波器Hk(z)的输出序列，yk-1(n)为输入序列，如图3-1所示。由（3.2）式可得到差分方程为：当k=1时，yk-1(n)=x(n)。所以H(z)对x(n)的总响应序列y(n)可以用顺序迭代算法得到。即依次对k=1，2，3，求解差分方程（3.3），最后得到y3(n)=y(n)。仿真程序就是实现上述求解差分方程和顺序迭代算法的通用程序。也可以直接调用MATLABfilter函数实现仿真。(3)在通用计算机上运行仿真滤波程序，并调用通用绘图子程序，完成实验内容（2）和（3）。(4)本实验中用到的心电图信号采用序列x(n)人体心电图信号在测量过程往往受到工业高频干扰，所以必须经过低通滤波处理后，才能作为判断心脏功能的有用信息。下面给出一实际心电图信号采样序列样本x(n)，其中存在高频干扰。在实验中，以x(n)作为输入序列，滤除其中的干扰成分。程序运行结果：直接运行程序，结果输出滤波器幅频特性曲线图，有噪声的心电图采集信号波形图和经过三级二阶滤波器滤波后的心电图信号波形图，见图3-2，对比分析就可以看出低通滤波器滤除信号中高频噪声的滤波效果。五.实验报告要求(1)简述实验原理及目的。目的：(1)熟悉用双线性变换法设计IIR数字滤波器的原理与方法。(2)掌握数字滤波器的计算机仿真方法。(3)通过观察对实际心电图信号的滤波作用，获得数字滤波的感性知识。原理：为了克服用脉冲响应不变法产生频谱混叠现象，可以采用非线性频率压缩方法（正切变换），从s平面映射到s1平面，再从s1平面映射到z平面,即实现了双线性变换。(2)由所打印的|H(ejw)|特性曲线及设计过程简述双线性变换法的特点。特点:靠频率的严重非线性关系得到S平面与Z平面的单值一一对应关系,整个jΩ轴单值对应于单位圆一周。双线性变换不存在混迭效应。(3)对比滤波前后的心电图信号波形，说明数字滤波器的滤波过程与滤波作用。答：滤波之后波形更加平滑，加大了取样间隔，排除了误差较大的点。(4)简要回答思考题。六、思考题