数字信号处理第七章-有限单位冲激响应数字滤波器的设计-课件.ppt

1、第七章 FIR滤波器的设计IIR数字滤波器：数字滤波器：可以利用模拟滤波器设计可以利用模拟滤波器设计但相位非线性但相位非线性FIR数字滤波器：数字滤波器：可以严格线性相位，又可任意幅度特性可以严格线性相位，又可任意幅度特性因果稳定系统因果稳定系统可用可用FFT计算计算但阶次比但阶次比IIR滤波器要高得多滤波器要高得多主要内容 线性相位FIR滤波器的特点 窗函数设计法 频率抽样设计法 IIR与FIR比较7.1 引言一、一、FIR滤波器的主要特点：滤波器的主要特点：q 单位冲激响应只有有限多项单位冲激响应只有有限多项q 可以设计成线性相位系统可以设计成线性相位系统q 只在原点处有极点，因此系统总是

2、稳定的只在原点处有极点，因此系统总是稳定的q 便于便于DSPDSP实现（并可用立即数乘加指令编程，节实现（并可用立即数乘加指令编程，节约存储器）约存储器）二、FIR与IIR相比较：q 首先在相频特性控制上可以做到线性相位，IIR而不能做到这一点，这一点在通信等领域中要求却很重要；q 其次，FIR不存在稳定性问题，其非递归结构不会产生极限环现象等有限精度问题；q 最后，FIR还可以FFT用来滤波。故FIR应用越来越多。三、线性相位设计的重要性1、系统的相移会造成信号波形的改变时间时间 t幅幅度度原始信号原始信号时间时间 t幅幅度度相移相移90o时间时间 t幅幅度度相移相移 180o2、系统非线性

3、相移造成输出信号失真dd)()(f1 f2f时时延延f1 f2f时时延延f1 f2f()f1 f2f()系统相位特性决定了信号不同频率的时延系统相位特性决定了信号不同频率的时延3、忽略相位信息的后果输入波形输入波形DFT变换变换忽略相忽略相位信息位信息IDFT变换变换输出波形输出波形4 4、要求线性相位的例子 通信系统：调制解调器、综合业务数据网（ISDN）等。希尔伯特变换器：要求输入输出信号正交。高保真音响系统：音乐的相位失真必须减到最小，尽可能逼真地重现原来的声音。理想微分器：线性相位要求：5 5、线性相位的FIRFIR滤波器设计基础constant)(ddgg)(NpgpTTph00)(

4、sin-系统的群延迟系统的群延迟7.2 线性相位FIR滤波器特点 FIR滤波器的单位冲激响应：滤波器的单位冲激响应：()01h nnN10()()NnnH zh n z系统函数：系统函数：在在 z 平面有平面有N 1 个零点个零点在在 z=0 处是处是N 1 阶极点阶极点 一线性相位条件一线性相位条件如果如果FIR DFFIR DF的单位抽样响应的单位抽样响应h h（n n）为实数，）为实数，而且满足偶对称而且满足偶对称h h（n n）=h=h（N-1-nN-1-n），或满），或满足奇对称足奇对称h h（n n）=-h=-h（N-1-nN-1-n），其对称中心），其对称中心在在 处，可证明处，

5、可证明filterfilter就具有就具有准确的线性相位。准确的线性相位。N N又分为偶数和奇数两种情况，所以有又分为偶数和奇数两种情况，所以有4 4种种线性相位线性相位FIR DFFIR DF，如下所述。，如下所述。12Nn1、N为奇数的偶对称例如 N=11，对称中心为n01234567891011 15,()(10)2nh nhn2、N为偶数时的偶对称例如 N=10，对称中心为n012345678910 14.5,()(9)2nh nhn3、N为奇数时的奇对称例如，N=11，对称中心为 n0123456789105,()(10)nh nhn 4、N为偶数时的奇对称例如，N=10，对称中心为

6、4.5，n012 3456789()(9)h nhn 有两类准确的线性相位，分别要求满足：有两类准确的线性相位，分别要求满足：()()其中其中、均为常数均为常数因此有：111()1()()()()()()Njj njjnNjj njjnH eh n eH eeH eh n eH ee令上两式实部虚部相等，则有：10()cos()()cos()NjnH eh nn10()sin()()sin()NjnH eh nn若将前页的实部虚部相除：1010()sin()sin()tan()cos()()cos()NnNnh nnh nn从而有：110010()cos()sin()()sin()cos()0

7、()sin()0NNnnNnh nnh nnh nn要使上式成立，必须：12()(1)Nh nh Nn 对另一相位形式，必须有：12()(1)Nh nh Nn 二线性相位频率响应特点命题：设FIR单位冲激响应h(n)为实序列，且满足偶对称（或奇对数）条件：)1()(nNhnhor)1()(nNhnh()()()()()jjjh nH eH eHe 1()()2N (1)()22Nor 则：则：)1()(nNhnh10)1(10)1(11010)()()1()()(NmmNNmmNnNmNnnNnnzmhzzmhznNhznhzH)()(1)1(zHzzHN证明：证明：1 1、偶对称时：、偶对称

8、时：即：即：)()(21)(1)1(zHzzHzHN10)1()(21NnnNnzzznh 所以有所以有：10)2)1()2)1(2)1(2121)(NnNnNnNZZnhz)()(1)1(zHzzHN)1()2)1(cos)()()(102)1()(HNneNjjNnnheeHj则则 为线性相位。为线性相位。2)1()(N其物理意义：其物理意义：该该FIR有有(N-1)/2个个 采样周期的群时延。采样周期的群时延。10)2)1()2)1(2)1(2121)()(NnNnNnNZZnhzzH)21N()1N()(,2)1N(0)(2)21N()(,2奇对称时)1()(nNhnh10)1(10)

9、1(11010)()()1()()(NmmNNmmNnNmNnnNnnzmhzzmhznNhznhzH)()(1)1(zHzzHN即即所以有：所以有：)()(21)(1)1(zHzzHzhN10)1()(21NnnNnzzznh)()(1)1(zHzzHN10)2)1()2)1(2)1(2121)(NnNnNnNZZnhz10)2)1()2)1(2)1(2121)()(NnNnNnNZZnhzzH()(1)120()(1)()()sin()(2)2jNNjjneHNH ejeh nn 1022)1()2)1(sin)()(NnjNjjNnnheeH或或(1)()22N 则则 为线性相位为线性相

10、位 可见，其相位特性是线性相位，而且还产生一个可见，其相位特性是线性相位，而且还产生一个90900 0相移，这样就使得通过相移，这样就使得通过filterfilter的所有频率都相移的所有频率都相移90900 0，因此称它为因此称它为正交变换网络正交变换网络。（相移。（相移90900 0的信号与原信的信号与原信号为正交的）。号为正交的）。)23N()(,2)12N()(,2)(,0)23(N)12(N)(202二幅度特点二幅度特点1、h(n)偶对称，偶对称，N为奇数为奇数对（对（1）式）式10)21(cos)()(NnNnnhH由于由于:)()1(nhnNh)()()(jjeHeH11cos(

11、1)cos()221cos(1)2NNNnnNNn (3)/201(1)/22111()()2()cos()2211()2()cos22NnNNmnmNNHhh nnNNhhmm 2/)1(0)cos()(Nnnna111(0)(),()2(),(1,2,.)222NNNaha nhnn由于由于得得其中其中:ncos2,0)(H2,0由于由于 对对 是偶对称的。是偶对称的。因此，因此，对对 为偶对称。为偶对称。线性相位滤波器的幅度特点线性相位滤波器的幅度特点其中，其中，/2 10/2211()2()cos()212()cos()22NnNNnmmNHh nnNhmm 2/1)21(cos)(N

12、nnnb(0)0,()2(),(1,2,.)22NNbb nhnn2、h(n)偶对称，偶对称，N为偶数为偶数)21(Nh对（对（1）式与如上合并项，注意到由于）式与如上合并项，注意到由于N为偶数，为偶数，项即为项即为0，则，则 0)21(cosn)0|)(1zzH)(H 由于由于 时，时，且对且对 呈奇对称。因此，呈奇对称。因此，对对 呈奇对称。呈奇对称。并有并有:111sin()sin()sin(1)222NNNnnNn )1()(nNhnh)21(|)1()21(21NhnNhNhNn0)21(Nh3、h(n)奇对称，奇对称，N为奇数为奇数所以有：所以有：1sin()2Nn为奇对称的，即：

13、10(3)/201()()sin()212()sin()2NnNnNHh nnNh nn(1)/21()sinNmc nn 代入式代入式:10(1)()()sin()2NjnNH eh nn111sin()sin()sin(1)222NNNnnNn 0)21(Nh其中其中,11(0)0,()2()(1,2,.)22NNcc nhnnnsin2,0由于由于 在在 均为均为0 0并对这些点呈奇对称。并对这些点呈奇对称。)1()(nNhnh10)21(sin)()(NnNnnhH/2 1012()sin()2NnNh nn/2112()sin()22NmNhmm(1)/211()sin()2Nmd

14、nm)1()12(2)(,0)0(nnNhndd其中：其中：对（对（2 2）式）式4、h(n)奇对称，奇对称，N为偶数为偶数线性相位滤波器的幅度特点2,02,0)21(sinn)(H由于由于 在在 处为处为0 0。因此，因此，对对 呈奇对称。呈奇对称。线性相位滤波器的幅度特点总结：（1）第1，2种一般为低通特性；第3，4种一般为高通、带通特性。（2）当N，h(n)均为偶（或奇）时，H(w)为奇对称。当N，h(n)为一奇一偶时，H(w)为偶对称。)2,0(:()2,0(:(四、系统函数四、系统函数H（Z）的零点分布情况）的零点分布情况 1、零点的分布原则、零点的分布原则所以，如果所以，如果 是零

15、点，则是零点，则 也一也一定是定是H（Z）的零点，的零点，h（n）为实数时，）为实数时，H（Z）的零点必成共轭对出现，即的零点必成共轭对出现，即 也一定是也一定是H（Z）的零点，）的零点，也一定是也一定是H（Z）的零）的零点。点。(1)1()()NH ZZH Z iZZ1/iZZ*iZZ*1/iZZ2、零点的位置（1）既不在实轴上，也不在单位圆上，则零 点是互为倒数的两组共轭对，22/1,22/1jZjZiiiZ*iZ/1ZjImZiiZiZ/1ZRe10*1111,4444iiZjZj(2)不在实轴上，但在单位圆上，共轭对的倒数就是它们本身，如iZ*iiZ/1Z i*iZ1Z 01*2222

16、,222222221/,1/2222iiiiZjZjZjZj（3）在实轴上，不在单位圆上，实数零点，没复共轭；只有倒数。例如，iZ/101122iZ1/2,1/2iiZZiZ（4）既在实轴上也在单位圆上。此时，只有一个零点，且有两种可能，或位于Z=1，或位于Z=-1。iZN为偶数时的偶对称为其零点；N为偶数奇对称H（0）=0，有Z=1零点；N为奇数奇对称有零点Z=1，和Z=-1。1Z,0)(H,0)(H)0(H1iZ 1iZ 7-3 窗函数设计法一、设计方法 1、设计思想 先给定理想filter的频响 ，所要求设计一个FIR的filter的频响为 ，使 逼近 2、设计过程 设计是在时域进行的，

17、先用傅氏反变换求出理想filter的单位抽样响应 ，然后加时间窗对 截断，以求得FIR filter的单位抽样响应h(n)。)e(Hjd)e(Hj)e(Hj)e(Hjd)n(hd)(nw)n(hd)()()()(21)(nhnwnhdeeHnhdnjjdd例如，低通filter)(Hd0cc 是矩形的，则 一定是无限长的且是非因果的。)e(Hjd)n(hd二、窗函数对频响的影响 1、理想LF的单位抽样响应理想低通filter的频响 为)(nhd)(jdeH)e(Hjdcccje,0,10)e(Hjdcc0)(为群延时因为其相位 ，所以 是偶对称，其对称中心为 ，这是因为 时，即 为其最大，故

18、为其对称中心。又是无限长的非因果序列)n(hdccc)n(j)n(jnjjjd1)n()nsin(e)n(2 j11de21dee21)e(HFccccc)n(hd)(n/)(cdh)(nhd)n(hd2/)1N(1Nn)n(RN1Nn0.12、加矩形窗 加窗就是实行乘操作，而矩形窗就是截断数据，这相当于通过窗口 看 ，称 为窗口函数。)()(nRnwNR)n(RN)n(hd)(nwR)()()(nwnhnhRd1Nn0),n(hd,0其他n值 因h（n）是偶对称的。长度为N，所以其对称中心应为 ，所以h（n）可写作2/)1N(h（n）=1Nn0,)21Nn()21Nnsin(ccc,0n为其

19、他值3、h（n）的频响 h（n）的频响 可通过傅式变换求得，为了便于与 的频响 相比较，利用卷积定理)e(Hj)()(nhFeHj)n(hd)e(HjddeWeHeHnwnhnhjRjdjRd)()(21)()()()()(1)对于矩形窗的频响10)()()(NnnjRRjRenwnwFeW2/sin2Nsinee1e1e)21N(jjNj1N0nnj)21()(NjReW 其中，为幅度函数，为相位函数。)2sin(/)2Nsin()(WR)21()(N（2）对于理想LF的频响)21N(jdjde)(H)e(H 其中，为幅度函数，为相位函数。)(Hdc,1c,0)21()(N（3）h（n）的频

20、响)e(Hjd)(W)(H21ede)(We)(H21Rd)21N(j)(21N(jR)21N(jd其中，为幅度函数，为相位函数。d)(W)(H21)(HRd)21()(N4、窗函数频响产生的影响从几个特殊频率点的卷积过程就可看出其影响:（1）时，0ccd)(W21d)(W121)0(HRR也就 在 到 全部面积的积分。因此，H（0）/H（0）=1（用H（0）归一化）。)(WR)(Hd0cc)(WRN/2N/20NWR/2)(的主瓣宽度的一半为注意：（2）时，正好与 的一半相重叠。这时有 。c)(WR)(Hd5.0)0(/)(HHc(3)时，的主瓣全部在的通带内，这时应出现正的肩峰。Nc2)(

21、RW)(dH（4）时，主瓣全部在通带外，时，主瓣全部在通带外，出现负的肩峰。出现负的肩峰。Nc/2（5）当 时，随 增加，左边 旁瓣的起伏部分扫过通带，卷积 也随着 的旁瓣在通带内的面积 变化而变化，故 将围绕着零值而波动。Nc2)(WR)(H)(WR)(H N/2c(6)当 时，的右边旁瓣将进入 的通带，右边旁瓣的起伏造成 值围绕 值而波动。N2c)(WR)(Hd)(H)0(H100.5)0(H/)(H 5、几点结论（1）加窗后，使频响产生一过渡带，其宽度正好等于窗的频响 的主瓣宽度（2）在 处出现肩峰，肩峰两侧形成起伏振荡，其振荡幅度取决于旁瓣的相对幅度，而振荡的多少则取决于旁瓣的多少。（

22、3）吉布斯（Gibbs）效应 因为窗函数的频响的幅度函数为这是一个很特殊的函数，分析表明，当改变N时仅能改变 的绝对值的大小，和主瓣的宽度 ，旁瓣的宽度 ，但不能改变主瓣与旁瓣的相对比例，也就是说，不会改变归一化频响 的肩峰的相对值。对于矩形窗最大相对肩峰为8.95%，不管N怎样改变，最大肩峰总是 8.95%，这种现象称作吉布斯效应。)(WRN4)(H N2c)2sin(/)2Nsin()(WR)(WR)N/4()N/2()(H 三、各种窗函数 1、基本概念（1）窗谱：窗函数的频响的幅度函数亦称作窗谱。（2）对窗函数要求 a）希望窗谱主瓣尽量窄，以获得较陡的过渡带，这 是因为过渡带等于主瓣宽度

23、。b）尽量减少窗谱最大旁瓣的相对幅度，这样可使肩峰 和波纹减少。2、矩形窗 时域表达式：频域表达式（频谱）：幅度函数：)()()(nwnRnwRN)21N(jRjRe)(W)e(W)2sin(/)2Nsin()(WR3、三角形（Bartlett）窗时域表达式：)(nw21Nn0,1Nn21Nn21N,1Nn22121210 1 2 3 4 频谱：)21N(j2je)2sin()41Nsin(1N2)e(W1N,e)2sin()4Nsin(1N2)21N(j2 第一对零点为 ，即 ，所以主瓣宽度 ，比矩形宽一倍。)e(Wj4NN4N/84、汉宁窗（升余弦窗）其窗谱可利用如下方法求出，将 变形为又

24、由于 其中又考虑到 ，这里)()12cos(1 21)(nRnNnwN)1Nn0()(nw)(41)(41)(21)()12()12(nRenRenRnwNnNjNnNjN)21()()()(NjRRjReWnwFeW)2sin(/)2Nsin()(WR)n(xe F)e(Xnj)(j001N20所以有)e(Wj)21()21()()12(41)12(41)(21)(NjNjRRReWeNWNWWnwF当 时，窗谱分析 可知，它等于三部分之和，旁瓣较大程度地互相抵消，但主瓣加宽一倍，即为1N N1N)N2(W41)N2(W41)(W21)(WRRR)(W N8)(W21R)N2(W41RN4N

25、2N2N4)(W N4N4汉宁窗是 时，特例2)n(R)1Nn(sin)n(N5、海明窗，又称作改进升余弦窗 其窗函数为仿照汉宁窗的分析方法可以得其频响的幅度函数为 其主瓣宽度仍为 ，（旁瓣峰值/主瓣峰值）1%有99.963%的能量集中在主瓣内。海明窗是下一类窗的特例)N2(W)N2(W23.0)(W54.0)1N2(W)1N2(W23.0)(W54.0RRRRRR)()12cos(46.054.0)(nRNnnWN)(W N8)54.0()()12cos()1()(nRNnnwN6、布拉克曼窗，又称二阶余弦窗 加上余弦的二次谐波分量，可以进一步抑制旁瓣相应的幅度函数为 其主瓣宽度为 ，是矩形

26、窗的三倍。)()14cos(08.0)12cos(5.042.0)(nRNnNnnwN)(W)1N4(W)1N4(W04.0)1N2(W)1N2(W25.0)(W42.0RRRRRN/127、五种窗函数的比较（1）时域窗布拉克曼三角矩形海明21N1Nn)n(（2）各个窗的幅度函数，如P.340，图7-11，注意图中 是dB表示的。（3）理想LF加窗后的幅度函数（响应）如P340，图7-12所示。四、窗函数法的设计 1、设计步骤（1）给定频响函数（2）求出单位抽样响应（3）根据过渡带宽度和阻带最小衰减，借助窗函数 基本参数表（P202表3）确定窗的形式及N的大小（4）最后求 及 2、设计举例)e

27、(Hjd)e(HF)n(hjd1d)()()(nwnhnhd)e(Hj例1：分别利用矩形窗与汉宁窗设计具有线性相位的 FIR 低通滤波器，具体要求：)e(Hjd,0,ecj其他;rad1,s12c并画出相应的频响特性解：（1）由于 是一理想LF，所以 可以得出 （2）确定N 由于相位函数 ，所以 呈 偶对称，其对称中心为 ，因此 )e(Hjd)n(hd)n()n(sin)n(hcccd)()n(hd2/)1N(2512N)12n()12nsin(1)n(hd（3）加矩形窗)()()()()(25nRnhnwnhnhdd24,2,1,0n),12n(/)12nsin(则有可以求出h（n）的数值，

28、注意偶对称，对称中心122/)1N(31831.0)12(14472.0)14()10(06022.0)16()8(01482.0)18()6(03936.0)20()4(01931.0)22()2(;01423.012/12sin)24()0(hhhhhhhhhhhhh26785.0)13(h)11(h01497.0)15(h)9(h06104.0)17(h)7(h02987.0)19(h)5(h01457.0)21(h)3(h02893.011/11sin)23(h)1(h)n(hn1224由于h（n）为偶对称，N=25为奇数，所以)(H 121n2/)1n(1n2/)1N(0n)ncos

29、()n12(h2)12(h)ncos()n21N(h2)21N(h)ncos()n(a例如 H（0）=0.94789，可以计算 的值，画如下图()H（4）加汉宁窗 由于 可以求出序列的各点值240),242cos(1 21)(nnnw1)12(9330.0)14()10(75.0)16()8(5.0)18()6(25.0)20()4(06698.0)22()2(0)24()0(wwwwwwwwwwwww9829.0)13()11(85355.0)15()9(62940.0)17()7(37059.0)19()5(1464.0)21()3(01903.0)23()1(wwwwwwwwwwww通过

30、 可求出加窗后的h（n）)()()(nwnhnhd31831.0)12()12()12(whhd13502.0)14(h)10(h04516.0)16(h)8(h00741.0)18(h)6(h00984.0)20(h)4(h00116.0)22(h)2(h0)24(h)0(h26326.0)13(h)11(h1277.0)15(h)9(h003841.0)17(h)7(h01107.0)19(h)5(h00213.0)21(h)3(h00049.0)23(h)1(h相应幅度函数可用下式求得：121n)ncos()n12(h2)12(h)(H如H（0）=0.98460，图如下例2.设计一个线性

31、相位FIR低通滤波器，给定抽样频率为 ，通带截止频率为 ，阻带起始频率为 ，阻带衰减不小于-50dB.幅度特性如图所示。421.5 10(/sec)srad 321.5 10(/sec)prad 323 10(/sec)strad 解：（1）求对应的数字频率：0dB50dB00.51.5stcp通带截止频率为：20.2pppssf阻带起始频率为：20.4stststssf（2）求 ：设 为理想线性相位滤波器()dh n()jdHe,()0,jcjdeHe其它31()22.25 10(/sec)2cpstrad 20.3cccssf()11()221sin(),1()2,jj njndcch n

32、eedednnNnn其中（3）求窗函数：由阻带衰减确定窗形状，由过渡带确定N。250dB 查表7-3知可选海明窗。所要求的过渡带宽（数字频域）：20.2stps 又海明窗的过渡带满足：6.6N由此求得：33N 所以：1162N（4）求h(n)：由海明窗的w(n)确定FIR滤波器的h(n)。2()0.540.64cos()()1Nnw nRnN1sin()2()1()2cdNnh nNn所以：1sin()22()()()0.540.64cos()()11()2cdNNnnh nw n h nRnNNn7-4、凯泽（、凯泽（Kaiser）窗及其滤波器设计）窗及其滤波器设计 上述几种窗函数：矩形窗、

33、汉宁窗、海明窗等，为了压制旁瓣，是以加宽主瓣为代价的。而且，每一种窗的主瓣和旁瓣之比是固定不变的，而凯泽窗可以在主瓣宽度与旁瓣衰减之间自由选择。一、凯泽窗 凯泽在1966（1974）发现，利用第一类零阶修正（变形）贝赛尔函数可以构成一种近似最佳的窗函数。凯泽窗定义为：1。定义其它,010,)()/)(1()(02/120NnInInW2/)1(N)(0I其中，为第一类零阶修正贝塞尔函数，是一个可自由选择的参数。2.特点0可同时调整主瓣宽度与旁瓣可同时调整主瓣宽度与旁瓣;越大，越大，窗越窄。频谱旁瓣越小，而主瓣窗越窄。频谱旁瓣越小，而主瓣相应增加；相应增加；相当于矩形窗相当于矩形窗;)(nW94

34、通常选择通常选择，它们相当于旁瓣与主，它们相当于旁瓣与主瓣幅度为瓣幅度为 3.1%-0.047%；凯泽窗随凯泽窗随 变化的曲线如下图：变化的曲线如下图：注：第一类零阶修正贝塞尔函数为.)!3()2/()!2()2/()2/(1!)2/(1)(26242120 xxxkxxIkk由图可以看出,为对称中心，且是偶对称,2/)1(Nn即即)1()(nNWnWkk1)()()2/1()(00IINWnWkk3.凯泽经验公式凯泽经验公式该公式可使该公式可使filter设计人员根据设计人员根据filter的设计指标的设计指标,估算出估算出值和值和 N 值。值。且，且，11sp1)(jeH：通带截止频率，由

35、：通带截止频率，由 定定;：止带截止频率，由：止带截止频率，由 定定.)(jeHps过渡带宽度过渡带宽度2/)(psc1285.2/)8(21,0.05021,07886.0)21(5842.050),7.8(1102.0lg20.4.0ANAAAAAAps4.设计举例设计举例利用凯泽窗设计一利用凯泽窗设计一FIR低通低通filter，要求，要求6.0,4.0,001.0sp2.04.06.0ps解：解：6010lg20lg203A65326.5)7.860(1102.0,22.3712.0285.2/)860(N取取38将将N=38,=5.653代入代入 表达式，得表达式，得)(nWk)()

36、()653.5()37(3065.0()(0000IxIInnInWknx)(nWk)()(00IxI)(0 xI0 37 0.0 1.000 0.0204 0.021 36 1.8336 2.030 0.0415 0.042 35 2.5568 3.345 0.0704 0.078 29 4.6548 19.96 0.4082 0.413 34 3.086 5.251 0.1074 0.11 4 33 3.5111 7.441 0.1522 0.155 32 3.8656 10.11 0.2067 0.216 31 4.1678 13.10 0.2679 0.297 30 4.4286 16

37、.44 0.3362 0.3417 20 5.6350 48.03 0.9822 0.98nx)(nWk)()(00IxI)(0 xI9 28 4.8512 23.83 0.4873 0.4910 27 5.0215 27.73 0.5671 0.5711 26 5.1682 31.72 0.6489 0.6512 25 5.2931 35.33 0.7225 0.72 13 24 5.3980 39.01 0.7978 0.8014 23 5.4838 41.93 0.8575 0.8615 22 5.5515 44.67 0.9135 0.9116 21 5.6017 46.74 0.95

38、58 0.9618 19 5.6515 48.90 1.0 1.000481216181925293337215.02/)4.06.0(2/)(psc)(2sin)()()()(sin)(00nWyyIxInnnhkcyy2sinyy2sin)(nWkn012345637363534333231-0.01220.01290.0139-0.01458-0.015590.016940.018480.020.040.070.110.150.210.27-0.000240.0005160.00096-0.0016-0.00230.00350.0049 yy2sin)(nh 78910111213143

39、029282726252423-0.01965-0.021520.02379-0.02659-0.03013-0.034770.041090.050220.340.410.490.570.650.720.800.86-0.0067-0.00880.0120.015-0.0196-0.0250.03290.0431516171822212019-0.06451-0.090400.15070.45200.910.960.981.00-0.059-0.0870.1480.45)(nh的图形如下所示的图形如下所示7-5、频率取样设计法、频率取样设计法一、设计思想一、设计思想窗函数设计法是从时域出发，把

40、理想的窗函数设计法是从时域出发，把理想的 用一定用一定形状的窗函数截取成有限长的形状的窗函数截取成有限长的 ，以，以 来近似来近似 )(nhd)(nhd)(nh)(nh)(jdeH)(jeH从而使频响从而使频响 近似理想频响近似理想频响 。频率取样法是从频域出发，对理想的频响频率取样法是从频域出发，对理想的频响 )(jdeH进行等间隔取样，以有限个频响采样去近似理想频响进行等间隔取样，以有限个频响采样去近似理想频响)(jdeH，即：，即：，)()(2kHeHdkNjd等间隔取样等间隔取样并且并且1,.,1,0),()(NkkHkHd二、利用二、利用N个频域采样值重构个频域采样值重构FIR的系统

41、函数与频响的系统函数与频响1.重构重构FIR的的单位抽样响应的的单位抽样响应h(n)根据频域抽样理论（根据频域抽样理论（p99），由），由N个频域采样点个频域采样点可以唯一确定可以唯一确定h(n),即对即对 H(k)进行进行IDFT1,.,1,0,)(1)(10/2NnekHNnhNkNnkj2.重构系统函数重构系统函数H(Z)1101/21010/21010/2101011)(111)(1)(1)(1)()(ZWZkHNZeZkHNZekHNZekHNZnhZHkNNNkNnkjNNknNnNnkjNknNkNnkjNnNnnNjNeW/23.FIR的频响的频响jeZ)(ZH将将 代入代入

42、表达式可得表达式可得)()(2/)/2sin()2/sin()(11)1)(1)(10)21(1010/2jkNkNkNjNkNkjNnkjNjjekHeNkNkHNeeekHNeH其中其中,)21(2/)/2sin()2/sin(1)(NkNjjkeNkNNe为大家所知的内插函数为大家所知的内插函数.分析分析 可知，当可知，当 时（采样点）时（采样点）)(jke1,.,1,0,2NiiN有：有：1,.,1,0,0,1)(2NikikieNijk这说明，重构的频响这说明，重构的频响 ，在采样上严格等于，在采样上严格等于H(k),)(jeH而在采样点之间，频响则由加权的内插函数延伸叠加而成。而在

43、采样点之间，频响则由加权的内插函数延伸叠加而成。三、线性相位的约束条件三、线性相位的约束条件以以h(n)为偶对称，为偶对称，N为奇数的情况进行分析为奇数的情况进行分析.1.FIR的频响具有线性相位的一般表达式的频响具有线性相位的一般表达式当当h(n)为偶对称，为偶对称，N为奇数时，则为奇数时，则)21()()(NjjeHeH（P191，表，表6-1）)(H而且幅度函数而且幅度函数 应为偶对称，即应为偶对称，即)2()(HH2.采样值采样值H(k)具有线性相位的约束具有线性相位的约束kkjkjkNjeHekNHeHkH)2()()(2其中，其中，表示采样值的模（纯标量），表示采样值的模（纯标量）

44、，表示表示)2(kNHHkk其相角。因此，在采样点上具有线性相位的条件应为：其相角。因此，在采样点上具有线性相位的条件应为：)11(221NkkNNk而且，而且，必须满足偶对称，即必须满足偶对称，即kNkHHkH四四.设计步骤设计步骤1.根据指标要求，根据指标要求，画出频率采样序列的图形；画出频率采样序列的图形；kH2.依据依据 的对称特点，可以使问题得以简化；的对称特点，可以使问题得以简化；3.根据线性相位的约束条件，求出根据线性相位的约束条件，求出 ；4.将将 代入代入FIR的频响表达式；的频响表达式；5.由由 的表达式画出实际的表达式画出实际 频响。频响。kkjkeHkH)(kH四四.设

45、计举例设计举例例例 试用频率采样法，设计一个具有线性相位试用频率采样法，设计一个具有线性相位的低通的低通FIR数字数字filter，其理想频率特性为：，其理想频率特性为：5.0,05.00,1)(jdeH5.0c已知已知 ，采样点，采样点N=33.由于由于h(n)为偶对称，且为偶对称，且 N=33为奇数，所以为奇数，所以 对于对于是偶对称。所以上图可画一半（到是偶对称。所以上图可画一半（到 ）kH截止频率截止频率5.0c，即，即33183316c解解：相位约束条件：相位约束条件：33/32)11(kNkk而而 为为kH258,03225,80,1kkkHk33/32)(kjjkeeHkHk将将 代入代入FIR的频响，得的频响，得16320)3316(3332320)21(102/)33/2sin()2/33sin(3312/)33/2sin()2/33sin(3312/)/2sin()2/sin()(1)(jkkjkjkNkNjNkjekkeekeNkNkHNeH考虑到考虑到 时，时，所以将负频部分加进去，所以将负频部分加进去258 k有：有：0kH1681)33/2/sin()2/33sin()33/2/sin()2/33sin()2/sin()2/33sin(331)(jkjekkkkeH的图形如下所示：的图形如下所示：)(jeH)(jeH5.00