第七章有限冲激响应滤波器的设计.ppt

第7章有限冲激响应滤波器（有限冲激响应滤波器（FIR）的设计）的设计 71 线性相位线性相位FIR滤波器的特点滤波器的特点 7.2 窗函数设计法窗函数设计法7.3 频率抽样设计法频率抽样设计法7.4 应用应用MATLAB设计设计FIR数字滤波器数字滤波器有限长单位冲激响应数字滤波器的特点有限长单位冲激响应数字滤波器的特点:有限长单位冲激响应（FIR）可以做成具有严格的线性相位，同时又可以具有任意的幅度特性。FIR滤波器的单位抽样响应是有限长的，因而FIR滤波器一定是稳定的。只要经过一定的延时，任何非因果有限长序列都能变成因果的有限长序列，总能用因果系统来实现。FIR滤波器由于单位冲激响应是有限长的，因而可以用快速傅里叶变换（FFT）算法来实现过滤信号，从而可大大提高运算效率。但是，要取得很好的衰减特性，FIR滤波器的阶次比HR滤波器的要高。本章主要讨论线性相位滤波器的设计。本章主要讨论线性相位滤波器的设计。7.1 线性相位线性相位FIR滤波器的特点滤波器的特点7.1.1 线性相位条件线性相位条件如果一个线性移不变系统的频率响应有如下形式：（7.1）则其具有线性相位。这里 是一个实数。因而，线性相位系统有一个恒定的群延时 （7.2）()()()|()|jjjjH eHeH ee 在实际应用中，有两类准确的线性相位，分别要求满足 （7.3）（7.4）FIR滤波器具有式（7.3）的线性相位的充分必要条件是：单位抽样响应 关于群延时 偶对称，即满足（7.5）（7.6）()()()h n12N()(1)01h nh NnnN 满足式（7.5）和式（7.6）的偶对称条件的FIR滤波器分别称为I型线性相位滤波器和型线性相位滤波器。FIR滤波器具有式（7.4）的线性相位的充分必要条件是：单位抽样响应 关于群延时 奇对称，即满足（7.7）（7.8）（7.9）()h n12N2()(1)01h nh NnnN 把满足式（7.7）、（7.8）和式（7.9）的奇对称条件的FIR滤波器分别称为型线性相位滤波器和型线性相位滤波器。1I型线性相位滤波器型线性相位滤波器 7.1.2 线性相位滤波器频率响应的特点线性相位滤波器频率响应的特点由于偶对称性，一个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,22型线性相位滤波器型线性相位滤波器 一个型线性相位滤波器，由于N是偶数，所以，的对称中心在半整数点 。其频率响应可以表示为：（7.13）其中()h n12N/2(1)/20()1()cos()2Njj NnH eeb kk()2()1,2,.,22NNb khkk幅度函数为 （7.14）相位函数为 （7.15）()(1)2N/20()1()cos()2NnHb kk型线性相位滤波器的幅度函数和相位函数的特点：幅度函数和相位函数的特点：幅度函数的特点：幅度函数的特点：（1）当 时，=0，也就是说 在 处必然有一个零点；（2）对 呈奇对称，对 呈偶对称。相位函数的特点：同相位函数的特点：同I I型线性相位滤波器。型线性相位滤波器。()H()Hz1z ()H0,23型线性相位滤波器型线性相位滤波器 由于型线性相位滤波器关于 奇对称，且 为整数，所以，其频率响应可以表示为 （7.16）其中 12N(1)/2(1)/21()()sin()Njj NnH ejec kk11()2()1,2,.,22NNc khkk幅度函数为 （7.17）相位函数为 （7.18）(1)/21()()sin()NnHc kk()(1)22N 型线性相位滤波器的幅度函数和相位函数的特点幅度函数和相位函数的特点：幅度函数的特点：幅度函数的特点：（1）当 时，=0，也就是说 在 处都为零点；（2）对 均呈奇对称。相位函数的特点：相位函数的特点：既是准确的线性相位，又包括 的相移，所以又称 移相器，或称正交变换网络。()H()Hz()H0,21z 0,2/2904型线性相位滤波器型线性相位滤波器 型线性相位滤波器关于 奇对称，且N为偶数，所以为非整数。其频率响应可以表示为 （7.19）其中 12N/2(1)/21()1()sin()2Njj NnH ejed kk()2()1,2,.,22NNd khkk幅度函数为 （7.20）相位函数为 （7.21）()(1)22N/21()1()sin()2NnHd kk型线性相位滤波器的幅度函数和相位函数的特点幅度函数和相位函数的特点：幅度函数的特点幅度函数的特点：（1）在 处必为零，也就是说 在 处为零点；（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一个线性相位滤波器零点的约束条件一个线性相位滤波器零点的约束条件 线性相位滤波器的级联结构实现 5N 122N在此情况下，111112212211()(1)(1)(1)(1)1 12(cos)12(cos)iiiijjjjiiiiiiiiiiiH zz rez rezezerrrzr zrzzr（1）零点 既不在实轴上，也不在单位圆上，零点是两组互为倒数的共轭对，其基本因子为 izz（7.22）,1,0ijiiiizrer在此情况下，3N 112N（2）零点 在单位圆上，但不在实轴上，即 ，此时零点的共轭值就是他的倒数，其基本因子为 izz（7.23）1112()(1)(1)12(cos)iijjiiH zz ez ezz 1,0,iiir在此情况下，3N 112N（3）零点 在实轴上，但不在单位圆上，即 ，此时零点是实数，他没有复共轭部分，只有倒数，倒数也在实轴上，其基本因子为 izz（7.24）1,0 iir或 111211()(1)(1)1()iiiiiH zrzzrzzrr 式中“”号相当于 ，零点在负实轴上，“”相当于零点在正实轴上。i0i（4）零点 既在单位圆上，但在实轴上，即 ，此时零点只有两种情况，即 ，这时零点既是自己的复共轭，又是倒数，其基本因子为 izz（7.25）1,0iir 或1()1iH zz 式中“”号相当于 ，零点在负实轴上，“”相当于零点在正实轴上。i0i在此情况下，即有半个抽样的延时。2N 11/22N线性相位滤波器只能由以上这几种因子的组合而成。11zz 或7.2 窗函数设计法窗函数设计法 7.2.1 设计方法设计方法 给出所要求的理想低滤波器频率响应()jdeH设计一个FIR滤波器频率响应 10()()Njj nneeHh n逼近()jdeH窗函数设计法设计是在时域进行（7.26）()jdeH由 的傅里叶反变换导出 sin()11()()22ccjj nj ncddeeenh nHddn()jdeH()dh n()h n由于 是矩形频率特性，故 一定是无限长的序列，而且是非因果的，而要设计的是FIR滤波器，必然 是有限长的。所以要用有限长的 来逼近无限长的 ，最有效的方法是截断 ，即用一个有限长度的窗函数序列 来截取 ，并将截短后的 移位，得()h n()dh n()dh n()n()dh n()dh n11()()()22dNNh nnh n（7.27）窗函数序列的形状及长度的选择很关键。取 ，即取矩形窗 例例7.1 设计一低通滤波器，所希望的频率响应截止频率 在 之间为1，在 之间为0，分别取N=11,21,41，观察其频谱响应的特点。()jdeH00.250.25解：解：(1)/2 00.25()0 0.25j NjdeeH1 01()0 nNn其它 由式（7.27）1sin0.25()12()()12()2dNnNh nh 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取的过小时，通频带过窄，且阻带内波纹较大，过渡带较宽，当N增大时，与 的近似程度越来越好。但当N增大时，通带内出现了波纹，而且随着N的继续增大，这些波纹并不消失，只是最大的尖峰处越来越接近于间断点，这种现象称作吉布斯现象。()jeH()jdeH吉布斯现象的产生是由于对 突然截短的结果。()dh n为了减少吉布斯现象，应选取旁瓣较小的窗函数。1矩形窗矩形窗 7.2.2 各种窗函数各种窗函数 窗函数为()()NnRn（7.28）幅度函数为（7.29）sin()2()|()|sin()2jRRNWW e主瓣宽度 ，过渡带宽 。4/2 2/NN0.9 2/N2汉宁（汉宁（Hanning）窗（又称升余弦窗）窗（又称升余弦窗）窗函数为 2()0.50.5cos()()1NnnRnN（7.30）幅度函数为（7.31）22()0.5()0.25()()11RRRnnWWWWNN主瓣宽度 ，过渡带宽 。4 2/8/NN3.1 2/N3海明（海明（Hamming）窗（又称改进的升余弦窗）窗（又称改进的升余弦窗）窗函数为（7.32）幅度函数为（7.33）主瓣宽度 ，过渡带宽 。4.4.凯泽（凯泽（Kaiser）窗）窗 窗函数为（7.34）2()0.540.46cos()()1NnnRnN22()0.54()0.23()()11RRRnnWWWWNN4 2/8/NN3.3 2/N20021(1)1(),01()nINnnNI其中 为第一类变形零阶贝塞尔函数，是一个可自由选择的参数，改变 值就可对主瓣宽度与旁瓣衰减进行选择，一般选择 。过渡带宽 。0I495 2/N 窗函数窗谱性能指标加窗后滤波器性能指标旁瓣峰值（dB）主瓣宽度()过渡带宽()阻带最小衰减（dB）矩形窗汉宁窗海明窗凯泽窗-13-31-41-57244 0.93.13.35-21-44-53-80最小阻带衰减只由窗形决定，不受N的影响，而过渡带宽则随N的增加而减小。表表7.1 几种窗函数的基本参数比较几种窗函数的基本参数比较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带阻数字滤波器的设计带阻数字滤波器的设计 令（7.39）则 求得 （7.40）(1)/2|,|()0 j NjlhdeeH其它 111()()()222111()222hllhNNNjnjnjndeeeh nddd111sinsinsin222()1()2lhdNNNnnnh nNn比较式（7.36）、（7.38）、（7.40）可知，一个高通滤