无限长单位冲激响应数字滤波器设计方法

1、信号处理方法及应用IIR数字滤波器的设计方法IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter721 数字滤波器的设计步骤按照任务要求确定滤波器的性能按照任务要求确定滤波器的性能指标。指标。用一个因果稳定的系统函数去逼用一个因果稳定的系统函数去逼近这一性能要求。近这一性能要求。利用有限精度算法实现这一系统利用有限精度算法实现这一系统函数。函数。IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter732 数字滤波器的性能指标低通滤波器幅度响应的容限图低通滤波器幅度响应的容限图通带允许的最大衰减及阻带应达到的通带允许的最大衰减及阻带应达到的最小衰减最小衰减 |,| )

2、(|, 1| )(|121stjpjeHeHmax1max2|()|20lg20lg|()|20lg(1)()|()|()|20lg20lg|()|20lg()|()|ppststjjpjjjsjH eRH edBH eH eAH edBH e IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter74通带波纹通带波纹阻带波纹阻带波纹理想低通滤波器逼近的误差容限理想低通滤波器逼近的误差容限1 2 IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter75表征数字滤波器频率响应特性的三个表征数字滤波器频率响应特性的三个参量参量l幅度平方响应幅度平方响应l相位响应相位响应l

3、群延迟响应群延迟响应2()jH e 2*1()()()( )()jjjjz eH eH eHeH z H z *11()1( )()lnln2()2()jjjjz eH eH zejHejH z IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter763 IIR数字滤波器的设计方法分类IIR DF设计方法设计方法间接法间接法直接法直接法冲激响应不变法冲激响应不变法阶跃响应不变法阶跃响应不变法双线性变换法双线性变换法频域最小均方误差逼近法频域最小均方误差逼近法及最小及最小P误差法误差法零极点配置累试法零极点配置累试法时域最小均方误差逼近法时域最小均方误差逼近法IIRIIR数字滤波器

4、的设计方法数字滤波器的设计方法Chapter774 模拟滤波器设计模拟低通滤波器的设计模拟低通滤波器的设计步骤l给定滤波器技术指标给定滤波器技术指标l选定滤波器原型选定滤波器原型l计算滤波器所需阶数计算滤波器所需阶数l通过查表或计算来确定归一化低通滤波通过查表或计算来确定归一化低通滤波器的系统函数器的系统函数l将将 转换为所需类型的低通滤波器系转换为所需类型的低通滤波器系统函数统函数( )anHs( )aHs( )anHsIIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter78模拟巴特沃思低通滤波器模拟巴特沃思低通滤波器l幅度平方函数幅度平方函数l巴特沃斯滤波器的特性巴特沃斯滤波

5、器的特性-3dB截止频率截止频率c通带内，频响具有最大平坦的幅度特性；通带内，频响具有最大平坦的幅度特性；N越大，通带内幅度特性越平坦越大，通带内幅度特性越平坦在过渡带及阻带中，在过渡带及阻带中，频响单调减小；频响单调减小；N越大，衰减得越快越大，衰减得越快 2211/aNcHj 1()0.7072acHj IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter79是全极点滤波器是全极点滤波器极点在极点在s平面是象限对称的平面是象限对称的l 的极点为的极点为l 的系统函数的系统函数( )aHs 21122,1,2,kNjkcsekN ( )aHs 211( )1()NcaNNckk

6、Hss jss IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter710l归一化低通滤波器的系统函数归一化低通滤波器的系统函数l实际的低通滤波器传输函数实际的低通滤波器传输函数Ha(s)是由归一化是由归一化原型低通滤波器去归一化后的结果。原型低通滤波器去归一化后的结果。 211( )1()NcaNNckkHss jss 0212( )1anNNdHsa sa sa s 当直流增益为当直流增益为1时，取时，取d0=a0=1( )( )()sccsaanansHsHsH IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter711l巴特沃思低通滤波器的设计参数的确定巴特

7、沃思低通滤波器的设计参数的确定通过指定通带截止频率、通带允许的最大通过指定通带截止频率、通带允许的最大衰减、阻带截止频率和阻带应达到的最衰减、阻带截止频率和阻带应达到的最小衰减来确定阶数和截止频率。小衰减来确定阶数和截止频率。20lg|()|20lg|()|ppsstRH jAH j 21()1aNcHj IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter712确定滤波器的阶数确定滤波器的阶数根据阶数，查表得到滤波器的归一化低通原根据阶数，查表得到滤波器的归一化低通原型型求截止频率，去归一化后得到实际滤波器求截止频率，去归一化后得到实际滤波器10.1210.12(101)(10

8、1)psRNccpstANccsst 0.10.1lg/ lg101/101spsARsstpNgg IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter713切比雪夫切比雪夫I型低通逼近型低通逼近l特点：特点：1、通带内具有等波纹；、通带内具有等波纹；2、阻带内单调下降、阻带内单调下降IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter714l幅度平方函数：幅度平方函数：通带波纹参数，小于：通带波纹参数，小于1的正数，其值越大，波的正数，其值越大，波纹也越大纹也越大归一化低通滤波器的系统函数归一化低通滤波器的系统函数2221()1(/)aNpHjC 1101cos

9、(cos),1,()ch(ch),1. ()2()(),1, ()1().NNNNNarcCNarcCCCNCC ，且且切比雪夫多项式切比雪夫多项式：迭代关系迭代关系 1201212( )NanNNHsaa sa sa s IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter715l设计参数：设计参数：求确定阶数11( )CPcharcchN 220lg|()| 10lg(1)ppRH j p/10 1210R 222221()1(/)111chch(/)asNspspHjCANarc 0.12(1 /)( 101 /)(/)(/)sAstpstparcchAarcchNarcc

10、harcchIIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter716切比雪夫切比雪夫II型低通逼近型低通逼近l特点：特点：1、通带内单调；、通带内单调；2、阻带内等波纹、阻带内等波纹IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter717椭圆函数低通滤波器（考尔滤波器）椭圆函数低通滤波器（考尔滤波器）简介简介l通带和阻带都是等波纹的通带和阻带都是等波纹的贝塞尔低通滤波器简介贝塞尔低通滤波器简介l在某一频带范围具有近似的线性相位在某一频带范围具有近似的线性相位l瞬态响应过冲小瞬态响应过冲小IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter718五类

11、模拟滤波器比较五类模拟滤波器比较l从幅度特性比较从幅度特性比较l从过渡带宽比较从过渡带宽比较l从阶次比较从阶次比较l从对参数量化的灵敏度比较从对参数量化的灵敏度比较l相位特性比较相位特性比较IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter719例例.给定模拟低通滤波器的性能指标为：通给定模拟低通滤波器的性能指标为：通带截止频率带截止频率内，通带最大衰内，通带最大衰减减 ；阻带截止频率；阻带截止频率阻带最小衰减阻带最小衰减 ，请分别用巴特沃思，请分别用巴特沃思和切比雪夫和切比雪夫I型逼近该滤波器。型逼近该滤波器。解：解：3000pfHz p2dBR 6000stfHz 30dB

12、sA 23000rad /sp 26000rad /sst IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter720巴特沃思逼近巴特沃思逼近l确定确定N取取N=6l查表得到归一化的传递函数查表得到归一化的传递函数l求求3dB截止频率截止频率l去归一化去归一化0.10.1lg(101)/(101)5.3692lg(/)psRAstpN 18637. 34641. 71416. 94641. 78637. 31)(23456 sssssssHansradNRpcp/1 .19711)110(211 . 0 ( )( )()sccsaanansHsHsH IIRIIR数字滤波器的设计

13、方法数字滤波器的设计方法切比雪夫切比雪夫I型逼近型逼近l求通带波纹参数求通带波纹参数l求阶数求阶数取取N=4l查表得到归一化的传递函数查表得到归一化的传递函数l去归一化去归一化Chapter721p/10 1100.76478R 0.1( 101 /)3.35(/)sAstparcchNarcch 143212( )0.7161.25650.51680.2058NanHsssss 164438212162.06336 10( )1.35 104.46 103.46 102.6 10H sssss IIRIIR数字滤波器的设计方法数字滤波器的设计方法用模拟频带变换法设计各种模拟滤波用模拟频带变换

14、法设计各种模拟滤波器器l变换函数是有理函数变换函数是有理函数l映射前后稳定性不变映射前后稳定性不变Chapter722归一化归一化 模拟低通滤波器模拟低通滤波器 1 p)(sHan)( jHan各类模拟各类模拟滤波器滤波器)(sHa)( jHa)(sGs )( jGjIIRIIR数字滤波器的设计方法数字滤波器的设计方法低通到高通的映射低通到高通的映射二者互为倒数关系二者互为倒数关系Chapter723 00:00: 11 pppp ppss,IIRIIR数字滤波器的设计方法数字滤波器的设计方法低通到带通的映射低通到带通的映射变换关系变换关系Chapter724 0 1p 1p ,0 00,pp

15、 0 , 21,pp 12,pp 22201221()ppppppsssB ss 22201221()ppppppB IIRIIR数字滤波器的设计方法数字滤波器的设计方法低通到带阻的映射低通到带阻的映射变换关系变换关系参数参数 可用可用来求解来求解 0st 0st 0 0, 21222012()stsststststststB sssss 21222012()stsststststststB st 1p Chapter725IIRIIR数字滤波器的设计方法数字滤波器的设计方法例设计一个模拟带通滤波器，分别采用巴特例设计一个模拟带通滤波器，分别采用巴特沃斯及切比雪夫沃斯及切比雪夫I型。给定指标为

16、型。给定指标为解：解：1）角频率指标为）角频率指标为1212200,300,2100,400,20pppststsfHz fHz RdBfHz fHz AdB Chapter72612122200,2300,22100,2400,20pppststsRdBAdBIIRIIR数字滤波器的设计方法数字滤波器的设计方法2）归一化低通的相关指标）归一化低通的相关指标3）确定滤波器的阶数）确定滤波器的阶数巴特沃斯：巴特沃斯：切比雪夫切比雪夫I：21012222220102121122100,2244.952.5,5min(,)2.5ppppppstpstpstststpstpstststBBB 0.10

17、.1lg(101)/(101)2.82lg(/)psRAstpN0.10.1( 101 /101)2.07(/)psRAstparcchNarcch 27Chapter7IIRIIR数字滤波器的设计方法数字滤波器的设计方法4）查归一化的低通原型）查归一化的低通原型巴特沃斯：巴特沃斯：切比雪夫切比雪夫I：5）计算归一化巴特沃斯）计算归一化巴特沃斯3dB截止频率截止频率321( )221anHssss 2 1212( )0.80381640.8230603anHsss 0.7647831 10.12(101)1.0935pRNcp Chapter728IIRIIR数字滤波器的设计方法数字滤波器的设

18、计方法6)去归一化去归一化巴特沃斯：巴特沃斯：切比雪夫切比雪夫I：7）转换为带通滤波器）转换为带通滤波器( )()csaanHsH ( )()( )psaananHsHHs 220( )( )ppsbpassBHsHs Chapter729IIRIIR数字滤波器的设计方法数字滤波器的设计方法Chapter730使用使用MATLAB进行模拟滤波器设计进行模拟滤波器设计l巴特沃斯滤波器巴特沃斯滤波器N,Wn=buttord(Wp,Ws,Rp,Rs,s,)z,p,k=buttap(N)num,den=butter(N,Wn,type,s)l切比雪夫切比雪夫I型滤波器型滤波器N,Wn=cheb1ord

19、(Wp,Ws,Rp,Rs,s,)num,den=cheby1(N,Rp,Wn,type,s)z,p,k=cheb1ap(N,Rp)IIRIIR数字滤波器的设计方法数字滤波器的设计方法6 间接法的IIR数字滤波器设计方案对转换关系的要求对转换关系的要求l因果稳定性不变。因果稳定性不变。l频响模仿。频响模仿。模拟频率模拟频率变换变换模拟归一模拟归一化低通原型化低通原型模拟低通、高通模拟低通、高通带通、带阻带通、带阻数字低通、高通数字低通、高通带通、带阻带通、带阻数字化数字化模拟归一模拟归一化低通原型化低通原型数字低通、高通数字低通、高通带通、带阻带通、带阻模拟归一模拟归一化低通原型化低通原型数字低

20、通数字低通数字低通、高通数字低通、高通带通、带阻带通、带阻数字化数字化数字频率数字频率变换变换Chapter731IIRIIR数字滤波器的设计方法数字滤波器的设计方法7 冲激响应不变法变换思路变换思路l使数字滤波器的单位脉冲响应模仿模拟使数字滤波器的单位脉冲响应模仿模拟滤波器的单位冲激响应。滤波器的单位冲激响应。变换过程变换过程l设模拟滤波器的系统函数只有单阶极点设模拟滤波器的系统函数只有单阶极点l其冲激响应其冲激响应)()(nThnha NkkkassAsH1)( NktskatueAthk1)()(Chapter732IIRIIR数字滤波器的设计方法数字滤波器的设计方法l抽样抽样lZ变换变

21、换l为使频率响应的增益保持不变为使频率响应的增益保持不变取取有有 NknTskanueAnThnhk1)()()()( NkTsknnzeAznhzHk111)()()nT(Th)n(ha NkTskzeTAzHk111)(Chapter733IIRIIR数字滤波器的设计方法数字滤波器的设计方法系统函数之间的关系系统函数之间的关系l单极点的映射单极点的映射l部分分式的系数不变部分分式的系数不变l稳定性不变稳定性不变l零点位置无映射关系零点位置无映射关系混叠失真混叠失真l原因：数字滤波器的频率响应是模拟滤原因：数字滤波器的频率响应是模拟滤波器频率响应的周期延拓，而不是简单波器频率响应的周期延拓，

22、而不是简单的重复。的重复。ks Tkse12()()jakkH eHjTT Chapter734IIRIIR数字滤波器的设计方法数字滤波器的设计方法频率关系频率关系l线性变换线性变换特点特点l频率之间是线性关系频率之间是线性关系l要求模拟滤波器严格带限于折叠频率，要求模拟滤波器严格带限于折叠频率，不能用于设计高通、带阻滤波器不能用于设计高通、带阻滤波器ls平面与平面与z平面之间没有简单的代数映射平面之间没有简单的代数映射关系，必须展成部分分式关系，必须展成部分分式T Chapter735IIRIIR数字滤波器的设计方法数字滤波器的设计方法8 阶跃响应不变法变换思路变换思路l使数字滤波器的阶跃响

23、应模仿模拟滤波使数字滤波器的阶跃响应模仿模拟滤波器的阶跃响应。器的阶跃响应。变换公式变换公式)()(nTgnga )(1)()(1)(zGzzzHzHzzzG )(1)()(11sHsLsGLtgaaa)()(nTgnga )(1)(ngZzzzH Chapter736IIRIIR数字滤波器的设计方法数字滤波器的设计方法混叠失真混叠失真混叠现象比冲激响应不变法小混叠现象比冲激响应不变法小12()()jakkH eHjTT )()(sGzGaezsT kaezkTjsGTzHzzsT 21)(1 kaezskTjsHTzHzzsT 21)(1Chapter737IIRIIR数字滤波器的设计方法数

24、字滤波器的设计方法9 双线性变换法基本思路基本思路l将整个将整个s平面平面压缩压缩变换到一变换到一中介平面中介平面的一的一条横带里，再将此横带变换到整个条横带里，再将此横带变换到整个z平面平面Chapter738IIRIIR数字滤波器的设计方法数字滤波器的设计方法变换过程变换过程压缩压缩扩展到整个复平面扩展到整个复平面再通过标准映射再通过标准映射变换常数变换常数c的选择：采用使模拟滤波器的选择：采用使模拟滤波器与数字滤波器在低频处有较确切的对与数字滤波器在低频处有较确切的对应关系。应关系。TsTsees1111 )(tanc2T1scsczz1z1cs11Tsez1 T2c Chapter73

25、9IIRIIR数字滤波器的设计方法数字滤波器的设计方法逼近的情况：逼近的情况：ls平面的虚轴与平面的虚轴与z平面的单位圆对应；平面的单位圆对应；l稳定的模拟滤波器经双线性变换后得到稳定的模拟滤波器经双线性变换后得到的数字滤波器也稳定。的数字滤波器也稳定。频率变换的关系是非线性的频率变换的关系是非线性的 2tan2 TChapter740IIRIIR数字滤波器的设计方法数字滤波器的设计方法特点特点l避免了频率的混叠避免了频率的混叠失真失真ls与与z之间是简单的之间是简单的代数关系，变换方代数关系，变换方便便l局限性：频率间的局限性：频率间的非线性非线性频率预畸频率预畸Chapter741IIRI

26、IR数字滤波器的设计方法数字滤波器的设计方法10 第一种间接法举例设计步骤设计步骤l将数字指标转换为模拟指标将数字指标转换为模拟指标衰减衰减不变不变频率指标按对应的方法转换频率指标按对应的方法转换l设计对应的模拟滤波器设计对应的模拟滤波器将模拟频率映射为归一化低通的频率指标将模拟频率映射为归一化低通的频率指标确定归一化低通滤波器的阶数确定归一化低通滤波器的阶数查表得到归一化的模拟低通原型滤波器查表得到归一化的模拟低通原型滤波器spAR 、Chapter742IIRIIR数字滤波器的设计方法数字滤波器的设计方法10 第一种间接法举例设计步骤设计步骤l将数字指标转换为模拟指标将数字指标转换为模拟指

27、标衰减衰减不变不变频率指标按对应的方法转换频率指标按对应的方法转换l设计对应的模拟滤波器设计对应的模拟滤波器将模拟频率映射为归一化低通的频率指标将模拟频率映射为归一化低通的频率指标确定归一化低通滤波器的阶数确定归一化低通滤波器的阶数查表得到归一化的模拟低通原型滤波器查表得到归一化的模拟低通原型滤波器spAR 、Chapter743IIRIIR数字滤波器的设计方法数字滤波器的设计方法低通滤波器设计低通滤波器设计例例.给定抽样频率为给定抽样频率为10kHz,要求在频率要求在频率小于小于1kHz的通带内，幅度特性下降小的通带内，幅度特性下降小于于1dB；在频率大于在频率大于1.5kHz的阻带内的阻带

28、内, 衰减大于衰减大于15dB。指定模拟滤波器采。指定模拟滤波器采用巴特沃斯低通滤波器。试分别用冲用巴特沃斯低通滤波器。试分别用冲激响应不变法和双线性变换法设计滤激响应不变法和双线性变换法设计滤波器。波器。Chapter744IIRIIR数字滤波器的设计方法数字滤波器的设计方法解：解：数字低通的技术指标数字低通的技术指标0.2,10.3,15ppssraddBraddB Chapter745IIRIIR数字滤波器的设计方法数字滤波器的设计方法（1 1）冲激冲激响应不变法响应不变法A.模拟低通的技术指标模拟低通的技术指标dBAsraddBRsradsTsspp15,/1031,/10210334

29、 B.设计巴特沃斯低通归一化原型滤波器设计巴特沃斯低通归一化原型滤波器0.10.10.11/23lg (101)/(101)5.886lg(/)(101)7.032 10/pspRAstpRNcpNrad s Chapter746IIRIIR数字滤波器的设计方法数字滤波器的设计方法C.查表得归一化传输函数查表得归一化传输函数D.去归一化，得到实际的传输函数去归一化，得到实际的传输函数E.用用冲激冲激响应不变法得到数字滤波器的传响应不变法得到数字滤波器的传递函数。递函数。654328637. 34641. 71416. 94641. 78637. 311)(sssssssHan 21121121

30、12570. 09972. 016304. 08558. 13699. 00691. 111454. 11428. 26949. 01297. 014466. 02871. 0)( zzzzzzzzzzHlpChapter747IIRIIR数字滤波器的设计方法数字滤波器的设计方法（2）双线性变换法双线性变换法A.模拟低通的技术指标模拟低通的技术指标4444421,1022 100.10.65 10/ ,12 100.151.019 10/ ,15ppststgTTtgradsdBtgradsdB B.设计巴特沃斯低通归一化原型滤波器设计巴特沃斯低通归一化原型滤波器0.10.10.11/23lg

31、(101)/(101)5.306lg(/)(101)7.662 10/pssRAstpANcstNrads Chapter748IIRIIR数字滤波器的设计方法数字滤波器的设计方法C.查表得归一化传输函数查表得归一化传输函数D.去归一化，得到实际的传输函数去归一化，得到实际的传输函数E.用双线性变换法得到数字滤波器的传用双线性变换法得到数字滤波器的传递函数递函数。Chapter749IIRIIR数字滤波器的设计方法数字滤波器的设计方法用双线性法设计一个用双线性法设计一个Butterworth数数字带通滤波器，技术指标为字带通滤波器，技术指标为解解：1）将数字频率预畸为模拟频率将数字频率预畸为模

32、拟频率dBAdBRsppstpst15,2,6 . 0,8 . 0,24. 0,12. 02211 376. 1)2/tan(396. 0)2/tan(0777. 3)2/tan(19. 0)2/tan(22112211 ppppstststst Chapter750IIRIIR数字滤波器的设计方法数字滤波器的设计方法2）计算归一化原型的相关指标）计算归一化原型的相关指标545. 098. 0212012 ppppppB96. 272. 2220222120211 pstpststpstpststBB 172. 2,min21 pststst0.10.1lg(101)/(101)1.978lg

33、(/)psRAstpN取N=2Chapter751IIRIIR数字滤波器的设计方法数字滤波器的设计方法3）归一化低通巴特沃斯模拟低通的系统函数）归一化低通巴特沃斯模拟低通的系统函数4）3dB截止频率截止频率5）去归一化，得到对应）去归一化，得到对应的低通的低通6）模拟频带变换，得到样本模拟带通滤波器）模拟频带变换，得到样本模拟带通滤波器7）双线性变换得到数字带通）双线性变换得到数字带通21( )1.41421LPHsss 1435. 1)110(211 . 0 NRpcp canasHsH)(ppsBssabpsHsH202)()( 1111)()( zzsbpsHzHChapter752II

34、RIIR数字滤波器的设计方法数字滤波器的设计方法11由样本模拟低通直接设计数字滤波器带通变换带通变换l变换关系变换关系0 0 0 0 s js0 jez 1 z22011cos21)1)(1()(00 zzzzzezezsjj Chapter753IIRIIR数字滤波器的设计方法数字滤波器的设计方法l频率关系频率关系l设计参量的关系设计参量的关系解得解得 sincoscos0 220sincoscosppp 110sincoscosppp 21p1p2021p1p2cos()sin()2cos=sinsincos()2pppp Chapter754IIRIIR数字滤波器的设计方法数字滤波器的设

35、计方法例设计一个抽样频率为例设计一个抽样频率为100KHz的数字的数字带通滤波器，其性能要求为：通带范带通滤波器，其性能要求为：通带范围从围从20kHz到到30kHz，在此两频率，在此两频率处衰减不大于处衰减不大于2dB，在，在10kHz和和45kHz处衰减不小于处衰减不小于20dB，采用巴，采用巴特沃思型滤波器。特沃思型滤波器。解：解：1.数字带通的技术指标数字带通的技术指标21210.6 ,0.40.9 ,0.2ppstst Chapter755IIRIIR数字滤波器的设计方法数字滤波器的设计方法2.计算带通的中心频率计算带通的中心频率3.计算样本低通的截止频率计算样本低通的截止频率取取4

36、.计算样本低通滤波器阶数计算样本低通滤波器阶数21021cos()2cos=0cos()2pppp 0111coscos1.376sinststst 0222coscos3.078sinststst 12min(,)1.376ststst 325. 0sincoscos220 ppp 取N=2Chapter7560.10.1lg(101)/(101)1.777lg(/)psRAstpNIIRIIR数字滤波器的设计方法数字滤波器的设计方法5.查表得二阶归一化巴特沃思低通滤波器原查表得二阶归一化巴特沃思低通滤波器原型型6.求样本低通的求样本低通的3dB截止频率截止频率7.去归一化，得到所需样本低通

37、去归一化，得到所需样本低通8.求所需带通滤波器的系统函数求所需带通滤波器的系统函数21( )1.41421anHsss 37. 0)110(411 . 0 pRpc canasHsH)(22011cos21)()( zzzsabpsHzH Chapter757IIRIIR数字滤波器的设计方法数字滤波器的设计方法带阻变换带阻变换l变换关系变换关系0 0 , 0 0 s js1 z2012cos211)()1)(1(00 zzzezezzzsjj 0 jez Chapter758IIRIIR数字滤波器的设计方法数字滤波器的设计方法l频率关系频率关系l设计参量的关系设计参量的关系解得解得0cosco

38、ssin 011coscossin ststst022coscossin stststChapter759211202112cos()sin()2cos=sinsincos()2stststststststst IIRIIR数字滤波器的设计方法数字滤波器的设计方法高通变换高通变换l变换关系变换关系l频率关系频率关系 0 0 0 s s1 z1 zzzs 11 2cot Chapter760IIRIIR数字滤波器的设计方法数字滤波器的设计方法12 数字频域频带变换数字频域频带变换的基本要求数字频域频带变换的基本要求l l 是全通函数是全通函数l因果稳定性不变因果稳定性不变低通低通DF的的Z平面平

39、面所需所需DF的的z平面平面)(11 zGZ)(1 zG jeZ jez Chapter761IIRIIR数字滤波器的设计方法数字滤波器的设计方法数字低通变换数字低通变换l频率映射频率映射l变换函数变换函数l变换函数中的极点变换函数中的极点11)(1111 实实数数zzzGZ 2sin2sinpppp 0 p 0 p Chapter762IIRIIR数字滤波器的设计方法数字滤波器的设计方法数字高通变换数字高通变换l频率映射频率映射l变换函数变换函数l变换函数中的极点变换函数中的极点1,1111111 实实数数zzzzZ 2cos2cospppp 0 p 0 p Chapter763IIRIIR数字滤波器的设计方法数字滤波器的设计方法数字带通变换数字带通变换l频率映射频率映射0 p 0 0 , 21,pp Chapter764IIRIIR数字滤波器的设计方法数字滤波器的设计方法l变换函数变换函数式中式中1,111)(112221121*11*111 zdzddzdzzzzzzGZ111221 kkdkkd 0121212cos2cos2cos2cot2tan pppppppkChapter765IIRIIR数字滤波器的设计方法数字滤波器的设计方法数字带阻变换数