(数字信号处理)第5章无限长单位脉冲响应（IIR）.ppt

需要特别强调的是，若模拟滤波器Ha(s)为低通滤波器，应用 变换得到的数字滤波器H(z)也是低通滤波器； 若Ha(s)为高通滤波器，应用 变换得到的数字滤波器H(z)也是高通滤波器; 若为带通、带阻滤波器也是如此。 在IIR数字滤波器的设计中，当强调模仿滤波器的瞬态响应时，采用脉冲响应不变法较好; 而在其余情况下，大多采用双线性变换法。 例 5-4 设计一个一阶数字低通滤波器，3 dB截止频率为ωc=0.25π，将双线性变换应用于模拟巴特沃思滤波器。 解 数字低通滤波器的截止频率为ωc=0.25π，相应的巴特沃思模拟滤波器的 3 dB截止频率是Ωc，就有 模拟滤波器的系统函数为 将双线性变换应用于模拟滤波器，有 由上题可知，T不参与设计，即双线性变换法中用 设计与用 设计得到的结果一致。 例5-5 用双线性变换法设计一个三阶巴特沃思数字低通滤波器，采样频率为fs=4 kHz（即采样周期为T=250 μs），其3dB截止频率为fc=1 kHz。 三阶模拟巴特沃思滤波器为 解 首先，确定数字域截止频率ωc=2πfcT=0.5π。 第二步，根据频率的非线性关系式（5-46），确定预畸变的模拟滤波器的截止频率 第三步，将Ωc代入三阶模拟巴特沃思滤波器Ha(s)，得 最后，将双线性变换关系代入就得到数字滤波器的系统函数 应该注意，这里所采用的模拟滤波器Ha(s)并不是数字滤波器所要模仿的截止频率fc=1 kHz的实际滤波器，它只是一个“样本”函数，是由低通模拟滤波器到数字滤波器的变换中的一个中间变换阶段。 图5-16给出了采用双线性变换法得到的三阶巴特沃思数字低通滤波器的幅频特性。由图可看出，由于频率的非线性变换， 使截止区的衰减越来越快。最后在折叠频率处形成一个三阶传输零点。这个三阶零点正是模拟滤波器在Ωc=∞处的三阶传输零点通过映射形成的。 图 5-16 用双线性变换法设计得到的三阶巴特沃思数字低通滤波器的频响 5.6 设计IIR滤波器的频率变换法 图 5-17 两种等效的设计方法 (a) 先模拟频率变换， 再数字化; (b) 将(a)的两步合成一步设计 对于第一种方案，重点是模拟域频率变换，即如何由模拟低通原型滤波器转换为截止频率不同的模拟低通、高通、带通、带阻滤波器，这里我们不作详细推导，仅在表5-3列出一些模拟到模拟的频率转换关系。一般直接用归一化原型转换，取Ωc=1, 可使设计过程简化。 表5-3 截止频率为Ωc的模拟低通滤波器到其它频率选择性滤波器的转换公式 第二种方法实际上是把第一种方法中的两步合成一步来实现，即把模拟低通原型变换到模拟低通、高通、带通、带阻等滤波器的公式与用双线性变换得到相应数字滤波器的公式合并，就可直接从模拟低通原型通过一定的频率变换关系， 一步完成各种类型数字滤波器的设计，因而简捷便利，得到普遍采用。 此外，对于高通、带阻滤波器，由于脉冲响应不变法不能直接采用，或者只能在加了保护滤波器以后使用，因此，脉冲响应不变法使用直接频率变换要有许多特殊考虑，故对于脉冲响应不变法来说，采用第一种方案有时更方便一些。 我们在下面只考虑双线性变换，实际使用中多数情况也正是这样。 5.6.1 模拟低通滤波器变换成数字低通滤波器 首先，把数字滤波器的性能要求转换为与之相应的作为“样本”的模拟滤波器的性能要求，根据此性能要求设计模拟滤波器， 这可以用查表的办法， 也可以用解析的方法。然后，通过脉冲响应不变法或双线性变换法，将此“样本”模拟低通滤波器数字化为所需的数字滤波器H(z)。例5-5已经说明了用双线性变换法设计低通滤波器的过程，这里再用脉冲响应不变法来讨论一下例5-5的低通滤波器设计问题。 例 5-6 用脉冲响应不变法设计一个三阶巴特沃思数字低通滤波器，采样频率为fs=4 kHz（即采样周期为T=250μs），其3 dB截止频率为fc=1 kHz。 解 查表可得归一化三阶模拟巴特沃思低通滤波器的传递函数 然后, 以s/Ωc代替其归一化频率，则可得三阶模拟巴特沃思低通滤波器的传递函数为 式中，Ωc=2πfc。上式也可由巴特沃思滤波器的幅度平方函数求得。  为了进行脉冲响应不变法变换，将上式进行因式分解并表示成如下的部分分式形式： 将此部分分式系数代入(5-40)式就得到 式中，ωc=ΩcT=2πfcT=0.5π是数字滤波器数