实验二-IIR数字滤波器的设计.docx

双线性 变换法时域 双线性 变换法 时域 逼近法 频域 逼近法 k T -- 实验二 I Ｉ R 数字滤波器的设计 1、 实验目的 （ 1） 掌握脉冲响应不变法和双线性变换法设计 IIR 数字滤波器的具体方法和 原理，熟悉双线性变换法和脉冲响应不变法设计低通、带通 IIR 数字滤波器的计 算机编程； （2） 观察双线性变换法和脉冲响应不变法设计的数字滤波器的频域特性，了 解双线性变换法和脉冲响应不变法的特点和区别； （3） 熟悉 Bu ｔ terｗ orth 滤波器、 Cheby ｓ heｖ滤波器和椭圆滤波器的频率特 性。 ２、实验原理与方法 II Ｒ数字滤波器的设计方法可以概括为如图所示， 本实验主要掌握Ｉ IR 滤波 器的第一种方法 ,即利用模拟滤波器设计 IIR 数字滤波器， 这是 IIR 数字滤波器设 计最常用的方法。 利用模拟滤波器设计 ,需要将模拟域的 Ha(s)转换为数字域Ｈ (z)， 最常用的转换方法为脉冲响应不变法和双线性变换法。 IIR 数字滤波器的设计方法 从模拟滤波器设计 IIR 数字滤波器 脉冲响应 不变法 （ 1）脉冲响应不变法 用数字滤波器的单位脉冲响应序列ｈ 让ｈ (n）正好等于 ha(t)的采样值，即 直线设计 IIR 数字滤波器 零、极点 累试法 (n)模仿模拟滤波器的冲激响应 ha( ｔ )， h(n) ha (nT ) 其中 T 为采样间隔。 如果以 Ha( ｓ )及 H(z) 分别表示 ha( ｔ )的拉氏变换及 h(ｎ) 的 Z 变换 ,则 H ( z) |z esT 1 T H a ( s j 2 k ) 在 MA Ｔ LAB 中 ,可用函数 iｍｐ invar 实现从模拟滤波器到数字滤波器的脉 冲响应不变映射。 -- T T 2 2 -- (2）双线性变换法 S 平面与ｚ平面之间满足下列映射关系 T 1 z 12 1 T 1 z 1 或 z 2 T 2 T Ｓ平面的虚轴单值地映射于 s s z平面的单位圆上 ,s平面的左半平面完全映射到 z平面的单位圆内。双线性变换不存在频率混叠问题。 在Ｍ ATLＡ B 中 ,可用函数 b ｉ li ｎｅａ r实现从模拟滤波器到数字滤波器的双 线性变换映射。 双线性变换是一种非线性变换，即 tan ，这种非线性引起的幅频特 性畸变可通过预畸变得到校正。 (3）设计步骤 Ｉ IR 数字滤波器的设计过程中，模拟滤波器的设计是关键。模拟滤波器的 设计一般是采用分布设计的方式， 这样设计原理非常清楚 ,具体步骤如前文所述。 MATLAB 信号处理工具箱也提供了模拟滤波器设计的完全工具函数 :ｂ ut ｔｅ r、 cheby1、ｃ heｂ y2、 elli ｐ、 besself。用户只需一次调用就可完成模拟滤波器的设 计，这样虽简化了模拟滤波器的设计过程，但设计原理却被屏蔽了。 模拟滤波器设计完成之后 ,利用 imp ｉｎｖａ r 或ｂ il ｉ near 函数将模拟滤波 器映射为数字滤波器 ,即完成了所需数字滤波器的设计。 下图给出了实际低通、 高通、 带通和带阻滤波器的幅频特性和各截止频率的 含义。另外 ,为了描述过渡带的形状，还引入了通带衰减和阻带衰减的概念。 -- | | H (ej p ) | dB -- 通带衰减 : p 图 实际滤波器的幅频特性和各截止频率的含义 20log | H (ej 0 ) | 阻带衰减： s 20log | H (ej 0 ) | | H (ej s ) | dＢ d 在 MA TLAB 信号处理工具箱中 ,通常用 Rp 和Ｒ ｓ 来表示 αp 和 αs。 3、实验内容 (1）参照教材 5．5 节所述滤波器设计步骤 ,利用双线性变换法设计一个 Cheby ｓ hev I 型数字高通滤波器，观察通带损耗和阻带衰减是否满足要求。已知滤波 器的指标为 fｐ=0．３ｋ Hz, αｐ =1.2dB,fs= ０．２ kHz, αs=20dＢ ,T＝１ ms。 (2）已知 fp=0.2kＨ z， αp ＝1 ｄＢ ,fs=0.3 ｋ H ｚ , αｓ=25dＢ， T=1ms，分别用脉 冲响应不变法和双线性变