第七章_有限脉冲响应数字滤波器的设计.pptVIP

第七章FIR数字滤波器的设计 学习目标 掌握线性相位FIR数字滤波器的特点 掌握窗函数设计法 了解频率采样设计法 掌握IIR与FIR数字滤波器的相应特点 IIR数字滤波器： 优良的幅度特性 可以利用模拟滤波器设计 但相位非线性 FIR数字滤波器 可以严格线性相位，又可任意幅度特性 因果稳定系统 可用FFT计算 但阶次比IIR滤波器要高很多 7.1 线性相位FIR滤波器的特点 FIR滤波器的单位冲击响应： 系统函数 在z平面上有N-1个零点 在z=0处有N-1阶极点 1、线性相位FIR数字滤波器 h(n)为实序列时，其频率响应： Hg(ω)：幅度特性 θ(ω)：相位特性。 线性相位是指：θ(ω)是ω的线性函数，即群时延是一个常数 第一类线性相位： 第二类线性相位： ，θ0是起始相位 2、线性相位FIR滤波器的时域约束条件 （1）第一类线性相位对h(n)的约束条件： 满足第一类线性相位的条件 （2）第二类线性相位对h(n)的约束条件： 满足第二类线性相位的条件 2、线性相位FIR滤波器幅度特性Hg(ω)的特点 情况1： h(n)=h(N－n－1), N为奇数。 情况2： h(n)=h(N－n－1), N为偶数。 情况3： h(n)=-h(N－n－1), N为奇数。 情况4： h(n)=-h(N－n－1), N为偶数。 cos［ω(n-τ)］关于ω=0, π, 2π三点偶对称 Hg(ω)关于ω=0, π, 2π三点偶对称。因此情况1可以实现各种（低通、高通、带通、带阻）滤波器。 cos［ω(n－τ)］关于过零点奇对称，关于ω=0和2π偶对称。 Hg(π)=0，Hg(ω)关于ω=π奇对称，关于ω=0和2π偶对称。因此，情况2不能实现高通和带阻滤波器。 情况3： h(n)=－h(N－n－1)，N为奇数。 当ω=0，π, 2π时，sin[ω(n－τ)]=0, 且sin[ω(n－τ)]关于过零点奇对称。因此Hg(ω)关于ω=0, π, 2π三点奇对称。由此可见，情况3只能实现带通滤波器。 情况4： h(n)=－h(N－n－1), N为偶数。 用情况3的推导过程可以得到:  　　　　　　　　　　　　　　　　　　　 N是偶数，τ=(N－1)/2=N/2－1/2。 当ω=0, 2π时，sin［ω(n－τ)］=0；当ω=π时，sin［ω(n－τ)］=(－1)n－N/2, 为峰值点。而且sin［ω(n－τ)］关于过零点ω=0和2π两点奇对称，关于峰值点ω=π偶对称。 Hg(ω)关于ω=0和2π两点奇对称，关于ω=π偶对称。由此可见，情况4不能实现低通和带阻滤波器。 h(n)为偶对称时，适用于选频滤波器 h(n)为奇对称时，有90度相移，适用于微分器和90度移相器 3、线性相位FIR滤波器的零点分布特点 特殊情况： 零点是实数：两个零点 零点是纯虚数且在单位圆上：两个零点 零点是实数且在单位圆上：只有一个零点 7.2 利用窗函数设计FIR滤波器 1、设计方法 以低通滤波器为例讨论： 线性相位理想低通滤波器的频率响应： 加窗处理后对频率响应的影响： 加窗函数的影响： 在理想特性不连续点ω=ωc附近形成过渡带。过渡带的宽度近似等于WRg(ω)主瓣宽度4π/N。 在 处出现肩峰值，通带与阻带中波纹的情况与窗函数的幅度谱有关， WRg(ω)旁瓣幅度的大小直接影响Hg(ω)波纹幅度的大小。 吉布斯效应 吉布斯效应由hd[n] 直接截断引起，也称为截断效应。 该效应引起通带和阻带内的波动性，尤其使阻带的衰减变小。影响滤波器的性能。 如何减少Gibbs效应？ 增大N？ 主瓣幅度和旁瓣幅度相对值不变 主瓣和旁瓣宽度变窄，波动频率加快 改变窗函数的形状 主瓣包含更多能量，旁瓣幅度更小 7.2.2 典型窗函数介绍 窗函数的要求： 窗谱主瓣尽可能窄以获得较陡的过渡带 尽量减少窗谱最大旁瓣的相对幅度以减小肩峰和波纹 两者之间取折中 矩形窗 三角形窗 汉宁窗——升余弦窗 哈明窗——改进的升余弦窗 布莱克曼窗——二阶升余弦窗 常用窗函数 凯塞——贝塞尔窗 7.2.3窗函数法设计FIR滤波器的步骤 4. 加窗得到设计结果：h(n)=hd(n)w(n) FIR滤波器设计实例 1、设计一个线性相位FIR低通滤波器， 7.5 IIR和FIR数字滤波器的比较 IIR