c巴特沃斯低通滤波器的设计方法实用.pptx

《c巴特沃斯低通滤波器的设计方法实用.pptx》由会员分享，可在线阅读，更多相关《c巴特沃斯低通滤波器的设计方法实用.pptx（71页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1数字滤波器：是指输入输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。高精度、稳定、体积小、重量轻、灵活，不要求阻抗匹配，可实现特殊滤波功能优点：5.1 数字滤波器的基本概念第1页/共71页2 1.数字滤波器的分类(1)总体分类：经典滤波器，信号与干扰分占不同频带：选频滤波器 现代滤波器，信号与干扰频带重叠：维纳滤波器、卡尔曼滤波器、自适应滤波器等(2)滤波特性：低通、高通、带通、带阻等(3)网络结构：IIR、FIR第2页/共71页3理想低通、高通、带通、带阻数字滤波器幅度特性 第3页/共71页4 2.数字滤波器的技术指标 我们通常用的数字滤波器

2、一般属于选频滤波器。假设数字滤波器的传输函数H(ej)用下式表示：数字低通滤波器的技术要求通带纹波幅度阻带纹波幅度通带截止频率3dB通带截止频率阻带截止频率第4页/共71页5 通带内和阻带内允许的衰减一般用dB数表示，通带内允许的最大衰减用ap表示，阻带内允许的最小衰减用as表示，ap和as分别定义为：(5.1.3)(5.1.4)如将|H(ej0)|归一化为1，(5.1.3)和(5.1.4)式则表示成：(5.1.5)(5.1.6)第5页/共71页6 3.数字滤波器设计方法概述 SpecificationsDesired IIRButterworth,Chebyshev,elliptic,Bes

3、sel,etc.脉冲响应不变法阶跃响应不变法双线性变换法 IIR滤波器和FIR滤波器的设计方法是不同，IIR滤波器设计需要借助于模拟滤波器来完成 回想控制器设计问题第6页/共71页75.2 模拟滤波器的设计 模拟滤波器的理论和设计方法已发展得相当成熟，且有若干典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Elliptic)滤波器、贝塞尔(Bessel)滤波器等，这些滤波器都有严格的设计公式、现成的曲线和图表供设计人员使用。第7页/共71页8 图5.2.1 各种理想模拟滤波器的幅频特性 第8页/共71页9 1.模拟低通滤波器的

4、设计指标及逼近方法 模拟低通滤波器的设计指标有ap,p,as和s。p：通带截止频率s：阻带截止频率 ap：通带最大衰减系数as：阻带最小衰减系数(5.2.1)(5.2.2)将=0处幅度已归一化到1，即|Ha(0)|=1,得到(5.2.3)(5.2.4)第9页/共71页10 技术指标给定后，设计一个传输函数Ha(s)，希望其幅度平方函数满足给定的指标ap和as。一般滤波器的单位冲激响应为实数，因此(5.2.5)注意：1.从数学上讲，有无限多种Ha(s)能满足指标2.实际中，针对结构给定的Ha(s)（如：Butterworth,Chebyshev,Elliptic等），选取合适的参数第10页/共7

5、1页11 幅度平方函数：(5.2.6)两个参数：N,c2.巴特沃斯低通滤波器的设计方法巴特沃斯低通滤波器的设计方法第11页/共71页12 将|Ha(j)|2写成s的函数：(5.2.7)N阶Butterworth滤波器，|Ha(j)|2=Ha(s)Ha(-s)有2N个极点，极点sk为：(5.2.8)三阶巴特沃斯滤波器极点分布|Ha(j)|2的2N个极点均匀分布在半径为c的圆上第12页/共71页13 为形成稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成Ha(s)，例：设N=3，极点有6个，它们分别为(5.2.10)取s平面左半平面的极点s0,s1,s2组成Ha(s)：第13页/共71页

6、14 归一化：由于各滤波器的幅频特性不同，为使设计统一，需要将所有的频率归一化这里采用对3dB截止频率c归一化，归一化后的Ha(s)表示为 令归一化复变量p=s/c，pk=sk/c，得到归一化巴特沃斯的传输函数(5.2.11)(5.2.12)第14页/共71页15 式中，pk为归一化极点，用下式表示：将 极 点 表 示 式(5.2.13)代 入(5.2.12)式，得 到 的Ga(p)的分母是p的N阶多项式，用下式表示：(5.2.13)(5.2.15)归一化后Ga(p)的N个极点均匀分布在左半单位圆上第15页/共71页表6.2.1 巴特沃斯归一化低通滤波器参数 第16页/共71页第17页/共71

7、页第18页/共71页19求参数N:通带指标:阻带指标:N可能有小数部分，应取大于等于N的最小整数(5.2.18)第19页/共71页20求参数 c:3dB截止频率c可由通带指标/阻带指标分别得到:通带指标阻带指标(5.2.19)(5.2.20)第20页/共71页21低通巴特沃斯滤波器的设计步骤总结：(2)根据N求出归一化极点pk，以及归一化传输函数Ga(p):(4)Ga(p)去归一化:(3)求c:将p=s/c,pk=sk/c代入Ga(p)，得到实际的传输函数Ha(s)(1)根据技术指标p,ap,s和as，求滤波器的阶数N:第21页/共71页22 例5.2.1 已知通带截止频率fp=5kHz，通带

8、最大衰减p=2dB，阻带截止频率fs=12kHz，阻带最小衰减s=30dB，按照以上技术指标设计巴特沃斯低通滤波器。解：(1)确定阶数N 第22页/共71页23 (2)其极点为归一化传输函数为第23页/共71页24 上式分母可以展开成为五阶多项式，或者将共轭极点放在一起，形成因式分解形式。这里不如直接查表5.2.1；由N=5，直接查表得到：极 点：-0.3090j0.9511,-0.8090j0.5878;-1.0000 其中，a0=1.0000,a1=3.2361,a2=5.2361,a3=5.2361,a4=3.2361第24页/共71页25 (3)求3dB截止频率c 按照(5.2.19)

9、式，得到：(4)Ga(p)去归一化，将p=s/c,pk=sk/c代入Ha(p)中得到：第25页/共71页26Matlab实现常用函数：buttap,buttord,butter%Butterworth_Lowwp=2*pi*5000;ws=2*pi*12000;Ap=2;As=30;N,wc=buttord(wp,ws,Ap,As,s);B,A=butter(N,wc,s);fk=0:14000/512:14000;wk=2*pi*fk;figure(1);freqs(B,A,wk);Hk=freqs(B,A,wk);figure(2);plot(fk/1000,20*log10(abs(Hk

10、);xlabel(Frequency(kHz);ylabel(Magnitude(dB);axis(0,14,-40,5);第26页/共71页27第27页/共71页283.切比雪夫低通滤波器的设计方法Motivation：Butterworth的频率特性曲线在通带和阻带内都是单调递减若通带（阻带）边界处满足指标，则通带（阻带）内会有较大余量更有效的方法在通带（或阻带、两者）之内等波纹(equiripple)，可降低滤波器阶数Chebyshev和Elliptic滤波器都是根据此原则设计的第28页/共71页29切比雪夫低通滤波器：I型通带等波纹，阻带单调递减II型通带单调递减，阻带等波纹第29页/

11、共71页30 0 1，表示通带内幅度波动的程度，愈大，波动幅度愈大；p为通带截止频率；CN(x)为N阶切比雪夫多项式切比雪夫I型低通滤波器幅度平方函数：(5.2.24)双曲函数：p给定，两个参数和N第30页/共71页31当：N=0时，C0(x)=1 N=1时，C1(x)=x N=2时，C2(x)=2x2-1 N=3时，C3(x)=4x3-3x可归纳出高阶切比雪夫多项式的递推公式为 CN+1(x)=2xCN(x)-CN-1(x)(5.2.25)CN(x)特性：(1)切比雪夫多项式的过零点在|x|1的范围内；(2)当|x|1时，|CN(x)|1,在|x|1时，CN(x)是双曲线函数，随x单调上升。

12、第31页/共71页32N=0,4,5切比雪夫多项式曲线Observations:CN(x)有N个过零点N为偶数时，在x正半轴有N/2个过零点N为奇数时，在x正半轴有(N+1)/2个过零点CN(x)在x正半轴的过零点对应着|Ha(j)|=1第32页/共71页33切比雪夫型滤波器幅频特性 N=3N=4第33页/共71页34|Ha(j)|的性质：x=0(或=0)x=1(或=p)0 x 1(或0p)|Ha(j)|在1和 之间振荡x 1(或 p)|Ha(j)|单调递减到0第34页/共71页35 切比雪夫型与巴特沃斯低通的A2()曲线N相同时，切比雪夫I型性能更好第35页/共71页36|Ha(s)|的极点

13、分布：N阶|Ha(s)|2有2N个极点，令分布在左平面的N个稳定极点为|Ha(s)|的极点。通过繁琐的计算可知极点sj为，极点位于椭圆之上短轴为ap，在实轴上长轴为bp，在虚轴上第36页/共71页37三阶切比雪夫滤波器的极点分布三个极点，i分别为:0=2/31=2=4/3第37页/共71页38归一化与去归一化(5.2.38)c为待定系数，根据幅度平方函数(5.2.24)式可得：c=2N-1(5.2.40)归一化：令p=s/p，pi=si/p，去归一化：(5.2.39)第38页/共71页39指标与参数、N之间的关系决定通带内最大衰减ap|Ha(j)|2最大值：|Ha(j)|2最小值：通带内最大衰

14、减ap：N影响过渡带宽度，以及通带内波动的疏密（N等于通带内极值的总个数）(5.2.27)(5.2.26)求得第39页/共71页40和N决定阻带内最大衰减as(5.2.44)结合求得N第40页/共71页41(1)给定技术指标ap,p,as和s切比雪夫型滤波器设计步骤(2)计算参数和N(3)求归一化极点pi及系统Ga(p)(4)去归一化得到实际系统Ha(s)第41页/共71页42 例5.2.2设计低通切比雪夫滤波器，要求通带截止频率fp=3kHz，通 带 最 大 衰 减ap=0.1dB，阻 带 截 止 频 率fs=12kHz，阻带最小衰减as=60dB。解：(1)滤波器的技术要求：(2)求阶数N

15、和：第42页/共71页43 (3)求Ga(p):可求出N=5时的极点pi，得到：(4)将Ga(p)去归一化，得到：第43页/共71页44Matlab实现常用函数：cheb1ap,cheb1ord,cheby1,cheb2ap,cheb2ord,cheby2%Chebyshev1_Lowwp=2*pi*3000;ws=2*pi*12000;Ap=0.1;As=60;N,wp=cheb1ord(wp,ws,Ap,As,s);B,A=cheby1(N,Ap,wp,s);fk=0:14000/512:14000;wk=2*pi*fk;figure(1);freqs(B,A,wk);Hk=freqs(B

16、,A,wk);figure(2);plot(fk/1000,20*log10(abs(Hk);xlabel(Frequency(kHz);ylabel(Magnitude(dB);axis(0,14,-70,5);第44页/共71页45第45页/共71页464.椭圆低通滤波器简介Motivation：Butterworth的频率特性曲线在通带和阻带内都是单调递减Chebyshev-I在通带等波纹，Chebyshev-II在阻带等波纹Elliptic在通带和阻带内都是等波纹几种滤波器比较：Butterworth滤波器是Chebyshev滤波器的特殊形式Chebyshev滤波器是Elliptic滤

17、波器的特殊形式Elliptic滤波器性能最好第46页/共71页47 UN(x)为N阶雅可比椭圆函数。椭圆低通滤波器幅度平方函数：p给定，两个参数和N。的求法与切比雪夫滤波器相似，N的求法如下：第47页/共71页48Observations：椭圆滤波器通带和阻带波纹幅度固定时，阶数越高，过渡带越窄；椭圆滤波器阶数固定时，通带和阻带波纹幅度越小，过渡带就越宽General Conclusions：阶数越高，性能越好，系统越复杂；波纹幅度和过渡带宽度是一对矛盾第48页/共71页49Matlab实现常用函数：ellipap,ellipord,ellip%Elliptic_Lowwp=2*pi*3000

18、;ws=2*pi*12000;Ap=0.1;As=60;N,wp=ellipord(wp,ws,Ap,As,s);B,A=ellip(N,Ap,As,wp,s);fk=0:14000/512:14000;wk=2*pi*fk;figure(1);freqs(B,A,wk);Hk=freqs(B,A,wk);figure(2);plot(fk/1000,20*log10(abs(Hk);xlabel(Frequency(kHz);ylabel(Magnitude(dB);axis(0,14,-80,5);第49页/共71页50第50页/共71页51几种滤波器比较：N,ap,as相同时的过渡带宽度

19、：Butterworth Chebyshev I=Chebyshev II EllipticButterworth：单调下降Chebyshev I/II：通带/阻带等波纹，过渡带和阻带/通带单调下降Elliptic：通带、阻带都等波纹相位逼近情况：Butterworth和Chebyshev在大约3/4通带接近线性相位特性；Elliptic在大约1/2通带接近线性相位特性复杂性：在满足相同幅频响应指标下的阶数N，Butterworth Chebyshev I/II Elliptic第51页/共71页525.频率变换与模拟高通、带通和带阻滤波器设计目标：设计指标为as、ap的模拟高通、带通或带阻滤

20、波器Hd(s)第52页/共71页53步骤：归一化低通滤波器G(p)的截止频率根据as,ap设计G(p)得到Hd(s)Hd(s)的截止频率频率变换公式Butterworth,Chebyshev,Elliptic频率变换公式：p=F(s)第53页/共71页54(1)低通到高通的频率变换公式(5.2.50)复频域：频域：(5.2.51)高通滤波器通带截止频率 归一化低通滤波器G(p)频率高通滤波器Hd(s)频率通带0,1ph,+阻带s,+0,sh高通滤波器阻带截止频率 由对称性可得第54页/共71页55低通与高通滤波器的幅度特性第55页/共71页56例：设计高通滤波器,ph=200Hz,sh=100

21、Hz，幅度特性单调下降，ph处最大衰减为3dB，阻带最小衰减as=15dB。解 （1）高通技术要求：ph=200Hz,ap=3dB;sh=100Hz,as=15dB （2）频率转换：p=1,第56页/共71页57 （3）设计归一化低通G(p)。采用巴特沃斯滤波器，故 （4）求模拟高通Hd(s)：第57页/共71页58Matlab实现常用函数：lp2hp%Butt_LP2HPwp=1;ws=2;Ap=3;As=15;N,wc=buttord(wp,ws,Ap,As,s);B,A=butter(N,wc,s);wph=2*pi*200;BH,AH=lp2hp(B,A,wph);figure(1);

22、freqs(B,A);figure(2);freqs(BH,AH);第58页/共71页59第59页/共71页60其实可以直接来%Butt_HPwp=2*pi*200;ws=2*pi*100;Ap=0.1;As=40;N,wc=buttord(wp,ws,Ap,As,s);B,A=butter(N,wc,high,s);freqs(B,A);第60页/共71页61(2)低通到带通的频率变换公式(5.2.53)复频域：频域：(5.2.54)带通滤波器通带的高、低截止频率 归一化低通滤波器G(p)频率高通滤波器Hd(s)频率通带0,1-1,0pl,00,pu阻带s,+-,-s0,slsu,对数意义下

23、的中心频率第61页/共71页62低通与带通滤波器的幅度特性 第62页/共71页63例：设计模拟带通滤波器，通带带宽Bw=2200rad/s，中心频率0=21000rad/s，通带内最大衰减ap=3dB，阻带sl=2830rad/s,su=21200rad/s，阻带最小衰减as=15dB。解 (1)模拟带通的技术要求：0=21000rad/s,ap=3dB sl=2830rad/s,su=21200rad/s as=15dB,Bw=2200rad/s;(2)频率转换：p=1,第63页/共71页64（3）设计模拟归一化低通滤波器G(p)，采用巴特沃斯型，有（4）求模拟带通Hd(s)：第64页/共7

24、1页65Matlab实现%Butt_BPwp=2*pi*905,1105;%905*1105 approximates 1000*1000ws=2*pi*830,1200;Ap=3;As=15;N,wc=buttord(wp,ws,Ap,As,s);B,A=butter(N,wc,s);fk=0:2000/512:2000;wk=2*pi*fk;Hk=freqs(B,A,wk);plot(fk/1000,20*log10(abs(Hk);xlabel(Frequency(kHz);ylabel(Magnitude(dB);axis(0,2,-20,5);第65页/共71页66(3)低通到带阻的

25、频率变换公式(5.2.58)复频域：频域：(5.2.59)归一化低通滤波器G(p)频率高通滤波器Hd(s)频率通带0,1-,00,pl pu,阻带s,+-,-ssl,0 0,sl第66页/共71页67低通与带阻滤波器的幅频特性第67页/共71页68例：设计模拟带阻滤波器，其技术要求为：pl=2905rad/s,sl=2980rad/s,su=21020rad/s,pu=21105rad/s,ap=3dB,as=25dB。试设计巴特沃斯带阻滤波器。解 （1）模拟带阻滤波器的技术要求：20=plpu=421000025,Bw=su-sl=2200;（2）频率转换：p=1,第68页/共71页69（3

26、）设计归一化低通滤波器G(p)：（4）带阻滤波器的Hd(s)为：第69页/共71页70Matlab实现%Butt_BSwp=2*pi*905,1105;ws=2*pi*980,1020;Ap=3;As=25;N,wc=buttord(wp,ws,Ap,As,s);B,A=butter(N,wc,stop,s);fk=0:2000/512:2000;wk=2*pi*fk;Hk=freqs(B,A,wk);plot(fk/1000,20*log10(abs(Hk);xlabel(Frequency(kHz);ylabel(Magnitude(dB);axis(0,2,-30,5);第70页/共71页71感谢您的欣赏！第71页/共71页