第三章 FIR数字滤波器设计

1、1第三章第三章 FIR数字滤波器的设计数字滤波器的设计3.1 数字滤波器设计概述数字滤波器设计概述3.2 窗口设计法（时间窗口法）窗口设计法（时间窗口法） 3.3 频率取样法频率取样法3.4 等波纹优化设计等波纹优化设计23.1 数字滤波器设计概述数字滤波器设计概述滤波器设计滤波器设计：根据给定滤波器的频率特性，求根据给定滤波器的频率特性，求得满足该特性的传输函数。得满足该特性的传输函数。数字滤波器数字滤波器0)()()(nnjezjenhzHeHj任何线性时不变任何线性时不变离散系统离散系统都是一个数字滤波器（简都是一个数字滤波器（简称称DF）。设计数字滤波器的）。设计数字滤波器的任务任务就

2、是寻求一个就是寻求一个因因果稳定的线性时不变系统果稳定的线性时不变系统，使其系统函数，使其系统函数H(z)具有具有指定的指定的频率特性频率特性。33.1 数字滤波器设计概述数字滤波器设计概述 DFDF的幅度特性性能要求（低通为例）的幅度特性性能要求（低通为例）)(H0psp1p1s通带通带过渡带过渡带阻带阻带p：通带波纹：通带波纹s：阻带波纹：阻带波纹：通带截止频率：通带截止频率p：阻带截止频率：阻带截止频率s：过渡带：过渡带43.1 数字滤波器设计概述数字滤波器设计概述 DFDF的相位特性性能要求的相位特性性能要求从信号不失真角度，其相位特性通常必须满足以下两个从信号不失真角度，其相位特性通

3、常必须满足以下两个要求之一：要求之一： 相位线性。即：相位线性。即： 具有群恒时延特性。即：具有群恒时延特性。即：)(常数)()(dd其中 称为时延常数时延常数)(53.1 数字滤波器设计概述数字滤波器设计概述数字滤波器设计采用角频率，如何与实际信号频率数字滤波器设计采用角频率，如何与实际信号频率对应？对应？角频率角频率w w，采样频率，采样频率f fs s , ,实际信号频率实际信号频率f f的转换关系的转换关系2/swff 采样频率的角频率为采样频率的角频率为263.1 数字滤波器设计概述数字滤波器设计概述 例如：设计一个带通滤波器，通带可通过例如：设计一个带通滤波器，通带可通过中心频率中

4、心频率为为3MHz3MHz，带宽为带宽为1.2MHz1.2MHz的带通信号的带通信号，系统采样频率为，系统采样频率为24MHz24MHz。带通滤波器通带。带通滤波器通带3dB3dB衰减处距离带通信号通带边缘衰减处距离带通信号通带边缘100kHz100kHz，阻带衰减，阻带衰减40dB40dB处距离带通信号通带边缘处距离带通信号通带边缘300kHz300kHz。根据数字带通滤波器的角频率技术指标：根据数字带通滤波器的角频率技术指标： 通带截止频率（以通带截止频率（以3dB3dB衰减为准）衰减为准） 阻带截止频率（以阻带截止频率（以40dB40dB衰减为准）衰减为准） 带通滤波器中心频率带通滤波器

5、中心频率画出数字带通滤波器的角频率幅频响应频谱示意图画出数字带通滤波器的角频率幅频响应频谱示意图 703fM1.2BMHz100pfKHZ3000.3cfKHZMHZ24sfMHZ102 *(/2)/0.1917ppswfBff202 *(/2)/0.31ppswfBff102 *(/2)/0.175ccswfBff202 *(/2)/0.325ccswfBff002 */0.25swff解：（解：（1）3.1 数字滤波器设计概述数字滤波器设计概述8（2 2）数字带通滤波器的角频率频谱示意图为：）数字带通滤波器的角频率频谱示意图为：3.1 数字滤波器设计概述数字滤波器设计概述93.1 数字滤波

6、器设计概述数字滤波器设计概述设计设计滤波器理想频响滤波器理想频响包含：低通，高通，带通，带阻包含：低通，高通，带通，带阻对应到数字滤波器的频谱为：对应到数字滤波器的频谱为：10低通滤波器低通滤波器1 ()0 cjwcwwH eww幅频响应幅频响应相频响应相频响应awcw称为截止频率称为截止频率wcw()jwH ecw3.1 数字滤波器设计概述数字滤波器设计概述 1( )sin()2 ()cccwjawjwncwwnah needww nanana11低通滤波器低通滤波器3.1 数字滤波器设计概述数字滤波器设计概述 1( )sin()2 ()cccwjawjwncwwnah neednw nan

7、anaa12高通滤波器高通滤波器0 ()1 cjwcwwH eww幅频响应幅频响应相频响应相频响应cw称为通带开始频率称为通带开始频率wcw()jwH ecw3.1 数字滤波器设计概述数字滤波器设计概述aw1 ( )sin() 0()ccwnah nw nanna133.1 数字滤波器设计概述数字滤波器设计概述高通高通低通低通时域除时域除n=a点外点外符号相符号相反反频域上频域上相当于相当于低通平低通平移移14带通滤波器带通滤波器11220 01 ()0 cjwcccwwwwwH eww幅频响应幅频响应相频响应相频响应12,ccww称为通带截止频率称为通带截止频率w1cw()jwH e1cw2

8、cw2cw3.1 数字滤波器设计概述数字滤波器设计概述aw2121 ( )sin()sin() ()ccccwwnah nwnawnanana可看成是两个低通滤波器时域相减得到可看成是两个低通滤波器时域相减得到153.1 数字滤波器设计概述数字滤波器设计概述低通低通wc1低通低通wc2带通带通wc1wc216带阻滤波器带阻滤波器11221 00 ()1 cjwcccwwwwwH eww幅频响应幅频响应相频响应相频响应cw称为通带开始频率称为通带开始频率w1cw()jwH e1cw2cw2cw3.1 数字滤波器设计概述数字滤波器设计概述aw21211 ( )sin()sin() ()ccccww

9、nah nwnawnanana与带通滤波器互补与带通滤波器互补17（二）带阻滤波器（二）带阻滤波器3.1 数字滤波器设计概述数字滤波器设计概述带通滤波器带通滤波器带阻滤波器带阻滤波器183.1 数字滤波器设计概述数字滤波器设计概述无论哪种数字滤波器，设计过程一般包括以下三个基本问题：无论哪种数字滤波器，设计过程一般包括以下三个基本问题： 根据实际要求根据实际要求确定数字滤波器性能指标确定数字滤波器性能指标； 用一个用一个因果稳定的系统函数去逼近因果稳定的系统函数去逼近这个理想性能指标；这个理想性能指标； 用一个用一个有限精度的运算去实现有限精度的运算去实现这个传输函数。这个传输函数。 现在要讨

10、论的是问题现在要讨论的是问题2 2数字滤波器数字滤波器IIR数字滤波器数字滤波器：其系统函数通常可表示：其系统函数通常可表示成关于成关于 的有理分式的有理分式1zFIR数字滤波器数字滤波器：其系统函数通常可表示：其系统函数通常可表示成关于成关于 的多项式的多项式1z19FIRFIR数字滤波器的特点：数字滤波器的特点： 优点优点 ：（：（1 1）很容易获得严格的线性相位，很容易获得严格的线性相位，避免被处理避免被处理 的信号产生相位失真，这一特点在的信号产生相位失真，这一特点在 宽频带信号处理、阵宽频带信号处理、阵 列信号列信号处理、数据传输等系统中非常重要；处理、数据传输等系统中非常重要； （

11、2 2）可得到多带幅频特性可得到多带幅频特性； （3 3）极点全部在原点（永远稳定），极点全部在原点（永远稳定），无稳定性问题；无稳定性问题； （4 4）任何一个非因果的有限长序列，总可以通过一）任何一个非因果的有限长序列，总可以通过一 定的延时，转变为因果序列，定的延时，转变为因果序列， 所以所以因果性总是满足；因果性总是满足； （5 5）无反馈运算，）无反馈运算，运算误差小运算误差小。3.1 数字滤波器设计概述数字滤波器设计概述20缺点：（缺点：（1 1）因为无极点，要获得好的过渡带特性，需以）因为无极点，要获得好的过渡带特性，需以较较 高的阶数高的阶数为代价；为代价； （2 2）无法利用

12、模拟滤波器的设计结果无法利用模拟滤波器的设计结果，一般无解析，一般无解析设计公式，要借助计算机辅助设计程序完成。设计公式，要借助计算机辅助设计程序完成。3.1 数字滤波器设计概述数字滤波器设计概述213.1 数字滤波器设计概述数字滤波器设计概述线性相位线性相位FIRFIR数字滤波器数字滤波器102121cos)(| )()(NnNjezjnNnhezHeHj11/220()( )|1( )sin2jjz eNNjjnH eH zNeh nn223.1 数字滤波器设计概述数字滤波器设计概述下表给出了上述4种类型的线性相位滤波器的相位响应相位响应、时时域幅度响应域幅度响应和频域幅度响应频域幅度响应

13、的示意图。23四种线性相位四种线性相位FIR DFFIR DF特性，特性，第一种情况第一种情况 ，偶对称、奇数长度，四种滤波器都可设计偶对称、奇数长度，四种滤波器都可设计。第二种情况，偶对称、偶数长度，可设计低、带通滤波器第二种情况，偶对称、偶数长度，可设计低、带通滤波器，不能设计高通和带阻。在，不能设计高通和带阻。在处有零点处有零点第三种情况，奇对称、奇数长度，只能设计带通滤波器，第三种情况，奇对称、奇数长度，只能设计带通滤波器，其它滤波器都不能设计。在其它滤波器都不能设计。在0 0，2 处都有零点处都有零点第四种情况，奇对称、偶数长度，可设计高通、带通滤波第四种情况，奇对称、偶数长度，可设

14、计高通、带通滤波器，不能设计低通和带阻。器，不能设计低通和带阻。0 0， 2 2 处有零点处有零点3.1 数字滤波器设计概述数字滤波器设计概述243.1 数字滤波器设计概述数字滤波器设计概述例例1 1 N N=5, =5, h h (0) = (0) = h h (1) = (1) = h h (3) = (3) = h h (4) = -1/2, (4) = -1/2, h h (2) = 2 (2) = 2，求频响函数，求频响函数H H () ()，画出幅频响应，画出幅频响应253.1 数字滤波器设计概述数字滤波器设计概述例例1 1 N N=5, =5, h h (0) = (0) = h

15、 h (1) = (1) = h h (3) = (3) = h h (4) = -1/2, (4) = -1/2, h h (2) = 2 (2) = 2，求频响函数，求频响函数H H () ()，画出幅频响应，画出幅频响应H () = e-j2w (2 - cos- cos2)26 四种四种FIRFIR数字滤波器的数字滤波器的相位特性相位特性只取决于只取决于h(n)h(n)的对称性的对称性，而与而与h(n)h(n)的值无关。的值无关。幅度特性取决于幅度特性取决于h(n)h(n)。设计设计FIRFIR数字滤波器数字滤波器时，在保证时，在保证h(n)h(n)对称的条件下，对称的条件下，只只要完

16、成幅度特性的逼近即可。要完成幅度特性的逼近即可。注意：当注意：当H()H()用用H()H()表示时，表示时，H()H()为为奇对称奇对称时，时，其相频特性中还应加一个其相频特性中还应加一个固定相移固定相移。3.1 数字滤波器设计概述数字滤波器设计概述结论：结论：273.2 窗口设计法（时域）窗口设计法（时域）FIRFIR滤波器的设计就在于寻找一个传递函数滤波器的设计就在于寻找一个传递函数 去逼近去逼近 ，逼近方法有三种：，逼近方法有三种： 窗口设计法（时域逼近）窗口设计法（时域逼近） 频率采样法（频域逼近）频率采样法（频域逼近） 最优化设计（等波纹逼近）最优化设计（等波纹逼近） 时间窗口设计法

17、时间窗口设计法是从单位脉冲响应序列着手，使是从单位脉冲响应序列着手，使h(n)h(n)逼近理想的单位脉冲响应序列逼近理想的单位脉冲响应序列h hd d(n)(n)。我们知道。我们知道h hd d(n)(n)可以从可以从理想频响通过付氏反变换获得理想频响通过付氏反变换获得 10)(NnjnjenheH)(jdeH221)(onjjdddeeHnh28理想频响理想频响 是分段恒定，在边界频率处有突变是分段恒定，在边界频率处有突变点点理想单位脉冲响应理想单位脉冲响应h hd d(n)(n)往往都是往往都是无限长序列无限长序列，而且，而且是非因果是非因果办法办法: :直接截取一段直接截取一段 h hd

18、 d(n) (n) 代替代替 h(n)h(n) )(jdeH 3.2 窗口设计法（时域）窗口设计法（时域）29 3.2 窗口设计法（时域）窗口设计法（时域）h(n)h(n)是通过一个是通过一个“窗口窗口”所看到所看到的一段的一段h hd d(n) (n) 表达为表达为h(n)h(n)和一和一个个“窗函数窗函数”的的乘积，即乘积，即 h(n)=w(n) hh(n)=w(n) hd d(n) (n) 在这里窗口函数在这里窗口函数就是矩形脉冲函就是矩形脉冲函数数R RN N（n n）30设计步骤：设计步骤：)()()()(nwnhnheHddjd)()(nheHj)()(nheHdjd设10)(Nn

19、jnjenheH由定义3.2 窗口设计法（时域）窗口设计法（时域）该式为有限项，该式为有限项，NN越大，误差越小。但对越大，误差越小。但对矩形窗截取还存在矩形窗截取还存在“吉吉布斯（布斯（GibbsGibbs）效应）效应”，这将使滤波器的特性很这将使滤波器的特性很差。差。31一一. .矩形窗口法矩形窗口法)(jdeH则 以一个截止频率为 c的线性相位理想低通滤波器线性相位理想低通滤波器为例,讨论FIR的设计问题。a. 对于给定的理想低通滤波器 ，计算：低通滤波器的延时：低通滤波器的延时)(nhd3.2 窗口设计法（时域）窗口设计法（时域） ( )sin() ()cdcwnah nw nanan

20、a为保证所得到的是线性相位FIR滤波器，延时延时 a应为应为h(n)长度长度N的一半的一半,即(1)/2N32为其它值nNnonhnwnhnhdRd01)()()()(其中)()(nRnwNRb.计算)(nh3.2 窗口设计法（时域）窗口设计法（时域）33c.计算 。设 为窗口函数的频谱: 用幅度函数和相位函数来表示，则有 其线性相位部分 则是表示延时一半长度 ， )(jeWnNnjjNnjnjRjeeeenweW1011)()()2/sin()2/sin(21NeNjjRjeWeW)()(je2/ ) 1( N)(jeH)(*)()(jRjdjeWeHeH3.2 窗口设计法（时域）窗口设计法

21、（时域）34两个信号时域的乘积对应于频域卷积，所以有deWeHeWeHeHjRjdjRjdj)(21)(*)()()(3.2 窗口设计法（时域）窗口设计法（时域）()1( )()21( )()2jjdRjdRHeWedeHWd 相频响应未发生变化相频响应未发生变化35则实际FIR滤波器的幅度函数H（）为正好是理想滤波器幅度函数与窗函数幅度函数的卷积理想滤波器幅度函数与窗函数幅度函数的卷积。jjeHeH)()(dWHHRd)()(21)(3.2 窗口设计法（时域）窗口设计法（时域）363.2 窗口设计法（时域）窗口设计法（时域）37)(RWNc2)(RW4 4个特殊频率点看卷积结果：个特殊频率点

22、看卷积结果：（1 1）=0=0时时,H(0),H(0)等于等于在在-c, c-c, c内的积分面积内的积分面积因一般因一般故故H(0) H(0) 近似为近似为在在-, -, 内的积分面积内的积分面积3.2 窗口设计法（时域）窗口设计法（时域）(2)=(2)=c c时，一半重叠，时，一半重叠，H(H(c c)=0.5 H(0);)=0.5 H(0);(3)=(3)=c c2/N2/N时，第一旁瓣（负数）在通带外，时，第一旁瓣（负数）在通带外，出现出现正肩峰正肩峰； (4)=(4)=c c+2/N +2/N 时，第一旁瓣（负数）在通带内，时，第一旁瓣（负数）在通带内，出现出现负肩峰负肩峰。38窗口

23、函数对理想特性的影响： 改变了理想频响的边沿特性，形成过渡带，宽为 ， 等于WR()的主瓣宽度。（决定于窗长） 过渡带两旁产生肩峰和余振（带内、带外起伏），取决于 WR()的旁瓣，旁瓣多，余振多；旁瓣相对值大，肩峰强 ，与 N无关。（决定于窗口形状） N增加,过渡带宽减小,肩峰值不变。 因主瓣附近 其中x=N/2,所以N的改变不能改变主瓣与旁瓣的比例关系，只能改变WR（）的绝对值大小和起伏的密度，当N增加时，幅值变大，频率轴变密，而最大肩峰永远为8.95%，这种现象称为吉布斯（Gibbs）效应。 xxNNNNNWRsin2/)2/sin()2/sin()2/sin()(N43900.250.5

24、0.751-40-30-21-100N=15N=31用矩形窗设计的c=/2 FIR滤波器的幅度响应3.2 窗口设计法（时域）窗口设计法（时域）40 改变窗函数的形状改变窗函数的形状，可改善滤波器的特性，窗函数有许多种，但要满足以下两点要求：窗谱主瓣宽度要窄，以获得较陡的过渡带；相对于主瓣幅度，旁瓣要尽可能小旁瓣要尽可能小，使能量尽量集中在主瓣中，这样就 可以减小肩峰和余振，以提高阻带衰减和通带平稳性。但实际上对同样长度的窗同样长度的窗这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑制。 肩峰值的大小肩峰值的大小决定了滤波器通带内的平稳程度和阻带内决定了滤波器通带内的平稳程度和阻带内的衰

25、减，所以对滤波器的性能有很大的影响。的衰减，所以对滤波器的性能有很大的影响。3.2 窗口设计法（时域）窗口设计法（时域）41几种常用的窗函数： 1. 矩形窗，矩形窗，上面已讲过，不再细述 2. 汉宁窗汉宁窗（升余弦窗） 利用付氏变换的移位特性，汉宁窗频谱的幅度函数W（）可用矩形窗的幅度函数表示为： )(12cos1 21)(nRNnnwN)(25.0)(5.01212nReenRNNnjNnjN 211221122121121225. 05 . 011225. 05 . 0NjRRRNNjRNNjRNjRjeNWNWWeNWeNWeWeW3.2 窗口设计法（时域）窗口设计法（时域）42)12(

26、)12(25. 0)(5 . 0)(NWNWWWRRR三部分矩形窗频谱相加，使旁瓣互相抵消，能量集中在主瓣，旁瓣大大减小，主瓣宽度增加1倍，为 。 N83.2 窗口设计法（时域）窗口设计法（时域）433. 汉明窗汉明窗（改进的升余弦窗） 它是对汉宁窗的改进，在主瓣宽度（对应第一零点的宽度）相同的情况下，旁瓣进一步减小，可使旁瓣进一步减小，可使99.96%的能量集的能量集中在窗谱的主瓣内中在窗谱的主瓣内。 4. 布莱克曼窗布莱克曼窗（三阶升余弦窗） 增加一个二次谐波余弦分量，可进一步降低旁瓣，但主瓣可进一步降低旁瓣，但主瓣宽度进一步增加宽度进一步增加，为 。增加N可减少过渡带。频谱的幅度函数为：

27、 )(12cos46.054.0)(nRNnnwN)(14cos08.012cos5.042.0)(nRNnNnnwN)12()12(25. 0)(42. 0)(NWNWWWRRR)14()14(04. 0NWNWRRN1244下图为N=31时，矩形窗矩形窗、三角窗三角窗、汉宁窗汉宁窗、汉明窗汉明窗及布莱布莱克曼克曼这5种窗口函数的包络曲线3.2 窗口设计法（时域）窗口设计法（时域）45下图为下图为N=51N=51时时矩形窗矩形窗、汉宁窗汉宁窗、汉明窗汉明窗及及布莱克曼布莱克曼4 4种窗口种窗口函数的幅度响应函数的幅度响应3.2 窗口设计法（时域）窗口设计法（时域）465 . 051cN下图为

28、下图为N=51N=51时用时用矩形窗矩形窗、汉宁窗汉宁窗、汉明窗汉明窗及及布莱克曼布莱克曼设计设计的低通滤波器的幅度响应的低通滤波器的幅度响应3.2 窗口设计法（时域）窗口设计法（时域）475.凯塞窗 以上四种窗函数，都是以增加主瓣宽度为代价来降低旁瓣。凯塞窗则可自由选择主瓣宽度和旁瓣衰减。 101/211)(2NnINnInwooI0(x)是零阶修正贝塞尔函数，参数可自由选择，决定主瓣宽度与旁瓣衰减。越大，w(n)窗越窄，其频谱的主瓣变宽，旁瓣变小。一般取 49。=5.44 接近汉明=8.5 接近布莱克曼=0 为矩形3.2 窗口设计法（时域）窗口设计法（时域）483.2 窗口设计法（时域）窗

29、口设计法（时域）2crc1）根据逼近滤波器的理想频率响应，计算截止频率）根据逼近滤波器的理想频率响应，计算截止频率设计过程设计过程2）确定通带和阻带波动，通带最大衰减为）确定通带和阻带波动，通带最大衰减为Ap，阻带，阻带最大衰减为最大衰减为Ar0.050.05101101ppApA3）由于通带和阻带等波动，因此选择最小的波动）由于通带和阻带等波动，因此选择最小的波动0.0510rAr4）用最小波动重新计算）用最小波动重新计算Ap， Ar5）计算过渡带）计算过渡带 =wr-wp6）计算）计算值值210!)2/(1)(kkkxxI49过渡带宽过渡带宽 Ar为最小衰减，此为经验公式为最小衰减，此为经

30、验公式0.40.1102(8.7),500.5842(21)0.07886(21),21500,2182.286rrrrrrrAAdBAAdBAdBAdBAN3.2 窗口设计法（时域）窗口设计法（时域）503.2 窗口设计法（时域）窗口设计法（时域） 过渡带 通带波纹（ dB） 阻带最小衰减(dB) 2.120 3.00 /N 0.27 -30 3.384 4.46 /N 0.08647 -40 4.538 5.86 /N 0.0274 -50 5.658 7.24 /N 0.00868 -60 6.764 8.64 /N 0.00275 -70 7.865 10.0 /N 0.000868

31、-80 8.960 11.4 /N 0.000275 -90 10.056 12.8 /N 0.000087 -100 513.2 窗口设计法（时域）窗口设计法（时域）窗函数设计滤波器步骤：窗函数设计滤波器步骤：1 1、根据要求求出理想滤波器的单位脉冲响应，通常无限长、根据要求求出理想滤波器的单位脉冲响应，通常无限长2 2、根据阻带衰减的要求，选择窗函数类型；根据过渡带宽、根据阻带衰减的要求，选择窗函数类型；根据过渡带宽要求，估计窗口的长度要求，估计窗口的长度3 3、计算、计算FIRFIR滤波器的窗函数滤波器的窗函数h h（n n）hdhd（n n）w w（n n）4 4、计算、计算h h（n

32、 n）的）的DFTDFT验证频谱响应是否满足技术要求，如不验证频谱响应是否满足技术要求，如不满足，重复满足，重复2 2，3 3，4 4H edj()H edj()02nkMjkMjdMkMeeHMnh2210)(1)(hnd( )若得不到 封闭式或不能用上式计算 时，可对 在 到 间等间隔采样M，用下式代替上式的积分523.2 窗口设计法（时域）窗口设计法（时域）例例1 1：根据下面的技术指标设计一个：根据下面的技术指标设计一个FIRFIR线性相位低通滤波器线性相位低通滤波器0.99()1.01, 00.19()0.01, 0.21jwjwH ewH ew例例2 2：确定幅度频响应为下图的理想

33、脉冲响应：确定幅度频响应为下图的理想脉冲响应21/6/3533.2 窗口设计法（时域）窗口设计法（时域）例：根据下面的技术指标设计一个例：根据下面的技术指标设计一个FIRFIR线性相位线性相位LPFLPF20log(0.01)40dB 阻带衰减：阻带衰减：过渡带：过渡带：0.21 0.190.02汉宁汉宁：0.026.6 / 6.6/0.02330N N汉明：汉明：0.026.2 / 6.2/0.02310N N布莱克曼：布莱克曼：0.0211 / 11/0.02550N N543.2 窗口设计法（时域）窗口设计法（时域）例：确定幅度频响应为下图的理想脉冲响应例：确定幅度频响应为下图的理想脉冲

34、响应21/6/3121/3 ( )( )( )sin ()/3sin ()/6,()dddnah nhnhnnananana553.2 窗口设计法（时域）窗口设计法（时域）56nnnndeenhccnjjdcc/,)()(sin21)(2 . 0cr302 . 0285. 2850N55. 4)7 . 850(1102. 03.2 窗口设计法（时域）窗口设计法（时域）573.2 窗口设计法（时域）窗口设计法（时域）58 wn=kaiser(30,4.55);nn=0:1:29;alfa=(30-1)/2;hd=sin(0.4*pi*(nn-alfa)./(pi*(nn-alfa);h=hd.*

35、wn;h1,w1=freqz(h,1);plot(w1/pi,20*log10(abs(h1);axis(0,1,-80,10);grid;xlabel(归一化频率/)ylabel(幅度/dB)3.2 窗口设计法（时域）窗口设计法（时域）593.3 频率采样法频率采样法一、基本思想一、基本思想 使所设计的FIR数字滤波器的频率特性在某些离散频率某些离散频率点上的值准确地等于所需滤波器在这些频率点处的值点上的值准确地等于所需滤波器在这些频率点处的值，在其它频率处的特性则有较好的逼近。jnhNIDFTNNkjdjdeHnhkHeHeHd )(2)()(不同于点点频率取样确定内插公式60二二.设计方

36、法设计方法1）确定频率采样值2）利用IDFT计算3）计算系统函数或频率响应,kkH(),()jwHeHZ)(nh,)()(2kjkNkjdeHkHeH0,1,1kN,)(1)(10/2NkNnkjekHNnh1, 1 , 0Nn1010()( )( )( )NjwjwnnNnnH eh n eH zh n z3.3 频率采样法频率采样法61三、 约束条件 为了设计线性相位的FIR滤波器，采样值 H（k）要满足一定的约束条件。 前已指出,具有线性相位的线性相位的FIR滤波器，其单位脉冲滤波器，其单位脉冲响应响应h(n)是实序列，且满足是实序列，且满足 ，由此得到的幅频和相频特性，就是对H(k)的

37、约束。)1()(nNhnh3.3 频率采样法频率采样法下面以偶对称条件下面以偶对称条件h(n)=h(N-n-1)为例说明设计公式为例说明设计公式。62kNjeHkH2)()(H ej()( )02 对 在 等间隔N点采样得H(k) ,k=0,1,N-1 ) 1(21)(N对线性相位滤波器，其相位响应为： )()()(jgjeHeH 幅度响应幅度响应)(nh)()()(kjgekHkH令kNggHkH2)()(则：kNk2)()(,3.3 频率采样法频率采样法63)(kHg对频率采样的幅度值 有：)()(kNHkHggN N为奇数为奇数 )()(kNHkHgg0)2()2(NHNHgg及N N为

38、偶数为偶数 按线性相位的偶对称条件: )(gH)2()(ggHH若若N N为偶数，为偶数，对对为奇对称为奇对称: : )(gH)2()(ggHH对对为偶对称为偶对称: : 若若N为奇数，为奇数，kNNNkNNNkN1) 1()(1)(kNNkNNk12) 1(21)( ，k =0,1,(N-1)/2按式 对相位进行频率采样得： 2/) 1()(N3.3 频率采样法频率采样法64NkNjjNkNjNjee) 1() 1() 1(对对N N为偶数为偶数，(N-1)为奇数 取kNNkN1)( ，k=0,1,N/2-1NkNjNkNjNjee) 1() 1() 1(对对N N为奇数为奇数，(N-1)为

39、偶数kNNkN1)( 取 ，k=0,1,(N-1)/2综合以上分析，可得到设计公式如下：综合以上分析，可得到设计公式如下：3.3 频率采样法频率采样法65综合以上分析，可得到设计公式如下：综合以上分析，可得到设计公式如下：)()(kNHkHggkNNk1)(kNNkN1)(,对对N N为奇数，设计公式为：为奇数，设计公式为：k=0,1,(N-1)/20)2(,NHg)()(kNHkHggkNNk1)(kNNkN1)(,，k=0,1,N/2-1对对N N为偶数，设计公式为：为偶数，设计公式为：3.3 频率采样法频率采样法66 滤波器的频率响应jez 将 代入频率采样公式频率采样公式得：)2()(

40、)()(10kNkHzHeHNkezjj21)2/sin()2/sin(1)(NjeNN其中Nk21, 2 , 1 , 0,Nk1)2(Nk在采样点在采样点与NkjkkeH/2)(无误差)()(2NkjdeHkH但在采样点之间在采样点之间， 两者误差与 特性的平滑程度有关： )(jdeH在 幅度曲线的平滑段，误差较小幅度曲线的平滑段，误差较小，但在曲线的间曲线的间断点附近，会产生较大的误差断点附近，会产生较大的误差，使得滤波器的阻带性能变坏 。误差还与采样点数N有关，N越大误差越小越大误差越小。)(jdeH3.3 频率采样法频率采样法67为提高阻带衰减，常用增加过渡带法如右图为提高阻带衰减，常

41、用增加过渡带法如右图3.3 频率采样法频率采样法68| )()(|jdjeHeH优化原则：优化原则：在通带内要求min| )()(|maxjdjeHeH在阻带内要求例例 用频率采样法设计一个低通滤波器，通带截止频率radp2 . 0的偶对称情况。，采样点数N=20，采用h(n)=h(N-n-1)2 , 0解解：N=20，在范围内等间隔采样，显然，在通带共有3个采样点，分别是k=0,1,2radN1 . 0/2采样间隔为：173019181201)(kkkkHg1918)20(95. 02095. 02025 . 9)(kkkkkk又N为偶数，可得：3.3 频率采样法频率采样法69)()()(k

42、jgekHkH)(kHg)(k将 和 代入)(kH)(zH)(nh求得 ，进而求得或其相位响应为线性相位。过渡带为：2/20 /10其幅度响应如下图，图中还给出了其单位脉冲响应3.3 频率采样法频率采样法70346035, 53904. 03936, 401)(kkkkkkHg3935)40(975. 050975. 04025 .19)(kkkkkk)(nh)(kHg)(k由上图 d）可见，所设计的滤波器的阻带衰减很小，只有-16dB。为了改进阻带衰减，在边界频率处增加一个过渡点，为了保证过渡带宽不变，将采样点数增加一倍，变为N=40，并将过渡点的采样值进行优化，取H1=0.3904，得到

43、和 分别为过渡带为（2/402） /10，求得的单位脉冲响应和幅度响应如下图所示。3.3 频率采样法频率采样法713.3 频率采样法频率采样法72)(kHg)(k由上图 d）可见，这时阻带衰减达到了-43dB。为了进一步增加阻带衰减，可再增加一个过渡采样点，并将采样点数增加到60，两个过渡样点值经优化分别为=0.5925和=0.1099, 相应的 和 分别为106,54590.59257,53( )0.10998,520851gkkkkHkkkk5952)60(983. 080983. 06025 .29)(kkkkkk)(nh过渡带为（2/603） /10所得单位脉冲响应和幅度响应如下图所示

44、。 3.3 频率采样法频率采样法733.3 频率采样法频率采样法74增大阻带衰减三种方法：增大阻带衰减三种方法：1）加宽过渡带宽，以牺牲过渡带换取阻带衰减的增加。加宽过渡带宽，以牺牲过渡带换取阻带衰减的增加。 2）过渡带的优化设计过渡带的优化设计利用线性最优化的方法确定过渡带采样点的值，得到要求的滤波器的最佳逼近（而不是盲目地设定一个过渡带值）。上例中H9=0.3904。对应的 H(ej)的幅频特性，比H9=0.5时 的阻带衰减大大改善。 3）增大增大N 如果要进一步增加阻带衰减，但又不增加过渡带宽，可增加采样点数N。代价是滤波器阶数增加，运算量增加。代价是滤波器阶数增加，运算量增加。75小结

45、：小结：频率采样设计法优点： 直接从频域进行设计，物理概念清楚，直观方便； 适合于窄带滤波器设计，这时频率响应只有少数几个非零值。典型应用：用一串窄带滤波器组成多卜勒雷达接收机，覆盖不同的频段，多卜勒频偏可反映被测目标的运动速度； 缺点：截止频率难以控制截止频率难以控制。 因频率取样点都局限在2/N的整数倍点上，所以在指定通带和阻带截止频率时，这种方法受到限制，比较死板。 充分加大N，可以接近任何给定的频率，但计算量和复杂性增加。76最优化设计的前提是最优准则的确定，在FIR滤波器最优化设计中，常用的准则有 最小均方误差准则 最大误差最小化准则。3.4 最大误差最小化设计最大误差最小化设计77

46、deEdeHeHjjjd22221211) 均方误差最小化准则， 若以E(ej)表示逼近误差，则 那么均方误差为)()(jjdjeHeHeE）（3.4 最大误差最小化设计最大误差最小化设计就是选择一组时域采样值，以使均方误差就是选择一组时域采样值，以使均方误差 ，这一方，这一方法注重的是在整个法注重的是在整个-频率区间内频率区间内总误差的全局最小总误差的全局最小，但，但不不能保证局部频率点的性能能保证局部频率点的性能，有些频率点可能会有较大的误差。，有些频率点可能会有较大的误差。min278对于窗口法窗口法FIR滤波器设计滤波器设计，因采用有限项的h(n)逼近理想的hd(n)，所以其逼近误差为

47、：如果采用矩形窗 则有ndnhnh22)()(其它01)()(Nnonhnhd1222| )()(| )()(|nNnddnhnhnhnh3.4 最大误差最小化设计最大误差最小化设计这是一个最小均方误差。矩形窗窗口设计法是一个最小均方矩形窗窗口设计法是一个最小均方误差误差FIR设计设计其优点是过渡带较窄过渡带较窄，缺点是局部点误差大局部点误差大,或者说误差分布不均匀。792) 最大误差最小化准则最大误差最小化准则（也叫最佳一致逼近准则）表示为其中F是根据要求预先给定的一个频率取值范围，可以是通带，也可以是阻带。最佳一致逼近即选择选择N个频率采样值个频率采样值 （ 或或时域时域 h(n) 值值

48、），在给定频带范围内使频响的最大逼近误差达），在给定频带范围内使频响的最大逼近误差达到最小。也叫等波纹逼近。到最小。也叫等波纹逼近。优点：可保证局部频率点的性能也是最优的，误差分布均匀，误差分布均匀， 相同指标下，可用相同指标下，可用最少的阶数最少的阶数达到最佳化达到最佳化。 min| )(|maxjeEF3.4 最大误差最小化设计最大误差最小化设计80切比雪夫最佳一致逼近切比雪夫最佳一致逼近用等波纹逼近法设计滤波器需要确定五个参数： M、c、r、1、2按上图所示的误差容限设计低通滤波器，就是说要在通带 0 p 内以最大误差 1 逼近1，在阻带r 内 以最大误差2逼近零。 3.4 最大误差最小

49、化设计最大误差最小化设计cr813.4 最大误差最小化设计最大误差最小化设计要同时确定上述五个参数较困难。常用的两种逼近方法两种逼近方法： 1）给定M、1、2，以c和r为变量。 缺点：边界频率不能精确确定。 2）给定M、c和r，以1和2为变量，通过迭代运算 ，使逼近误差1和2 最小，并确定h(n)切比雪切比雪 夫最佳一致逼近夫最佳一致逼近。 特点：能准确地指定通带和阻带边界频率。 82一.误差函数误差函数 定义逼近误差函数： )()(HHWEd E 为所设计的滤波器与理想滤波器的幅频特性在通带和阻带内的误差值， 是已知的权函数权函数，在不同频带可取不同的值， 所要设计的滤波器的幅频特性 理想滤

50、波器的幅频特性 W H dH3.4 最大误差最小化设计最大误差最小化设计83 rcdH001 rckW101例如，希望在固定 M, c, r 的情况下逼近一个低通滤波器，这时有21k21cos)()(0NMnnaHMn21, 2 , 1,212)(,21)0(NnnNhnaNha对于第一种线性相位滤波器（偶对称，奇数长度）3.4 最大误差最小化设计最大误差最小化设计840( )( )( )( )cosMdnEWHa nn切比雪夫逼近问题变为，寻求一组系数使逼近误差的最大切比雪夫逼近问题变为，寻求一组系数使逼近误差的最大值达到最小，即值达到最小，即, 1 , 0),(Mnna2minmaxE 给定后等效于求 最小。c0r21/k3.4 最大误差最小化设计最大误差最小化设计85二二.交替定理（最佳逼近定理）交替定理（最佳逼近定理） 令F表示闭区间 的任意闭子集，为了使 在 F 上唯一最佳地逼近于 ,其充分必要条件是误差函数 在 F 上至少应