第四章有限长单位脉冲响应滤波器的设计方法课件.ppt

1、第四章 有限长单位脉冲响应（FIR）滤波器的设计方法n本章主要内容引言引言4.1 线性相位FIR数字滤波器的特性4.2 窗口设计法（时间窗口法）4.3 频率采样法4.4 IIR与FIR数字滤器的比较引言n为何要设计为何要设计FIR滤波器？滤波器？n一、一、IIR滤波器的优缺点（回顾）滤波器的优缺点（回顾）IIR数字滤波器的优点：可以利用模拟滤波器设计的结果，而模拟滤波器的设计有大量图表可查，方便简单。IIR数字滤波器的缺点：相位的非线性，将引起频率的色散，若须线性相位，则要采用全通网络进行相位校正,使滤波器设计变得复杂，成本也高。NiiiMiiizbzazHIIR101)(滤波器的系统函数：二

2、、FIR滤波器的优点设FIR滤波器单位冲激响应h(n)长度为N，其系统函数H(z)为：10)(NnnizazH10)()(Niiinxany FIR数字滤波器的差分方程描述 例差分方程yn=xn+0.5xn-1-0.4xn-2a.单位脉冲响应有限 b.单位脉冲响应无限 答案：a系统函数说明H(z)是z-1的N-1次多项式，它在z平面上原点z=0是N-1阶重极点。即除原点外在Z平面上没有极点，H(z)总是稳定的。稳定和线性相位特性是FIR滤波器突出的优点，而且允许设计多通带（或多阻带）滤波器。其中线性相位和多通带滤波器设计都是IIR系统不易实现的FIR滤波器在保证幅度特性满足技术要求的同时，很容

3、易做到有严格的线性相位特性。三、FIR的缺点n1、由于FIR系统只有零点(只在原点有极点)，因此FIR系统不像IIR系统那样易取得比较好的通带与阻带衰减特性。要取得好的衰减特性，一般要H(z)的阶次要高，也即N大。n2、无法利用模拟滤波器的设计结果，一般无解析设计 公式，要借助计算机辅助设计程序完成。四、四、FIR滤波器应用滤波器应用（1）语音处理，图象处理以及数据传输要求线性相位，任意幅度。（即要求信道具有线性相位特性）而FIR数字滤波器具有严格的线性相位，而且同时可以具有任意的幅度特性。（2)另外FIR数字滤波器的单位抽样响应是有限长的，因而滤波器一定是稳定的只要经过一定的延时，任何非因果

4、有限长序列都变成因果的有限序列。（3)FIR可以用FFT算法来实现过滤信号。五、FIR DF设计思路FIR滤波器的设计方法和IIR滤波器的设计方法有很大不同。FIR DF设计的含义是：根据设计指标，求解所选运算结构要求的h(n)或H(z)：1）线性卷积和快速卷积型结构，求FIR DF的h(n).2）级联和频率采样型结构，求FIR DF 的H(z).4.1 线性相位FIR滤波器的特点4.1.1 线性相位的条件线性相位的条件对于长度为N的h(n),传输函数为：总是正值。负值，而称为幅度函数，可以取的实函数，为不同于注意，这里称为相位特性。称为幅度特性，式中，)()(,)()()()()()()()

5、()(10jwjwwjjwNnjwnjweHwwHeHwHwwHewHeHenheH1、H(ejw)线性相位线性相位意味着一个系统的相频特性是频率的线性函数，即式中为常数，此时通过这一系统的各频率分量的时延为一相同的常数，系统的群时延为 ddg)(，第二类线性相位或，第一类线性相位)()(补充定义n1、时延：所谓时延是指信号通过传输通道所需要的传输时间。n2、群时延：它是滤波器平均延迟的一个度量，定义为相频特性对角频率w的一阶导数的负值。即：的线性函数。是即相频特性具有线性相位，常数时，当)(arg)()()(arg)(jjeHDFdddeHd幅频特性和相频特性输出信号与输入信号的幅值比是的非

6、线性函数，称其输出信号与输入信号的幅值比是的非线性函数，称其为系统的幅频特性，记为为系统的幅频特性，记为A(A()它描述了在稳态情况下，当系统输入不同频率的谐波信号时，反映幅值比随频率而变化的规律，其幅值的衰减（A1)特性 输出信号与输入信号的相位差输出信号与输入信号的相位差(或称相移或称相移)也是的非也是的非线性函数，称为系统的相频特性线性函数，称为系统的相频特性它描述了在稳态情况下，当系统输入不同频率的谐波信号时，反映相位差随频率而变化的规律，其相位产生超前（)0或滞后（)0的特性对于物理系统，相位一般是滞后的，即一般是负值.2、FIR滤波器具有线性相位的条件 1010101010)()s

7、in()sin()cos()cos()sin()cos()sin()cos(-NnNnNnNnNnnjjjjnnhHnnhHnnhjnnhjHHenheHeHeH，）（若：明第一类线性相位条件证式中 H()是正或负的实函数。等式中间和等式右边的实部与虚部应当各自相等，同样实部与虚部的比值应当相等:0sincossinsincoscossincossin1010101010NnNnNnNnNnnnhnnhnnhnnhnnh )(21-N)221)()1:10,1210sin10阶数的一半为群时延处在为偶对称，其对称中心要条件即第一类线性相位的充利用数字归纳法可证得有解，此解必为唯一按付氏级数的性

8、质，若nhNnhNnnNhnhNnnhNn1）第一类线性相位线性相位条件证明)()()()()()()(1)1()()(:)1()()(1)1(10110)1(101(1010zHzzHzmhzHzmhzzmhzHmnNznNhznhzHnNhnhnhNNmmNmmNNmmNNnnNnn）令即有第一类线性相位特性为偶对称，该滤波器具另外，证明方法代入上式，得到：将表示为则可将jwNnNnNnNNnnNnNnnNNnnNezzznhzzzznhznhzznhzHzzHzHzH102)1(212110)1(10)1(101)1(2)()(21)()(21)()(21)(:)(有第一类线性相位。那么

9、该滤波器就一定具为偶对称，是实序列，且看出：只要其群时延为相位函数为幅度函数为)()()21()()21()(:)21(cos)()(:)21(cos)(2)()()(1010)21(10)21()21()21(nhnhNwwNwwNnnhwHwNnnheeenheeHNnNnwNjNnwNnjwNnjwNjjw2、第二类线性相位条件证明21)()()1()(221-N)()()()()(NanhnNhnheHeHddajjg为奇对称，对称中心为的充要条件即第二类线性相位类似可以得出，此时群延时若：与前一种不同之处在于增加了 的相移2 20)1(N 20)5.0(N2 偶对称)(nh 奇对称)

10、(nh线性相位特性注意从第二类线性相位看出：零频率w=0有2的截距，说明不仅有：也就是：h(n)为奇对称时，FIR滤波器是一个具有准确的线性相位的理想正交变换网络。个抽样间隔的延时，而且还产生一个90的相移，这种使频率皆为90的网络，称为正交变换网络，它具有重要的理论和实际意义。21N4.1.2 幅度特性由于h(n)的长度N取奇数还是偶数，对H(w)的特性有影响，因此，对于两类线性相位，下面我们分四种情况讨论其幅度特性的特点：（1）h(n)=h(N-1-n)，即h(n)为偶对称偶对称,N=奇数（2）h(n)=h(N-1-n),即h(n)为偶对称偶对称，N=偶数（3）h(n)=-h(N-1-n)

11、，即h(n)为奇对称奇对称,N=奇数（4）h(n)=-h(N-1-n),即h(n)为奇对称奇对称，N=偶数1、第一种情况：h(n)=h(N-1-n)偶对称，N=奇数 212301121211210102121)(NjNnnNjnjNNnnjNjNnnjNnnjjjeNheenhenheNhenhenheHeH 2121cos221)()(23021212302121NhNnnheNheenheeHNnNjNnjNnNnjNjj2/)3(021cos)(221)(NnNnnhNhH21)(N2102/)1(1cos)(212)(,21)0(cos)21(221)(21NnNmnnaHnNhnaN

12、hammNhNhHNnm整理后得：令令 看出：cos(nw)对于w=0,2皆为偶对称，所以幅度函数H(w)也对 w=0,2皆为偶对称。且H(0)、H(/2),H(),H(2)都可不为零。(只要h(N-1)/2)不为零。所以w从从0 2 范围内，无任何约束，可以设计成任何一种滤范围内，无任何约束，可以设计成任何一种滤波器。低通、高通、带通、带阻）波器。低通、高通、带通、带阻）n对称中心N=7)(wHw022关于w=0及w=偶对称)(nh可以设计任何一种滤波器可以设计任何一种滤波器2、第二种情况h(n)=h(N-1-n)偶对称，N=偶数 120211201120112021cos21NnNjNnn

13、NjnjNnnNjNnnjjNnnheeenhenNhenheH21)(N nNhnbnnbHmmNhHmNnNnnhHNnNmNn122)(21cos)(21cos1221221cos)(22/12/112/0令阻滤波器的滤波器，如高通、带）（时况不能用于设计呈奇对称；所以这种情对所以是奇对称，对且由于处，必然有一个零点。在即，（时，）当从上看出：（其中：上面式子可以表示成：0H)()21(cos)2(1)(,0)(0)21(cos12,2,1),2(2)()21(cos()()(21wwHwnwzzHHnwwNnnNhnbnwnbwHNnw=奇对称，H()=0(总是)只能设计低通和带通滤波

14、器。能设计低通和带通滤波器。n对称中心N=6)(wHw022)(nh03、第三种情况 h(n)=-h(N-1-n)奇对称，N=奇数 230221230112123021sin2NnNjNnnNjnjNNnnjNnnjjNnnheeenhenhenheH 2/)1(1211230sin212sin21221)21(sin)(2)(NmNmNnmmNhHmmNhHNmnNnnhH令2)21(-)(N nNhncnncHNn212)(sin)(211所以：看出：sin(nw)对于w=0,2处皆为0即H(w)在w=0,2处必为零。也即H(z)在z=1处都为零。(2)sin(nw)对w=0,2呈奇对称形

15、式不能用于 的滤波器设计，故不能用作低通、高通和带阻滤波器的设计。00)0(HH和n对称中心h(n)=-h(N-1-n)，N=7奇数)(wHw022关于w=0、w=奇对称H(0)=0、H()=0(总是)只能设计带通滤波器。能设计带通滤波器。)(nh04、第四种情况h(n)=-h(N-1-n)，N=偶数)21(sin)12(2)(1221sin221120221NmNnNjjmmNhHNnmNnnheeH令2)21()(N呈偶对称。处呈奇对称，对，对：）由于（处有一零点。在处为零。即在即处为零，在）由于由此看出：（其中：可设由式子：wwHwnzzHwwHwwnNnnNhndwnndwHwnnNh

16、wHNnNn20)()21(sin21)(2,0)(20)21(sin12,3,2,1)12(2)()21(sin)()()21(sin)12(2)(2/12/1关于w=0奇对称、w=偶对称H(0)=0(总是)只能设计带通、高通滤波器。能设计带通、高通滤波器。n对称中心h(n)=-h(N-1-n)，N=6偶数)(wHw022)(nh0总结第一种情况，偶、奇，四种滤波器都可设计。第二种情况，偶、偶，可设计低、带通滤波器，不能设计 高通和带阻。第三种情况，奇、奇，只能设计带通滤波器，其它滤波器 都不能设计。第四种情况，奇、偶，可设计高通、带通滤波器，不能设计低通和带阻四种FIR数字滤波器的相位特性

17、只取决于h(n)的对称性，而与h(n)的值无关。幅度特性取决于h(n)。设计FIR数字滤波器时，在保证h(n)对称的条件下，只要完成幅度特性的逼近即可。例例1 N=5,h(0)=h(1)=h(3)=h(4)=-1/2,h(2)=2，求幅度函数H()。解 为奇数并且h(n)满足偶对称关系a(0)=h(2)=2a(1)=2 h(3)=-1a(2)=2 h(4)=-1H()=2-cos-cos2 =2-(cos+cos2)4.1.3 零点特性 )()k(n-1-N,)n-1-N()1()(FIR1110110)1(1010zHzzkhzzhkzhzHznhzHnNhnhNNkkNNkkNNnnNnn

18、令又即：具有对称性滤波器的单位脉冲响应由于线性相位0)1()1(0)()()()()()(1)()()20)()()()()()(,0)()(11010111)1(111)1(1iiiiNnnNnniiiNiiNiiNiiizHzHzHzHzHznhznhzHzzzzzHnhzHzzHzzHzHzzHzHzzzHzz同理：也必是零点及所以，的零点还必须共轭成对是实数，由于也是零点的零点，则是）若讨论：所以线性相位滤波器的零点必须是互为倒数的共轭对，这种共轭对共有四种1）既不在单位园上，也不在实轴上，有四个互为倒数的两组共轭对，如图 zi z*i 1/zi 1/z*i2）在单位圆上，但不在实轴上

19、，因倒数就是自己的共轭，所以有一对共轭零点，zi,z*i3）不在单位圆上，但在实轴上，是实数,共轭就是自己，所以有一对互为倒数的零点,zi,1/zi4）又在单位圆上，又在实轴上，共轭和倒数都合为一点，所以只有一个零点，只有两种可能，zi=1或zi=-1总结1)h(n)偶对称，N为偶数H()=0 z=-1是H（z）的单根；2）h(n)奇对称，N为奇数，因 H(0)=0,H()=0 所以z=1，z=-1都是H（z）的单根；2)h(n)奇对称，N为偶数，H（0）=0，所以z=1是H（z）的单根。线性相位滤波器是FIR滤波器中最重要的一种，应用最广。实际使用时应根据需用选择其合适类型，并在设计时遵循其

20、约束条件。4.2 窗口设计法n4.2.1 设计思路设计思路n(1)先给定所要求设计的理想滤波器的频率响应Hd(ejw).n(2)设计一个可实现的FIR滤波器频率响应H(ejw)。n(3)由于设计是在时域中进行，使所设计滤波器的h(n)去逼近理想单位取样响应hd(n)如 果 希 望 得 到 的 滤 波 器 的 理 想 频 率 响 应为 ，那么 FIR滤波器的设计就在于寻找一个传递函数 去逼近 ，逼近方法有三种：窗口设计法（时域逼近）频率采样法（频域逼近）最优化设计（等波纹逼近）时间窗口设计法是从单位脉冲响应序列着手，使h(n)逼近理想的单位脉冲响应序列hd(n)。我们知道hd(n)可以从理想频响

21、通过付氏反变换获得 10)(NnjnjenheH)(jdeH221)(onjjdddeeHnh)(jdeH 但一般来说，理想频响 是矩形频率特性，所以，这样得到的理想单位脉冲响应hd(n)往往都是无限长序列，而且是非因果的。但FIR的h(n)是有限长的，问题是怎样用一个有限长的序列去近似无限长的hd(n)。最简单的办法是直接截取一段 hd(n)代替 h(n)。这种截取等效于在hd(n)上施加了一个长度为N的矩形窗，h(n)是通过一个“窗口”所看到的一段，因此，h(n)也可表达为h(n)和一个“窗函数”的乘积，即 h(n)=w(n)hd(n)这一方法通常称为窗口设计法。)(jdeH设计步骤：)(

22、)()()(nwnhnheHddjd)()(nheHj)()(nheHdjd设10)(NnjnjenheH1）由定义)()()2jeHnhDFT3）卷积插值给FIR滤波器加窗的目的在于a.使得单位脉冲响应有限长 b.获得理想滤波器特性 c.产生一正弦变化 练习答案：a4.2.2 矩形窗口法以一个截止频率为 c的线性相位理想低通滤波器为例,讨论FIR的设计问题。a.对于给定的理想低通滤波器 ，计算)()(sin(2121)(01)(nndeedeeHnhaeeHcnjjnjjddccjjdcc为低通滤波器的延时，)(jdeH)(nhd 这是一个以N-1/2为中心的偶对称的无限长非因果序列，如果截

23、取一段n=0N-1的hd(n)作为h(n)，则为保证所得到的是线性相位FIR滤波器，延时 a 应为h(n)长度N的一半,即 2/)1(N)()(01)()()()()(.nRnwnNnonhnwnhnhnhbNRdRd其中：为其它值计算)2/sin()2/sin()(21)()()2/sin()2/sin(11)()()()(*)()()(.R2110RNWNeWeWNeeeeenweWeWeWeHeHeHcRjRjNjnNnjjNnjnjRjjRjRjdjj，其中来表示用幅度函数和相位函数为窗口函数的频谱设计算 2/sin2/sinNWR对频响起作用的是它的幅度函数 理想频响也可以写成幅度函

24、数和相位函数的表示形式 Hd(ej)=Hd()e-j其中幅度函数为 两个信号时域的乘积对应于频域卷积，所以有|0|1)(ccdHdeWeHeWeHeHjRjdjRjdj)(21)(*)()()(deWeHjRjd)()()(21dWHeRdj)()(21如果也以幅度函数 和相位函数来表示 H（ej）,则实际FIR滤波器的幅度函数H（）为正好是理想滤波器幅度函数与窗函数幅度函数的卷积。jjeHeH)()(dWHHRd)()(21)()(H 矩形窗的卷积过程（P95的图4.5来说明）4个特殊频率点看卷积结果：（1）=0时,H(0)等于 在-c,c内的积分面积，（2）=c时，一半重叠，H(c)=0.

25、5 H(0);（3）=c 2/N时，第一旁瓣（负数）在通带外，出现正肩峰；（4）=c+2/N 时，第一旁瓣（负数）在通带内，出现负肩峰。)(RW()H随，绕零值波动()(0)HH随，绕波动2cN2cN窗口函数对理想特性的影响：改变了理想频响的边沿特性，形成过渡带，其宽度取决于窗函数的主瓣宽度，WR()的主瓣宽度为 。过渡带两旁产生肩峰和余振（带内、带外起伏），取决于 WR()的旁瓣，旁瓣所包围的面积越大，通带波动增加，阻带衰减减少。N增加,过渡带宽减小,肩峰值不变。因主瓣附近 其中x=N/2,所以N的改变只能改变窗谱的主瓣宽度不能改变主瓣与旁瓣的比例关系，只能改变WR（）的绝对值大小和起伏的密

26、度，当N增加时，幅值变大，频率轴变密，而最大肩峰永远为8.95%，这种现象称为吉布斯（Gibbs）效应。xxNNNNNWRsin2/)2/sin()2/sin()2/sin()(N4 肩峰值的大小决定了滤波器通带内的平稳程度和阻带内的衰减，所以对滤波器的性能有很大的影响改变窗函数的形状，可改善滤波器的特性，窗函数有许多种，但要满足以下两点要求：窗谱主瓣宽度要窄，以获得较陡的过渡带；相对于主瓣幅度，旁瓣要尽可能小，使能量尽量集中在主瓣中，这样就 可以减小肩峰和余振，以提高阻带衰减和通带平稳性。但实际上这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑制。例 用矩形窗设计一个FIR线性相位低

27、通数字滤波器。已知求出并h(n)的表达式。21,5.0Nc为其他故：其中，nnnnnwnhnhNnndeedeeHnheeHdccccnjjnjjddccjjdcc,0200,)10(2sin)()()(5.0 102/)1()()(sin21)(21)(01)(00.250.50.751-40-30-21-100N=15N=31用矩形窗设计的c=/2 FIR滤波器的幅度响应几种常用的窗函数：1.矩形窗，上面已讲过，不再细述2.汉宁窗（升余弦窗）)(12cos1 21)(nRNnnwN)(25.0)(5.01212nReenRNNnjNnjN3.汉明窗（改进的升余弦窗）4.布莱克曼窗（三阶升余

28、弦窗）)(12cos46.054.0)(nRNnnwN)(14cos08.012cos5.042.0)(nRNnNnnwN5.051cN窗函数主瓣宽度过渡带宽旁瓣峰值衰减(dB)阻带最小衰减(dB)矩形N/4N/8.1-13-21汉宁N/8N/2.6-31-44汉明N/8N/6.6-41-53布莱克曼N/12N/11-57-74对窗函数的总的要求n希望它频谱的主瓣尽量地窄，旁瓣尽量地小，使频域的能量能主要集中在主瓣内。n采用窗函数法，设计简单，方便，也实用。但要求用计算机，且边界频率不易控制.长度N也不易一次确定，要反复几次才能求得满意结果。nFIR DF设计的窗函数法不但可以用来设计普通的L

29、P，HP，BP及BS滤波器，也可以用来设计一些特殊的滤波器，例如差分滤波器，希尔伯特滤波器。3.矩形窗设计的FIR滤波器过渡带最窄。错 对 工程上，常给定频域上的技术指标，。因此，采用频率采样法更为直接，尤其对于Hd(ejw)公式较复杂，或Hd(ejw)不能用封闭公式表示而用一些离散值表示时，频率采样设计法更为方便，有效。jnhNIDFTNNkjdjdeHnhkHeHeHd)(2)()(不同于点点频率取样确定 内插公式4.3 频率采样法窗函数法与频率采样法的区别n1、窗 函 数 法 是 从 时 域 出 发，把 理 想 的 hd(n)用 一 定 形 状 的 窗 函 数 截 取 成 有 限 长 的

30、 h(n)，以 此 h(n)求 近 似 理 想 的hd(n)，这 样 得 到 的 频 率 响 应H(ejw)逼 近 于 所 求 的 理 想 频 率 响 应Hd(ejw)。n2、频 率 抽 样 法 则 从 频 域 出 发，把 给 定 的 理 想 频 率 响 应Hd(ejw)加 以 等 间 隔 抽 样 然 后 以 此 Hd(k)作 为 实 际 FIR 滤 波 器 的 频 率 特 性 的 抽 样 值 H(k)，知 道 H(k)后，由 DFT 定 义，可 用 频 域 的 这N 个 抽 样 值 H(k)求 唯 一 确 定 的 有 限 长 序 列 h(n)，利 用 这 N 个 频 域 抽 样 值 H(k)

31、同 样 可 得 FIR 滤 波 器 的 系 统 函 数 H(z)及 频 率 响 应。4.3.1频率抽样法基本原理n频 率 抽 样 法 从 频频 域域 出出 发发，把 给 定 的 理 想 频 率 响 应Hd(ejw)加 以 等 间 隔 抽 样得到Hd(k)，即10221,2,1,0,)(1)().()(1,2,1,0,)()(NkknNjddkNwjwddNnekHNnhnhIDFTkHNNkeHkH，得到进行点再对。型）式适合频率采样结构（滤波器）式可用来设计直接型（。形成滤波器的系统函数采样值上式就是直接利用频率看出：令（其系统函数为FIRFIRkHzWkHNzzHeWzekHNzzHznh

32、zHdNkkNNjNkkNjdNNnn21)(1)(1)(,)2(1)(1)1()()101/2101210单位圆上的频响为：1,1,001)(,1,1,0,22/2sin)2/sin(1)(2/2sin2/sin)(11)(122110102110/2NikikieNiiNeNkNNeekHeNkNkHNeekHNeeHiNjkNkNjjkNkjkNkNkNjNkjNkjNjj则令式中这是一个内插公式。内插公式表明：在每个采样点上，逼近误差为零，频响 严格地与理想频响的采样值 H(k)相等；在采样点之间，频响由各采样点的内插函数延伸迭加而形成，因而有一定的逼近误差，误差大小与理想频率响应的曲

33、线形状有关，理想特性平滑，则误差小；反之，误差大。在理想频率响应的不连续点附近，会产生肩峰和波纹。N增大，则采样点变密，逼近误差减小。)()(kHeHkj)(jeH)(jeHn抽样点上，频率响应严格相等n抽样点之间，加权内插函数的延伸叠加n变化越平缓，内插越接近理想值，逼近误差较小4.3.2.用频率采样法设计线性相位滤波器的约束条件 为了设计线性相位的FIR滤波器，单位脉冲响应序列，要满足一定的约束条件。另外，前已指出,具有线性相位的FIR滤波器，其单位脉冲响应h(n)是因果，有限长，实序列，由此得到的幅频和相频特性，就是对H(k)的约束。)1()(nNhnh1、第一种情况线性相位FIR滤波器

34、kNkkkjkjkNjkNjwNjjwHHHNkkNNeHekNHeHkHeHkHwHwHwHewHeHNnhkk必须满足偶对称：得出：（也用幅值与相角表示：（如果抽样值（即应为偶对称。其中幅度函数（为奇数；是偶对称，)11(2)21()2()()()2()()()(222)1(2、第二种情况线性相位FIR滤波器kNkkkjkjkNjkNjwNjjwHHHNkkNNeHekNHeHkHeHkHwHwHwHewHeHNnhkk必须满足奇对称：得出：（也用幅值与相角表示：（如果抽样值（即应为奇对称。其中幅度函数（为偶数；是偶对称，)11(2)21()2()()()2()()()(222)1(3、第

35、三种情况线性相位FIR滤波器kNkkkjkjkNjkNjwNjjwHHHNkkNNeHekNHeHkHeHkHwHwHwHewHeHNnhkk必须满足偶对称：得出：（也用幅值与相角表示：（如果抽样值（即应为偶对称。其中幅度函数（为奇数；是奇对称，)11(2)21()2()()()2()()()(222)1(4、第四种情况线性相位FIR滤波器kNkkkjkjkNjkNjwNjjwHHHNkkNNeHekNHeHkHeHkHwHwHwHewHeHNnhkk必须满足奇对称：得出：（也用幅值与相角表示：（如果抽样值（即应为奇对称。其中幅度函数（为偶数；是奇对称，)11(2)21()2()()()2()

36、()()(222)1(例：设计一个FIR数字 LP 滤波器，其理想特性为 采样点数 N=33，要求线性相位。解：能设计低通线性相位数字滤波器的只有单位脉冲响应偶对称两种情况，因N为奇数，所以只能选择第一种情况。即 h(n)=h(N-1-n)，幅频特性关于偶对称，也即 HK 偶对称。利用 HK 的对称性,求频响采样值。5.005.001jdeH根据指标要求，共有33个取样点，所以第k点对应频率为 而截止频率 0.5位于 之间，所以，k=08时，取样值为1；根据对称性，故 k=2532时，取样值也为1，因 k=33 为下一周期，所以0区间有9个值为 1的采样点，2区间有8个值为 1 的采样点，因此

37、：k33293328332和258,312321330.,HHHHHHHH所以代入内插公式，求将258,02/33/2sin33233sin3312/2sin2/sin1)()(32033322124903225;801163201016322kHekkHeeNkNHNeHeHeHkHkkNkkHkjkkNkNkjNkjkjjjkkNkkk818181813225332sin33233sin332sin33233sin2sin233sin331)(332sin33233sin332sin33233sin33/)33(2sin33)33(233sin2/33/2sin33233sinkjnnnk

38、kkkkkeHnnknnnkkH小结频率采样设计法优点：直接从频域进行设计，物理概念清楚，直观方便；适合于窄带滤波器设计，这时频率响应只有少数几个非零值。典型应用：用一串窄带滤波器组成多卜勒雷达接收机，覆盖不同的频段，多卜勒频偏可反映被测目标的运动速度；缺点：截止频率难以控制。因频率取样点都局限在2/N的整数倍点上，所以在指定通带和阻带截止频率时，这种方法受到限制，比较死板。充分加大N，可以接近任何给定的频率，但计算量和复杂性增加。4.5 IIR和FIR滤波器的比较n前面讨论了IIR和FIR两种滤波器传输函数的设计方法。n这两种滤波器究竟各自有什么特点？在实际运用时应该怎样去选择它们呢？n为此

39、对这两种滤波器作一简单的比较。n1.从性能上比较n2.从结构上比较n3.从设计工具上比较1、从性能上进行比较n从性能上来说，IIR滤波器传输函数的极点可位于单位圆内的任何地方，因此可用较低的阶数获得高的选择性，所用的存贮单元少，所以经济而效率高。但是这个高效率是以相位的非线性为代价的。选择性越好，则相位非线性越严重。相反，FIR滤波器却可以得到严格的线性相位，然而由于FIR滤波器传输函数的极点固定在原点，所以只能用较高的阶数达到高的选择性；对于同样的滤波器设计指标，FIR滤波器所要求的阶数可以比IIR滤波器高510倍，结果，成本较高，信号延时也较大；如果按相同的选择性和相同的线性要求来说，则I

40、IR滤波器就必须加全通网络进行相位较正，同样要大增加滤波器的节数和复杂性。从结构上看nIIR滤波器必须采用递归结构，极点位置必须在单位圆内，否则系统将不稳定。相反，FIR滤波器主要采用非递归结构，不论在理论上还是在实际的有限精度运算中都不存在稳定性问题，运算误差也较小。此外，FIR滤波器可以采用快速付里叶变换算法，在相同阶数的条件下，运算速度可以快得多。从设计工具看nIIR滤波器可以借助于模拟滤波器的成果，因此一般都有有效的封闭形式的设计公式可供准确计算，计算工作量比较小，对计算工具的要求不高。FIR滤波器设计则一般没有封闭形式的设计公式。窗口法虽然仅仅对窗口函数可以给出计算公式，但计算通带阻

41、带衰减等仍无显式表达式。一般，FIR滤波器的设计只有计算程序可循，因此对计算工具要求较高。另外nIIR滤波器虽然设计简单，但主要是用于设计具有片段常数特性的滤波器，如低通、高通、带通及带阻等，往往脱离不了模拟滤波器的格局。而FIR滤波器则要灵活得多，尤其它能易于适应某些特殊的应用，如构成微分器或积分器，或用于Butterworth、Chebyshev等逼近不可能达到预定指标的情况，例如，由于某些原因要求三角形振幅响应或一些更复杂的幅频响应，因而有更大的适应性和更广阔的天地。总结n从上面的简单比较我们可以看到IIR与FIR滤波器各有所长，所以在实际应用时应该从多方面考虑来加以选择。例如，从使用要求上来看，在对相位要求不敏感的场合，如语言通讯等，选用IIR较为合适，这样可以充分发挥其经济高效的特点，而对于图像信号处理，数据传输等以波形携带信息的系统，则对线性相位要求较高，如果有条件，采用FIR滤波器较好，当然，在实际应用中应考虑经济上的要求以及计算工具的条件等多方面的因素。重点习题n4.1，4.6，4.8