ChIIR数字滤波器的设计方法优秀讲义.pptVIP

第六章 IIR数字滤波器的设计方法 6.1 引言 6.4 用模拟滤波器设计IIR数字滤波器 6.5 冲激响应不变法 6.7 双线性变换法 6.8 常用模拟低通滤波器特性 6.9 设计IIR滤波器的频率变换法 6.10 先用模拟域频带变换法，再利用数字化法设计数字各型滤波器 6.11 先将模拟归一化低通原型数字化为数字低通，再利用数字域频带变换法设计数字各型滤波器 下周一做实验（11月28日，6,7或8,9节） 3、数字滤波器设计方法概述 利用模拟滤波器来设计数字滤波器，就是从已知的模拟滤波器传递函数Ha(s)设计数字滤波器的系统函数H(z)。因此，它归根结底是一个由S平面映射到Z平面的变换，这个变换通常是复变函数的映射变换，为了保证转换后的H(z)稳定且满足技术要求，这个映射变换必须满足以下两条基本要求：  通过对连续函数 等间隔采样得到离散序列 使 ，T为采样间隔。它是一种时域上的转换方法。 转换步骤： 设模拟滤波器的传输函数为Ha(s)，相应的单位冲激响应是ha(t) 对ha(t)进行等间隔采样，采样间隔为T，得到： 上式表明将模拟信号ha(t)的拉氏变换在s平面上沿虚轴按照周期Ωs=2π/T延拓后，再按照 式映射关系，映射到z平面上，就得到H(z)。下面进一步分析这种映射关系。 那么 σ=0, r=1 σ0, r1 σ0, r1 优点： 冲激响应不变法使得数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，也就是时域逼近良好，而且模拟频率Ω和数字频率ω之间呈线性关系ω=ΩT。 因而，一个线性相位的模拟滤波器（例如贝塞尔滤波器）通过冲激响应不变法得到的仍然是一个线性相位的数字滤波器。 缺点： 有频率响应的混叠效应。 所以， 冲激响应不变法只适用于限带的模拟滤波器(例如， 衰减特性很好的低通或带通滤波器)，而且高频衰减越快，混叠效应越小。 如果要对高通和带阻滤波器采用冲激响应不变法， 就必须先对高通和带阻滤波器加一保护滤波器，滤掉高于频率一半以上的频率，然后再使用冲激响应不变法转换为数字滤波器。 但这样会进一步增加设计复杂性和滤波器的阶数。 双线性变换的由来 1、低频附近程线性 2、分段常数的幅频特性要求 模拟滤波器研究较早，理论已经十分成熟，且有若干典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Ellipse)滤波器、贝塞尔(Bessel)滤波器等，这些滤波器都有严格的设计公式、现成的曲线和图表供设计人员使用，利用这些现有技术来解决数字滤波器的设计问题； 采用这种方法时，先要设计一个合适的模拟滤波器，然后将它转换成满足给定指标的数字滤波器； 这种方法适合于设计幅频特性比较规则的滤波器，例如低通、高通、带通、带阻等。 寻找一个恰当的近似函数来逼近理想特性——所谓逼近问题。 最常用的具有优良性能的滤波器： 巴特沃思(Butterworth)滤波器 切比雪夫(Chebyshev)滤波器 椭圆(elliptic)函数或考尔（Cauer）滤波器 实现线性相位的贝塞尔滤波器 (1) 巴特沃思低通滤波器的幅度平方函数|Ha(jΩ)|2用下式表示： 将幅度平方函数|Ha(jΩ)|2写成s的函数： 由于各滤波器的幅频特性不同，为使设计统一，将所有的频率归一化。这里采用对3dB截止频率Ωc归一化，归一化后的Ha(s)表示为 式中，s/Ωc=jΩ/Ωc。 阶数N的大小主要影响幅度特性下降的速度，它应该由技术指标确定。将　　　代入幅度平方函数 中，再将幅度平方函数式代入 得： 用上式求出的N可能有小数部分，应取大于等于N的最小整数。关于3dB截止频率Ωc，如果技术指标中没有给出，可以按照 求出 低通巴特沃思滤波器的设计步骤如下： 1、根据技术指标Ωp,αp,Ωs和αs，用 求出滤波器的阶数N； 2、按照 ，求出归一化极点pk，将pk代入 式，得到归一化传输函数Ha(p)； 3、将Ha(p)去归一化。将p=s/Ωc代入Ha(p)，得到实际的滤波器传输函数Ha(s)