第七章FIR滤波器设计ppt课件

1、 第7章有限冲激响应滤波器FIR的设计 71 线性相位线性相位FIR滤波器的特点滤波器的特点 7. 2 窗函数设计法窗函数设计法7. 3 频率抽样设计法频率抽样设计法7. 4 应用应用MATLAB设计设计FIR数字滤波器数字滤波器有限长单位冲激响应数字滤波器的特点有限长单位冲激响应数字滤波器的特点:有限长单位冲激响应有限长单位冲激响应FIR可以做成具有严格的可以做成具有严格的线性相位，同时又可以具有任意的幅度特性。线性相位，同时又可以具有任意的幅度特性。FIR滤波器的单位抽样响应是有限长的，因而滤波器的单位抽样响应是有限长的，因而FIR滤波器一定是稳定的。滤波器一定是稳定的。只要经过一定的延时

2、，任何非因果有限长序列都只要经过一定的延时，任何非因果有限长序列都能变成因果的有限长序列，总能用因果系统来能变成因果的有限长序列，总能用因果系统来实现。实现。FIR滤波器由于单位冲激响应是有限长的，因而可滤波器由于单位冲激响应是有限长的，因而可以用快速傅里叶变换以用快速傅里叶变换FFT算法来实现过滤算法来实现过滤信号，从而可大大提高运算效率。信号，从而可大大提高运算效率。但是，要取得很好的衰减特性，但是，要取得很好的衰减特性，FIR滤波器的阶次滤波器的阶次比比IIR滤波器的要高。滤波器的要高。本章主要讨论线性相位滤波器的设计。本章主要讨论线性相位滤波器的设计。7. 1 线性相位线性相位FIR滤

3、波器的特点滤波器的特点7.1.1 线性相位条件线性相位条件如果一个线性移不变系统的频率响应有如下形式： （7.1）则其具有线性相位。这里 是一个实数。因而，线性相位系统有一个恒定的群延时 （7.2）()()()|() |jjjjH eHeH ee 位函数。分别称为幅度函数与相与)()(H在实际应用中，有两类准确的线性相位，分别要求满足 （7.3） （7.4）FIR滤波器具有式7.3的线性相位的充分必要条件是：单位抽样响应 关于群延时 偶对称，即满足 （7.5）（7.6）( ) ( ) ( )h n12N( )(1) 01h nh NnnN 满足式7.5和式7.6的偶对称条件的FIR滤波器分别称

4、为I型N为奇数线性相位滤波器和型N为偶数线性相位滤波器。 (7.3)(7.5)(7.6)的证明：将7.6式代入1010) 1()( then ,)()(NnnNnnznNhzHznhzH令令 m=N-n-1,m=N-n-1,则有则有)()()()()(1)1(10)1(10)1(zHzzHzmhzzmhzHNNmmNNmmN于是于是1021212110)1(1)1(21)( )(21)()(21)(NnNnNnNNnnNnNzznhzzzznhzHzzHzH将将 代入上式代入上式jez1021jj)21cos()(e)e (NnNNnnhH幅度函数与相位函数分别为10)21cos()()(Nn

5、NnnhH21)(N12N可见群延时 ，只要h(n)是实序列，且满足（7.5式，该滤波器具有第一类线性相位。FIR滤波器具有式7.4的线性相位的充分必要条件是：单位抽样响应 关于群延时 奇对称，即满足 （7.7）（7.8）（7.9）( )h n12N2 ( )(1) 01h nh NnnN 把满足式7.7）、（7.8和式7.9的奇对称条件的FIR滤波器分别称为型线性相位滤波器和型线性相位滤波器 。(7.7)-(7.9)的证明：将7.9式代入1010) 1()( then ,)()(NnnNnnznNhzHznhzH令令 m=N-n-1,m=N-n-1,则有则有)()()()()(1)1(10)

6、1(10)1(zHzzHzmhzzmhzHNNmmNNmmN于是于是1021212110)1(1)1(21)( )(21)()(21)(NnNnNnNNnnNnNzznhzzzznhzHzzHzH将将 代入上式代入上式jez10212j10212j1021jj)21sin()(e )21sin()(e )21sin()(je)e (NnNNnNNnNNnnhNnnhNnnhH幅度函数与相位函数分别为101( )( )sin() ,2NnNHh nn212)(N12N可见 ，21I型线性相位滤波器型线性相位滤波器 7.1.2 线性相位滤波器频率响应的特点线性相位滤波器频率响应的特点由于偶对称性，

7、一个I型线性相位滤波器的频率响应可表示为 （7.10）其中(1)/2(1)/20()( )cos()Njj NnH eea kk11( )2 () 1,2,.,22NNa khkk1(0)()2Nah幅度函数为 （7.11）相位函数为 （7.12）(1)/20()( )cos()NnHa kk()(1)2N I型线性相位滤波器的幅度函数和相位函数的特点：幅度函数对 偶对称，同时对 也呈偶对称；相位函数为准确的线性相位。 12N0, ,2证明：1021jj)21cos()(e)e (NnNNnnhH而幅度函数10)21cos()()(NnNnnhH21)(N) 1()(nNhnh相位函数为)21

8、cos()21cos( )211cos(NnNnNnN因为故)21cos()211cos(NnNnN2/ )3(0)21cos()(2)21()(NnNnnhNhH于是当N为奇数时令 m=(N-1)/2-n,则有2/ )1(1cos)21(2)21()(NmmmNhNhH或2/ )1(02/ )1(1cos)( cos)21(2)21()(NkNkkkakkNhNhH式中2/)1(,2, 1 ),21(2)()21()0(NkkNhkaNha2型线性相位滤波器型线性相位滤波器 一个型线性相位滤波器，由于N是偶数，所以， 的对称中心在半整数点 。其频率响应可以表示为： （7.13）其中 ( )h

9、 n12N/2(1)/21()1( )cos() 2Njj NkH eeb kk( )2 () 1,2,.,22NNb khkk幅度函数为 （7.14）相位函数为 （7.15）()(1)2N /21()1( )cos() 2NkHb kk型线性相位滤波器的幅度函数和相位函数的特点：幅度函数的特点：（1当 时， =0，也就是说 在 处必然有一个零点；（2） 对 呈奇对称，对 呈偶对称。相位函数的特点：同I型线性相位滤波器。 ()H( )Hz1z ()H0,2120)21cos()(2)(NnNnnhH当N为偶数时令 m=N/2-n,则有2/1)21cos()2(2)(NmmmNhH或2/12/1

10、)21cos()( )21cos()2(2)(NkNkkkbkkNhH式中2/,2, 1 ),2(2)(NkkNhkb3型线性相位滤波器型线性相位滤波器 由于型线性相位滤波器关于 奇对称，且 为整数，所以，其频率响应可以表示为 （7.16）其中 12N(1)/2(1)/21()( )sin()Njj NkH ejec kk11( )2 () 1,2,.,22NNc khkk幅度函数为 （7.17）相位函数为 （7.18）(1)/21()( )sin()NkHc kk()(1)22N 型线性相位滤波器的幅度函数和相位函数的特点：幅度函数的特点：（1当 时， =0，也就是说 在 处都为零点 ；（2

11、） 对 均呈奇对称。相位函数的特点：既是准确的线性相位，又包括 的相移，所以又称 移相器，或称正交变换网络。 ()H( )Hz()H0, ,21z 0, ,2/2904型线性相位滤波器型线性相位滤波器 型线性相位滤波器关于 奇对称，且N为偶数，所以为非整数。其频率响应可以表示为 （7.19）其中 12N/2(1)/21()1( )sin() 2Njj NkH ejed kk( )2 () 1,2,.,22NNd khkk幅度函数为 （7.20）相位函数为 （7.21）()(1)22N /21()1( )sin() 2NkHd kk型线性相位滤波器的幅度函数和相位函数的特点：幅度函数的特点：（1

12、） 在 处必为零，也就是说 在 处为零点 ；（2） 在 处呈奇对称，在 处呈偶对称 相位函数的特点：同型线性相位滤波器。 ( )Hz()H0,2()H1z 0,27.1.3 零点位置零点位置对于I型或型线性相位滤波器，意味着 ( )(1)h nh Nn (1)1( )()NH zzH z对于型或型线性相位滤波器，意味着 ( )(1)h nh Nn (1)1( )()NH zzH z 在上述两种情况下，假如 在 处等于零，则在处也一定等于零。所以 的零点呈倒数对出现。另外，假设 是实值的，则复零点呈共轭倒数对出现，或者说是共轭镜像的。 ( )H z0zz01/zz( )H z( )h n一个线性

13、相位滤波器零点的约束条件一个线性相位滤波器零点的约束条件 线性相位滤波器的级联结构实现 5N 122N在此情况下， 1111122212211( )(1)(1)(1)(1)1 12 (cos)2 (cos)iiiijjjjiiiiiiiiiiiiH zz rez rezezerrrzr zrrzzr（1零点 既不在实轴上，也不在单位圆上， 零点是两组互为倒数的共轭对，其基本因子为 izz（7.22） ,1,0ijiiiizrer注：此时注：此时H(z)H(z)为为N-1=4N-1=4次多项式。次多项式。在此情况下， 3N 112N（2零点 在单位圆上，但不在实轴上，即 ，此时倒数是其本身，只有

14、共轭对出现，基本因子为 izz（7.23） 1112( )(1)(1)1 2(cos)iijjiiH zz ez ezz 1,0,iiir在此情况下， 3N 112N（3零点 在实轴上，但不在单位圆上，即 ，此时零点是实数，没有复共轭部分，只有倒数，倒数也在实轴上，其基本因子为 （7.24） 111211( )(1)(1)1()iiiiiH zrzzrzzrr 式中“ ”号相当于 ，零点在负实轴上；“ ”相当于零点在正实轴上。 i0iizz1 ; 0,iirizz（4零点 既在单位圆上，又在实轴上，即 ，此时零点只有两种情况，即 ，这时零点既是自己的复共轭，又是倒数，其基本因子为 izz（7.

15、25） 1,0iir 或1( )1iH zz 式中“ ”号相当于 ，零点在负实轴上；“ ”相当于零点在正实轴上。 i0i在此情况下， ， ，即有半个抽样的延时。 2N 11/22N线性相位滤波器只能由以上这几种因子的组合而成。 11zz 或7. 2 窗函数设计法窗函数设计法 7. 2.1 设计方法设计方法 给出所要求的理想低通滤波器频率响应()jdeH设计一个FIR滤波器频率响应 10()( )Njj nneeHh n迫近 ()jdeH1 ()0 jjcdceHe窗函数设计法设计是在时域进行 （7.26） ()jdeH由 的傅里叶反变换导出 ()sin()11( )()22()ccjj njn

16、cddeeenhnHddn( )( ) ( )dh nh n w n（7.27） 窗函数序列的形状及长度的选择很关键。 为得到长度为N的因果且线性相位FIR滤波器h(n),必须有( ), 011( ) 20, ndh nnNNh n和 =其余上式是对hd(n)加了一个矩形窗：( )( )( )dNh nh n Rn设w(n)是关于(N-1)/2对称的长度为N的一个窗函数，一般有取 ，即取矩形窗 例例7.1 设计一低通滤波器，所希望的频率响应截止频率设计一低通滤波器，所希望的频率响应截止频率 在在 之间为之间为1，在，在 之间为之间为0，分别取，分别取N= 11,21,41，观察其频谱响应的特点

17、。，观察其频谱响应的特点。 ()jdeH00.250.25解：解： (1)/2 00.25()0 0.25j NjdeeH1 01( )0 nNn其它 由式7.26),(7.27） 1sin0.25()2( )( ) ( )1()2dNnh nh n w nNn（线性相位）当N=11时，求得 (0)(10)0.045, (1)(9)0, (2)(8)0.075(3)(7)0.1592, (4)(6)0.2251, (5)0.25hhhhhhhhhhh 显然 ，满足对称关系。152N根据序列 ，分别求得N=11,21,41时的幅频特性 ( )h n|()|jeH由图可以看出，当N取的过小时，通频

18、带过窄，且阻带内波纹较大，过渡带较宽，当N增大时， 与 的近似程度越来越好。但当N增大时，通带内出现了波纹，而且随着N的继续增大，这些波纹并不消失，只是最大的尖峰处越来越接近于间断点，这种现象称作吉布斯现象。 ()jeH()jdeH吉布斯现象的产生是由于对 突然截短的结果。 ( )dh n为了减少吉布斯现象，应选取旁瓣较小的窗函数。 1矩形窗矩形窗 7.2.2 各种窗函数各种窗函数 窗函数为 ( )( )Nw nRn（7.28） 幅度函数为 （7.29） sin()2( ) |()|sin()2jRRNWW e主瓣宽度 ，过渡带宽 。 4 /2 2 /NN0.9 2 / N2汉宁汉宁Hanni

19、ng窗又称升余弦窗）窗又称升余弦窗） 窗函数为 2( )0.50.5cos()( )1Nnw nRnN（7.30） 幅度函数为 （7.31） 22( )0.5( )0.25()()11RRRnnWWWWNN主瓣宽度 ，过渡带宽 。 4 2 /8 /NN3.1 2 / N3海明海明Hamming窗又称改进的升余弦窗）窗又称改进的升余弦窗） 窗函数为 （7.32） 幅度函数为 （7.33） 主瓣宽度 ，过渡带宽 。 4. 4. 凯泽凯泽KaiserKaiser窗窗 窗函数为 （7.34） 2( )0.540.46cos()( )1Nnw nRnN22( )0.54( )0.23()()11RRRn

20、nWWWWNN4 2 /8 /NN3.3 2 / N20021 (1)1( ), 01( )nINw nnNI其中 为第一类变形零阶贝塞尔函数， 是一个可自由选择的参数，改变 值就可对主瓣宽度与旁瓣衰减进行选择，一般选择 。过渡带宽 。 0I495 2 / N 窗函数窗函数窗谱性能指标加窗后滤波器性能指标旁瓣峰值（dB）主瓣宽度( )过渡带宽( )阻带最小衰减（dB）矩形窗汉宁窗海明窗凯泽窗-13-31-41-57244 5-21-44-53-80最小阻带衰减只由窗形决定，不受N的影响，而过渡带宽则随N的增加而减小。 表表7.1 7.1 几种窗函数的基本参数比较几种窗函数的基

21、本参数比较2/ N2/ N1高通数字滤波器的设计高通数字滤波器的设计 7.2.3 其他各型其他各型FIR滤波器的设计方法滤波器的设计方法 令 （7.35） (1)/2 ()0 0j NjcdceeH那么 11()()2211( )22ccNNjnjndeeh ndd求得 11sinsin22( )1()2cdNNnnh nNn（7.36） 2带通数字滤波器的设计带通数字滤波器的设计 令 （7.37） 那么 求得 （7.38） (1)/2 ()0 j NjlhdeeH其它11()()2211( )22lhhlNNjnjndeeh ndd11sinsin22( )1()2hldNNnnh nNn3

22、带阻数字滤波器的设计带阻数字滤波器的设计 令 （7.39） 那么 求得 （7.40） (1)/2 |, | |()0 j NjlhdeeH其它 111()()()222111( )222hllhNNNjnjnjndeeeh nddd111sinsinsin222( )1()2lhdNNNnnnh nNn比较式7.36）、（7.38）、（7.40可知，一个高通滤波器相当于用一个全通滤波器减去一个低通滤波器；一个带通滤波器相当于两个低通滤波器相减，其中一个截止频率在 ，另一个在 ；一个带阻滤波器相当于一个低通滤波器加上一个高通滤波器，低通滤波器的截止频率 ，高通滤波器在 。 hlhl选取一个满意的

23、窗函数，令 ( )( )( ), 0,1,.,1dh nw n h nnN（7.41） 那么 即为要设计的滤波器的单位抽样响应。 ( )h n按上述方法设计的滤波器，由于满足了 的对称关系，因此都具有线性相位。 ( )(1)h nh Nn 7. 3 频率抽样设计法频率抽样设计法 频率抽样法是从频域出发，把给定的理想频率响应 加以等间隔抽样，即 ()jdeH2()( )jddkNeHHk令 2( )( )(), 0,1,.,1jddkNeH kHkHkN（7.42） 由DFT定义，得2101( )( ), 0,1,.,1NjnkNdkh nHk ekNN（7.43） 可求得滤波器的系统函数（7.

24、44） 21110002110012011( )( )( )1 ( )11 ( )1NNNjnknnNdnnkNNjnknNdknNNdjkkNH zh n zHk ezNHkezNzHkNez该系统的频率响应为12011()( )( )1jNNjdjjkz ekjNeH eH zHkNee（7.45） 2jkNNNNezz经过推导，有 （7.46） 由式7.46可知， 是由内插函数（7.47） 1(1)/2(1)/0sin(2/)/2()( )sin (2/)/2Njj Nj NkNdkNk NH eeHk eNk N()jeH(1)2sin2( )sin2NjNeN的插值所决定的，即 10

25、2()( ) ()NjdkH eHkkN 由内插公式7.47可知，在各频率抽样点上，滤波器的实际频率响应严格地和理想频率响应值相等。但是在抽样点之间的频率响应则是由N个离散值 作为权重和插值函数 线性组合的结果。显然抽样点N取得越大，近似程度越好，N的选取要视在通带和阻带内的技术要求而定。 ( )dHk( )（7.47） 102()( ) ()NjdkH eHkkN 的指定原则的指定原则( )dHk（3由 求出的 应具有线性相位。（1在通带内可令| |=1，阻带内| |=0，且在通带内赋给 一相位函数；( )dHk( )dHk( )dHk（2指定的 应保证由式7.43求出的 是实序列；( )d

26、Hk( )h n( )h n()jeH 的指定的指定 ( )dHk由式7.46知，若保证 (1)/( )j NkNdHk e实数那么 就具有线性相位， 。 ()jeH1( )2N 并考虑| |=1，等效地指定 ( )dHk(1)/( ), 0,1,.,1j NkNdHkekN（7.48） 根据DFT的性质可知，为保证 是实序列，应满足下列对称关系 ( )h n*( )()()dddHkHkHNk（7.49） （7.46） 1(1)/2(1)/0sin(2/)/2()( )sin (2/)/2Njj Nj NkNdkNk NH eeHk eNk N2*( )( )( )1dddHkHk Hk由于

27、 (1)()/(1)(1)/(1)*()( )j NNkNj Nj NkNj NddeeeeHNkHk（7.50） 当N为偶数时， ；当N为奇数时， 。这样当N为偶数时，若按式7.48对 赋值，就不能满足式7.49的对应关系。由此，按如下原则对 赋值。 (1)1j Ne (1)1j Ne( )dHk( )dHkN为偶数时 (1)/(1)/ 0,1,.,/2 1( )0 /2 /2 1,.,1j NkNdj NkNeekNHkkNkNNN为奇数时 (1)/( ) 0,1,.,1j NkNdeHkkN（7.51） （7.52） *( )()()dddHkHkHNk（7.49） 用频率抽样法设计用频

28、率抽样法设计FIRFIR数字滤波器的步骤：数字滤波器的步骤： （1根据所设计的滤波器的通带与阻带的要求，根据N为偶数还是奇数，按式7.51）、（7.52指定 ，在阻带内， =0； ( )dHk( )dHk（2由指定的 构成所设计的滤波器的转移函数7.44） ，也可由式7.47求得频率响应 。 ( )dHk()jeH以低通滤波器设计为例：设采样点数为N,截止频率为c , c 对应的样本点序号为kc ， 那么 kc 的计算依据如下：设抽样间隔为t, 抽样频率为FC=1/ t (Hz), 频率分辨率为f=1/(N t),角频率分辨率为= 2f =2 /(N t) ， kc= c / = N t c

29、/（ 2）， 表示取整。采样间隔缺省时为t=1。N为偶数时 (1)/(1)/ 0,1,.,( )0 +1, 1 ,.,1j NkNcdccj NkNceekkHkkkNkkNkNN为奇数时 (1)/ 0,1,.,; ,1( )0 j NkNccdekkkNkNHkk为其他值（7.53） （7.54） 例例7.2 用频率抽样法设计一个低通滤波器，其截止频率是抽样频用频率抽样法设计一个低通滤波器，其截止频率是抽样频率的率的1/10，取，取N=20。 解：解： 抽样频率为FC=1/t=1 (Hz)，截止频率为fc=0.1FC =1/10，kc= fc/f=fcNt= fcN=2此处N为偶数，且在通带

30、内对 抽样时，仅得两个点(1/10*N=2)，由式7.53），有 ()jeH(0)1dH19/20(1)jdeH*19/20(19)(20 1)(1)jdddeHHH在其它点处， ( )0dHk *38/20(18)(202)(2)jdddeHHH38/20(2)jdeH将 代入式7.43）得 序列如下 ( )dHk( )h nh(0)=h(19)=0.046337h(1)=h(18)=0.019678h(2)=h(17)=-0.020711h(3)=h(16)=-0.054178h(4)=h(15)=-0.060749h(5)=h(14)=-0.029462h(6)=h(13)=0.0366

31、21h(7)=h(12)=0.12071h(8)=h(11)=0.19788h(9)=h(10)=0.243870.1Hz , 20.2cccff7. 4 应用应用MATLAB设计设计FIR数字滤波数字滤波器器 1窗函数窗函数 （1bartlett.m （三角窗）（三角窗）（2blackman.m （布莱克曼窗）（布莱克曼窗）（3boxcar.m （矩形窗）（矩形窗） （4hamming.m （海明窗）（海明窗）（5hanning.m （汉宁窗）（汉宁窗）（6triang.m （三角窗）（三角窗）（7chebwin.m （切比雪夫窗）（切比雪夫窗）（8kaiser.m （凯泽窗）（凯泽窗）7.

32、4.1 与本章内容有关的与本章内容有关的MATLAB文件文件 （1fir1.m本文件采用窗函数法设计本文件采用窗函数法设计FIR数字滤波器，其调用格式是数字滤波器，其调用格式是1b=fir1N,Wn）2b=fir1N,Wn,high）3b=fir1N,Wn,stop） 2FIR数字滤波器的文件数字滤波器的文件 式中N为滤波器的阶次，因此滤波器的长度为N+1；Wn是通带截止频率，其值在01之间，1对应抽样频率的一半；b是设计好的滤波器系数。 对于格式1），若Wn是一标量，则可用来设计低通滤波器；若Wn是 的向量，则用来设计带通滤波器；若Wn是 的向量，则可用来设计带滤波器，此时，格式将变为：b=

33、fir1N,Wn,DC-1）或b=fir1N,Wn,DC-0） 格式2用来设计高通滤波器；格式3用来设计带阻滤波器。 21L1（2 2fir2.mfir2.m本文件采用窗函数法设计具有任意幅频特性的本文件采用窗函数法设计具有任意幅频特性的FIRFIR数字滤波器。其调用数字滤波器。其调用格式是格式是b=fir1b=fir1N,F,MN,F,M）其中其中F F是频率向量，其值在是频率向量，其值在0101之间，之间，M M是与是与F F相对应的所希望的幅频响相对应的所希望的幅频响应。不指定窗函数的类型，则自动选择汉明窗。应。不指定窗函数的类型，则自动选择汉明窗。 （3remez.m本文件用来设计采用

34、切比雪夫最佳一致逼近本文件用来设计采用切比雪夫最佳一致逼近FIR数字滤波器。同时，数字滤波器。同时，还可以用来设计希尔伯特变换器和差分器。其调用格式是还可以用来设计希尔伯特变换器和差分器。其调用格式是1b=remezN,F,A）2b=remezN,F,A,W）3b=remezN,F,A,W,hilbert）4b=remezN,F,A,W,differentiator）其中，其中，N是给定的滤波器的阶次；是给定的滤波器的阶次；b是设计的滤波器的系数，其长度为是设计的滤波器的系数，其长度为N+1；F是频率向量，其值在是频率向量，其值在01之间；之间；A是对应是对应F的各频段上的理想幅的各频段上的理

35、想幅频响应；频响应；W是各频段上的加权向量。是各频段上的加权向量。 注意：若注意：若b b的长度为偶数，设计高通和带阻滤波器时有可能出现错误，的长度为偶数，设计高通和带阻滤波器时有可能出现错误，因此最好保证因此最好保证b b的长度为奇数，即的长度为奇数，即N N应为偶数。应为偶数。 （4 4remexord.m remexord.m 本文件采用切比雪夫一致逼近设计本文件采用切比雪夫一致逼近设计FIRFIR数字滤波器时所需要的滤波器阶数字滤波器时所需要的滤波器阶次。其调用格式是次。其调用格式是N,Fo,Ao,W=remexordN,Fo,Ao,W=remexordF,A,DEV,FsF,A,DE

36、V,Fs）式中，式中，F F、A A的含义同文件的含义同文件3 3），是通带和阻带上的偏差；该文件输出），是通带和阻带上的偏差；该文件输出的是符合要求的滤波器阶次的是符合要求的滤波器阶次N N、频率向量、频率向量FoFo、幅度向量、幅度向量AoAo和加权向量和加权向量W W。若设计者事先不能确定自己要设计的滤波器的阶次，那么，调用若设计者事先不能确定自己要设计的滤波器的阶次，那么，调用remexordremexord后，就可利用这一族参数再调用后，就可利用这一族参数再调用remezremez，即，即b=remezb=remezN,Fo,Ao,WN,Fo,Ao,W），从而设计出所需要的滤波器。因

37、而，通常），从而设计出所需要的滤波器。因而，通常remezremez和和remexordremexord结合使用。结合使用。 说明：说明：remexordremexord给出的阶次给出的阶次N N有可能偏低，这时适当增加有可能偏低，这时适当增加N N即可；另外，即可；另外，若若N N为奇数，就可令其加为奇数，就可令其加1 1，使其变为偶数，这样，使其变为偶数，这样b b的长度为奇数。的长度为奇数。 （5sgolay.m本文件用来设计本文件用来设计Savitzky-Golay平滑滤波器。其调用格式是平滑滤波器。其调用格式是b=sgolayk,f）式中式中k是多项式的阶次，是多项式的阶次，f是拟合

38、的双边点数。要求是拟合的双边点数。要求 kf ，且，且f为奇数。为奇数。 （6firls.m本文件用最小平方法设计线性相位本文件用最小平方法设计线性相位FIR数字滤波器。可设计任意给定数字滤波器。可设计任意给定的理想幅频特性。的理想幅频特性。（7fircls.m用带约束的最小平方法设计线性相位用带约束的最小平方法设计线性相位FIR数字滤波器。可设计任意给数字滤波器。可设计任意给定的理想幅频特性。定的理想幅频特性。（8fircls1.m用带约束最小平方法设计线性相位用带约束最小平方法设计线性相位FIR低通和高通滤波器。可设计任低通和高通滤波器。可设计任意给定的理想幅频特性。意给定的理想幅频特性。

39、（9firrcos.m用来设计低通线性相位用来设计低通线性相位FIR数字滤波器，其过渡带为余弦函数形状。数字滤波器，其过渡带为余弦函数形状。7.4. 2 应用应用MATLAB设计设计FIR数字滤波器数字滤波器 例例7.3 令令N=10，分别用矩形窗和海明窗重复例，分别用矩形窗和海明窗重复例7.1。 解：解： 根据要求编制MATLAB程序如下：clear all;N=10;b1=fir1(N,0.25,boxcar(N+1); % 用矩形窗作为冲激响应的窗函数b2=fir1(N,0.25,hamming(N+1); % 用Hamming窗作为冲激响应的窗函数%M=128;h1=freqz(b1,1,M);h2=freqz(b2,1,M);% 分别求两个滤波器的频率响应；t=0:10;subplot(221)stem(t,b2,.);hold on;plot(t,zeros(1,11);grid;f=0:0.5/M:0.5-0.5/M;M1=M/4;for