数字信号处理第6章无限脉冲响应数字滤波器的设计_.ppt

第6章 无限脉冲响应数字滤波器的设计 6.1 数字滤波器的基本概念 6.2 模拟滤波器的设计 6.3 用脉冲响应不变法设计IIR数字低通滤波器 6.4 用双线性变换法设计IIR数字低通滤波器 6.5 数字高通、带通和带阻滤波器的设计 6.6 IIR 数字滤波器的直接设计法 6.1 数字滤波器的基本概念 1. 数字滤波器的分类 从实现的网络结构或者从单位脉冲响应分类 无限脉冲响应(IIR)滤波器 有限脉冲响应(FIR)滤波器 通带内和阻带内允许的衰减一般用dB数表示，通带内允许的最大衰减用αp表示，阻带内允许的最小衰减用αs表示，αp和αs分别定义为： 6.2 模拟滤波器的设计 模拟滤波器的理论和设计方法已发展得相当成熟，且有若干典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些滤波器都有严格的设计公式、现成的曲线和图表供设计人员使用。 1.模拟低通滤波器的设计指标及逼近方法 滤波器的技术指标给定后，需要设计一个传输函数Ha(s)，希望其幅度平方函数满足给定的指标αp和αs，一般滤波器的单位冲激响应为实数，因此 要求Ha(s)必须稳定 2 巴特沃斯（ Butterworth ）低通滤波器的设计方法 2 Butterworth模拟低通滤波器设计 将幅度平方函数|Ha(jΩ)|2写成s的函数： 由于各滤波器的幅频特性不同，为使设计统一，将所有的频率归一化。这里采用对3dB截止频率Ωc归一化，归一化后的Ha(s)表示为 式中，s/Ωc=jΩ/Ωc。 令λ=Ω/Ωc，λ称为归一化频率；令p=jλ，p称为归一化复变量，这样归一化巴特沃斯的传输函数为 式中，pk为归一化极点，用下式表示： 将极点表示式(6.2.12)代入(6.2.11)式，得到的Ha(p)的分母是p的N阶多项式，用下式表示： 由(6.2.14)和(6.2.15)式得到： 用上式求出的N可能有小数部分，应取大于等于N的最小整数。关于3dB截止频率Ωc，如果技术指标中没有给出，可以按照(6.2.14)式或(6.2.15)式求出，由(6.2.14)式得到： Butterworth模拟低通滤波器设计 Butterworth模拟低通滤波器设计 常用归一化(wc =1) Butterworth模拟滤波器的系统函数 例6.2.1 已知通带截止频率fp=5kHz，通带最大衰减αp=2dB，阻带截止频率fs=12kHz，阻带最小衰减αs=30dB，按照以上技术指标设计巴特沃斯低通滤波器。 解 (1) 确定阶数N。 (2) 按照(6.2.12)式，其极点为 上式分母可以展开成为五阶多项式，或者将共轭极点放在一起，形成因式分解形式。这里不如直接查表6.2.1简单，由N=5，直接查表得到： 极点：-0.3090±j0.9511,-0.8090±j0.5878; -1.0000 (3) 为将Ha(p)去归一化，先求3dB截止频率Ωc。 按照(6.2.17)式，得到： 阶数N为奇数与偶数时的切比雪夫Ⅰ型滤波器幅频特性。幅度平方函数用A2(Ω)表示： 阶数N=0,4,5时的切比雪夫多项式特性 (1)过零点在|x|≤1的范围内； (2)当|x|1时，|CN(x)|≤1, 具有等波纹性； (3)当|x|1时，CN(x)是双曲线函数，随x单调上升。 按照(6.2.19)式，平方幅度函数与三个参数即ε,Ωp和N有关。其中ε与通带内允许的波动大小有关，定义允许的通带波纹δ用下式表示： 设阻带的起始点频率(阻带截止频率)用Ωs表示，在Ωs处的A2(Ωs)用(6.2.19)式确定： 3dB截止频率用Ωc表示， 以上Ωp,ε和N确定后，可以求出滤波器的极点，并确定Ha(p)，p=s/Ωp。求解的过程请参考有关资料。下面仅介绍一些有用的结果。 ? 设Ha(s)的极点为si=σi+jΩi，可以证明： 上式中仅取正号，得