测试信号分析与处理实验指导书.docxVIP

《测试信号分析与处理》 实验指导书 翟任何王明武编写 适用专业：测控专业 陕西理工学院机械工程学院 二零壹零年十月 目 录 实验一：数字滤波器的设计 第2页 实验二：数字振荡器的实现 第14页 实验一：数字滤波器的设计 实验学时：4 实验类型：综合 实验要求：必修 一、实验目的: 1、 熟悉Matlab界回并进行操作。 2、 掌握数字滤波器的计算机仿真方法。 2、 掌握用双线性变换法设计I1R数字滤波器的原理少方法，通过观察对实际心电图 信号的滤波作川，获得数字滤波的感性知识。 3、 掌握用窗函数法设计FIR数字滤波器的原理与方法，了解各种窗函数对滤波特性 的影响。 二、实验内容： （一） IIR数字滤波器设计 （1） 川双线性变换法设计一个巴特沃斯低通IIR数字滤波器，设计指标参数为：在 通带内频率低于0.2pi吋,最大衰减小于ldB;在阻带内［0.3pi, pi］频率区间上，，最小衰减大 于 15dB； （2） 以0.02pi为采样间隔，打印出数字滤波器在频率区间［0,0.5pi］上的幅频响应特 性曲线； （3） 川所设计的滤波器对实际心电图信号采样序列（存本实验后面给出）进行仿真滤 波处理，并分别打印出滤波前后的心电图波形图，观察总结滤波作用与效杲。 （二） FIR数字滤波器设计 （1） 川以种窗函数设计线性相位低通FIR数字滤波器，截止频率K/4rad, N=256。 （2） 绘制相应的幅频特性曲线，观察3dB带宽和20dB带宽以及阻带最小衰减。 （3） 比较四种窗函数对滤波器特性的影响。 三、实验原理： （一）用双线形变换法设计IIR数字低通滤波器 脉冲响应不变法的主要缺点是会产牛频谱混叠现彖，使数字滤波器的频响偏离模拟滤 波器的频响特性，产牛的原因是模拟低通滤波器不是带限丁?折叠频率31 AT,在数字化后产 生了频谱混盏，再通过映射关系2 =幺”，使数字滤波器在0 = 71附近形成频谱混證现象。 为了克服这一缺点，可以采用非线性频率压缩方法，将整个模拟频率轴压缩到土 h/TZ 唤，再用乙=尹转换到Z平而上。设Hgs = jQ,然厉经过非线性频率轴压缩用 朮（$i）,$i= jQi 表示。 用正切变换实现频率压缩 Q = 7tan（i°ir） （2-1） 式中T仍是采样间隔，当Q］从n/T经过0变化到:n/T时，Q则由8经过0变化到+ 实现了 s平甸上整个虚轴完全压缩到si平面上虚轴的土 H/T之间的转换。 这样便有 (2-2)s = -r/iC-QJ) = -l~e 7. (2-2) T 2 1 :Tl + r 再通过转换到Z 再通过 转换到Z平而上, 得到: (2-3) 川窗函数设计FIR数字低通滤波器 如果所希望的滤波器的理想频率响应函数为比(宀,则其对应的单位脉冲响应为 hd(n) =——「H2Wdco (2-4) 9 1T J 〃 窗函数设计法的基本原理是用有限长单位脉冲响应序列/2(〃)逼近hn)。山于 打⑺)往往是无限长序列，且是非因果的，所以用窗函数⑵⑺)将打⑺)截断，并进行加 权处理，得到： h(n) = hd (/2)69(n) (2-5) 加仍作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数H(l(ej(0)为: N-l 工 (2-6) n =0 式屮，N为所选窗函数血⑺)的长度。 窗函数法设计的滤波器性能取决于窗函数⑵⑺)的类型及窗口长度N的取值。设计过 程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。 表2?1各种窗函数的基本参： 窗换数 旁瓣峰值幅度/dB 过渡带宽 阻带最小衰减/dB 矩形窗 -13 4咒/N -12 三角形窗 -25 8弭/N -25 汉宁窗 -31 8h/N -44 哈明窗 -41 8n/N -53 不莱克曼窗 -57 12 ji/N -74 凯塞窗(a二7. 865) -57 10n/N -80 H(严)是否满足要求，要进行验算。一般在/2(7?)尾部加零使长度满足于2的整数 次慕，以便用FFT计算H(e) 如果要观察细节，补零点数增多即可。如果H(严)不 满足耍求，则要重新选择窗函数类型和长度N,再次验算，直至满足要求。 如果要求线性相位特性，则力⑺)还必须满足 h(n) = ±h(N -1-n) (2-4) 根据上式中的正负号和长度N的奇偶性又将线性相位FTR滤波器分成四类。要根据设计的 滤波特性正确选择其中一类。 实验条件： 微型计算机matlab软件 五、实验步羽 五、 实验步羽 1、 复习有关巴特沃斯模拟滤波器设计和用双线性变换法设计1IR数字滤波器的内容。 根据设计指标,调用MATLAB信号处理工具箱函数buttord和butter,可得到H (z)。编 写滤波器仿真程序，计算H(z)对心电图信号采样序列x