FIR数字滤波器的设计方法

《FIR数字滤波器的设计方法》由会员分享，可在线阅读，更多相关《FIR数字滤波器的设计方法（97页珍藏版）》请在装配图网上搜索。

1、第第7章学习目标章学习目标 掌握线性相位掌握线性相位FIR数字滤波器的特点数字滤波器的特点 掌握窗函数设计法掌握窗函数设计法 理解频率抽样设计法理解频率抽样设计法 理解理解IIR与与FIR数字滤波器的比较数字滤波器的比较本章作业练习本章作业练习 P370:P370:1 1 7 7 9 9第第7章章 FIR数字滤波器的设计方法数字滤波器的设计方法IIR数字滤波器：数字滤波器：可以利用模拟滤波器设计可以利用模拟滤波器设计但相位非线性但相位非线性FIR数字滤波器：数字滤波器：可以严格线性相位，又可任意幅度特性可以严格线性相位，又可任意幅度特性因果稳定系统因果稳定系统可用可用FFT计算计算但阶次比但阶

2、次比IIR滤波器要高得多滤波器要高得多7.1 线性相位线性相位FIR滤波器的特点滤波器的特点 FIR滤波器的单位冲激响应：滤波器的单位冲激响应：()01h nnN10()()NnnH zh n z系统函数：系统函数：在在 z 平面有平面有N 1 个零点个零点在在 z=0 处是处是N 1 阶极点阶极点 h(n)为实序列时，其频率响应：为实序列时，其频率响应：1、线性相位条件、线性相位条件()()jjH ee 即群延时 是常数()dd 0()第二类线性相位：第二类线性相位：()第一类线性相位：第一类线性相位：10()()Njj nnH eh n e()()jHe 线性相位是指线性相位是指 是是 的

3、线性函数的线性函数 ()()jjH ee 10()()Njj nnH eh n e()第一类线性相位：第一类线性相位：()jjH ee 10()coscosNjnH eh nn 10()sinsinNjnH eh nn 1010sinsincoscosNnNnh nntgh nn 1100sincoscossin0NNnnh nnh nn 10sin0Nnh nn 第一类线性相位第一类线性相位 的充要条件：的充要条件：()()(1)01h nh NnnN 12Nn=(N 1)/2 为为h(n)的偶对称中的偶对称中心心 10sin0Nnh nn 第二类线性相位第二类线性相位 的充要条件：的充要条

4、件：0()()(1)01h nh NnnN 12N0/2 n=(N 1)/2 为为h(n)的奇对称中心的奇对称中心2、线性相位、线性相位FIR滤波器频率响应的特点滤波器频率响应的特点1100()()(1)NNnnnnH zh n zh Nn z 1(1)0()NNmmh m z(1)1()NzH z 1mNn 令系统函数：系统函数：()(1)01h nh NnnN 由由1(1)0()NNmmzh m z(1)11()()()2NH zH zzH z得11(1)001()()2NNnNnnnh n zzh n z1(1)01()2NnNnnh nzzz11122120()2NNnnNNnzzzh

5、 n(1)1()NH zzH z 由11221cos 221sin 2jNNnnz eNnzzNjn 11122120()2NNnnNNnzzH zzh n112011201()cos2()()1()sin2jNNjnjz eNNjnNeh nnH eH zNjeh nn cos2jxjxeex 频率响应：频率响应：()(1)h nh Nn 11201()()()cos2jNNjjz enNH eH zeh nn12N1）h(n)偶对称偶对称为线性相位，称第一类线性相位为线性相位，称第一类线性相位1()2N 相位函数：相位函数：频率响应：频率响应：()(1)h nh Nn 11201()()(

6、)sin2jNNjjz enNH eH zjeh nn12N112201()sin2NNjjnNeh nn0/22）h(n)奇对称奇对称1()22N 相位函数：相位函数：为线性相位，称第二类线性相位为线性相位，称第二类线性相位3 3、幅度函数的特点、幅度函数的特点1）h(n)偶对称，偶对称，N为奇数为奇数11cos(1)cos22NNNnn 11cos22NNn对呈偶对称101()()cos2NnNHh nn幅度函数：幅度函数：1cos2Nn-32011()2()cos22NnNNHhh nn121112cos()22NmNNhhmm12Nnm令120()()cos()NnHa nn1(0)2

7、Nah其中：其中：11,.,2Nn1()22Na nhn120()()cos()NnHa nn1(0)2Nah11,.,2Nn其中：其中：1()22Na nhn()0,2 H对呈偶对称cos()0,2 n对，呈偶对称120()()cos()NnHa nn可以设计低通、高通、带通、带阻滤波器可以设计低通、高通、带通、带阻滤波器2）h(n)偶对称，偶对称，N为偶数为偶数12012()cos2NnNh nn101()()cos2NnNHh nn幅度函数：幅度函数：2112cos22NmNhmm2Nnm令/211()()cos2NnHb nn()22Nb nhn1,.,2Nn 其中：其中：1201()

8、2()cos2NnNHh nn/211()()cos2NnHb nn()22Nb nhn1,.,2Nn 其中：其中：()H对呈奇对称()01Hz 则是零点1 cos02n时 1z 为零点 故不能设计成高通、带阻滤波器故不能设计成高通、带阻滤波器 ()0,2H对呈偶对称/211()()cos2NnHb nn3）h(n)奇对称，奇对称，N为奇数为奇数11sin(1)sin22NNNnn 11sin22NNn对呈奇对称101()()sin2NnNHh nn幅度函数：幅度函数：1sin2Nn-3201()2()sin2NnNHh nn12112sin()2NmNhmm12Nnm令121()()sin(

9、)NnHc nn1()22Nc nhn11,.,2Nn其中：其中：1()02Nh nNh奇对称且 为奇数121()()sin()NnHc nn1()22Nc nhn11,.,2Nn其中：其中：()0,2H故对，呈奇对称()01Hz 则是零点0,2 sin()0n时 121()()sin()NnHc nn sin()0,2 n因对，呈奇对称只能设计带通滤波器只能设计带通滤波器4）h(n)奇对称，奇对称，N为偶数为偶数101()()sin2NnNHh nn幅度函数：幅度函数：12012()sin2NnNh nn1201()2()sin2NnNHh nn2112sin22NmNhmm2Nnm令/21

10、1()()sin2NnHd nn()22Nd nhn1,.,2Nn 其中：其中：/211()()sin2NnHd nn()22Nd nhn1,.,2Nn 其中：其中：()01Hz 则是零点10,2 sin02n时 ()0,2H对呈奇对称 h(n)为奇对称时，有为奇对称时，有900相移，适用于微分器和相移，适用于微分器和900移相器。不能设计低通、带阻滤波器。移相器。不能设计低通、带阻滤波器。()H对呈偶对称/211()()sin2NnHd nn4、零点位置、零点位置()0iH z*,1/iizz即 也是零点(1)1()()NH zzH z 得：得：由由1）若）若 z=zi 是是H(z)的零点，

11、则的零点，则 z=zi-1 也是零点也是零点2）h(n)为实数，则零点共轭成对为实数，则零点共轭成对线性相位滤波器的零点是互为倒数的共轭对线性相位滤波器的零点是互为倒数的共轭对 即共轭成对且镜像成对即共轭成对且镜像成对(1)1()()0NiiiH zzH z 11()(1)(1)iijjiiiH zrezrez111111iijjiiezezrr1222112 cosiiiirzr zr2122 cosiiirrzz1522NN10ijiiiizrer 或1）11iiiijjjjiiiirereeerr零点：零点：12221()12 cosiiiiiH zrzr zr2122 cosiiirr

12、zz11()11iijjiH zezez1212 cosirzz 1312NN10ijiiiizrer 或2），即零点在单位圆上，即零点在单位圆上iijjee零点：零点：111()11iiiH zrzzr1211iirzzr 1312NN i负实轴上 0 i正实轴上10ijiiiizrer 或3），即零点在实轴上，即零点在实轴上1iirr零点：零点：1()(1)iH zz11222NN 1 01iizz 10ijiiiizrer 或4）即零点既在实轴上，又在单位圆上即零点既在实轴上，又在单位圆上1零点：零点：7.2 窗函数设计法窗函数设计法 7.2.1 设计方法设计方法10()()()Njj

13、njdnH eh n eHe1()2jj nddh nHeed()()()dh nw n h nw(n)：窗函数序列：窗函数序列要选择合适的形状和长度要选择合适的形状和长度2、设计过程、设计过程 )()()()(21)(nhnwnhdeeHnhdnjjdd以低通滤波器为例讨论：以低通滤波器为例讨论：线性相位理想低通滤波器的频率响应：线性相位理想低通滤波器的频率响应：()0,jjccdcceHe 1sin()()2()ccjj nccdcnh needn其理想单位抽样响应：其理想单位抽样响应：中心点为中心点为 的偶对称无限长非因果序列的偶对称无限长非因果序列 取矩形窗：取矩形窗：()()Nw n

14、Rn 1sin201 120cccNnnNh nNnn其它()01()()()0ddh nnNh nh n w nn其它则则FIR滤波器的单位抽样响应：滤波器的单位抽样响应：12N按第一类线性相位条件，得按第一类线性相位条件，得sin()()()ccdcnh nn加窗处理后对频率响应的影响：加窗处理后对频率响应的影响：时域乘积相当于频域卷积时域乘积相当于频域卷积1()2jjjdH eHeW ed 1120sin2()()sin2NNjjj nRnNW ew n ee而矩形窗的频率响应：而矩形窗的频率响应：()()()dh nhn w nsin2()sin2RNW其幅度函数：理想滤波器的频率响应

15、：理想滤波器的频率响应：12()()NjjddHeHe11()221()()()2NNjjjdRH eHeWed 121()()2NjdReHWd1()0cdcH其幅度函数：其幅度函数：则则FIR滤波器的频率响应：滤波器的频率响应：1()()2dRHHWd其幅度函数：(0)()RHW近似于的全部积分面积()0.5(0)cHH2cHN为最大值，正肩峰2cHN为最小值，负肩峰()H随，绕零值波动()(0)HH随，绕波动0 c 2cN 2cN 2cN 2cN 1()()()2dRHHWd幅度函数：幅度函数：加窗函数的影响：加窗函数的影响：不连续点处边沿加宽形成过渡带，其宽度（两肩不连续点处边沿加宽形

16、成过渡带，其宽度（两肩峰之间的宽度）等于窗函数频率响应的主瓣宽度。峰之间的宽度）等于窗函数频率响应的主瓣宽度。在在 处出现肩峰值，两侧形成起伏振处出现肩峰值，两侧形成起伏振 荡，振荡的幅度和多少取决于旁瓣的幅度和多少荡，振荡的幅度和多少取决于旁瓣的幅度和多少2cN 改变改变N只能改变窗谱的主瓣宽度，但不能改变主只能改变窗谱的主瓣宽度，但不能改变主瓣与旁瓣的相对比例。其相对比例由窗函数形状瓣与旁瓣的相对比例。其相对比例由窗函数形状决定，称为决定，称为吉布斯吉布斯(Gibbs)效应。效应。7.2.2 各种窗函数各种窗函数 窗函数的要求：窗函数的要求：窗谱主瓣尽可能窄以获得较陡的过渡带窗谱主瓣尽可能

17、窄以获得较陡的过渡带 尽量减少窗谱最大旁瓣的相对幅度尽量减少窗谱最大旁瓣的相对幅度以减小肩峰和波纹以减小肩峰和波纹()()Nw nRn 矩形窗矩形窗4N主瓣宽度最窄：主瓣宽度最窄：旁瓣幅度大旁瓣幅度大 1120()()NNjjj nRRnWew n eWe窗谱：窗谱：sin2()sin2RNW幅度函数：幅度函数：三角形（三角形（Bartlett）窗）窗21012()212112nNnNw nnNnNN8N主瓣宽度宽：主瓣宽度宽：旁瓣幅度较小旁瓣幅度较小 12()NjjW eWe窗谱：窗谱：2sin24()sin2NWN幅度函数：幅度函数：1N 汉宁（汉宁（Hanning）窗）窗（升余弦窗）（升

18、余弦窗）12()1cos()21Nnw nRnN8N主瓣宽度宽：主瓣宽度宽：旁瓣幅度小旁瓣幅度小22()0.5()0.25RRRWWWWNN幅度函数：幅度函数：1N 海明（海明（Hamming）窗）窗（改进的升余弦窗）（改进的升余弦窗）2()0.540.46cos()1Nnw nRnN8N主瓣宽度宽：主瓣宽度宽：旁瓣幅度更小旁瓣幅度更小22()0.54()0.23RRRWWWWNN幅度函数：幅度函数：1N 布莱克曼（布莱克曼（Blackman）窗）窗（二阶升余弦窗）（二阶升余弦窗）24()0.420.5cos0.08cos()11Nnnw nRnNN12N主瓣宽度最宽：主瓣宽度最宽：旁瓣幅度最

19、小旁瓣幅度最小22()0.42()0.25RRRWWWWNN440.04RRWWNN幅度函数：幅度函数：1N 凯泽（凯泽（Kaiser）窗）窗 ：第一类变形零阶：第一类变形零阶 贝塞尔函数贝塞尔函数0()I 2002111()()nINw nI01nN改变 可同时调整主瓣 宽度和旁瓣幅度 旁瓣幅度 但主瓣宽度窗函数窗函数窗谱性能指标窗谱性能指标加窗后滤波器性能指标加窗后滤波器性能指标旁瓣峰值旁瓣峰值/dB主瓣宽度主瓣宽度过渡带宽过渡带宽阻带最小衰减阻带最小衰减/dB矩形窗矩形窗三角形窗三角形窗汉宁窗汉宁窗海明窗海明窗布拉克曼窗布拉克曼窗凯泽窗凯泽窗-13-25-31-41-57-5724446

20、0.92.13.13.35.55-21-25-44-53-74-80/2/N/2/N7.865阻带最小衰减只由窗形状决定阻带最小衰减只由窗形状决定过渡带宽则与窗形状和窗宽过渡带宽则与窗形状和窗宽N都有关都有关7.2.3 窗函数法的设计窗函数法的设计 1、设计步骤、设计步骤（1）给定频响函数）给定频响函数（2）求出单位抽样响应）求出单位抽样响应（3）根据过渡带宽度和阻带最小衰减，借助窗函数）根据过渡带宽度和阻带最小衰减，借助窗函数 基本参数表（基本参数表（P342表表7-3）确定窗的形式及）确定窗的形式及N的的 大小，一般大小，一般N要通过试探而最后确定；要通过试探而最后确定；（4）最后求）最后

21、求 及及 (5)求求H(e j)=DTFTh(n)，检验是否满足设计要检验是否满足设计要 求，如不满足，则需重新设计。求，如不满足，则需重新设计。)e(Hjd)e(HF)n(hjd1d)()()(nwnhnhd)e(Hj 2、设计举例、设计举例例：分别利用矩形窗与汉宁窗设计具有线性相位的例：分别利用矩形窗与汉宁窗设计具有线性相位的 FIR 低通滤波器，具体要求：低通滤波器，具体要求：)e(Hjd,0,ecj其他其他;rad1,s12c并画出相应的频响特性并画出相应的频响特性解：（解：（1）由于）由于 是一理想是一理想LF，所以，所以 可以得出可以得出 （2）确定）确定N 由于相位函数由于相位函

22、数 ，所以，所以 呈呈 偶对称，其对称中心为偶对称，其对称中心为 ，因此，因此 )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)的数值，注意偶对称，对称中心的数值，注意偶对称，对称中心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(;

23、01423.012/12sin)24()0(hhhhhhhhhhhhh26785.0)13()11(01497.0)15()9(06104.0)17()7(02987.0)19()5(01457.0)21()3(02893.011/11sin)23()1(hhhhhhhhhhhh)n(hn1224由于由于h(n)为偶对称，为偶对称，N=25为奇数，所以为奇数，所以)(H1212/)1(12/)1(0)cos()12(2)12()cos()21(2)21()cos()(nnnNnnnhhnnNhNhnna例如例如 ，可以计算，可以计算H()的值，的值，如下图如下图（4）加汉宁窗）加汉宁窗 由于由

24、于 可以求出序列的各点值可以求出序列的各点值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通过通过 可求出加窗后的可求出加窗后的h(n)()()(nwnhnhd31831.0)12()12()12(whhd13502.0)14(h)10(h

25、04516.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如，如下图如，如下图7.3 凯泽（凯泽（Kaiser）窗及其滤波器设计）窗及其滤波器设计 上述几种窗函数：矩形窗、汉宁窗、海明窗等，上述几种窗函数：矩形窗、

26、汉宁窗、海明窗等，为了压制旁瓣，是以加宽主瓣为代价的。而且，为了压制旁瓣，是以加宽主瓣为代价的。而且，每每一种窗的主瓣和旁瓣之比是固定不变的一种窗的主瓣和旁瓣之比是固定不变的，而凯泽窗，而凯泽窗可以在主瓣宽度与旁瓣衰减之间自由选择。可以在主瓣宽度与旁瓣衰减之间自由选择。凯泽在凯泽在1966发现，利用第一类零阶修正发现，利用第一类零阶修正（变形）贝塞尔函数可以构成一种近似最佳的（变形）贝塞尔函数可以构成一种近似最佳的窗函数。凯泽窗定义为：窗函数。凯泽窗定义为：1、定义定义其它,010,)()/)(1()(02/120NnInInW2/)1(N其中其中 是一个可自由选择的参数。是一个可自由选择的参

27、数。)(0I为第一类零阶修正贝塞尔函数，为第一类零阶修正贝塞尔函数，可同时调整主瓣宽度与旁瓣可同时调整主瓣宽度与旁瓣;越大，频谱旁瓣越小，而主瓣相应增加；越大，频谱旁瓣越小，而主瓣相应增加；0相当于矩形窗相当于矩形窗;凯泽窗随凯泽窗随 变化的曲线如下图：变化的曲线如下图：2.特点特点通常选择通常选择，它们相当于旁瓣与主，它们相当于旁瓣与主瓣幅度比为瓣幅度比为 3.1%-0.047%；94注：第一类零阶修正贝塞尔函数为注：第一类零阶修正贝塞尔函数为.)!3()2/()!2()2/()2/(1!)2/(1)(26242120 xxxkxxIkk由图可以看出由图可以看出,为对称中心，且是偶对称为对称

28、中心，且是偶对称,2/)1(Nn即即)1()(nNWnWkk1)()()21()(0021IINWnWkNnk3.凯泽经验公式凯泽经验公式 该公式可使该公式可使filter设计人员根据设计人员根据filter的设计指的设计指标，估算出标，估算出 值和值和 N 值。值。且且2222224.0222lg20286.295.721,05021),21(07886.0)21(5842.050),7.8(1102.0psNdBdBdBdB112sp1)(jeH：通带截止频率，由：通带截止频率，由 定定;：阻带截止频率，由：阻带截止频率，由 定定.2)(jeHps过渡带宽度过渡带宽度2/)(psc例：利用

29、凯泽窗设计一例：利用凯泽窗设计一FIR低通低通filter，要求，要求6.0,4.0,001.02sp2.04.06.0ps解：解：6010lg20lg2032265326.5)7.860(1102.04设计举例设计举例,22.362.0286.295.760N取取37将将N=37,=5.653代入代入 表达式，得表达式，得)(nWk)()()653.5()37(3065.0()(0000IxIInnInWknx)(nWk)()(00IxI)(0 xI3 34 3.086 5.251 0.1074 0.11 nx)(nWk)()(00IxI)(0 xI12 25 5.2931 35.33 0.

30、7225 0.72 0481216181925293336215.02/)4.06.0(2/)(psc)(2sin)()()()(sin)(00nWyyIxInnnhkc)(nWk n012345637363534333231-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.0049yy2sin)(nh 78910111213143029282726252423-.01965-.021520.02379

31、-.02659-.03013-.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的图形如下所示的图形如下所示 凯塞凯塞贝塞尔窗参数对滤波器的性能影贝塞尔窗参数对滤波器的性能影响响过渡带宽过渡带宽(/(/N N)通带波纹通带波纹(dB)(dB)阻带最小衰阻带最小衰减减(dB)(dB

32、)2.1202.1203.003.000.270.27-30-303.3843.3844.464.460.08640.0864-40-404.5384.5385.865.860.02740.0274-50-505.5685.5687.247.24 0.008680.00868-60-606.7646.7648.648.64 0.002750.00275-70-707.8657.865 10.0 10.0 0.0008680.000868-80-808.9608.96011.411.4 0.0002750.000275-90-9010.05610.05610.810.8 0.0000870.00

33、0087-100-1005、线性相位、线性相位FIR低通滤波器的设计低通滤波器的设计/2/0.2ppspsf /2/0.4ststsstsf 250dB解解：1）求数字频率）求数字频率例：设计一个线性相位例：设计一个线性相位FIR低通滤波器，低通滤波器，给定抽样频率给定抽样频率 ，421.5 10(/sec)srad 321.5 10(/sec)prad 通带截止频率通带截止频率 ，323 10(/sec)strad 阻带起始频率为阻带起始频率为 ，阻带衰减不小于阻带衰减不小于-50dB，幅度特性如图所示，幅度特性如图所示2）求hd(n)()0,jjccdcceHe ccsf()11()22c

34、cjj njndh needed1sin()()ccnnnn12N1/220.3psts 2()0.540.46cos()1Nnw nRnN20.2stps 1162N4）确定）确定N 值值250dB3）选择窗函数：由）选择窗函数：由 确定海明窗（确定海明窗（-53dB）6.6N海明窗带宽：33,6.62.0NN5）确定）确定FIR滤波器的滤波器的h(n)()()()dh nh n w n33sin 0.3160.540.46cos()1616nnRnn6）求）求 ，验证，验证()jH e若不满足，则改变若不满足，则改变N或窗形状重新设计或窗形状重新设计6、线性相位、线性相位FIR高通滤波器的

35、设计高通滤波器的设计()()1()2ccjnjndh neded1sinsin()1()ccnnnnn()=()cc高通滤波器全通滤波器低通滤波器其单位抽样响应：其单位抽样响应：12N()0jjcdeHe其它理想高通的频响：理想高通的频响：7、线性相位、线性相位FIR带通滤波器的设计带通滤波器的设计1221()()1()2jnjndh neded21211sinsin1nnnnn1221(,)=()()带通滤波器低通滤波器低通滤波器其单位抽样响应：其单位抽样响应：120()0jjdeHe其它理想带通的频响：理想带通的频响：12N8、线性相位、线性相位FIR带阻滤波器的设计带阻滤波器的设计211

36、2()()()1()2jnjnjndh nededed12121sinsinsin1nnnnnn1221()=()()带阻滤波器，高通滤波器+低通滤波器其单位抽样响应：其单位抽样响应：12N120,()0jjdeHe 其它理想带阻的频响：理想带阻的频响：7.4 频率抽样设计法频率抽样设计法 1、设计方法、设计方法 对理想频率响应等间隔抽样对理想频率响应等间隔抽样 作为实际作为实际FIR数字滤波器的频率特性的抽样值数字滤波器的频率特性的抽样值2()()()jddkNH kHkHe0,1,.,1kN()h n()H z()jH e10()()()Njj njdnH eh n eHe1()()2jj

37、 nddh nHeed()()()dh nw n h n窗函数设计法：窗函数设计法：窗函数设计法是从时域出发，把理想的窗函数设计法是从时域出发，把理想的 用一定用一定形状的窗函数截取成有限长的形状的窗函数截取成有限长的 ，以，以 来近似来近似 )(nhd)(nhd)(nh)(nh)(jdeH)(jeH从而使频响从而使频响 近似理想频响近似理想频响 。频率取样法是从频域出发，对理想的频响频率取样法是从频域出发，对理想的频响 进行等间隔取样，以有限个频响采样去近似理想频响进行等间隔取样，以有限个频响采样去近似理想频响)(jdeH，)(jdeH1101()()1NNkkNzH kH zNWz 102

38、NjkH eH kkN12sin12()sin2NjNeN 11(1)20sin21sin2kNNjNjjNkkNNH eeH kekNN内插公式：内插公式：110Int824()10Int191622ccNNkH kNNk 例：利用频率抽样法设计一个频率特性为矩形的理想例：利用频率抽样法设计一个频率特性为矩形的理想低通滤波器，截止频率为低通滤波器，截止频率为，抽样点数为，抽样点数为N=33，要，要求滤波器具有线性相位。求滤波器具有线性相位。解：解：100cjdHe理想低通频率特性其它按频率抽样方式，按频率抽样方式，N=33，得抽样点，得抽样点得线性相位得线性相位FIR滤波器的频率响应：滤波器

39、的频率响应：816133sin 33sin 33sin233233233sin33sin33sin2233233jjkkkH eekk2/33过渡带宽：过渡带宽：阻带衰减：阻带衰减：-20dB 过渡带抽样的优化设计过渡带抽样的优化设计 增加一点过渡带抽样点增加一点过渡带抽样点令令H4/33过渡带宽：过渡带宽：阻带衰减：阻带衰减：-40dB 增加两点过渡带抽样点增加两点过渡带抽样点且增加抽样点数为且增加抽样点数为N=65令令H H6/65过渡带宽：过渡带宽：阻带衰减：阻带衰减：-60dB 增加过渡带抽样点，可加大阻带衰减，但导致增加过渡带抽样点，可加大阻带衰减，但导致过渡带变宽过渡带变宽 增加增

40、加N，使抽样点变密，减小过渡带宽度，但，使抽样点变密，减小过渡带宽度，但增加了计算量增加了计算量7.6 IIR和和FIR数字滤波器的比较数字滤波器的比较IIR滤波器滤波器FIR滤波器滤波器 h(n)无限长无限长 h(n)有限长有限长 极点位于极点位于z平面任意位置平面任意位置 滤波器阶次低滤波器阶次低 非线性相位非线性相位 递归结构递归结构 不能用不能用FFT计算计算 可用模拟滤波器设计可用模拟滤波器设计 用于设计规格化的选频用于设计规格化的选频滤波器滤波器 极点固定在原点极点固定在原点 滤波器阶次高得多滤波器阶次高得多 可严格的线性相位可严格的线性相位 一般采用非递归结构一般采用非递归结构 可用可用FFT计算计算 设计借助于计算机设计借助于计算机 可设计各种幅频特性和可设计各种幅频特性和相频特性的滤波器相频特性的滤波器