第四有限长单位脉冲响应fir滤波器设计方法.pptxVIP

;;缺点：（1）因为无极点，要获得好的过渡带特性，需以较 高的阶数为代价； （2）无法利用模拟滤波器的设计结果，一般无解 析设计公式，要借助计算机辅助设计程序完成。;§4.1 线性相位FIR数字滤波器的特性 ;FIR滤波器的DTFT为;将上式两边交叉相乘，再将等式右边各项移到左边，应用三角函数的恒等关系 ;另外一种情况是，除了上述的线性相位外，还有一附加的相位，即 ; 偶对称;;;第11页/共82页;令 ， 则;;;3. h(n)奇对称，N为奇数，h(n)=-h(N-1-n);;;4.h(n)奇对称，N为偶数 ;;四种线性相位FIR滤波器;四种线性相位FIR DF特性，参考 P91 表 第一种情况 ，偶、奇，四种滤波器都可设计。;例1 N=5, h (0) = h (1) = h (3) = h (4) = -1/2, h (2) = 2，求幅度函数H (ω)。 解 Ｎ为奇数并且h(n)满足偶对称关系 a (0) = h (2) = 2 a (1) = 2 h (3) = -1 a (2) = 2 h (4) = -1 H (ω) = 2 - cosω- cos2ω = 2- (cosω+cos2ω) ;小结：;;;第26页/共82页;;§4.2 窗口设计法（时域）; 但一般来说，理想频响 是分段恒定，在边界频率处有突变点，所以，这样得到的理想单位脉冲响应hd(n)往往都是无限长序列，而且是非因果的。但FIR的h(n)是有限长的，问题是怎样用一个有限长的序列去近似无限长的hd(n)。最简单的办法是直接截取一段 hd(n) 代替 h(n) 。这种截取可以形象地想象为h(n)是通过一个“窗口”所看到的一段hd(n)，因此 ，h(n)也可表达为h(n)和一个“窗函数”的乘积，即 h(n)=w(n) hd(n) 在这里窗口函数就可以是矩形脉冲函数RN（n），当然以后我们还可看到，为了改善设计滤波器的特性，窗函数还可以有其它的形式，相当于在矩形窗内对hd(n)作一定的加权处理。;窗函数法设计流程：;一.矩形窗口法;;;;矩形窗函数及其幅度函数（见P94图）;;; 矩形窗的卷积过程（P95的图来说明）;4个特殊频率点看卷积结果：;（2）ω=ωc时，一半重叠， H(ωc)=0.5 H(0); （3） ω=ωc –2π/N时，第一旁瓣（负数）在通带外，出现正肩峰； （ 4） ω=ωc +2π/N 时，第一旁瓣（负数）在通带内，出现负肩峰。 ;窗口函数对理想特性的影响： ①改变了理想频响的边沿特性，形成过渡带，宽为 ， 等于WR(ω)的主瓣宽度。（决定于窗长） ②过渡带两旁产生肩峰和余振（带内、带外起伏），取决于 WR(ω)的旁瓣，旁瓣多，余振多；旁瓣相对值大，肩峰强 ，与 N无关。（决定于窗口形状） ③N增加,过渡带宽减小,肩峰值不变。 因主瓣附近 其中x=Nω/2,所以N的改变不能改变主瓣与旁瓣的比例关系，只能改变WR（ω）的绝对值大小和起伏的密度，当N增加时，幅值变大，频率轴变密，而最大肩峰永远为8.95%，这种现象称为吉布斯（Gibbs）效应。 ;; 改变窗函数的形状，可改善滤波器的特性，窗函数有 许多种，但要满足以下两点要求：;;三部分矩形窗频谱相加，使旁瓣互相抵消，能量集中在主瓣，旁瓣大大减小，主瓣宽度增加1倍，为 。 ?;第46页/共82页;;第48页/共82页;窗口函数的频谱 N=51，A=20lg|W(ω)/W(0)|;第50页/共82页;;例题：用矩形窗设计一个线性相位低通滤波器 ;解： （1） 因此 ;（2）因为 偶对称，所以若N为奇数，则属于第一种线性相位滤波器，若N为偶数，则属于第二种线性相位滤波器。 （3）若为汉宁窗，则 ;§4.3 频率采样法;;二.设计方法;三、?约束条件;;因此，Hk 也必须满足奇对称性： 相位关系同上， 其它两种线性相位FIR数字滤波器的设计，同样也要满足幅度与相位的约束条件。?;四、逼近误差;;;式中