数字信号处理程佩青第三版课件_第六章__IIR滤波器的设计方法-1

整理文档很辛苦,赏杯茶钱您下走!

免费阅读已结束,点击下载阅读编辑剩下 ...

阅读已结束,您可以下载文档离线阅读编辑

资源描述

第六章IIR滤波器的设计主要内容理解数字滤波器的基本概念了解最小相位延时系统理解全通系统的特点及应用掌握冲激响应不变法掌握双线性变换法掌握Butterworth、Chebyshev低通滤波器的特点了解利用模拟滤波器设计IIR数字滤波器的设计过程了解利用频带变换法设计各种类型数字滤波器的方法6.1引言数字滤波器:是指输入输出均为数字信号,通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。高精度、稳定、体积小、重量轻、灵活,不要求阻抗匹配,可实现特殊滤波功能优点:1、滤波器的基本概念(1)滤波器的功能滤波器的功能是对输入信号进行滤波以增强所需信号部分,抑制不要的部分。a)时域说明b)频域说明(2)四种基本的滤波器四种基本滤波器为低通(LP)、高通(HP)、带通(BP)和带阻滤波器(BRF):(3)四种基本滤波器的数字表示低通高通带通带阻)(jeH)(jeH)(jeH)(jeHω2ππππ2ω2πππ2πω2πππ2πω2πππ2π2、LP到其他滤波器的变换由LP实现的HPLP实现的BPLP实现的BRF3、滤波器的性能指标带宽:当幅度降低到0.707时的宽度称为滤波器的带宽(3dB带宽)通带、阻带与过渡带:信号允许通过的频带为通带,完全不允许通过的频带为阻带,通带与阻带之间为过渡带。滚降与滚降率:滤波器幅频特性在过渡带的衰减和衰减速度称为滚降与滚降率。阻带衰减:输入信号在阻带的衰减量带内平坦度:通带和阻带内的平坦程度4、数字滤波器的设计步骤数字滤波器的设计三个步骤:(1)按要求确定滤波器的性能参数;(2)用一个因果稳定的离散线性移不变系统的系统函数去逼近去逼近这一性能要求;(3)用有限精度的运算实现;实现可以采用通用计算机,也可以采用DSP。5、数字滤波器的技术要求选频滤波器的频率响应:()()()jjjjHeHee为幅频特性:表示信号通过该滤波器后各频率成分的衰减情况()jHe为相频特性:反映各频率成分通过滤波器后在时间上的延时情况()j:通带截止频率c:阻带截止频率st:通带容限1:阻带容限2st2()jHe阻带:cst过渡带:c11()1jHe通带:理想滤波器不可实现,只能以实际滤波器逼近通带最大衰减:1011()20lg20lg()20lg(1)()ccjjjHeHeHe阻带最小衰减:2022()20lg20lg()20lg()ststjjjHeHeHe其中:0()1jHe当时,()2/20.707cjHe称为3dB通带截止频率13dBc6、表征滤波器频率响应的特征参量幅度平方响应2*()()()jjjHeHeHe1()()()()jjjzeHeHeHzHz的极点既是共轭的,又是以单位圆成镜像对称的1()()HzHzRe[]zIm[]jz01aa*a*1/aH(z)的极点:单位圆内的极点相位响应()()()Re()Im()jjjjejjHeHeeHejHe*()()()jjjjeHeHee*1()()ln2()jjjHeejHe11()ln2()jzeHzjHz()jHe2()*()()jjjejHeeHeIm[()]()arctanRe[()]jjjHeeHe相位响应:()()jjdeed()1Re()jzedHzzdzHz群延迟响应相位对角频率的导数的负值()je若滤波器通带内=常数,则为线性相位滤波器7、IIR数字滤波器的设计方法先设计模拟滤波器,再转换为数字滤波器01()1MkkkNkkkbzHzaz用一因果稳定的离散LSI系统逼近给定的性能要求:,kkab即为求滤波器的各系数计算机辅助设计法s平面逼近:模拟滤波器z平面逼近:数字滤波器6.2最小与最大相位延时系统、最小与最大相位超前系统LSI系统的系统函数:1()11111(1)()()(1)()MMmmNMmmNNkkkkczzcHzKKzdzzd()arg[()]11()()()()jMjmjjNMjjHemNjkkecHeKeHeeed频率响应:模:11()MjjmmNjkkecHeKed各零矢量模的连乘积各极矢量模的连乘积11()argarg[]arg[]()jMNjjmkmkHeecedNMK相角:02,2当位于单位圆内的零/极矢量角度变化为2位于单位圆外的零/极矢量角度变化为0Re[]zIm[]jz011()argarg[]arg[]()jMNjjmkmkHeecedNMKiommMioppN令:单位圆内零点数为mi单位圆外的零点数为mo单位圆内的极点数为pi单位圆外的极点数为po2()arg2()22jiiHeNMmpK则:11()argarg[]arg[]()jMNjjmkmkHeecedNMK全部极点在单位圆内:po=0,pi=N因果稳定系统2()arg222()jiiHempNMK22imM1)全部零点在单位圆内:,0iomMmarg[]02)全部零点在单位圆外:0,iommMarg[]2M20om为最小相位延时系统为最大相位延时系统,1zrrn0时,h(n)=0相位延时系统逆因果稳定系统1)全部零点在单位圆内:,0iomMmarg[]2N2()arg222()jiiHempNMK22()0imNM2)全部零点在单位圆外:0,iommMarg[]2()NM全部极点在单位圆外:po=N,pi=0为最大相位超前系统为最小相位超前系统相位超前系统,1zrrn0时,h(n)=0最小相位延时系统的性质1)在相同的系统中,具有最小的相位滞后()jHe2)最小相位延时系统的能量集中在n=0附近,而总能量相同5)级联一个全通系统,可以将一最小相位系统转变成一相同幅度响应的非最小相位延时系统4)在相同的系统中,唯一()jHemin()hn3)最小相位序列的最大:min(0)hmin(0)(0)hh22min00()()mmnnhnhn1mN1122min00()()NNnnhnhn6.3全通系统对所有,满足:()1japHe称该系统为全通系统一阶全通系统:11()1apzaHzaaz为实数01aza极点:1/za零点:零极点以单位圆为镜像对称za极点:1*1()1apzaHzaaz为复数01a*1/za零点:实系数二阶全通系统1*111()11*apzazaHzazaz1a两个零点(极点)共轭对称*zaa,极点:*1/zaa,1/零点:零点与极点以单位圆为镜像对称N阶数字全通滤波器1*11()1NkkkzaHzaz1(1)111(1)11...1...NNNNNNNNddzdzzdzdzdz1()()NzDzDz极点:的根()Dz1jpzrer零点:的根1()Dz11jozerr全通系统的应用min()()()apHzHzHz1)任一因果稳定系统H(z)都可以表示成全通系统Hap(z)和最小相位系统Hmin(z)的级联11*100()()()()HzHzzzzz令:其中:H1(z)为最小相位延时系统,为单位圆外的一对共轭零点*0001/1zzz,1/,*1111*00100*110011()()11zzzzHzHzzzzzzzzz11**1100100*1100()1111zzzzHzzzzzzzzzmin()()apHzHz把H(z)单位圆外的零点:映射到单位圆内的镜像位置:构成Hmin(z)的零点。*0001/,1/,1zzzz*00,zzz而幅度响应不变:minmin()()()()jjjjapHeHeHeHeP231图6-62)级联一个全通系统可以使非稳定滤波器变成一个稳定滤波器把非稳定系统的单位圆外的极点映射到单位圆内单位圆外极点:11jzerr,1111()11jjapjjzrezreHzrezrez()()()dapHzHzHz0()()()()dapdd3)作为相位均衡器,校正系统的非线性相位,而不改变系统的幅度特性()()()jjjdapHeHeHe()()()()dapjjjdapHeHee22200()()()apde利用均方误差最小准则求均衡器Hap(z)的有关参数6.4用模拟滤波器设计IIR数字滤波器设计思想:s平面z平面模拟系统数字系统()()aHsHzH(z)的频率响应要能模仿Ha(s)的频率响应,即s平面的虚轴映射到z平面的单位圆因果稳定的Ha(s)映射到因果稳定的H(z),即s平面的左半平面Re[s]0映射到z平面的单位圆内|z|1设计方法:-冲激响应不变法-阶跃响应不变法-双线性变换法6.5冲激响应不变法一、变换原理数字滤波器的单位冲激响应h(n)模仿模拟滤波器的单位冲激响应ha(t))()(nThnhaT—抽样周期)()(nTxnxa)()(nTynya)()(nThnha)(txa)(tya)(tha设)()(sHthaa)()(zHnh则:masTezTjmsHTzH)2(1|)(从频率响应来看:kakajTkjHTkTjjHTeH)2(1)2(1)(数字滤波器的频率响应是模拟滤波器频率响应的周期延拓,周期为2/T只有当模拟滤波器的频率响应是带限的,且带限于折叠频率以内时,即2,0)(saTjH才能使数字滤波器的频响在折叠频率以内重现模拟滤波器的频响而不产生混迭失真)(),(1)(jHTeHaj,sfTTTTccTT2ssfTT混迭实际系统不可能严格限带,都会混迭失真,在||s/2处衰减越快,失真越小当滤波器的设计指标以数字域频率c给定时,不能通过提高抽样频率来改善混迭现象二、模拟滤波器的数字化1()NkakkAHsss111kNksTkAez()()()()()aaaHshthnThnHz11()[()]()kNstaakkhtLHsAeut1()()()kNsnTakkhnhnTAeunT()()nnHzhnz01kNnsTnknkAez110kNnsTkknAez系数相同:kA1()NkakkAHsss11()1kNksTkAHzez极点:s平面z平面kssksTze稳定性不变:s平面z平面Re[]0ks1ksTe1()jaHeHjTT11()1kNksTkTAHzez当T很小时,数字滤波器增益很大,易溢出,需修正()()ahnThnT令:2()jakkHeH

1 / 104
下载文档,编辑使用

©2015-2020 m.777doc.com 三七文档.

备案号:鲁ICP备2024069028号-1 客服联系 QQ:2149211541

×
保存成功