数字信号处理第7章课件

第七章有限长单位脉冲响应（FIR） 滤波器的设计方法

序言§7.1线性相位FIR数字滤波器的特性§7.2窗口设计法（时间窗口法）

§7.3频率取样法

§7.4FIR数字滤波器的最优化设计§7.5IIR与FIR数字滤器的比较1第七章有限长单位脉冲响应（FIR）

FIR数字滤波器的差分方程描述①

对应的系统函数

因为它是一种线性时不变系统，可用卷积和形式表示

③

比较①、③得：

序言

2

序言2

FIR数字滤波器的特点(与IIR数字滤波器比较)：

优点：（1）很容易获得严格的线性相位，避免被处理的信号产生相位失真，这一特点在宽频带信号处理、阵列信号处理、数据传输等系统中非常重要（2）永远稳定，无稳定性问题

3

FIR数字滤波器的特点(与IIR数字滤波器比较)：3（3）任何一个非因果的有限长序列，总可以通过一定的延时，转变为因果序列，所以因果性总是满足；（4）无反馈运算，运算误差小。4（3）任何一个非因果的有限长序列，总可以通过一定的延时，转缺点：（1）因为无极点，要获得好的过渡带特性，需以较高的阶数为代价；（2）无法利用模拟滤波器的设计结果，一般无解析设计公式，要借助计算机辅助设计程序完成。5缺点：5§7.1线性相位FIR数字滤波器的特性

1线性相位的条件线性相位意味着一个系统的相频特性是频率的线性函数，即式中为常数，此时通过这一系统的各频率分量的时延为一相同的常数，系统的群时延为

6§7.1线性相位FIR数字滤波器的特性1线性相位的条FIR滤波器的DTFT为式中H(ω)是正或负的实函数。等式中间和等式右边的实部与虚部应当各自相等，同样实部与虚部的比值应当相等:7FIR滤波器的DTFT为式中H(ω)是正或负的实函数。等式将上式两边交叉相乘，再将等式右边各项移到左边，应用三角函数的恒等关系满足上式的条件是8将上式两边交叉相乘，再将等式右边各项移到左边，应用三角函数的另外一种情况是，除了上述的线性相位外，还有一附加的相位，即利用类似的关系，可以得出新的解答为9另外一种情况是，除了上述的线性相位外，还有一附加的相位，即

偶对称

奇对称10偶对称奇对称10

分四种情况(1)．h(n)偶对称，N为奇数h(n)=h(N-1-n)2线性相位FIR滤波器的幅度特性11

分四种情况2线性相位FIR滤波器的幅度特性111212令，则令则由于偶对称，因此对这些频率也呈偶对称。13令，则令则

(2)．h(n)偶对称，N为偶数h(n)=h(N-1-n)14(2)．h(n)偶对称，N为偶数h(n)=h(N令，则15令，则15

或写为：

由于奇对称，所以对也为奇对称，且由于时，处必有一零点，因此这种情况不能用于设计时的滤波器，如高通、带阻滤波器。16

或写为：由于(3).h(n)奇对称，N为奇数，h(n)=-h(N-1-n)17(3).h(n)奇对称，N为奇数，h(n)=-h(N-1

令n=m+(N-1)/2，得：

18

18

所以

由于点呈奇对称，所以对这些点也奇对称。由于时，相当于H（z）在处有两个零点，不能用于的滤波器设计，故不能用作低通、高通和带阻滤波器的设计。19

所以由于(4).h(n)奇对称，N为偶数令20(4).h(n)奇对称，N为偶数令20

由于在ω=0,２π处为零，所以H(ω)在ω=0,2π处为零，即H(z)在z=1上有零点，并对ω=0,2π呈奇对称。对ω=π呈偶对称。

21

由于四种线性相位FIR滤波器关于0、π、2π偶对称关于π奇对称H(π)=0关于0、π、2π奇对称关于0、2π奇对称π偶对称22四种线性相位FIR滤波器关于0、π、2π偶对称关于π奇对称H四种线性相位FIRDF特性:第一种情况，偶、奇，四种滤波器都可设计。第二种情况，偶、偶，可设计低、带通滤波器，不能设计

高通和带阻。第三种情况，奇、奇，只能设计带通滤波器，其它滤波器都不能设计。第四种情况，奇、偶，可设计高通、带通滤波器，不能设计低通和带阻。23四种线性相位FIRDF特性:第二种情况，偶、偶，可设计低、例1N=5,h(0)=h(1)=h(3)=h(4)=-1/2,h(2)=2，求幅度函数H(ω)。解Ｎ为奇数并且h(n)满足偶对称关系a(0)=h(2)=2a(1)=2h(1)=-1a(2)=2h(0)=-1H(ω)=2

-cosω-cos2ω=2-(cosω+cos2ω)24例1N=5,h(0)=h(1)=h(3)2525小结：

四种FIR数字滤波器的相位特性只取决于h(n)的对称性，而与h(n)的值无关。幅度特性取决于h(n)。设计FIR数字滤波器时，在保证h(n)对称的条件下，只要完成幅度特性的逼近即可。注意：当H(ω)用│H(ω)│表示时，当H(ω)为奇对称时，其相频特性中还应加一个固定相移π。26小结：26

3线性相位FIR滤波器的零点特性

27

3线性相位FIR滤波器的零点特性27由该式可看出，若z=zi是H(z)的零点，则z=zi-1也一定是H（z）的零点。由于h(n)是实数，H(z)的零点还必须共轭或对，所以z=zi*及z=1/z*也必是零点。所以线性相位滤波器的零点必须是互为倒数的共轭对，即成四出现，这种共轭对共有四种28由该式可看出，若z=zi是H(z)的零点，则z=zi-1也一可能的情况：①既不在单位园上，也不在实轴上，有四个互为倒数的两组共轭对，ziz*i1/zi1/z*i图4.2(a)②在单位圆上，但不在实轴上，因倒数就是自己的共轭，所以有一对共轭零点，zi,z*i图4.2（b）③不在单位圆上，但在实轴上，是实数,共轭就是自己，所以有一对互为倒数的零点,zi,1/zi图4.2（c）④又在单位圆上，又在实轴上，共轭和倒数都合为一点，所以成单出现，只有两种可能，zi=1或zi=-1图4.2(d),p9229可能的情况：293030我们从幅度响应的讨论中已经知道，对于第二种FIR滤波器（h(n)偶对称，N为偶数），，

即是的零点，既在单位圆，又在实轴，所以必有单根同样道理，对于第三种FIR滤波器，h(n)奇对称，N为奇数，因所以z=1，z=-1都是H（z）的单根；对于第四种滤波器，h(n)奇对称，N为偶数，H（O）=0，所以z=1是H（z）的单根。所以，h(n)奇对称→H(0)=0N为偶数→H(π)=031我们从幅度响应的讨论中已经知道，对于第二种FIR滤波器（h(线性相位滤波器是FIR滤波器中最重要的一种，应用最广。实际使用时应根据需用选择其合适类型，并在设计时遵循其约束条件。32线性相位滤波器是FIR滤波器中最重要的一种，应用最广。实际使3333§7.2窗口设计法（时域）

如果希望得到的滤波器的理想频率响应为，那么FIR滤波器的设计就在于寻找一个传递函数去逼近逼近方法有三种：窗口设计法（时域逼近）频率采样法（频域逼近）最优化设计（等波纹逼近）

34§7.2窗口设计法（时域）如果希望得到的滤波时间窗口设计法是从单位脉冲响应序列着手，使h(n)逼近理想的单位脉冲响应序列hd(n)

。hd(n)可以从理想频响的傅立叶反变换获得35时间窗口设计法是从单位脉冲响应序列着手，使h(n)逼近理想的但一般来说，理想频响是分段恒定，在边界频率处有突变点，所以，这样得到的理想单位脉冲响应hd(n)往往都是无限长序列，而且是非因果的。但FIR的h(n)是有限长的，问题是怎样用一个有限长的序列去近似无限长的hd(n)。

36但一般来说，理想频响是分段恒定，在最简单的办法是直接截取一段hd(n)代替h(n)。这种截取可以形象地想象为h(n)是通过一个“窗口”所看到的一段hd(n)，因此，h(n)也可表达为hd(n)和一个“窗函数”的乘积，即h(n)=w(n)hd(n)在这里窗口函数就是矩形脉冲函数RN(n)，当然以后我们还可看到，为了改善设计滤波器的特性，窗函数还可以有其它的形式，相当于在矩形窗内对hd(n)作一定的加权处理。37最简单的办法是直接截取一段hd(n)代替h(n)。这种3838

设计步骤：1）由定义3）卷积插值39

设计步骤：1）由定义3）卷积插值394040

以一个截止频率为ωc的线性相位理想低通滤波器为例,讨论FIR的设计问题。a.对于给定的理想低通滤波器，计算一.矩形窗口法：低通滤波器的延时41以一个截止频率为ωc的线性相位理想低通滤波器则42则42

这是一个以为中心的偶对称的无限长非因果序列，如果截取一段n=0～N-1的hd(n)作为h(n)，则为保证所得到的是线性相位FIR滤波器，延时应为h(n)长度N的一半,即

43这是一个以为中心的偶对称的无限长非其中b.计算44其中b.计算44

c.计算设为窗口函数的频谱:

用幅度函数和相位函数来表示，则有

其线性相位部分则是表示延时一半长度，

45

c.计算45对频响起作用的是它的幅度函数46对频响起作用的是它的幅度函数46

理想频响也可以写成幅度函数和相位函数的表示形式

Hd(ejω)=Hd(ω)e-jωα其中幅度函数为

47理想频响也可以写成幅度函数和相位函数的表示形

两个信号时域的乘积对应于频域卷积，所以有48

两个信号时域的乘积对应于频域卷积，所以有48如果也以幅度函数和相位函数来表示H（ejω）,则实际FIR滤波器的幅度函数H（ω）为正好是理想滤波器幅度函数与窗函数幅度函数的卷积。49如果也以幅度函数和相位函数

矩形窗的卷积过程（P95的图4.5来说明）50矩形窗的卷积过程（P95的图4.5来说明）504个特殊频率点看卷积结果：（1）ω=0时,H(0)等于在[-ωc,ωc]内的积分面积因一般故H(0)近似为:在[-π,π]内的积分面积514个特殊频率点看卷积结果：（1）ω=0时,H(0)等于在（2）ω=ωc时，一半重叠，H(ωc)=0.5H(0);（

3）ω=ωc–2π/N时，第一旁瓣（负数）在通带外，出现正肩峰；

（4）ω=ωc+2π/N时，第一旁瓣（负数）在通带内，出现负肩峰。52（2）ω=ωc时，一半重叠，52窗口函数对理想特性的影响：①改变了理想频响的边沿特性，形成过渡带，宽为，等于WR(ω)的主瓣宽度。（决定于窗长）②过渡带两旁产生肩峰和余振（带内、带外起伏），取决于WR(ω)的旁瓣，旁瓣多，余振多；旁瓣相对值大，肩峰强，与N无关。（决定于窗口形状）

53窗口函数对理想特性的影响：53③N增加,过渡带宽减小,肩峰值不变。

因主瓣附近

其中x=Nω/2,所以N的改变不能改变主瓣与旁瓣的比例关系，只能改变WR（ω）的绝对值大小和起伏的密度，当N增加时，幅值变大，频率轴变密，而最大肩峰永远为8.95%，这种现象称为吉布斯（Gibbs）效应。

肩峰值的大小决定了滤波器通带内的平稳程度和阻带内的衰减，所以对滤波器的性能有很大的影响。54③N增加,过渡带宽减小,肩峰值不变。肩峰值的大

改变窗函数的形状，可改善滤波器的特性，窗函数有许多种，但要满足以下两点要求：①窗谱主瓣宽度要窄，以获得较陡的过渡带；②相对于主瓣幅度，旁瓣要尽可能小，使能量尽量集中在主瓣中，这样就可以减小肩峰和余振，以提高阻带衰减和通带平稳性。但实际上这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑制。55改变窗函数的形状，可改善滤波器的特性，窗函

几种常用的窗函数：1.矩形窗2.汉宁窗（升余弦窗）

56

几种常用的窗函数：56三部分矩形窗频谱相加，使旁瓣互相抵消，能量集中在主瓣，旁瓣大大减小，主瓣宽度增加1倍，为。

利用付氏变换的移位特性，汉宁窗频谱的幅度函数W（ω）可用矩形窗的幅度函数表示为：57三部分矩形窗频谱相加，使旁瓣互相抵消，能量集中在主瓣，旁瓣大5858

3.汉明窗（改进的升余弦窗）

它是对汉宁窗的改进，在主瓣宽度（对应第一零点的宽度）相同的情况下，旁瓣进一步减小，可使99.96%的能量集中在窗谱的主瓣内。

4.布莱克曼窗（三阶升余弦窗）

增加一个二次谐波余弦分量，可进一步降低旁瓣，但主瓣宽度进一步增加，为。增加N可减少过渡带。频谱的幅度函数为：

59

3.汉明窗（改进的升余弦窗）596060窗口函数的频谱N=51，A=20lg|W(ω)/W(0)|四种窗函数的比较61窗口函数的频谱N=51，A=20lg|W(ω)/W(06262

63

635.凯塞窗以上四种窗函数，都是以增加主瓣宽度为代价来降低旁瓣。凯塞窗则可自由选择主瓣宽度和旁瓣衰减。I0(x)是零阶修正贝塞尔函数，参数β可自由选择，决定主瓣宽度与旁瓣衰减。β越大，w(n)窗越窄，其频谱的主瓣变宽，旁瓣变小。一般取4<β<9。β=5.44接近汉明β=8.5接近布莱克曼β=0为矩形645.凯塞窗I0(x)是零阶修正贝塞尔函数，参数β可自由选择，

零阶贝塞尔函数

:过渡带宽

At:阻带最小衰减65

零阶贝塞尔函数:过渡带宽

At:阻带最图2凯塞窗函数图1零阶修正贝塞尔函数I0(x)x0166图2凯塞窗函数图1零阶修正贝塞尔函数I0(x)x016767用窗函数法设计FIR－DF的步骤：1）根据技术要求确定待求滤波器的单位取样响应hd(n)当较为复杂时，不容易由反付里叶变换求得。这时一般可用离散解里叶变换代替连续付里叶变换，求得近似值：令

hd(n)而当M>>N时，hM(n)≈hd(n)68用窗函数法设计FIR－DF的步骤：1）根据技术要求确定待求滤2）根据对过渡带及阻带衰减的要求，选择窗函数的形式，并估计窗口长度N。原则是：在保证阻带衰减满足要求的情况下，尽量选择主瓣窄的窗函数。3）计算DF的单位取样响应h（n），h（n）＝w（n）hd(n)4）验算技术指标是否满足要求。1）由定义3）卷积若不满足，要根据具体情况重复2）3）4）步，直到满足要求。692）根据对过渡带及阻带衰减的要求，选择窗函数的形式，并估计窗设计一个线性相位FIR低通DF，给定抽样频率通带截止频率为阻带截止频率为阻带衰减不小于－50DB，幅度特性如图所示。例：70设计一个线性相位FIR低通DF，例：7071717272737374747575幅度幅度归一化频率/归一化频率/（a）（b）凯塞窗设计举例76幅度幅度归一化频率/归一化频率/（a）（b）凯塞窗设计举7777787879798080§7.3频率采样法

工程上，常给定频域上的技术指标，所以采用频域设计更直接。一、基本思想使所设计的FIR数字滤波器的频率特性在某些离散频率点上的值准确地等于所需滤波器在这些频率点处的值，在其它频率处的特性则有较好的逼近。内插公式81§7.3频率采样法

工程上，常给定频域二.设计方法1）确定2）计算3）计算82二.设计方法1）确定82三、

约束条件

为了设计线性相位的FIR滤波器，采样值H（k）要满足一定的约束条件。前已指出,具有线性相位的FIR滤波器，其单位脉冲响应h(n)是实序列，且满足，由此得到的幅频和相频特性，就是对H(k)的约束。（表7.1.1）。例如，要设计第一类线性相位FIR滤波器，令N为奇数，h(n)偶对称，则幅度函数H（ω）应具有偶对称性：83三、

约束条件

为了设计线性相位的FIR滤波

令则必须满足偶对称性：而必须取为：

同样，若要设计第二种线性相位FIR滤波器，N为偶数，h(n)偶对称，由于幅度特性是奇对称的，84

令同样，若要设计第二种线性相位FIR滤波因此，Hk也必须满足奇对称性：

相位关系同上，

其它两种线性相位FIR数字滤波器的设计，同样也要满足幅度与相位的约束条件。

85因此，Hk也必须满足奇对称性：

相位关系同上，设用理想低通作为希望设计的滤波器，截止频率为，采样点数N，和的计算公式如下：N=奇数时，N=偶数时，，另外，对于高通和带阻，N只能为奇数86设用理想低通作为希望设计的滤波器，截止频率为，N1）时域分析：由频域采样定理知道，四、误差分析871）时域分析：由频域采样定理知道，四、误差分析872）频域分析：

由或H（z）。由上述设计过程得到的与的逼近程度，以及与H（k）的关系？由882）频域分析：由

令，则

89

89

单位圆上的频响为：这是一个内插公式。90

单位圆上的频响为：这是一个内插公式。90式中

为内插函数令则91式中91

内插公式表明：在每个采样点上，逼近误差为零，频响严格地与理想频响的采样值H(k)相等；在采样点之间，频响由各采样点的内插函数延伸迭加而形成，因而有一定的逼近误差，误差大小与理想频率响应的曲线形状有关，理想特性平滑，则误差小；反之，误差大。在理想频率响应的不连续点附近，会产生肩峰和波纹。N增大，则采样点变密，逼近误差减小。

92

内插公式表明：92图频率采样的响应93图频率采样的响应93例：设计一个FIR数字LP滤波器，其理想特性为

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

将代入内插公式，求H(ejω)：考虑到8<k<25时Hk=0，而其它k时，Hk=1,令k=33-n，则96

将代

97

979898

从图上可以看出，其过渡带宽为一个频率采样间隔2π/33，而最小阻带衰减略小于20dB。对大多数应用场合，阻带衰减如此小的滤波器是不能令人满意的。增大阻带衰减三种方法：1）加宽过渡带宽，以牺牲过渡带换取阻带衰减的增加。例如在本例中可在k=9和k=24处各增加一个过渡带采样点H9=H24=0.5，使过渡带宽增加到二个频率采样间隔4π/33，重新计算的H(ejω)见图4.12(c)，其阻带衰减增加到约-40dB。99从图上可以看出，其过渡带宽为一个频率采样间隔

2）过渡带的优化设计根据H(ejω)的表达式，H(ejω)是Hk的线性函数，因此还可以利用线性最优化的方法确定过渡带采样点的值，得到要求的滤波器的最佳逼近（而不是盲目地设定一个过渡带值）。例如，本例中可以用简单的梯度搜索法来选择H9、H24,使通带或阻带内的最大绝对误差最小化。

要求使阻带内最大绝对误差达到最小（也即最小衰减达到最大），可计算得H9=0.3904。对应的H(ejω)的幅频特性，比H9=0.5时的阻带衰减大大改善,衰减约-50dB。如果还要进一步改善阻带衰减，可以进一步加宽过渡区，添上第二个甚至第三个不等于0的频率取样值，当然也可用线性最优化求取这些取样值。100

2）过渡带的优化设计要求使阻带内最大绝

3）增大N如果要进一步增加阻带衰减，但又不增加过渡带宽，可增加采样点数N。例如，同样截止频率ωc=0.5π,以N=65采样，并在k=17和k=48插入由阻带衰减最优化计算得到的采样值H17=H48=0.5886,在k=18、47处插入经阻带衰减最优化计算获得的采样值H17=H48=0.1065,这时得到的H(ejω)，过渡带为6π/65，小于33点采样时插入一个过渡带采样点的过渡带宽，而阻带衰减增加了20多分贝，达-60dB以上，当然，代价是滤波器阶数增加，运算量增加。

101

3）增大N101小结：频率采样设计法优点：①

直接从频域进行设计，物理概念清楚，直观方便；②

适合于窄带滤波器设计，这时频率响应只有少数几个非零值。典型应用：用一串窄带滤波器组成多卜勒雷达接收机，覆盖不同的频段，多卜勒频偏可反映被测目标的运动速度；缺点：截止频率难以控制。因频率取样点都局限在2π/N的整数倍点上，所以在指定通带和阻带截止频率时，这种方法受到限制，比较死板。充分加大N，可以接近任何给定的频率，但计算量和复杂性增加。102小结：102§7.4FIR数字滤波器的最优化设计

前面介绍了FIR数字滤波器的两种逼近设计方法，即窗口法（时域逼近法）和频率采样法（频域逼近法），用这两种方法设计出的滤波器的频率特性都是在不同意义上对给定理想频率特性Hd(ejω)的逼近。说到逼近，就有一个逼近得好坏的问题，对“好”“坏”的恒量标准不同，也会得出不同的结论，我们前面讲过的窗口法和频率采样法都是先给出逼近方法，所需变量，然后再讨论其逼近特性，如果反过来要求在某种准则下设计滤波器各参数，以获取最优的结果，这就引出了最优化设计的概念，最优化设计一般需要大量的计算，所以一般需要依靠计算机进行辅助设计。

103§7.4FIR数字滤波器的最优化设计前

最优化设计的前提是最优准则的确定，在FIR滤波器最优化设计中，常用的准则有①最小均方误差准则②最大误差最小化准则。1)

均方误差最小化准则，若以E(ejω)表示逼近误差，则

那么均方误差为104最优化设计的前提是最优准则的确定，在FIR滤

均方误差最小准则就是选择一组时域采样值，以使均方误差，这一方法注重的是在整个-π～π频率区间内总误差的全局最小，但不能保证局部频率点的性能，有些频率点可能会有较大的误差，对于窗口法FIR滤波器设计，因采用有限项的h(n)逼近理想的hd(n)，所以其逼近误差为：如果采用矩形窗则有105

均方误差最小准则就是选择一组时域采样值，以可以证明，这是一个最小均方误差。所以，矩形窗窗口设计法是一个最小均方误差FIR设计，根据前面的讨论，我们知道其优点是过渡带较窄，缺点是局部点误差大,或者说误差分布不均匀。2)

最大误差最小化准则（也叫最佳一致逼近准则）表示为

其中F是根据要求预先给定的一个频率取值范围，可以是通带，也可以是阻带。最佳一致逼近即选择N个频率采样值（或时域h(n)值），在给定频带范围内使频响的最大逼近误差达到最小。也叫等波纹逼近。106可以证明，这是一个最小均方误差。表示为106

优点：可保证局部频率点的性能也是最优的，误差分布均匀，相同指标下，可用最少的阶数达到最佳化。

例如，我们提到的频率采样最优化设计，它是从已知的采样点数N、预定的一组频率取样和已知的一组可变的频率取样（即过渡带取样）出发，利用迭代法（或解析法）得到具有最小的阻带最大逼近误差（即最大的阻带最小衰减）的FIR滤波器。但它只是通过改变过渡带的一个或几个采样值来调整滤波器特性。如果所有频率采样值（或FIR时域序列h(m)）都可调整，显然，滤波器的性能可得到进一步提高。

107

优点：可保证局部频率点的性能也是最优的，误差分布均匀，相同窗函数法和频率采样法为使整个频域满足要求，平坦区域必将超过技术要求。所以，要引入一种新的设计方法，切比雪夫逼近法，它是一种等波纹逼近法，能使误差在整个频带均匀分布，对同样的技术指标，这种逼近法所需的滤波器阶数要低，而对于同样的滤波器阶数，这种逼近法的最大误差最小。108窗函数法和频率采样法为使整个频域满足要求，平坦区域108低通滤波器的误差分配109低通滤波器的误差分配109

切比雪夫最佳一致逼近如图，用等波纹逼近法设计滤波器需要确定五个参数：M、ωp、ωr、δ1、δ2按上图所示的误差容限设计低通滤波器，就是说要在通带0ωωp内以最大误差δ1逼近1，在阻带ωr

ω内以最大误差δ2逼近零。要同时确定上述五个参数较困难。常用的两种逼近方法：1）给定M、δ1、δ2，以ωp和ωr为变量。缺点：边界频率不能精确确定。2）给定M、ωp和ωr，以δ1和δ2为变量，通过迭代运算，使逼近误差δ1和δ2最小，并确定h(n)——切比雪夫最佳一致逼近。110

切比雪夫最佳一致逼近110

等波动逼近的低通滤波器pr

特点：能准确地指定通带和阻带边界频率。

111

等波动逼近的低通滤波器pr特点：能准确地指定通带和阻

一.误差函数定义逼近误差函数：

为所设计的滤波器与理想滤波器的幅频特性在通带和阻带内的误差值，是已知的权函数，在不同频带可取不同的值，所要设计的滤波器的幅频特性理想滤波器的幅频特性112

一.误差函数为所设计的滤波器与理想滤波器的幅例如，希望在固定M,p,r的情况下逼近一个低通滤波器，这时有对于表4.1中的第一种滤波器，113例如，希望在固定M,p,r的情况下逼近一个低通于是切比雪夫逼近问题变为，寻求一组系数使逼近误差的最大值达到最小，即

给定后等效于求最小。114于是切比雪夫逼近问题变为，寻求一组系数

二.交替定理（最佳逼近定理）令F表示闭区间的任意闭子集，为了使在F上唯一最佳地逼近于,其充分必要条件是误差函数在F上至少应有（M+2）次“交替”，即其中，且属于F。

1）至少有M+2个极值，且极值正负相间，具有等波纹的性质，2）由于是常数，所以的极值也就是的极值。115

二.交替定理（最佳逼近定理）115

借助于低通滤波器的设计，可以直观地解释这个定理。这时，闭子集F包括区间和。因为滤波器频响是逐段恒定的，所以对应于误差函数各峰值点的频率同样也对应于恰好满足误差容限时的频率。根据前面的讨论，在开区间内至多有M-1个极值，此外，根据通带和阻带的定义，令的约束条件为，，再加上和π处的极值，误差曲线最多有M+1个极值频率（交替）满足定理。116借助于低通滤波器的设计，可以直观地解释这个定理。逼近方法：固定k、M、和，以作为参变量。按照交错定理，如果F上的M+2个极值点频率已知，则由（1）式可得到M+2个方程：为极值点频率对应的误差函数值117逼近方法：固定k、M、和，以

注意：极值点频率必须位于和区间内。由于和固定，因而和必为这些极值频率中的一个，设，则应有求解上述方程组可得到全部系数。问题：1）实际情况下，M+2个极值点频率未知；2）直接求解上述非线性方程组比较困难。雷米兹（Remez）算法给出了求解切比雪夫最佳一致逼近问题的方法。118

注意：极值点频率必须位于雷米兹交替算法119雷米兹交替算法119三.雷米兹（Remez）算法1）在频率子集F上均匀等间隔地选取M+2个极值点频率

120三.雷米兹（Remez）算法1）在频率子集F上均匀等间隔2）由求和利用重心形式的拉格朗日插值公式，其中如在频带F上，对所有频率都有，则为所求，即为极值点频率。1212）由3）对上次确定的极值点频率中的每一点，在其附近检查是否在某一频率处有，如有，则以该频率点作为新的局部极值点。对M+2个极值点频率依次进行检查，得到一组新的极值点频率。重复步骤1）、2），求出，完成一次迭代。重复上述步骤，直到的值改变很小，迭代结束，这个即为所求的最小值。由最后一组极值点频率求出，反变换得到,完成设计。优点：可准确确定；逼近误差均匀分布，相同指标下，滤波器所需阶数低。

1223）对上次确定的极值点频率有一些估算公式可用于决定最佳滤波器长度N：对于窄带低通滤波器，对滤波器长度N起主要作用：123有一些估算公式可用于决定最佳滤波器长度N：对于窄带低通滤波器1241241251250200400600800100012001400160018002000-80-70-60-50-40-30-20-10010频率/Hz·幅度/dBRemez交替法设计举例126020040060080010001200140016001§7.5IIR与FIR数字滤器的比较

FIRIIR设计方法一般无解析的设计公式，要借助计算机程序完成利用AF的成果，可简单、有效地完成设计设计结果可得到幅频特性和线性相位（最大优点）只能得到幅频特性，相频特性未知（一大缺点），如需要线性相位，须用全通网络校准，但增加滤波器阶数和复杂性稳定性极点全部在原点（永远稳定）无稳定性问题有稳定性问题

阶数高

结构非递归递归系统运算误差一般无反馈，运算误差小有反馈，由于运算中的四舍五入会产生极限环快速算法可用FFT实现，减少运算量无快速运算方法低127§7.5IIR与FIR数字滤器的比较

FIRIIR设计作业P.2351，2，4，13，15，16128作业P.235128第七章有限长单位脉冲响应（FIR） 滤波器的设计方法

序言§7.1线性相位FIR数字滤波器的特性§7.2窗口设计法（时间窗口法）

§7.3频率取样法

§7.4FIR数字滤波器的最优化设计§7.5IIR与FIR数字滤器的比较129第七章有限长单位脉冲响应（FIR）

FIR数字滤波器的差分方程描述①

对应的系统函数

因为它是一种线性时不变系统，可用卷积和形式表示

③

比较①、③得：

序言

130

序言2

FIR数字滤波器的特点(与IIR数字滤波器比较)：

优点：（1）很容易获得严格的线性相位，避免被处理的信号产生相位失真，这一特点在宽频带信号处理、阵列信号处理、数据传输等系统中非常重要（2）永远稳定，无稳定性问题

131

FIR数字滤波器的特点(与IIR数字滤波器比较)：3（3）任何一个非因果的有限长序列，总可以通过一定的延时，转变为因果序列，所以因果性总是满足；（4）无反馈运算，运算误差小。132（3）任何一个非因果的有限长序列，总可以通过一定的延时，转缺点：（1）因为无极点，要获得好的过渡带特性，需以较高的阶数为代价；（2）无法利用模拟滤波器的设计结果，一般无解析设计公式，要借助计算机辅助设计程序完成。133缺点：5§7.1线性相位FIR数字滤波器的特性

1线性相位的条件线性相位意味着一个系统的相频特性是频率的线性函数，即式中为常数，此时通过这一系统的各频率分量的时延为一相同的常数，系统的群时延为

134§7.1线性相位FIR数字滤波器的特性1线性相位的条FIR滤波器的DTFT为式中H(ω)是正或负的实函数。等式中间和等式右边的实部与虚部应当各自相等，同样实部与虚部的比值应当相等:135FIR滤波器的DTFT为式中H(ω)是正或负的实函数。等式将上式两边交叉相乘，再将等式右边各项移到左边，应用三角函数的恒等关系满足上式的条件是136将上式两边交叉相乘，再将等式右边各项移到左边，应用三角函数的另外一种情况是，除了上述的线性相位外，还有一附加的相位，即利用类似的关系，可以得出新的解答为137另外一种情况是，除了上述的线性相位外，还有一附加的相位，即

偶对称

奇对称138偶对称奇对称10

分四种情况(1)．h(n)偶对称，N为奇数h(n)=h(N-1-n)2线性相位FIR滤波器的幅度特性139

分四种情况2线性相位FIR滤波器的幅度特性1114012令，则令则由于偶对称，因此对这些频率也呈偶对称。141令，则令则

(2)．h(n)偶对称，N为偶数h(n)=h(N-1-n)142(2)．h(n)偶对称，N为偶数h(n)=h(N令，则143令，则15

或写为：

由于奇对称，所以对也为奇对称，且由于时，处必有一零点，因此这种情况不能用于设计时的滤波器，如高通、带阻滤波器。144

或写为：由于(3).h(n)奇对称，N为奇数，h(n)=-h(N-1-n)145(3).h(n)奇对称，N为奇数，h(n)=-h(N-1

令n=m+(N-1)/2，得：

146

18

所以

由于点呈奇对称，所以对这些点也奇对称。由于时，相当于H（z）在处有两个零点，不能用于的滤波器设计，故不能用作低通、高通和带阻滤波器的设计。147

所以由于(4).h(n)奇对称，N为偶数令148(4).h(n)奇对称，N为偶数令20

由于在ω=0,２π处为零，所以H(ω)在ω=0,2π处为零，即H(z)在z=1上有零点，并对ω=0,2π呈奇对称。对ω=π呈偶对称。

149

由于四种线性相位FIR滤波器关于0、π、2π偶对称关于π奇对称H(π)=0关于0、π、2π奇对称关于0、2π奇对称π偶对称150四种线性相位FIR滤波器关于0、π、2π偶对称关于π奇对称H四种线性相位FIRDF特性:第一种情况，偶、奇，四种滤波器都可设计。第二种情况，偶、偶，可设计低、带通滤波器，不能设计

高通和带阻。第三种情况，奇、奇，只能设计带通滤波器，其它滤波器都不能设计。第四种情况，奇、偶，可设计高通、带通滤波器，不能设计低通和带阻。151四种线性相位FIRDF特性:第二种情况，偶、偶，可设计低、例1N=5,h(0)=h(1)=h(3)=h(4)=-1/2,h(2)=2，求幅度函数H(ω)。解Ｎ为奇数并且h(n)满足偶对称关系a(0)=h(2)=2a(1)=2h(1)=-1a(2)=2h(0)=-1H(ω)=2

-cosω-cos2ω=2-(cosω+cos2ω)152例1N=5,h(0)=h(1)=h(3)15325小结：

四种FIR数字滤波器的相位特性只取决于h(n)的对称性，而与h(n)的值无关。幅度特性取决于h(n)。设计FIR数字滤波器时，在保证h(n)对称的条件下，只要完成幅度特性的逼近即可。注意：当H(ω)用│H(ω)│表示时，当H(ω)为奇对称时，其相频特性中还应加一个固定相移π。154小结：26

3线性相位FIR滤波器的零点特性

155

3线性相位FIR滤波器的零点特性27由该式可看出，若z=zi是H(z)的零点，则z=zi-1也一定是H（z）的零点。由于h(n)是实数，H(z)的零点还必须共轭或对，所以z=zi*及z=1/z*也必是零点。所以线性相位滤波器的零点必须是互为倒数的共轭对，即成四出现，这种共轭对共有四种156由该式可看出，若z=zi是H(z)的零点，则z=zi-1也一可能的情况：①既不在单位园上，也不在实轴上，有四个互为倒数的两组共轭对，ziz*i1/zi1/z*i图4.2(a)②在单位圆上，但不在实轴上，因倒数就是自己的共轭，所以有一对共轭零点，zi,z*i图4.2（b）③不在单位圆上，但在实轴上，是实数,共轭就是自己，所以有一对互为倒数的零点,zi,1/zi图4.2（c）④又在单位圆上，又在实轴上，共轭和倒数都合为一点，所以成单出现，只有两种可能，zi=1或zi=-1图4.2(d),p92157可能的情况：2915830我们从幅度响应的讨论中已经知道，对于第二种FIR滤波器（h(n)偶对称，N为偶数），，

即是的零点，既在单位圆，又在实轴，所以必有单根同样道理，对于第三种FIR滤波器，h(n)奇对称，N为奇数，因所以z=1，z=-1都是H（z）的单根；对于第四种滤波器，h(n)奇对称，N为偶数，H（O）=0，所以z=1是H（z）的单根。所以，h(n)奇对称→H(0)=0N为偶数→H(π)=0159我们从幅度响应的讨论中已经知道，对于第二种FIR滤波器（h(线性相位滤波器是FIR滤波器中最重要的一种，应用最广。实际使用时应根据需用选择其合适类型，并在设计时遵循其约束条件。160线性相位滤波器是FIR滤波器中最重要的一种，应用最广。实际使16133§7.2窗口设计法（时域）

如果希望得到的滤波器的理想频率响应为，那么FIR滤波器的设计就在于寻找一个传递函数去逼近逼近方法有三种：窗口设计法（时域逼近）频率采样法（频域逼近）最优化设计（等波纹逼近）

162§7.2窗口设计法（时域）如果希望得到的滤波时间窗口设计法是从单位脉冲响应序列着手，使h(n)逼近理想的单位脉冲响应序列hd(n)

。hd(n)可以从理想频响的傅立叶反变换获得163时间窗口设计法是从单位脉冲响应序列着手，使h(n)逼近理想的但一般来说，理想频响是分段恒定，在边界频率处有突变点，所以，这样得到的理想单位脉冲响应hd(n)往往都是无限长序列，而且是非因果的。但FIR的h(n)是有限长的，问题是怎样用一个有限长的序列去近似无限长的hd(n)。

164但一般来说，理想频响是分段恒定，在最简单的办法是直接截取一段hd(n)代替h(n)。这种截取可以形象地想象为h(n)是通过一个“窗口”所看到的一段hd(n)，因此，h(n)也可表达为hd(n)和一个“窗函数”的乘积，即h(n)=w(n)hd(n)在这里窗口函数就是矩形脉冲函数RN(n)，当然以后我们还可看到，为了改善设计滤波器的特性，窗函数还可以有其它的形式，相当于在矩形窗内对hd(n)作一定的加权处理。165最简单的办法是直接截取一段hd(n)代替h(n)。这种16638

设计步骤：1）由定义3）卷积插值167

设计步骤：1）由定义3）卷积插值3916840

以一个截止频率为ωc的线性相位理想低通滤波器为例,讨论FIR的设计问题。a.对于给定的理想低通滤波器，计算一.矩形窗口法：低通滤波器的延时169以一个截止频率为ωc的线性相位理想低通滤波器则170则42

这是一个以为中心的偶对称的无限长非因果序列，如果截取一段n=0～N-1的hd(n)作为h(n)，则为保证所得到的是线性相位FIR滤波器，延时应为h(n)长度N的一半,即

171这是一个以为中心的偶对称的无限长非其中b.计算172其中b.计算44

c.计算设为窗口函数的频谱:

用幅度函数和相位函数来表示，则有

其线性相位部分则是表示延时一半长度，

173

c.计算45对频响起作用的是它的幅度函数174对频响起作用的是它的幅度函数46

理想频响也可以写成幅度函数和相位函数的表示形式

Hd(ejω)=Hd(ω)e-jωα其中幅度函数为

175理想频响也可以写成幅度函数和相位函数的表示形

两个信号时域的乘积对应于频域卷积，所以有176

两个信号时域的乘积对应于频域卷积，所以有48如果也以幅度函数和相位函数来表示H（ejω）,则实际FIR滤波器的幅度函数H（ω）为正好是理想滤波器幅度函数与窗函数幅度函数的卷积。177如果也以幅度函数和相位函数

矩形窗的卷积过程（P95的图4.5来说明）178矩形窗的卷积过程（P95的图4.5来说明）504个特殊频率点看卷积结果：（1）ω=0时,H(0)等于在[-ωc,ωc]内的积分面积因一般故H(0)近似为:在[-π,π]内的积分面积1794个特殊频率点看卷积结果：（1）ω=0时,H(0)等于在（2）ω=ωc时，一半重叠，H(ωc)=0.5H(0);（

3）ω=ωc–2π/N时，第一旁瓣（负数）在通带外，出现正肩峰；

（4）ω=ωc+2π/N时，第一旁瓣（负数）在通带内，出现负肩峰。180（2）ω=ωc时，一半重叠，52窗口函数对理想特性的影响：①改变了理想频响的边沿特性，形成过渡带，宽为，等于WR(ω)的主瓣宽度。（决定于窗长）②过渡带两旁产生肩峰和余振（带内、带外起伏），取决于WR(ω)的旁瓣，旁瓣多，余振多；旁瓣相对值大，肩峰强，与N无关。（决定于窗口形状）

181窗口函数对理想特性的影响：53③N增加,过渡带宽减小,肩峰值不变。

因主瓣附近

其中x=Nω/2,所以N的改变不能改变主瓣与旁瓣的比例关系，只能改变WR（ω）的绝对值大小和起伏的密度，当N增加时，幅值变大，频率轴变密，而最大肩峰永远为8.95%，这种现象称为吉布斯（Gibbs）效应。

肩峰值的大小决定了滤波器通带内的平稳程度和阻带内的衰减，所以对滤波器的性能有很大的影响。182③N增加,过渡带宽减小,肩峰值不变。肩峰值的大

改变窗函数的形状，可改善滤波器的特性，窗函数有许多种，但要满足以下两点要求：①窗谱主瓣宽度要窄，以获得较陡的过渡带；②相对于主瓣幅度，旁瓣要尽可能小，使能量尽量集中在主瓣中，这样就可以减小肩峰和余振，以提高阻带衰减和通带平稳性。但实际上这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑制。183改变窗函数的形状，可改善滤波器的特性，窗函

几种常用的窗函数：1.矩形窗2.汉宁窗（升余弦窗）

184

几种常用的窗函数：56三部分矩形窗频谱相加，使旁瓣互相抵消，能量集中在主瓣，旁瓣大大减小，主瓣宽度增加1倍，为。

利用付氏变换的移位特性，汉宁窗频谱的幅度函数W（ω）可用矩形窗的幅度函数表示为：185三部分矩形窗频谱相加，使旁瓣互相抵消，能量集中在主瓣，旁瓣大18658

3.汉明窗（改进的升余弦窗）

它是对汉宁窗的改进，在主瓣宽度（对应第一零点的宽度）相同的情况下，旁瓣进一步减小，可使99.96%的能量集中在窗谱的主瓣内。

4.布莱克曼窗（三阶升余弦窗）

增加一个二次谐波余弦分量，可进一步降低旁瓣，但主瓣宽度进一步增加，为。增加N可减少过渡带。频谱的幅度函数为：

187

3.汉明窗（改进的升余弦窗）5918860窗口函数的频谱N=51，A=20lg|W(ω)/W(0)|四种窗函数的比较189窗口函数的频谱N=51，A=20lg|W(ω)/W(019062

191

635.凯塞窗以上四种窗函数，都是以增加主瓣宽度为代价来降低旁瓣。凯塞窗则可自由选择主瓣宽度和旁瓣衰减。I0(x)是零阶修正贝塞尔函数，参数β可自由选择，决定主瓣宽度与旁瓣衰减。β越大，w(n)窗越窄，其频谱的主瓣变宽，旁瓣变小。一般取4<β<9。β=5.44接近汉明β=8.5接近布莱克曼β=0为矩形1925.凯塞窗I0(x)是零阶修正贝塞尔函数，参数β可自由选择，

零阶贝塞尔函数

:过渡带宽

At:阻带最小衰减193

零阶贝塞尔函数:过渡带宽

At:阻带最图2凯塞窗函数图1零阶修正贝塞尔函数I0(x)x01194图2凯塞窗函数图1零阶修正贝塞尔函数I0(x)x0119567用窗函数法设计FIR－DF的步骤：1）根据技术要求确定待求滤波器的单位取样响应hd(n)当较为复杂时，不容易由反付里叶变换求得。这时一般可用离散解里叶变换代替连续付里叶变换，求得近似值：令

hd(n)而当M>>N时，hM(n)≈hd(n)196用窗函数法设计FIR－DF的步骤：1）根据技术要求确定待求滤2）根据对过渡带及阻带衰减的要求，选择窗函数的形式，并估计窗口长度N。原则是：在保证阻带衰减满足要求的情况下，尽量选择主瓣窄的窗函数。3）计算DF的单位取样响应h（n），h（n）＝w（n）hd(n)4）验算技术指标是否满足要求。1）由定义3）卷积若不满足，要根据具体情况重复2）3）4）步，直到满足要求。1972）根据对过渡带及阻带衰减的要求，选择窗函数的形式，并估计窗设计一个线性相位FIR低通DF，给定抽样频率通带截止频率为阻带截止频率为阻带衰减不小于－50DB，幅度特性如图所示。例：198