c巴特沃斯低通滤波器的设计方法课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《c巴特沃斯低通滤波器的设计方法课件.ppt》由用户(三亚风情)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 巴特沃斯低通 滤波器 设计 方法 课件
- 资源描述:
-
1、1第第5章章 数字滤波器的设计数字滤波器的设计 5.1 数字滤波器的基本概念数字滤波器的基本概念5.2 模拟滤波器的设计模拟滤波器的设计5.3 用脉冲响应不变法设计用脉冲响应不变法设计IIR数字低通滤波器数字低通滤波器5.4 用双线性变换法设计用双线性变换法设计IIR数字低通滤波器数字低通滤波器5.5 数字高通、带通和带阻滤波器的设计数字高通、带通和带阻滤波器的设计2数字滤波器:是指输入输出均为数字信号,通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。 高精度、稳定、体积小、重量轻、灵活,不要求阻抗匹配,可实现特殊滤波功能优点:5.1 数字滤波器的基本概念3 1.
2、数字滤波器的分类(1)总体分类: 经典滤波器,信号与干扰分占不同频带:选频滤波器 现代滤波器,信号与干扰频带重叠:维纳滤波器、卡尔曼滤波器、自适应滤波器等(2)滤波特性:低通、高通、带通、带阻等(3)网络结构:IIR、FIR1001( )( )( )1MrrNnrNknkkb zH zH zh n za z4理想低通、高通、带通、带阻数字滤波器幅度特性 )(ejH)(ejH)(ejH)(ejH0低通0高通0带通0带阻222222225 2.数字滤波器的技术指标 我们通常用的数字滤波器一般属于选频滤波器。假设数字滤波器的传输函数H(ej)用下式表示: 数字低通滤波器的技术要求()()()jjjH
3、 eH ee 通带纹波幅度阻带纹波幅度通带截止频率3dB通带截止频率阻带截止频率6 通带内和阻带内允许的衰减一般用dB数表示,通带内允许的最大衰减用ap表示,阻带内允许的最小衰减用as表示,ap和as分别定义为:00()20lgdB()()20lgdB()psjpjjsjH eaH eH eaH e(5.1.3) (5.1.4) 如将|H(ej0)|归一化为1,(5.1.3)和(5.1.4)式则表示成:20lg() dB20lg() dBpsjpjsaH eaH e (5.1.5) (5.1.6) 7 3. 数字滤波器设计方法概述 SpecificationsDesired IIRButter
4、worth, Chebyshev, elliptic, Bessel, etc.脉冲响应不变法阶跃响应不变法双线性变换法 IIR滤波器和FIR滤波器的设计方法是不同,IIR滤波器设计需要借助于模拟滤波器来完成 回想控制器设计问题85.2 模拟滤波器的设计模拟滤波器的设计 模拟滤波器的理论和设计方法已发展得相当成熟,且有若干典型的模拟滤波器供我们选择,如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Elliptic)滤波器、贝塞尔(Bessel)滤波器等,这些滤波器都有严格的设计公式、现成的曲线和图表供设计人员使用。 9 图5.2.1 各种理想模拟滤波器的
5、幅频特性 )(jaH低通带通带阻高通)(jaH)(jaH)(jaH000c10 1.模拟低通滤波器的设计指标及逼近方法模拟低通滤波器的设计指标及逼近方法 模拟低通滤波器的设计指标有ap, p,as和s。 p:通带截止频率s:阻带截止频率 ap:通带最大衰减系数as:阻带最小衰减系数22()10lg()apapHjaHj(5.2.1) (5.2.2) 22()10lg()asasHjaHj 将=0处幅度已归一化到1,即|Ha(0)|=1,得到210lg()papaHj 210lg()sasaHj (5.2.3) (5.2.4) 11 技术指标给定后,设计一个传输函数Ha(s),希望其幅度平方函数
6、满足给定的指标ap和as。 一般滤波器的单位冲激响应为实数,因此2()( )()()()aaasjaaHjHs HsHjHj (5.2.5) 注意:1. 从数学上讲,有无限多种Ha(s)能满足指标2. 实际中,针对结构给定的Ha(s) (如:Butterworth, Chebyshev, Elliptic等),选取合适的参数12 幅度平方函数:221()1()aNcHj(5.2.6) 两个参数:N, c2.巴特沃斯低通滤波器的设计方法巴特沃斯低通滤波器的设计方法13 将|Ha(j)|2写成s的函数: 21( )()1()aaNcHs Hssj(5.2.7) N阶Butterworth滤波器,|
7、Ha(j)|2=Ha(s)Ha(-s)有2N个极点,极点sk为:12121()22( 1)()NkckjNcsje (5.2.8) 三阶巴特沃斯滤波器极点分布|Ha(j)|2的2N个极点均匀分布在半径为c的圆上14 为形成稳定的滤波器,2N个极点中只取s平面左半平面的N个极点构成Ha(s),LHP( )()NcakHsss例:设N=3,极点有6个,它们分别为(5.2.10) 22330121133345jjcccjjcccsessesesse 取s平面左半平面的极点s0, s1, s2组成Ha(s): 32233( )()()()cajjcccHssss15 归一化:由于各滤波器的幅频特性不同
8、,为使设计统一,需要将所有的频率归一化这里采用对3dB截止频率c归一化,归一化后的Ha(s)表示为 令归一化复变量p=s/c,pk=sk/c,得到归一化巴特沃斯的传输函数101( )()aNkkccHsss(5.2.11) 101( )()aNkkHppp(5.2.12) 16 式中,pk为归一化极点,用下式表示: 将极点表示式(5.2.13)代入(5.2.12)式,得到的Ga(p)的分母是p的N阶多项式,用下式表示: 121()22,0,1,1kjkNkcspekN(5.2.13)(5.2.15) 11101( )aNNNGppbpb pb归一化后Ga(p)的N个极点均匀分布在左半单位圆上表
9、6.2.1 巴特沃斯归一化低通滤波器参数 20求参数N: 通带指标: 阻带指标:/10/10101()101psapNas/10/10101lglg101psapasN/1021 ()10papNc N可能有小数部分,应取大于等于N的最小整数/1021 ()10saNsc/10/10101lglg101psapasN(5.2.18)21求参数 c: 3dB截止频率c可由通带指标/阻带指标分别得到: 通带指标10.12(101)paNcp 阻带指标(5.2.19)(5.2.20) 10.12(101)saNcs 22低通巴特沃斯滤波器的设计步骤总结: (2) 根据N求出归一化极点pk,以及归一化
10、传输函数Ga(p):/10/10101lglg101psapasN (4) Ga(p)去归一化:121()22kjNkpe101( )()aNkkGppp(3) 求c:110.10.122(101)or(101)psaaNNcpcs 将p=s/c, pk=sk/c代入Ga(p),得到实际的传输函数Ha(s) (1) 根据技术指标p, ap, s和as,求滤波器的阶数N:23 例5.2.1 已知通带截止频率fp=5kHz,通带最大衰减p=2dB,阻带截止频率fs=12kHz,阻带最小衰减s=30dB,按照以上技术指标设计巴特沃斯低通滤波器。 解: (1) 确定阶数N 0.10.11010.024
11、210122.42lg0.02424.255lg2.4psaaspffN 24 (2) 其极点为3455016523754,jjjjjpepepepepe归一化传输函数为401( )()akkGppp25 上式分母可以展开成为五阶多项式,或者将共轭极点放在一起,形成因式分解形式。这里不如直接查表5.2.1;由N=5,直接查表得到: 极点:-0.3090j0.9511,-0.8090j0.5878; -1.00005432432101( )aHppa pa pa pa pa 其中,a0=1.0000, a1=3.2361, a2=5.2361, a3=5.2361, a4=3.236126 (3
12、) 求3dB截止频率c 按照(5.2.19)式,得到:10.12(101)25.2755/paNcpkrad s (4) Ga(p)去归一化,将p=s/c, pk=sk/c代入Ha(p)中得到:55423324543210( )accccccHssasasasasa27Matlab实现常用函数: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
13、*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,-40,5);28103104-200-1000100200Frequency (rad/s)Phase (degrees)10310410-210-1100Frequency (rad/s)Magnitude02468101214-40-35-30-25-20-15-10-505Frequency(kHz)Magnit
14、ude(dB)293.切比雪夫低通滤波器的设计方法Motivation:Butterworth的频率特性曲线在通带和阻带内都是单调递减若通带(阻带)边界处满足指标,则通带(阻带)内会有较大余量更有效的方法在通带(或阻带、两者)之内等波纹(equiripple),可降低滤波器阶数Chebyshev和Elliptic滤波器都是根据此原则设计的30切比雪夫低通滤波器:I型通带等波纹,阻带单调递减II型通带单调递减,阻带等波纹31 0 1,表示通带内幅度波动的程度,愈大,波动幅度愈大;p为通带截止频率;CN(x)为N阶切比雪夫多项式cos(arccos )1( )cosh(arcosh )1NNxxC
15、xNxx 22221()1()aNpAHjC 切比雪夫I型低通滤波器幅度平方函数:(5.2.24) cosh2xxeex双曲函数:p给定,两个参数和N32当: 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单调上升。33N=0,4,5切比雪夫多项式曲线Observat
16、ions:CN(x)有N个过零点N为偶数时,在x正半轴有N/2个过零点N为奇数时,在x正半轴有(N+1)/2个过零点CN(x)在x正半轴的过零点对应着|Ha(j)|=134切比雪夫型滤波器幅频特性 N=3N=435|Ha(j)|的性质:x = 0(或 = 0)x = 1(或 = p)0 x 1(或0p) |Ha(j)|在1和 之间振荡x 1 (或 p) |Ha(j)|单调递减到02( 0)1for odd1( 0)for even1aaHjNHjN21( 1)for all 1aHjN21136 切比雪夫型与巴特沃斯低通的A2()曲线N相同时,切比雪夫I型性能更好37|Ha(s)|的极点分布:
17、N阶|Ha(s)|2有2N个极点,令分布在左平面的N个稳定极点为|Ha(s)|的极点。通过繁琐的计算可知极点sj为,,0,.,1cos,sin2122iiiipiipiisjiNabiN 极点位于椭圆之上11112,22111NNNNab短轴为ap ,在实轴上长轴为bp,在虚轴上38三阶切比雪夫滤波器的极点分布三个极点,i分别为:0 = 2/31 = 2 = 4/339归一化与去归一化11( )()aNiiGpcpp(5.2.38)c为待定系数,根据幅度平方函数(5.2.24)式可得:c = 2N-1111( )2()aNNiiGppp(5.2.40) 归一化:令p = s/p, pi = s
展开阅读全文