上机验.pptVIP

} K = LH; while(J=K) { J=J-K; K=(int)(K/2+0.5); } J=J+K; } } 10.5 实验三： 用双线性变换法设计IIR数字滤波器 1. 实验目的 (1) 熟悉用双线性变换法设计IIR数字滤波器的原理与方法。 (2) 掌握数字滤波器的计算机仿真方法。 (3) 通过观察对实际心电图信号的滤波作用， 获得数字滤波的感性知识。 2. 实验内容 (1) 用双线性变换法设计一个巴特沃斯低通IIR数字滤波器。 (2) 以 0.02π为采样间隔， 打印出数字滤波器在频率区间［0, π/2］上的幅频响应特性曲线。 (3) 用所设计的滤波器对实际心电图信号采样序列(在本实验后面给出)进行仿真滤波处理，并分别打印出滤波前后的心电图信号波形图， 观察总结滤波作用与效果。 3. 实验步骤 (1) 复习有关巴特沃斯模拟滤波器设计和用双线性变换法设计IIR数字滤波器的内容， 按照 例 6.4.2， 用双线性变换法设计数字滤波器系统函数H(z)。 例 6.4.2 中已求出满足本实验要求的数字滤波器系统函数： (10.5.1) A=0.090 36 B1=1.2686,C1=-0.7051 B2=1.0106,C2=-0.3583 B3=0.9044,C3=-0.2155 (10.5.2) 由(10.5.1)式和(10.5.2)式可见， 滤波器H(z)由三 个二阶滤波器H1(z), H2(z)和H3(z)级联组成， 如图 10.5.1 所示。 图10.5.1 滤波器H(z)的组成 (2) 编写滤波器仿真程序, 计算H(z)对心电图信号采样序列x(n)的响应序列y(n)。设yk(n)为第k级二阶滤波器Hk(z)的输出序列， y k-1 (n)为输入序列， 如图10.5.1 所示。 由(10.5.2)式可得到差分方程： yk(n)=Ay k-1 (n)+2Ay k-1 (n-1)+Ay k-1 (n-2)+Bkyk(n-1)+Ckyk(n-2) (10.5.3) (3) 在通用计算机上运行仿真滤波程序， 并调用通用绘图子程序， 完成实验内容(2)和(3)。 4. 思考题 用双线性变换法设计数字滤波器过程中， 变换公式 中T的取值， 对设计结果有无影响? 为什么? 5. 实验报告要求 (1) 简述实验目的及原理。 (2) 由所打印的|H(ejω)|特性曲线及设计过程简述双线性变换法的特点。 (3) 对比滤波前后的心电图信号波形， 说明数字滤波器的滤波过程与滤波作用。 (4) 简要回答思考题。 第10章 上 机实验 第10章 上 机实验 10.1 引言 10.2 关于实验用计算机语言 10.3 实验一： 信号、 系统及系统响应 10.4 实验二： 用FFT作谱分析 10.5 实验三： 用双线性变换法设计IIR数字滤波器 10.6 实验四： 用窗函数法设计FIR数字滤波器 10.1 引言 第一、 二章是全书的基础内容。 第三、 四章DFT、 FFT是数字信号处理的重要数学工具， 它有广泛的使用内容。 10.2 关于实验用计算机语言 MATLAB是一种交互式的以矩阵为基本数据结构的系统。 在生成矩阵对象时， 不要求明确的维数说明。 所谓交互式， 是指MATLAB的草稿纸编程环境。 与C语言或FORTRON语言作科学数值计算的程序设计相比较， 利用MATLAB可节省大量的编程时间。 10.3 实验一： 信号、 系统及系统响应 1. 实验目的 (1) 熟悉连续信号经理想采样前后的频谱变化关系， 加深对时