第7章 有限脉冲响应数字滤波器的设计C

1、第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 7.1 线性相位线性相位FIR数字滤波器的条件和特点数字滤波器的条件和特点 7.2 利用窗函数法设计利用窗函数法设计FIR滤波器滤波器7.3 利用频率采样法设计利用频率采样法设计FIR滤波器滤波器7.4 利用等波纹最佳逼近法设计利用等波纹最佳逼近法设计FIR滤波器滤波器7.5 IIR和和FIR数字滤波器的比较数字滤波器的比较第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计)e (Hjd )n(hd)n()n(h)n(hd )e(Hj de )e(H21)n(hnjjdd截截断断逼逼进进一、设计思想一、

2、设计思想nj1N0nje )n(h)e(H 7.2窗函数法设计数字滤波器窗函数法设计数字滤波器第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计（1） 在理想特性不连续点在理想特性不连续点=c附近形成附近形成过渡带过渡带。过。过渡带的宽度近似等于渡带的宽度近似等于WRg()主瓣宽度主瓣宽度4/N。（2） 通带内产生了通带内产生了波纹波纹，最大的峰值在，最大的峰值在c2/N处。处。阻带内产生了阻带内产生了余振余振，最大的负峰在，最大的负峰在c+2/N处。通带与阻处。通带与阻带中波纹的情况与窗函数的幅度谱有关，带中波纹的情况与窗函数的幅度谱有关， WRg()旁瓣幅旁瓣幅度的大小直

3、接影响度的大小直接影响Hg()波纹幅度的大小波纹幅度的大小。以上两点就是对以上两点就是对hd(n)用矩形窗用矩形窗截断截断后，在频域的反映，后，在频域的反映，称为称为吉布斯效应吉布斯效应。这种效应直接影响滤波器的性能。这种效应直接影响滤波器的性能。加窗影响：加窗影响：第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计表7.2.2 6种窗函数的基本参数 第第7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计(1)(1)根据对过渡带及阻带衰减的指标要求，选择窗根据对过渡带及阻带衰减的指标要求，选择窗函数的函数的类型类型，并估计窗口长度，并估计窗口长度N N。 (2)

4、 (2) 构造希望逼近的频率响应函数构造希望逼近的频率响应函数Hd(eHd(ejj) ) (3) (3) 计算计算h hd d(n(n) )7.2.3 用窗函数法设计用窗函数法设计FIR滤波器的步骤滤波器的步骤(4) (4) 加窗得到设计结果：加窗得到设计结果：h(n)=hh(n)=hd d(n)w(n(n)w(n) ) (5) (5) 由由h(nh(n) )求出系统函数求出系统函数H(zH(z) )(1) fir1 fir1用窗函数法设计用窗函数法设计线性相位线性相位FIR数字滤波器的数字滤波器的工具箱工具箱函数函数，以实现，以实现线性相位线性相位FIR数字滤波器的数字滤波器的标准窗函数法标

5、准窗函数法设计。设计。这里的所谓这里的所谓“标准标准”，是指在设计，是指在设计低通、高通、带通和带阻低通、高通、带通和带阻FIR滤波器滤波器时，时，Hd(ej)分别表示相应的线性相位分别表示相应的线性相位理想理想低通、低通、高通、带通和带阻滤波器的频率响应函数。因而将所设计的高通、带通和带阻滤波器的频率响应函数。因而将所设计的滤波器的滤波器的频率响应称为标准频率响应频率响应称为标准频率响应。7.2.4 窗函数法的窗函数法的MATLAB设计函数简介设计函数简介第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 hn=fir1(M, wc) 返回返回6 dB截止频率为截止频率

6、为wc的的M阶（单位脉冲响应阶（单位脉冲响应h(n)长长度度N=M+1）FIR低通（低通（wc为标量）滤波器系数向量为标量）滤波器系数向量hn，默，默认选用认选用哈明窗哈明窗。 滤波器单位脉冲响应滤波器单位脉冲响应h(n)与向量与向量hn的关系为的关系为h(n)=hn(n+1) n=0, 1, 2, , M满足线性相位条件满足线性相位条件: h(n)=h(N1n)。wc为对为对归一化的数字频率，归一化的数字频率，0wc1。当当wc=wcl, wcu时，得到的是时，得到的是带通带通滤波器，其滤波器，其6 dB通通带为带为wclwcu。调用格式及功能：调用格式及功能：第第7 7章章 有限脉冲响应数

7、字滤波器的设计有限脉冲响应数字滤波器的设计 hn=fir1(M, wc, ftype) 可设计高通和带阻FIR滤波器。当ftype=high时，设计高通FIR滤波器；当ftype=stop，且wc=wcl, wcu时，设计带阻FIR滤波器。 应当注意，在设计高通和带阻FIR滤波器时，阶数M只能取偶数（h(n)长度N=M+1为奇数）。不过，当用户将M设置为奇数时，fir1会自动对M加1。hn=fir1(M, wc, window)，可以指定窗函数向量window。如果缺省window参数，则fir1默认为哈明窗。例如:第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 hn

8、=fir1(M, wc, bartlett(M+1)，使用Bartlett窗设计；hn =fir1(M, wc, blackman(M+1)，使用blackman窗设计；hn=fir1(M, wc, ftype, window) 通过选择wc、ftype和 window参数（含义同上），可以设计各种加窗滤波器。(2) fir2 fir2为任意形状幅度特性的窗函数法设计函数，用fir2设计时，可以指定任意形状的Hd(ej)，它实质是一种频率采样法与窗函数法的综合设计函数。主要用于设计幅度特性形状特殊的滤波器（如数字微分器和多带滤波器等）。用help命令查阅其调用格式及调用参数的含义。第第7 7章

9、章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 例7.2.2 的设计程序ep721.m如下: ep721.m: 例7.2.2 用窗函数法设计线性相位高通FIR数字滤波器wp=pi/2; ws=pi/4;Bt=wp-ws; 计算过渡带宽度N0=ceil(6.2*pi/Bt); 根据表7.2.2汉宁窗计算所需h(n)长度N0，ceil(x)取大于等于x的最小整数N=N0+mod(N0+1, 2); 确保h(n)长度N是奇数wc=(wp+ws)/2/pi; 计算理想高通滤波器通带截止频率(关于归一化)hn=fir1(N-1, wc, high, hanning(N); 调用fir1计算

10、高通FIR数字滤波器的h(n)略去绘图部分第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 运行程序得到h(n)的25个值:h(n)= 0.0004 0.0006 0.0028 0.0071 0.0000 0.0185 0.0210 0.0165 0.0624 0.0355 0.10610.2898 0.6249 0.2898 0.1061 0.0355 0.0624 0.0165 0.0210 0.0185 0.0000 0.0071 0.00280.00060.0004高通FIR数字滤波器的h(n)及损耗函数如图7.2.9所示。第第7 7章章 有限脉冲响应数字滤波器

11、的设计有限脉冲响应数字滤波器的设计 例7.2.2 高通FIR数字滤波器的h(n)波形及损耗函数曲线第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 【例例7.2.3】 对模拟信号进行低通滤波处理，要求通带0f 1.5kHz内衰减小于1 dB， 阻带2.5kHzf 上衰减大于40 dB。希望对模拟信号采样后用线性相位FIR数字滤波器实现上述滤波，采样频率Fs=10 kHz。用窗函数法设计满足要求的FIR数字低通滤波器，求出h(n)，并画出损耗函数曲线。为了降低运算量，希望滤波器阶数尽量低。解解 （1） 确定相应的数字滤波器指标: 通带截止频率为3 . 0100001500

12、22sppFf第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 阻带截止频率为 阻带最小衰减为 s =40dB（2） 用窗函数法设计FIR数字低通滤波器，为了降低阶数选择凯塞窗。根据式（7.2.16）计算凯塞窗的控制参数为5 . 010000250022sssFf0.4ss0.5842(21)0.07886(21)3.3953第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 指标要求过渡带宽度Bt=sp=0.2，根据式（7.2.17）计算滤波器阶数为取满足要求的最小整数M=23。所以h(n)长度为N=M+1=24。如果用汉宁窗，h(n)长度为N=

13、40。理想低通滤波器的通带截止频率c=(s+p)/2=0.4，所以得到:式中，w(n)是长度为24（ =3.395）的凯塞窗函数。st840822.28872.2852.2850.2MBdsin0.4 ()1( )( ) ( )( ) , 11.5()2nNh nh n w nw nn第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 实现本例设计的MATLAB程序为ep722.m。ep722.m: 例7.2.3 用凯塞窗函数设计线性相位低通FIR数字滤波器 fp=1500;fs=2500;rs=40; wp=2*pi*fp/Fs;ws=2*pi*fs/Fs; Bt=ws

14、-wp; %计算过渡带宽度 alph=0.5842*(rs-21)0.4+0.07886*(rs-21); %根据(7.2.16)式计算kaiser窗的控制参数 N=ceil(rs-8)/2.285/Bt); %根据(7.2.17)式计算kaiser窗所需阶数N第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 wc=(wp+ws)/2/pi; %计算理想高通滤波器通带截止频率(关于归一化） hn=fir1(N,wc,kaiser(N+1,alph); %调用kaiser计算低通FIRDF的h(n) 以下绘图部分省去运行程序得到h(n)的24个值:h(n)=0.00390

15、.00410.00620.01470.00000.02860.0242 0.03320.07550.00000.1966 0.37240.37240.1966 0.00000.07550.0332 0.0242 0.02860.00000.0147 0.00620.00410.0039第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 例7.2.3 低通FIR数字滤波器的h(n)波形及损耗函数曲线第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 【例例7.2.4】 窗函数法设计一个线性相位FIR带阻滤波器。要求通带下截止频率lp =0.2，阻带下截

16、止频率ls=0.35，阻通带上截止频率us=0.65，通带上截止频率up=0.8， 通带最大衰减 p=1 dB，阻带最小衰减 s=60 dB。解解 本例直接调用fir1函数设计。因为阻带最小衰减 s=60 dB，所以选择布莱克曼窗，再根据过渡带宽度选择滤波器长度N，布莱克曼窗的过渡带宽度Bt=12/N，所以lpls120.350.20.15N 第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 解之得N=80。调用参数lplsusupc,22设计程序为ep723.m，参数计算也由程序完成。ep723.m: 例7.2.4 用窗函数法设计线性相位带阻FIR数字滤波器wlp=0

17、.2*pi;wls=0.35*pi;wus=0.65*pi;wup=0.8*pi; %设计指标参数赋值B=wls-wlp; %过渡带宽度N=ceil(12*pi/B); %计算阶数N,ceil(x)为大于等于x的最小整数第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 wp=(wls+wlp)/2/pi,(wus+wup)/2/pi; %设置理想带通截止频率hn=fir1(N,wp,stop,blackman(N+1); %带阻滤波器要求h(n)长度为奇数，所以取N+1省略绘图部分程序运行结果: N=81由于h(n)数据量太大，因而仅给出h(n)的波形及损耗函数曲线，如

18、图7.2.11所示。第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 例7.2.4 带阻FIR数字滤波器的h(n)波形及损耗函数曲线第第7 7章章 有限脉冲响应数字滤波器的设计有限脉冲响应数字滤波器的设计 )( jdeHnkNNndWkHNnh 101)()(kNeHkHjdd 2 )()( 一、设计思想一、设计思想7.37.3频率采样法设计频率采样法设计FIRFIR滤波器滤波器 10Nnnn)zhzH()(kN2采样 10111NkkNdNZWkHNZzH)()()( jeH第7章 有限脉冲响应数字滤波器的设计 逼逼进进1）FIR滤波器具有线性相位的条件是h(n)为实

19、序列，且满足h(n)=h(Nn1)，其频响函数应满足的条件是：jj ()ddg(e )( )eHH （7.3.5）（7.3.6）21)(N(7.3.7)dgdg( )(2)HHN = 奇数 (7.3.8)dgdg( )(2)HH N = 偶数 二、设计线性相位滤波器时对二、设计线性相位滤波器时对Hd(k)的约束条件的约束条件呈呈偶偶对对称称对对为为奇奇数数：,)( 20 dgHN呈呈偶偶对对称称对对呈呈奇奇对对称称，对对为为偶偶数数：,)( 20 HNdg第7章 有限脉冲响应数字滤波器的设计 在=02区间上N个等间隔的采样频点为将=k代入（7.3.5）（7.3.8）式中，并写成k的函数：kNk

20、2k = 0，1，2，N 1 （7.3.9）)(jgde )()(kkHkH（7.3.10）kNNkNNk1221 )( （7.3.11）)()(kNHkHgg N = 奇数奇数 （7.3.12） )()(kNHkHgg N = 偶数偶数 H Hd d(k(k) )的约的约束条件束条件：第7章 有限脉冲响应数字滤波器的设计 设用理想低通作为设用理想低通作为希望逼近的滤波器希望逼近的滤波器Hd(ej)，截止频率为，截止频率为c，采样点数为采样点数为N，Hg(k)和和(k)用下列公式计算：用下列公式计算： （7.3.13） 12101)(1210)(2101)()(cccgcggNkkNNkkNk

21、kkkHkkkNHkH，ggcccgc()1 0,1, 2,()0 1,2,1()1 1, 2,1() 0,1, 2,1cHkkkHkkkkNkHNkkkNkkkNN 式中式中kc是是通带内最后一个采样点的序号通带内最后一个采样点的序号，kc值取不大于值取不大于cN/(2)的最大整数。另外，的最大整数。另外，对于高通和带阻滤波器，这里对于高通和带阻滤波器，这里N只能取奇数。只能取奇数。 N为偶为偶数时，数时， N为奇为奇数时，数时，第7章 有限脉冲响应数字滤波器的设计 2）FIR滤波器具有线性相位的条件是h(n)为实序列，且满足h(n)=-h(Nn1)，其频响函数应满足的条件是：)()(kNH

22、kHgg N = 奇数奇数 )()(kNHkHgg N = 偶数偶数 H Hd d(k(k) )的约的约束条件束条件：2112212 )()(NkNkNk第7章 有限脉冲响应数字滤波器的设计 3 逼近误差及其改进措施逼近误差及其改进措施如果待逼近的滤波器为Hd(ej)，对应的单位脉冲响应为hd(n)，则由频率域采样定理知道，在频域02范围等间隔采样N点，利用IDFT得到的h(n)应是hd(n)以N为周期的周期延拓的主值区序列，即如果Hd(ej)有间断点，那么相应的单位脉冲响应hd(n)应是无限长的。这样，由于时域混叠及截断，使h(n)与hd(n)有偏差。所以，频域的采样点数N愈大，时域混叠愈小

23、，设计出的滤波器频响特性愈逼近Hd(ej)。d( )()( )Nmh nh nmN Rn第7章 有限脉冲响应数字滤波器的设计 频域采样定理表明，频率域等间隔采样H(k)，经过IDFT得到h(n)，得到H(ej)=FTh(n)的内插表示形式：10j2)()e (NkkNkHH21je)2/sin()2/sin(1)(NNNl 在采样在采样点上点上实际实际频率响应频率响应严格等于理想频率响应严格等于理想频率响应l 采样点采样点之间之间实际频响由实际频响由各采样点加权内插函数延伸而成各采样点加权内插函数延伸而成，因，因而有误差，误差大小和理想频率响应的而有误差，误差大小和理想频率响应的曲线形状有关曲

24、线形状有关改进：在理想频率响应的改进：在理想频率响应的不连续的边缘点不连续的边缘点加上一些过渡的采样加上一些过渡的采样点点第7章 有限脉冲响应数字滤波器的设计 图7.3.1 频域幅度采样序列Hg(k)及其内插波形Hg()Hdg()特性愈平滑的区域，误差愈小；特性曲线间断点处，误差最大。特性愈平滑的区域，误差愈小；特性曲线间断点处，误差最大。第7章 有限脉冲响应数字滤波器的设计 表现形式为间断点变成倾斜下降的表现形式为间断点变成倾斜下降的过渡带曲线过渡带曲线，过，过渡带宽度近似为渡带宽度近似为2/N。通带和阻带内产生。通带和阻带内产生震荡波纹震荡波纹，且，且间断点附近振荡幅度最大，使阻带衰减减小

25、，往往不能间断点附近振荡幅度最大，使阻带衰减减小，往往不能满足技术要求。满足技术要求。 当然，增加当然，增加N可以使过渡带变窄可以使过渡带变窄，但是通带最大衰，但是通带最大衰减和阻带减和阻带最小衰减随最小衰减随N的增大并无明显改善的增大并无明显改善。且。且N太大，太大，会增加滤波器的阶数，即增加了运算量和成本。会增加滤波器的阶数，即增加了运算量和成本。 N=15和N=75两种情况下的幅度内插波形Hg()如图7.3.2所示，图中的空心圆和实心圆点分别表示N=15和N=75时的频域幅度采样。第7章 有限脉冲响应数字滤波器的设计 图7.3.2 N=15和N=75的幅度内插波形Hg()第7章 有限脉冲

26、响应数字滤波器的设计 在窗函数设计法中，通过加大过在窗函数设计法中，通过加大过渡带宽度换取阻带渡带宽度换取阻带衰减的增加衰减的增加。频率采样法同样满足这一规律。频率采样法同样满足这一规律。提高阻带衰减的具体方法：提高阻带衰减的具体方法：在频响在频响间断点附近区间内插间断点附近区间内插一个或几个过渡采样点一个或几个过渡采样点，使不，使不连续点变成缓慢过渡带连续点变成缓慢过渡带，这样，虽然这样，虽然加大了过渡带加大了过渡带，但阻带中相邻内插函数的旁，但阻带中相邻内插函数的旁瓣正负对消，明显瓣正负对消，明显增大了阻带衰减增大了阻带衰减。过渡带采样点的个数与阻带最小衰减 s的关系以及使阻带最小衰减 s

27、最大化的每个过渡带采样值求解都要用优化算法解决。其基本思想是将过渡带采样值设为自由量，用一种优化算法（如线性规划算法）改变它们，最终使阻带最小衰减 s最大。第7章 有限脉冲响应数字滤波器的设计 将过渡带采样点的个数m与滤波器阻带最小衰减 s的经验数据列于表7.3.1中，我们可以根据给定的阻带最小衰减 s选择过渡带采样点的个数m。表表7.3.1 过渡带采样点的个数过渡带采样点的个数m与滤波器阻带最小衰减与滤波器阻带最小衰减 s的经验数据的经验数据 第7章 有限脉冲响应数字滤波器的设计 四四 频率采样法设计步骤频率采样法设计步骤 （1） 根据阻带最小衰减 s选择过渡带采样点的个数m。 （2） 确定

28、过渡带宽度Bt，估算频域采样点数（即滤波器长度）N。 如果增加m个过渡带采样点，则过渡带宽度近似变成(m+1)2/N。当N确定时，m越大，过渡带越宽。如果给定过渡带宽度Bt，则要求(m+1)2/NBt ，滤波器长度N必须满足如下估算公式:（7.3.15） t2(1)NmB第7章 有限脉冲响应数字滤波器的设计 （3） 构造一个希望逼近的频率响应函数: 设计标准型片断常数特性的FIR数字滤波器时，一般构造幅度特性函数Hdg()为相应的理想频响特性，且满足表7.1.1要求的对称性。 （4） 按照（7.3.1）式进行频域采样: jj(1)/2ddg(e)( )eNHH（7.3.16） (1)jjd2(

29、 )(e) ( )e, 0, 1, 2, , 1NkNgkNH kHHkkN（7.3.17） dg2( ), 0, 1, 2, , 1gHkHkkNN第7章 有限脉冲响应数字滤波器的设计 并加入过渡带采样。过渡带采样值可以设置为经验值，或用累试法确定，也可以采用优化算法估算。 （5） 对H(k)进行N点IDFT，得到第一类线性相位FIR数字滤波器的单位脉冲响应：1 01( )IDFT( )( ), 0,1,2,1Nk nNkh nH kH k WnNN（7.3.18） （6） 检验设计结果。如果阻带最小衰减未达到指标要求，则要改变过渡带采样值，直到满足指标要求为止。如果滤波器边界频率未达到指标

30、要求，则要微调Hdg()的边界频率。 上述设计过程中的计算相当繁琐，所以通常借助计算机设计。第7章 有限脉冲响应数字滤波器的设计 【例例7.3.1】 用频率采样法设计第一类线性相位低通FIR数字滤波器，要求通带截止频率p=/3，阻带最小衰减大于40 dB，过渡带宽度Bt/16。解解 查表7.3.1， s=40 dB时，过渡带采样点数m=1。根据m=1和Bt/16估算滤波器长度: N(m+1)2/Bt=64留一点富余量，取N=65。构造Hd(ej)=Hdg()ej(N1)/2为理想低通特性，其幅度响应函数Hdg()如图7.3.3(a)中实线所示。第7章 有限脉冲响应数字滤波器的设计 图7.3.3

31、 一个过渡点的设计结果(T=0.38)第7章 有限脉冲响应数字滤波器的设计 设计由以下程序ep731.m完成：ep732.m: 用频率采样法设计FIR低通滤波器 T=input(T= ) 输入过渡带采样值T Bt=pi/16;wp=pi/3; 过渡带宽度为pi/16，通带截止频率为pi/3 m=1;N=ceil(m+1)*2*pi/Bt)+1; 按式（7.3.15）估算采样点数N Np=fix(wp/(2*pi/N); Np+1为通带0, wp上采样点数 Ns=N2*Np1; Ns为阻带wp, 2*piwp上采样点数 第7章 有限脉冲响应数字滤波器的设计 Hk=ones(1, Np+1), z

32、eros(1, Ns), ones(1, Np);N为奇数，幅度采样向量偶对称A(k)=A(Nk) Hk(Np+2)=T;Ak(NNp)=T; 加一个过渡采样 thetak=pi*(N1)*(0:N1)/N; 相位采样向量(k)=(N1)k/N, 0kN1 Hdk=Hk.*exp(j*thetak); 构造频域采样向量Hd(k) hn=real(ifft(Hdk); h(n)=IDFTH(k) Hw=fft(hn, 1024); 计算频率响应函数:DFTh(n)第7章 有限脉冲响应数字滤波器的设计 wk=2*pi* 0:1023/1024; Hgw=Hw.*exp(j*wk*(N1)/2);

33、计算幅度响应函数Hg() 计算通带最大衰减Rp和阻带最小衰减Rs Rp=max(20*log10(abs(Hgw) hgmin=min(real(Hgw);Rs=20*log10(abs(hgmin) 以下绘图部分略去第7章 有限脉冲响应数字滤波器的设计 运行程序，输入T=0.38，得到设计结果如图7.3.3所示，并输出通带最大衰减 p=0.4767 dB， 阻带最小衰减 s=43.4411 dB。 但是，如果过渡带采样值T=0.5和0.6，则得到阻带最小衰减 s=29.6896 dB和25.0690 dB。 由此可见，当过渡带采样点数给定时，过渡带采当过渡带采样点数给定时，过渡带采样值不同，

34、则逼近误差不同样值不同，则逼近误差不同。所以，对过渡带采样值进行优化设计才是有效的方法。第7章 有限脉冲响应数字滤波器的设计 hn= fir2(M,F,A,window(M+1) 设计一个M阶线性相位FIR数字滤波器，返回长度为N=M+1的单位脉冲响应序列向量hn。window表示窗函数名，缺省时默认选用Hamming窗。可供选择的窗函数有Boxcar、 Bartlett、Hann、Hamm、 Blackman、Kaiser和 Chebwin,当window=boxcar时，fir2就是纯粹的频率采样设计法。 希望逼近的幅度特性由边界频率向量F和相应的幅度向量A确定，plot(F, A)画出的

35、就是希望逼近的幅度特性曲线。MATLAB信号处理工具箱函数fir2是一种频率采样法与窗函数法相结合的FIR数字滤波器设计函数。第7章 有限脉冲响应数字滤波器的设计 图7.3.4 例7.3.1希望逼近的幅度特性（T=0.38）第7章 有限脉冲响应数字滤波器的设计 例例7.3.2 试用频率采样法，设计一个具有线性相位试用频率采样法，设计一个具有线性相位的低通的低通FIR数字数字filter，其，其理想频率特性理想频率特性为：为： 5 . 0 , 05 . 00 , 1)(jdeH5 . 0c已知已知 ，采样点，采样点N=33.H(0)=1,N=33，所以，所以 为第一类线性相位滤波器为第一类线性相

36、位滤波器kNkHH解解：32333211k0 /k)N(kk第7章 有限脉冲响应数字滤波器的设计 截止频率截止频率5 . 0c2582.kkNc 258 , 03225, 80 , 1kkkHk323332k0 eeH)k(H/kjjkk第7章 有限脉冲响应数字滤波器的设计 例例7.3.3 7.3.3 用用频率采样法频率采样法，设计情况，设计情况FIRFIR低通滤波器，采样点数低通滤波器，采样点数N=21N=21，技术指标为：技术指标为：p=0.3, p=5dB s=0.4, s=40dB解：解：根据指标画出理想低通滤波器的幅度特性及采样点根据指标画出理想低通滤波器的幅度特性及采样点根据上图确

37、定频率采样法的幅度约束条件根据上图确定频率采样法的幅度约束条件011111111742019183210HHHHHHHHHN=21时时 频率样本频率样本Hk第7章 有限脉冲响应数字滤波器的设计 相位约束条件相位约束条件)11(Nkk 21400213802136021602140212002019183210 对对H(k)进行进行IDFT得到得到h(n) 102)(1)(NkknNjekHNnh 1021NkknNjjkeeHNk 第7章 有限脉冲响应数字滤波器的设计 200)220(213cos2)220(212cos2)220(21cos21 211)( nnnnnh 计算计算振幅响应振幅

38、响应H()和幅度响应和幅度响应第7章 有限脉冲响应数字滤波器的设计 从幅度响应曲线可见，从幅度响应曲线可见，N=21时，最小的阻带衰时，最小的阻带衰减为减为14dB，没有达到设计要求。，没有达到设计要求。验证指标验证指标 根据设计出的根据设计出的h(n)计算其振幅响应计算其振幅响应H()和幅度响应和幅度响应第7章 有限脉冲响应数字滤波器的设计 N=61时时 频率样本频率样本Hk增大采样点数，能否使阻带衰减增大增大采样点数，能否使阻带衰减增大? ?0150116051100 HHHHHH 613600,613540,61180,61120,6160, 060593210 第7章 有限脉冲响应数字

39、滤波器的设计 N=61时时，最小的阻带衰减为，最小的阻带衰减为18dB，仍不满足指标。，仍不满足指标。结论：结论：频率采样法设计频率采样法设计FIRFIR滤波器滤波器，不能用增大采，不能用增大采样点数的方法提高阻带的衰减。样点数的方法提高阻带的衰减。第7章 有限脉冲响应数字滤波器的设计 在在H Hd d(e(ejj) )的间断点附近区间加的间断点附近区间加入若干个过渡采样点入若干个过渡采样点，能否使能否使H(eH(ejj) )与与 H Hd d(e(ejj) )的误差减小？的误差减小？ 例如：例如：增加增加两两个采样点个采样点T1=0.7和和T2=0.2。 2 . 07 . 001501151

40、104912605290 HHHHHHHHHH 613600,613540,61180,61120,6160, 060593210 第7章 有限脉冲响应数字滤波器的设计 从幅度响应曲线可见，从幅度响应曲线可见，N=61时，最小的阻带衰时，最小的阻带衰减为减为45dB，达到设计要求。，达到设计要求。结论：结论：为改善滤波器的特性，为改善滤波器的特性，可以在频响间断点附近可以在频响间断点附近插入一个或几个过渡采样点插入一个或几个过渡采样点，适当增加采样点数提高，适当增加采样点数提高阻带的衰减。阻带的衰减。N=61N=61，两，两个过渡点个过渡点第7章 有限脉冲响应数字滤波器的设计 窗函数设计法和频

41、率采样法简单方便，易于实现。但它们存在以下缺点缺点: 滤波器边界滤波器边界频率不易精确控制频率不易精确控制。 窗函数设计法窗函数设计法总使通带和阻带波纹幅度相等，总使通带和阻带波纹幅度相等，频率采样法频率采样法只只能依靠优化过渡带采样点的取值控制阻带波纹幅度，所以两种能依靠优化过渡带采样点的取值控制阻带波纹幅度，所以两种方法都不能分别控制通带和阻带波纹幅度。方法都不能分别控制通带和阻带波纹幅度。 所设计的滤波器在阻带边界频率附近的衰减最小，距阻带边所设计的滤波器在阻带边界频率附近的衰减最小，距阻带边界频率越远，衰减越大。所以，如果在阻带边界频率附近的衰界频率越远，衰减越大。所以，如果在阻带边界频率附近的衰减刚好达到设计指标要求，则阻带中其他频段的衰减就有很大减刚好达到设