数字信号处理曹成茂有限长单位脉冲响应fir滤波器的设计方法.pptxVIP

数字信号处理曹成茂有限长单位脉冲响应fir滤波器的设计方法;第五章 有限长单位脉冲响应（FIR）滤波器的设计方法;;;缺点：（1）因为无极点，要获得好的过渡带特性，需以较 高的阶数为代价； （2）无法利用模拟滤波器的设计结果，一般无解 析设计 公式，要借助计算机辅助设计程序完成。;第6页/共159页;第7页/共159页;第8页/共159页;第9页/共159页;第10页/共159页;第11页/共159页;第12页/共159页;第13页/共159页;第14页/共159页;第15页/共159页;第16页/共159页;第17页/共159页;第18页/共159页;第19页/共159页;第20页/共159页;第21页/共159页;第22页/共159页;第23页/共159页;第24页/共159页;第25页/共159页;第26页/共159页;第27页/共159页;第28页/共159页;第29页/共159页;第30页/共159页;第31页/共159页;第32页/共159页;第33页/共159页;第34页/共159页;第35页/共159页;§5.1 线性相位FIR数字滤波器的特性 ;线性相位FIR滤波器的DTFT为;将上式两边交叉相乘，再将等式右边各项移到左边，应用三角函数的恒等关系 ;另外一种情况是，除了上述的线性相位外，还有一附加的相位，即 ;;;;第43页/共159页;令 ，则;第45页/共159页;;;第48页/共159页;3. h(n)奇对称，N为奇数，h(n)=-h(N-1-n);;;第52页/共159页;4.h(n)奇对称，N为偶数 ;;四种线性相位FIR滤波器;四种线性相位FIR滤波器;四种线性相位FIR DF特性，参考表5.1 第一种情况 ，偶、奇，四种滤波器都可设计。;例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ω) ;小结：;;;第62页/共159页;;§5.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)作一定的加权处理。;设计步骤：;一.矩形窗口法;第68页/共159页;;;对频响起作用的是它的幅度函数 ; 理想频响也可以写成幅度函数和相位函数的表示形式 Hd(ejω)=Hd(ω)e-jωα 其中幅度函数为 两个信号时域的乘积对应于频域卷积，所以有;;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的改变不能改变