第7章有限脉冲响应数字滤波器的设计幻灯片资料.pptVIP

第7章 有限脉冲响应数字滤波器的设计 ;7.1 线性相位FIR数字滤波器的条件和特点 ; 式中，Hg(ω)称为幅度特性，θ(ω)称为相位特性。注意，这里Hg(ω)不同于|H(ejω)|, Hg(ω)为ω的实函数，可能取负值，而|H(ejω)|总是正值。H(ejω)线性相位是指θ(ω)是ω的线性函数，即 θ(ω)=-τω, τ为常数 (7.1.3) ? 如果θ(ω)满足下式也称为线性相位， θ(ω)=θ0-τω, θ0是起始相位 (7.1.4) 严格地说，此时θ(ω)不具有线性相位。但以上两种情况都满足群时延是一个常数，即 ; 所以两种情况都称为线性相位。一般称满足(7.1.3)式是第一类线性相位；满足(7.1.4)式为第二类线性相位。 FIR滤波器具有第一类线性相位的条件是：h(n)是实序列且对(N-1)/2偶对称，即 h(n)=h(N-n-1) ? FIR滤波器具有第二类线性相位的条件是：h(n)是实序列且对(N-1)/2奇对称，即 h(n)=-h(N-n-1) 2. 线性相位FIR滤波器幅度特性Hg(ω)的特点;; 3. 线性相位FIR滤波器零点分布特点 第一类和第二类线性相位的系统函数分别满足(7.1.7)式和(7.1.10)式，综合起来用下式表示： ; 4. 线性相位FIR滤波器网络结构 设N为偶数，则有 ;如果N为奇数，则将中间项h［(N-1)/2］单独列出， 从上面的公式可以看出，线性相位FIR滤波器比FIR滤波器的直接型结构节省乘法器近一半。线性相位滤波器的网络结构图如下： ;图7.1.2 第一类线性相位网络结构;图7.1.3 第二类线性相位网络结构;7.2 利用窗函数法设计FIR滤波器 ; 相应的单位取样响应hd(n)为 ; 我们实际实现的滤波器的单位取样响应为h(n),长度为N，其系统函数为H(z)，; 以上就是用窗函数法设计FIR滤波器的思路。另外，我们知道Hd(e jω)是一个以2π为周期的函数，可以展为傅氏级数，即 ;RN(ω)称为矩形窗的幅度函数；将Hd(ejω)写成下式：; 将H(ejω)写成下式： 其卷积结果用图表示比较直观。 ; 图7.2.2 矩形窗对理想低通 幅度特性的影响 ; 通过以上分析可知，对hd(n)加矩形窗处理后，H(ω)和原理想低通Hd(ω)差别有以下两点： (1)在理想特性不连续点ω=ωc附近H(ω)有过渡带。过渡带的宽度，近似等于RN(ω)主瓣宽度，即4π/N。 (2)通带内H(ω)有波动，最大的峰值在ωc-2π/N处。阻带内H(ω)有余振，最大的负峰在ωc+2π/N处。 ; 下面介绍几种??用的窗函数。设 h(n)=hd(n)w(n) ? 式中w(n)表示窗函数。 1. 矩形窗(Rectangle Window) 2. 三角形窗(Bartlett Window) 3. 汉宁(Hanning)窗——升余弦窗 4. 哈明(Hamming)窗——改进的升余弦窗 5. 布莱克曼(Blackman)窗;图7.2.4 常用的窗函数; 图7.2.5 常用窗函数的幅度特性 (a)矩形窗；(b)巴特利特窗(三角形窗)；(c)汉宁窗；(d)哈明窗；(e)布莱克曼窗 ; 图7.2.6 理想低通加窗后的幅度特性(N=51,ωc=0.5π) (a)矩形窗；(b)巴特利特窗(三角形窗)；(c)汉宁窗； (d)哈明窗；(e)布莱克曼窗; 下面介绍用窗函数设计FIR滤波器的步骤。 (1)根据技术要求确定待求滤波器的单位取样响应hd(n)。如果给出待求滤波器的频响为Hd(ejω)，那么单位取样响应用下式求出： ; (2)根据对过渡带及阻带衰减的要求，选择窗函数的形式，并估计窗口长度N。设待求滤波器的过渡带用Δω表示，它近似等于窗函数主瓣宽度。