有限脉冲响应数字滤波器设计原理.pptVIP

第6章 有限脉冲响应数字滤波器的设计 6.1 线性相位FIR数字滤波器的条件和特点 6.2 利用窗函数法设计FIR滤波器 6.3 利用频率采样法设计FIR滤波器 6.4 利用切比雪夫逼近法设计FIR滤波器 6.5 IIR和FIR数字滤波器的比较 设FIR滤波器单位脉冲响应h(n)长度为N，其系统函数H(z)为： H(z)是 的N-1次多项式，它在z平面上有N-1个零点，原点z=0是N-1阶重极点。因此H(z)永远稳定。 稳定和线性相位是FIR滤波器突出的优点。 6.1 线性相位FIR数字滤波器的条件和特点 本节主要介绍FIR滤波器具有线性相位的条件及幅度特性。 1. 线性相位条件 对于长度为N的h(n)，传输函数定义为 式中，Hg(ω)称为幅度特性，θ(ω)称为相位特性。注意，这里Hg(ω)不同于|H(ejω)|,Hg(ω)为ω的实函数，可能取负值，而|H(ejω)|总是正值。H(ejω)线性相位是指θ(ω)是ω的线性函数，即 θ(ω)=-τω, τ为常数 (6.1.3) ? 如果θ(ω)满足下式： θ(ω)=θ0-τω, θ0是起始相位 (6.1.4) 严格地说，此时θ(ω)不具有线性相位，但以上两种情况都满足群时延是一个常数，即 一般称满足(6.1.3)式是第一类线性相位；满足(6.1.4)式为第二类线性相位。 下面推导与证明满足第一类线性相位的条件是： h(n)是实序列且对(N-1)/2偶对称，即 h(n)=h(N-n-1) (6.1.5) 满足第二类线性相位的条件是： h(n)是实序列且对(N-1)/2奇对称，即 h(n)=-h(N-n-1) (6.1.6) (1) 第一类线性相位条件证明： 按照上式可以将H(z)表示为 (2) 第二类线性相位条件证明： 6.2 利用窗函数法设计FIR滤波器 设希望设计的滤波器传输函数为Hd(ejω),hd(n)是与其对应的单位脉冲响应，因此 如果已知 则可求出hd(n)，经过z变换可得到滤波器的系统函数。 但是一般 逐段恒定，在边界频率处有不连续点，故致使hd(n)是无限时宽的，且是非因果序列。 例如：理想的低通滤波器的传输函数 我们知道Hd(ejω)是一个以2π为周期的函数，可以展开为傅氏级数，即 将H(ejω)写成下式： 通过以上分析可知，对hd(n)加矩形窗处理后，H(ω)和原理想低通Hd(ω)差别有以下两点： (1)在理想特性不连续点ω=ωc附近形成过渡带。过渡带的宽度，近似等于RN(ω)主瓣宽度，即4π/N。 (2)通带内增加了波动，最大的峰值在ωc-2π/N处。阻带内产生了余振，最大的负峰在ωc+2π/N处。通带与阻带中的波动与窗函数的幅度谱有关系。RN(ω)波动越快(N加大时)，通带、阻带波动越快，RN(ω)旁瓣的大小直接影响H(ω)波动的大小。 通过实验分析：调整窗口的长度N可以有效的控制过渡带的宽度。减少带内波动以及加大阻带的衰减只能从窗函数的形状上找解决办法。 思路：能否找到窗函数的形状，使其谱函数的主瓣包含更多的能量，相应旁瓣幅度就减小了，旁瓣的减小可使通带、阻带波动减小，从而加大阻带衰减。 下面介绍几种常用的窗函数。设 h(n)=hd(n)w(n) ? 式中w(n)表示窗函数。 1. 矩形窗(Rectangle Window) wR(n)=RN(n) 前面已分析过，按照(6.2.5)式，其频率响应为 2. 三角形窗(Bartlett Window) 3. 汉宁(Hanning)窗——升余弦窗 4. 哈明(Hamming)窗——改进的升余弦窗 5. 布莱克曼(Blackman)窗 6. 凯塞－贝塞尔窗(Kaiser-Basel Window) MATLAB提供了窗函数的子程序： 1)w=boxcar(M)在数组w中产生M点的矩形窗函数； 2)w=tr