巴特沃斯滤波器设计.pdf-资料库

鞍山科技大学电信学院电子教研室第五章 IIR 数字滤波器设计电子教案数字信号处理电子教案第五章 IIR 数字滤波器设计若您在使用过程中遇到什么问题，请通过 EMAIL:askdqq@163.com 或或 EMAIL: askdqq@sina.com 鞍山科技大学电信学院网上答疑留言板留言网址: http://www.asust.edu.cn/dxxy/index.asp 鞍山科技大学电子教研室 2006 年 5 月 10 日 - 1 -

鞍山科技大学电信学院电子教研室第五章 IIR 数字滤波器设计电子教案 5.1 引言数字滤波器是指完成信号滤波处理功能的，用有限精度算法实现的离散时间线性非时变系统，其输入是一组数字量，其输出是经过变换的另一组数字量。因此，它本身即可以是用数字硬件装配成的一台完成给定运算的专用数字计算机，也可以是将所需运算编成程序，让通用计算机来执行。数字滤波器具有稳定性高、精度高、灵活性大等优点。随着数字技术的发展，用数字技术实现滤波器的功能越来越受到人们的注意和广泛的应用。一、常用滤波器的性能指标滤波器性能一般用系统频率特性 ( ωjeH ) 来说明，常用的性能指标主要有以下三个参数： 1. 幅度平方函数 H e ( j ω ) 2 = 该性能指标主要用来说明系统的幅频特性。 2. 相位函数 = = j ω ) j ω ) ⋅ j ω ) H e H e ( *( j ω− H e H e ( ( ) H z H z ( ) ( ) * j ω z e = H e ( j ω ) Re[ = H e ( j ω )] + j Im[ H e ( j ω )] = H e ( j ω ) e j ( β e j ω ) 其中： j ωβ e ( ) = arctg    Im[ Re[ eH ( eH ( j ω )] j ω )]    该指标主要用来说明系统的相位特性。 3. 群延时 ) ωτ ( −= d j )] ω e ( [ β d ω 定义为相位对角频率导数的负值，说明了滤波器对不同的频率成分的平均延时。当要求在通带内的群延迟是常数时，滤波器相位响应特性应该是线性的。二、实际滤波器的频率特性实际设计中所能得到的滤波器的频率特性与理想滤波器的频率特性之间存在着一些显著的差别，现以低通滤波器的频率特性为例进行说明。 1. 理想滤波器的特性：设滤波器输入信号为 )(tx ，信号中混入噪音 )(tu ，它们有不同的频率成分。滤波器的单位脉冲响应为 )(th 。则理想滤波器输出为： y t ( ) = x t [ ( ) + u t ( )] ∗ h t ( ) = K x t τ ) − ( ⋅ （5-1）即噪音信号被滤除 th )( ) Ω = 假定噪音信号被滤除，即 tu )( ∗ Y j ( ) 0 = ，而信号无失真只有延时和线性放大。对(5-1)式作傅里叶变换得： X j ( Ω （5-2） Ω = Ω + H j ( H j ( U j ( X j ( ) Ω ⋅ ) Ω ⋅ ) ) τ− Ω j Ke U j ( ) Ω ⋅ H j ( Ω = （5-3） ) 0 将（5-3）式代入（5-2）式整理得： - 1 -

鞍山科技大学电信学院电子教研室第五章 IIR 数字滤波器设计电子教案 H j ( Ω = ) Y j ( X j ( ) Ω ) Ω = Ke τ− Ω j 假定信号频率成分为： Ω ≤ Ω ，噪音频率成分为 c Ω > Ω 。则完成滤波的理想低通滤波器特性 c 是： H jΩ ( ) h(t) K 0 cΩ Ω 图 5-1 理想低通滤波器频率特性 H j ( Ω = ) 图 5-2 理性滤波器的单位脉冲响应( 0τ= ) t Y j ( X j ( j τ− Ω  ⋅ K e | ) Ω = Ω ) 0 |  Ω ≤ Ω Ω > Ω | | c c （5-4）即 | H j (  Ω =   ) | K | 0 | Ω ≤ Ω c Ω > Ω | | c arg( H j ( Ω = −Ω τ )) 系统的单位脉冲响应为： h t ( ) = 1 2 π Ω c −Ω ∫ − Ω j τ Ke c ⋅ e j t Ω d Ω = K t ) sin[( − Ω τ t ) ( − π τ ] c （5-5）理性低通滤波器的频率特性如图 5-1 所示，单位脉冲响应的波形如图 5-2 所示。 2. 实际滤波器特性理想滤波器具有非因果、无限长的单位脉冲响应和不连续的频率特性，要用稳定的线性时不变（LTI）系统来实现这样的特性是不可能的。工程上是用脉冲响应为有限长的、因果的、稳定的线性时不变系统或具有连续频率特性的线性时不变系统来逼近理想特性。在满足一定的误差要求的情况下来实现理想滤波特性。因此实际的滤波器的频率特性如图 5-3 所示。其中 )jH e ω ( 1 δ+ 1 1 σ− 1 2σ 通带 0 cω sω 过渡带阻带 ω cω ——截止频率 sω ——阻带起始频率 s ωω − ——过渡带宽 c 在通带内幅度响应以 1σ± 的误差接近于 1，即 1 − σ 1 ≤ j H e ω ) ( 1 ≤ + σ ω ω c ≤ 1 图 5-3 实际滤波器的频率特性 sω 为阻带起始频率，在阻带内幅度响应以小于 2σ 的误差接近于零，即 H e ω σ ω ω π ≤ )j ≤ ≤ ( 2 s - 2 -

鞍山科技大学电信学院电子教研室第五章 IIR 数字滤波器设计电子教案为了使逼近理想低通滤波器的方法成为可能，还必须提供一带宽为 s cω ω− 的不为零的过渡带。在这个频带内，幅度响应从通带平滑的下落到阻带。三、设计无限长单位脉冲响应数字滤波器的常用方法常用的 IIR 滤波器设计方法主要有以下几种： 1. 以模拟滤波器函数为基础的变换法；即先设计一满足指定条件的模拟滤波器 H(s)，再将该模拟滤波器转化为数字滤波器 H(z)。 2. 直接设计法：在 z 平面内，根据零、极点对系统特性的影响，调整零极点位置得 H(z)。 3. 最优化设计法(计算机辅助设计)，在某种最小化误差准则下，建立差分方程系数 a k、b i 对理想特性的逼近方程，使用迭代方法解方程组得到最佳逼近系统。由于此方法计算量大，需要借助于计算机进行设计。在本书中着重研究第一种方法，即由模拟滤波器设计数字滤波器的方法。 5.2 常用模拟滤波器设计为了从模拟滤波器设计 IIR 数字滤波器，必须先设计一个满足技术指标的模拟原型滤波器。设计“模拟原型”滤波器有多种方法，如模拟低通逼近有巴特沃斯（Butter worth ）型，切比雪夫（Chebyshev）型或椭圆（Elliptic）型。低通滤波器是最基本的，至于高通、带通、带阻等滤波器可以用频率变换的方法由低通滤波器变换得到。模拟滤波器设计就是将一组规定的设计要求，转化为相应的模拟系统函数 ( ) aH s ，使其逼近某个理想滤波器的特性。一、根据幅度平方函数确定系统函数模拟滤波器幅度响应常采用“幅度平方函数” 2( A Ω 表示。 ) A 2( Ω = ) H j ( a 2 Ω = ) H s H ( ) a ( − s ) a s = Ω j (5-6) 式中 ( ) aH s 是模拟滤波器的系统函数，它是 s 的有理函数。 ( aH jΩ 是其稳态响应，又称为滤波器 ) 的频率响应。 ( aH jΩ 是滤波器的稳态振幅特性。 ) 从模拟滤波器变换为数字滤波器是从 ( ) aH s 开始的，为此必须由已知的 2( A Ω 求得 ( ) aH s 。这 ) 就要将（5-6）式与 s 平面的解释联系起来。设 ( ) aH s 有一临界频率（极点或零点）位于 s s= ，则 ( aH 0 s− 必有一相应的临界 ) 频率落在 s s= − 的位置，即当 ( ) aH s 的临界频率是落在 0 aH a − ± 位置时，则 ( jb s− 相应的临界频率必落在 a +  的 jb ) jΩ 2 2 S 平面 δ 位置。应该指出，纯虚数的临界频率必然是二阶的。在 s 平面上，上述临界频率的特性如图 5-4 所示。所得到的对称形式称为象限对称。图中在 jΩ 轴上零点处所表示的数代表零点的阶次是二阶图 5-4 s 平面的零极点分布的。 - 3 -

鞍山科技大学电信学院电子教研室第五章 IIR 数字滤波器设计电子教案任何实际的滤波器都是稳定的，因此极点必落在 s 平面的左半平面。所以落于 s 左半平面的极点都属于 ( ) aH s ，落于 s 右半平面的极点都属于 ( aH s− 。 ) 零点的分布与滤波器的相位特性有关。如要求最小相位特性，则应选 s 平面左半平面的零点为 aH s 的零点；若对相位有特殊要求，则可以以各种不同的组合来分配左半平面和右半平面的零点。 ( ) 综上所述，可归纳出由 2( A Ω 确定 ( ) aH s 的方法是： ) （1）根据（5-6）式，代入 2 Ω = − 或 2s Ω = − 到 2( A Ω ，得到一个 s 平面的函数； js ) （2）求出第一步中所得 s 函数的所有零极点，将左半平面的极点分配给 ( ) aH s ，右半平面的极 aH 点分配给 ( s− ，如要求最小相位特性，则应选 s 平面左半平面的零点为 ( ) ) aH s 的零点；若对相位没有特殊要求，则可以各种不同的组合来分配左半平面和右半平面的零点。（3）根据具体情况，对比 ( )A Ω 与 ( ) aH s 的低频或高频特性就可以确定出增益常数 k。例 5.1 根据以下幅度平方函数 2( A Ω 确定滤波器的系统函数 ( ) ) aH s 。 2 A ( Ω = ) 解： 16(25 2 + Ω − Ω )(36 (49 2 2 ) + Ω 2 ) （1）根据（5-6）式，代入 2 Ω = − 到 2( A Ω ，得到一个 s 平面的函数； 2s ) H s H ( ) a ( s − = ) 2 A ( Ω ) a = 2 Ω =− 2 s 16(25 2 s − s + )(36 2 2 ) s − (49 2 ) （2）求出 ( ) H s H ( s− 中所有的零极点，将左半平面极点分配给 ( ) aH s ， ) a a s = ± ， s = ± 6 s j = ± （二阶） 7 5 极点为：零点为：选 s = − 6, s = − 及一对虚轴零点 7 s = ± 为的 ( ) aH s 极点和零点，即 5 j H s ( ) a = ( 2 k s 0( s + + s 6)( 25) 7) + （5-7）（3）根据 ( )A Ω 低频特性确定出增益常数。  A 2( Ω = ) H j ( a Ω ) 2 ∴ A 2(0) = H a (0) 2 代入（5-7）式得 2 k 0 = ⇒ = ，即 16 4 k 0 - 4 -

鞍山科技大学电信学院电子教研室第五章 IIR 数字滤波器设计电子教案 H s ( ) a = 2 s 25) 4( + s s 6)( 7) + + ( = s 4 2 + 2 100 + s 13 42 + s 二、巴特沃斯（Butterworth）低通滤波器的设计巴特沃斯低通滤波器的幅度平方函数定义为 H j ( 2 ) | λ = | 1 2 2 C λ 1 + 其中 C 为一常数参数，N 为滤波器阶数，λ为归一化低通截止频率, （5-8） N λ= Ω Ω 。 / p 式中 N 为整数，是滤波器的阶次。巴特沃斯低通滤波器在通带内具有最大平坦的振幅特性，这就是说，N 阶低通滤波器在 0Ω = 处幅度平方函数的前 2N-1 阶导数等于零，在阻带内的逼近是单调变化的。巴特沃斯低通滤波器的振幅特性如图 5-5 所示。滤波器的特性完全由其阶数 N 决定。当 N 增加时，滤波器的特性曲线变得更陡峭，这时虽然由（5-8）式决定了在 Ω = Ω 处的幅度函数总是 p 1 0.707 H jΩ ( ) 2 N=2 N=4 N=10 衰减 3dB，但是它们将在通带的更大范围内接近于 1，在阻带内更迅速的接近于零，因而振幅特性更接近于理想的矩形频率特性。滤波器的振幅特性对参数 N 的依赖关系如图 5-5 所示。设归一化巴特沃斯低通滤波器的归一化频率为λ，归一化传递函数为 ( pΩ Ω 图 5-5 巴特沃斯低通滤波器的振幅特性 )H p ，其中 p jλ= ，则由（5-6）式和（5-8）式得：由于 H j ( ) λ 2 = λ = p j 1 + C 2 1 ( 1)N − 2 N p H s H ( ) a ( s − = ) 2 A ( Ω ) a = Ω=− js 1 s Ω j c 2 N ) 1 ( + （5-9）所以巴特沃斯滤波器属于全极点滤波器。 1、常用设计巴特沃斯低通滤波器指标 pλ ：通带截止频率； pα ：通带衰减，单位：dB； sλ：阻带起始频率； sα ：阻带衰减，单位：dB。 H jΩ ( ) 1 pα sα pΩ sΩ Ω 说明：图 5-6 巴特沃斯低通滤波器指标（1）衰减在这里以分贝（dB）为单位；即 - 5 -

鞍山科技大学电信学院电子教研室第五章 IIR 数字滤波器设计电子教案 α = 10lg 1 H j ) ( λ 2 =  10lg 1 + C 2 2 λ N   （2）当 3dBα= 时 p Ω = Ω 为通常意义上的截止频率。 C （3）在滤波器设计中常选用归一化的频率 λ= Ω Ω ，即 / C λ p = Ω Ω p p = 1, λ s = Ω Ω s p 2、巴特沃斯低通滤波器设计实质根据设计指标要求 pλ ， pα ， sλ， sα 确定归一化巴特沃斯低通滤波器幅度平方函数中的待定系数 C 及滤波器的阶数 N；然后再根据幅度平方函数确定巴特沃斯低通滤波器的传递函数 H(s)。（1）将实际频率 Ω 归一化得 λ p = Ω Ω p p = 1 ， λ s = Ω Ω s p ，再根据已知的 pα ， sα ，幅度平方函 | H j ( 2 ) | λ 1 2 2 C λ （2）求 C 和 N = + 1 N 确定 C 和 N。由 αλ ( ) 10lg = 1 H j ) ( λ 2 = ( 10lg 1 + C 2 2 λ N ) 并带入 pλ ， pα ， sλ， sα 得  10 lg    10 lg ( 1 ( 1 + + C C N = 2 p 2 αλ p 2 2 λ α s s = N ) ) 数即      C C 2 2 λ p 2 2 λ s N N α p 10 α s 10 = = 10 10 − − 1 1 因为 pλ = ，所以 1 α p 1010 2 C = 1 − 由 N 2 λ s = α s 10 10 C 1 − 2 = α s 10 10 α p 10 10 − 1 − 1 两边取对数得：这样可以求出 C 和 N。 N = a lg lg s λ 其中 a = α s 10 10 α p 10 10 − 1 − 1 - 6 -

鞍山科技大学电信学院电子教研室第五章 IIR 数字滤波器设计电子教案注意：当 α = p dB 3 时， 2 C = α p 1010 1 10 − = 0.3 1 1 − = ，即 C=1，此时巴特沃斯滤波器只剩下一个参数 N。（3）确定巴特沃斯滤波器的传递函数 H(p)。由于 H p H p − ( ) ( ) = G j ( ) λ 2 = λ = p j 1 ( + 1 p j 2 N ) = 1 1 ( 1) + − N 2 N p 由 1 ( 1) + − N Np 2 = ，解得极点为： 0 kp = e 2 j 1 k N π+ − , N 2 k 将 p 左半平面的极点赋予 ( )H p 即 N =  1,2, ,2 H p ( ) = ( p − )( p 1 p − 1 p 2 ) 其中 kp = e 2 j 1 k N π+ − , N 2 k N =  1,2, , − p ( p N ) 为了便于设计，工程上已将当 pλ = 时，各阶巴特沃斯低通滤波器系统函数设计成表格供查阅， 1 该表如表 5-1 所示。在表 5-1 中的函数被称为归一化巴特沃斯原型低通滤波器系统函数。表 5-1 归一化巴特沃斯模拟低通滤波器系统函数表阶次归一化系统函数 1 1 +s 1 2 1 2 1 3.4 1 3 + + 2 s + s + 1 3 s + 2 s 2 s + 1 4 s + 3 2.6 s + 2 + 2.6 s + 1 s 1 2 3 4 5 5 s + s 3.2361 4 + s 5.2361 2 + s 3.2361 s 5.2361 1 pλ = ，而实际中截止频率为 pΩ ，所 + 1 （4）去掉归一化影响上面设计中采用归一化的频率即以要进行如下的变量代换： p = j jλ = Ω Ω p = s Ω p 即 H s H p = ( ) = ( ) p s Ω p 综上，归纳出设计巴特沃斯低通滤波器的方法如下： - 7 -

资料库

巴特沃斯滤波器设计.pdf

相关推荐

行业

热门标签

最新资料