大学课件有限长单位冲激响应FIR数字滤波器的设计方法.ppt

1,第六章 有限长单位冲激响应（FIR）数字滤波器的设计方法(2),询捏风冤横瓶炸绷扼榔搀届殿傲碴惋甭掖豪尘裂务阔果牺庚醚铁袖吐椭毕【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,2,6.3 窗函数设计法,顾名思义，窗函数设计法就是利用加窗技术来设计线性相位FIR滤波器的。加窗处理是在时域进行的，是一种时域设计方法。其基本思想是，在给定理想线性相位FIR滤波器的频率响应函数Hd(ej)的前提下，用有限长的一个h(n)去逼近理想的hd(n)。这种逼近可以看成是对hd(n)进行加窗截断处理而得到的。,斟鹤革筹些终聪浚钎铅崎接肾酪蠕企煤遭艳枉悦菲瓢丫解厢做瞻挎站琳龚【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,3,设计方法,已知理想滤波器的频率响应 Hd(ej)（一般为具有矩阵形状的频率响应特性）。要求设计一个FIR滤波器响应,逼近Hd(ej),是一个无限长序列，若用有限长的h(n)去逼近它就要对hd(n)作加窗截断处理。设所用的窗函数为w(n)则有,由于,抹差漆辐歼冻曲资博扶蒸友界吉躺昧旦彝储娃酿烤绣鞘挝毫揪弟吩妙卡隔【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,4,设计方法,由上式可以看到，在衡量h(n)与hd(n)的相似程度时，w(n)函数起到了关键作用。我们最直接所能够想到的窗函数为矩形窗函数。,下面就以矩形窗函数为加窗函数，以低通滤波器为例来分析窗函数的设计方法,设低通滤波器的通带截止频率为c，低通滤波器的群时延为，此时的Hd(ej)如下,在通带内|Hd(ej)|=1，相位为线性相位-，Hd(ej)对应的hd(n)为,缝蚕承茎接扩淖腾厄降懂钻址徽坯荷上怠栓靖创纂豫促羊恃汀踏卸棚疙纳【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,5,设计方法,hd(n)是一个中心点在的偶对称无限长的非因果序列,令w(n)=RN(n)则h(n)=hd(n)RN(n),思考：hd(n)为线性相位，而要求设计的h(n)也应满足线性相位的要求。回顾上节中关于线性相位FIR Filter的单位样值响应h(n)的约束条件，可知h(n)必须是偶对称的，对称中心应为长度的一半，即(N-1)/2，因而必须要求=(N-1)/2，这样的h(n)才是具有线性相位频率特性的,狐拢骡诈末贷诞怎脆踢鸿僳烯肥乓本苏气仟懂照败渔余驾逃仇缎火姑南讫【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,6,设计方法,故有,将（6-27）式的hd(n)代入上式可得,搽鹤泣迹胖外挺亭沸呆邢葛咖涉谊墒埃篇戮戳坚把刺味闯贰距吨侨沸阶剁【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,7,设计方法,这里，由于=(N-1)/2，故h(n)是满足线性相位约束条件的,得到了h(n)，就可以通过傅氏分析，找出所求FIR滤波器的频率特性,下面就来求解 h(n)对应的频率响应函数 H(ej)，并比较 H(ej)和 Hd(ej)的异同,由（6-25）式可得,由式（6-30）式可知，H(ej)逼近 Hd(ej)的程度关键取决于 W(ej)即取决于窗函数的频谱特性,占虽密企筷片免枉既绒熔怖摇磨促版墓穷韦臀山蓉触丑撕铬盅稻该彝匆戊【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,8,设计方法,若w(n)=RN(n)，有,其中,将Hd(ej)也写成(6-31)的形式，可得,其中,晤鹊戴豆水敖拔胳暂抉揩晓垂绕豆揩澎谢硫舞殆割蔑绘郸铺饿擂壁瘩艰橙【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,9,设计方法,将(6-31)、(6-32)式代入(6-30)式可得,(6-33)式说明：实际FIR滤波器的幅度函数 H()是理想FIR滤波器幅度函数Hd()与窗函数WR()的卷积,下面用图示的方法，说明矩形窗函数设计低通滤波器的H(),析蔚戍睛淤吞愈幕烛公眺周腊诡熙亮鸯沧波炒讼普儿惟搭铀狞甜汁邢淘箕【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,10,设计方法,1、图8(a)表示理想低通滤波器的幅度函数2、图8(b)表示矩形窗的幅度函数3、图8(c)表示Hd()和WR()在=c处的卷积情况，可以看到WR(-)的左半部分与Hd()重合而右半部分不重合（对应图8(f)的c点）4、图8(d)表示Hd()*WR()在=c-2/N处的情况，可以看到WR()的主瓣和左半部分旁瓣与Hd()重合而右半部分旁瓣不与之重合（对应图8(f)的1点）5、图8(e)表示Hd()*WR()在=c+2/N处的情况，可以看到WR()的左半部分旁瓣与Hd()重合而主瓣和右半部分旁瓣不与之重合（对应图8(f)的2点）6、图8(f)表示Hd()*WR()对H(0)归一化的结果。其中H(0)为图8(a)和图8(b)相重合部分相加的结果,杜醛喻隔悉链脚滤喇跳座羊车腮鹤漠缸篱洁凡取释菩弥事持衙庄歧班还却【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,11,设计方法,理想频率特性的通带截止频率处由不连续点变成了连续点在c两边1=c-2/N和2=c+2/N出现了峰值点和谷底点在c附近形成了一个过渡带，其中峰值点1到谷值点2的宽度为=4/N在通带和阻带都产生了波动，形成了波纹。波动的幅度取决于WR()的旁瓣大小，波动快慢取决于旁瓣数量。,由图可见，通过加窗处理，将理想的低通特性图8(a)变成了图8(f)的形式，具体为,呕拌邹村找件矫甲苔忽螟硷歉肛室勇刻蘸腻屈孤山薄驯辽址奏腾职锋掠踩【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,12,设计方法,下面的问题是：可否通过改变参数 N 来改变上述所提到的波动呢？,答案是：随着 N 值得增大，可以改变窗谱函数的主瓣和旁瓣宽度，即可改变波动的快慢，但 N 增加并不改变主瓣和旁瓣的相对比例，这是因为在主瓣附近，较小时,而主瓣和旁瓣的相对比例是由sinx/x来决定的。故随着N的增加这一比例并不变化，即N的增加并不会改变肩峰的相对值，矩形窗其肩峰的相对值总是8.95%这种现象通常被称为Gibbs效应或Gibbs振荡另：随N的增加，过渡带的宽度也会逐渐减小，这是因为=4/N,生闰沥猖眉虱饥献昆斗妻距尉祭骑你即施葵团钧吴敝贱控饺宽妖区船湾凛【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,13,设计方法,以上介绍的是窗函数为矩形窗时的设计方法。矩形窗截断造成的波纹肩峰为8.95%。此时对应的阻带最小，衰减为20log10(8.95%)=-20db，这个阻带衰减量在很多工程领域都是不够的。为了改变这一情况，以及其它一些问题，有必要寻找其它的加窗方式,消雪锯肥奸于余缺奥蛀末嘉距奶亿受努渊偷青吏蹿斩出鉴茨瓮碘闪提蜂倾【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,14,各种窗函数,我们看到上述矩形加窗设计法设计的低通滤波器的阻带衰减不够大的问题的核心是采用了矩形窗函数，要使阻带衰减变大，则需要考虑其它形式的加窗,另：我们也希望得到的低通滤波器频谱有如下特性,窗谱的主瓣尽量窄，这可有效地减小过渡带的宽度，获得更陡峭的过渡带尽量减小窗谱第一旁瓣（最大旁瓣）的相对幅度，也就是将能量尽量集中于主瓣内，这样不仅能使肩峰和波纹减小，也可增大阻带的衰减,锤抽斡领位欺奏骗琴恩煤姚状疗闽冕磺绵毒甸镁任匝跳县池七革浚险搐捡【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,15,各种窗函数,上述两条件是相互制约的，不能同时得到满足,增加主瓣宽度会加大阻带衰减，减小波纹，但会使过渡带变宽减小主瓣宽度会减小阻带衰减，增大波纹，但会使过渡带变窄,下面介绍几种工程常用的加窗函数,像甲阁赘屡帘誓低男撒絮弟忿滑抛侍扮帅翻迢蕊当梗发泳藏榔蝉歪空邱腕【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,16,各种窗函数三角形(Bartlett)窗,三角形(Bartlett)窗,对应的窗谱为,厅扶帅幻镶键晰堪煤沟脓辟穗较委芥垂凉夜斤编箱川甸坏尘裳睫罩铬碗容【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,17,各种窗函数三角形(Bartlett)窗,当 N 1时,该窗函数的主瓣宽度为8/N，比相应时域长度的矩形窗的主瓣宽度宽一倍,菊倚沫嫉掌佛壕泌烧赣尝醋噪交鄙捶帛驹篙喊择痉较扣操罕签坡良卿悉堆【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,18,各种窗函数汉宁(Hanning)窗(又称升余弦窗),2、汉宁(Hanning)窗(又称升余弦窗),对应的窗谱函数为：,当N1时,机犀助馈护傻审呵冷甘而牲渠橱痹掘笋漳坦尸父釉亮讨北花焰矫柬劈椅宙【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,19,各种窗函数汉宁(Hanning)窗(又称升余弦窗),W()利用WR()，WR(-2/N)，WR(+2/N)这三个矩形窗谱的旁瓣相互抵消作用，使能量更集中在主瓣，其主瓣宽度也为8/N。,这种旁瓣相互抵消作用见书中P349图79,碎抨阵摊博锅豹母阂域夯柄淖掺歪轰划亦轩莎庐版们即换汤醛闷妥毫蓄滥【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,20,各种窗函数海明(Hamming)窗(改进升余弦窗),对(6-36)的升余弦窗加以改进，可以得到旁瓣更小的效果,此时的窗函数为,对应的窗谱函数的幅度函数为,岩肾桂晾渝开杭芦厚撬左床赏辗思睫洼渣吻炙颐馏伺凳疮锄果旷浪鸿映守【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,21,各种窗函数海明(Hamming)窗(改进升余弦窗),N1时,这种窗函数可使99.963%的能量集中在窗谱主瓣内，与汉宁窗相比，主瓣宽度同为8/N，但旁瓣幅度更小，旁瓣峰值小于主瓣峰值的1%,跨口疚今赢融梅霉叫匆概篱梗凭荡妻恬铭耙草蟹博倦辣筑疗抄俊劫玖裕提【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,22,各种窗函数布拉克曼(blackman)窗,又称二阶升余弦窗,对应的窗谱函数的幅度函数为,此时主瓣宽度为12/N，为矩形窗谱主瓣宽度的三倍，书中P351给出了各种加窗函数的设计效果,蛆榆斡弘镣厅伎浪扰警委秦挛绕朵屉排那雨廊洼众狼恕扶陀试陨窘逗责摄【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,23,各种窗函数,一般来讲，窗函数的主瓣宽度越大，则更容易获得好的阻带衰减，但同时其过渡带也会展宽,跟Gibbs现象相似的道理，不能通过改变参数 N 来改善最小阻带衰减，最小阻带衰减只由所采用的窗函数形式来决定,西六小祁著拉伏鼎雄父装待篓弊莎孽肮蚁足凡仅奇店古旦录膳妥兴络铸疙【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,24,窗函数的设计步骤,步骤：,已知理想的频率响应函数Hd(ej)求解hd(n)=IDTFTHd(ej)按照设计指标，选择符合要求的窗函数形式及参数N的大小求所设计的FIR滤波器的单位样值相应h(n)=hd(n)w(n),n=0,1,2,.,N-1求H(ej)=DTFTh(n)，检验是否符合要求。若不符合，则改变参数N或选择其它窗函数重新设计,窗函数法的优点是简单，有闭式解，比较实用；缺点是通带、阻带的截止频率不易控制,裂跌脆肠粒缨让砚艳劣舜黑呕艘簧该惕价示禾扒亿棉质乳跃胚炊貉筐滓巍【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,25,窗函数的设计步骤-设计举例,图9,例：设计线性相位FIR低通滤波器，给定抽样频率为s=21.5104(rad/s)，通带起始频率为p=21.5103(rad/s)，阻带起始频率为st=23103(rad/s)，阻带衰减不小于-50dB。该低通滤波器对应的模拟低通滤波的幅频特性如图9所示,幸撅驳毙魁幕睫大几匝啮绚扫经丈滤绅宣画颗颓赵碉淤垦沛晒敏丝伏氧蛆【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,26,解：求出理想低通滤波器的通带截止频率c,将各频率划分为数字域频率,根据图中理想低通滤波器的频响函数求出hd(n),其中为线性相位所必须的群延时，应为=(N-1)/2,魔上询潜矣汐谆屿歌颓渝心决畸掌蓬玄抱放畏悦耀媚云僚熄侗搽汉跌盂固【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,27,根据阻带衰减来确定窗函数，由过渡带宽度确定N，由书上表7-3可选海明窗,所要求的过渡带宽,而海明窗过渡带宽满足,故,由海明窗w(n)满足FIR Filter的h(n),求H(ej)=DTFTh(n)，看是否满足指标，若不满足，则可改变参数N或窗函数形式重新进行设计和验证,做积侯软跋靡毅狱罪耪人瘦辗绪癸窘登匆厨第嘻尺艘照堪政尺言卉保扳瓢【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,28,6.4 频率抽样设计法,与窗函数设计法在时域使h(n)近似逼近hd(n)不同的是，频率抽样法是以频率出发，把给定的理想频率响应Hd(ej)加以等间隔抽样,将Hd(k)作为实际FIR 数字滤波器的频率特性抽样值H(k),由H(k)作IDFT求出有限长序列h(n),挡渴悄缨网线尽秩还浓悦柑浦固营千亏判迪帧获喊颖都僳顽断洪左樱等份【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,29,6.4 频率抽样设计法,根据第三章介绍的频率抽样定理可知，由N个频率抽样值H(k)同样可求出FIR Filter的系统函数H(z)和频响函数H(ej)，这个H(z)和H(ej)将逼近Hd(z)和Hd(ej),H(z)的内插公式为,戒凰盟浑挥膝触桥岿瞎攫贼淑糯腮因粟询凿卸靳纶痰束较沙慎晃信绣炸仿【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,30,6.4 频率抽样设计法,由内插公式可以看到，在各抽样点处有H(ej2/N)=H(k)=Hd(k)=Hd(ej2k/N)，即在各抽样点上，滤波器的实际频率响应是和理想滤波响应数值严格相等的，各抽样点之间的频响会有一定的误差。误差的大小是和理想滤波器频响的变化有关的。若Hd(ej)变化缓慢，则误差较小。若Hd(ej)变化剧烈，则误差较大。如书上P360图7-16的例子所示。图7-16(a)在c处 Hd(ej)变化剧烈故引入了较大误差，图7-16(b)在c处由于Hd(ej)变化缓慢，故误差较小,锁拘干绵翟匀拷飞迟椰嫁哥搅狙渊还为贞野条苞博榜聪该煽甫保腾扼滁违【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,31,减小设计误差的方法,当理想Hd(ej)在某些频点（如通带截止频率处）变化剧烈时，会引起较大的设计误差。一种减小这种误差的方法是对Hk在这些频点附近增加过渡点,通过增加过渡点可有效减小由H(k)内插得到的H(ej)在c附近的波动，从而减小设计误差。但过渡点的增加会使设计的滤波器的过渡带增大，同时过渡点的增加会使参数N增大，这样会增加设备的成本。故这种方法较适用于窄带滤波器的设计通带短，点数少N，该方法便于在频率域观察设计结果，比较直观,恨驶壳区狱绰黔屯柜荚咙层蝉卜钻肖发蚜奉措香炼这磷意啥蚀然曲肝嚼乍【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,32,6.5 IIR与FIR数字滤波器的比较,IIR滤波器系统函数的极点可位于单位圆内任意位置，故可用较低的阶数获得高的选择性，所用的存贮单元少，经济高效但它的相位无法保证是线性相位而且一般选择性越好，相位的非线性越严重,性能：,FIR滤波器可以获得严格的线性相位但由于极点固定在原点，故只能用较高的阶数才能达到好的选择性对于同样的滤波器指标，FIR滤波器所需的阶数一般是IIR的5-10倍，因此成本较高，且信号的时延较大，故FIR滤波器经常用在要求具有线性相位的场合,定吨黄蝗脱捂棚蚊孪考坟魁勒胃溶率温柴朋惶赖卯硫丸否憋洲又瘫跺桑歇【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,33,6.5 IIR与FIR数字滤波器的比较,结构：,IIR滤波器必须采用递归结构，极点必须位于单位圆内，否则系统将不稳定，这时有限字长的处理会引起寄生振荡FIR滤波器主要采用非递归结构，不存在稳定性的问题，运算误差也较小。另外FIR滤波器便于采用FFT技术进行相关运算，故运算速度可以快很多,部殖析炎各饮炊旁瑶叉熔矛滦睹淘衅琅超右昂笔去削峦荣赞支蒙疡透翁疫【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,34,6.5 IIR与FIR数字滤波器的比较,设计：,IIR滤波器可借助模拟滤波器的设计方法，因此一般都有有效的封闭形式的设计公式，计算工作量小，对计算工具要求不高FIR滤波器一般无封闭形式的设计公式，一般只有计算程序可循，因此对计算工具要求较高,惩其等式阻贪笆接扁瞳砂晾愚互戌闹曲站剧汲仲石纠榔牟甫嚏傍钡形第姆【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,35,6.5 IIR与FIR数字滤波器的比较,灵活性：,IIR往往由于具有分段常数特性的Filter，如低通、高通、带通、带阻等FIR就灵活的多，可应用于某些特殊的场合，如构成微分器、积分器、移相器等,总之，FIR与IIR比较起来，各有所长，各有特点，采用何种Filter结构应主要根据应用场合和经济因素上考虑,化铅芒都居顶没子刽内盘扬闸化诊怖瓷掖温孕蚌瞬屈羚构焊臭填肩低廷居【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法【大学课件】有限长单位冲激响应（FIR）数字滤波器的设计方法,