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

第7章 有限脉冲响应数字滤波器的设计 7.1 线性相位FIR数字滤波器的条件和特点 7.2 利用窗函数法设计FIR滤波器 7.3 利用频率采样法设计FIR滤波器 7.5 IIR和FIR数字滤波器的比较 7.1 线性相位FIR数字滤波器的条件和特点 特点： 1. 具有(N-1)个零点和一个(N-1)重极点z=0； 2. 收敛域包括单位圆，系统总是稳定的； 3. 易做到严格的相位特性。 1. 线性相位 FIR 数字滤波器 对于长度为 N 的h(n)，传输函数为 群时延： (2) 证明第二类线性相位条件： 2. 线性相位FIR滤波器零点分布特点 第一类和第二类线性相位的系统函数综合起来表示为： 3. 线性相位FIR滤波器网络结构 7.2 利用窗函数法设计FIR滤波器 滤波器传输函数 单位脉冲响应 怎样得到 h(n) 的频谱？ 时域相乘，频域卷积 H(ω)和理想低通Hd(ω)差别： (1)在理想特性不连续点ω=ωc附近形成过渡带。过渡带的宽度，近似等于RN(ω)主瓣宽度，即4π/N。 (2)通带内增加了波动，最大的峰值在ωc-2π/N处。阻带内产生了余振，最大的负峰在ωc+2π/N处。通带与阻带中的波纹与窗函数的幅度谱有关。 1. 矩形窗(Rectangle Window) wR(n)=RN(n) 3. 汉宁(Hanning)窗——升余弦窗 4. 哈明(Hamming)窗——改进的升余弦窗 5. 布莱克曼(Blackman)窗 6. 凯塞—贝塞尔窗(Kaiser-Basel Window) 7.2.3 窗函数设计FIR滤波器的步骤 (1)根据对阻带衰减及过渡带的指标，选择窗函数类型、并估计窗口长度 N ； (2)构造希望逼近的频率响应函数Hd(ejω) 。理想截止频率一般取ωc=(ωp+ωs )/2 ； (3)计算理想滤波器的单位取样响应hd(n)； (4) 计算 FIR 滤波器的单位脉冲响应： h(n)=hd(n)w(n) (5)求H (ejω)=FT[h(n)]检验是否满足设计要求。 例：用窗函数法设计线性相位FIR高通滤波器，要求通带截止频率ωp=0.5πrad ，阻带截止频率ωs=0.25πrad ，通带最大衰减αp=1dB，阻带最小衰减αs=40dB。 解： (1)根据阻带衰减αs及选择汉宁窗； 过渡带B =ωp-ωs=0.25πrad，估计窗口长度 N≥24.8 ； 根据P198表7.1.1选择情况I，确定N=25。 7.3 利用频率采样法设计FIR滤波器 在ω=0到2π之间，对待设计的滤波器的传输函数Hd(ejω)等间隔采样N点，得到Hd(k)， 1. 用频率采样法设计线性相位滤波器的条件 FIR滤波器具有线性相位的条件： h(n)是实序列，且满足h(n)=h(N-n-1) 传输函数应满足的条件是： 希望设计理想低通滤波器，截止频率为ωc，采样点数N，Hg(k)和θ(k)的计算公式： (1) N=奇数； (2) N=偶数。 (2) N=偶数时， 例：利用频率采样法设计线性相位低通滤波器，要求截止频率ωc=π/2rad，采样点数N=33，选用h(n)=h(N-1-n)情况。 解：用理想低通作为逼近滤波器， 7.5 IIR和FIR数字滤波器的比较 在相同的技术指标下，IIR滤波器存在着输出对输入的反馈，可以用比FIR滤波器较少的阶数来满足指标，所以存储单元少，运算次数少，较为经济。 FIR滤波器可以得到严格的线性相位；而IIR滤波器的频率选择性越好，其相位非线性越严重，需要加全通网络进行相位校正。 (3) FIR采用非递归结构，肯定是稳定的系统；IIR滤波器采用递归结构，极点必须在单位圆内，运算过程中对序列的舍入处理会引起寄生振荡。 (4) FIR滤波器冲击响应有限长，可以使用快速傅里叶变换算法，运算速度加快，但是对计算工具要求较高。 (5) IIR滤波器可以利用模拟滤波器设计的现成公式、数据和表格，因而计算量较小，对计算工具要求不高。 (6) FIR滤波器设计较灵活，可适应各种幅度特性和相位特性的要求，适应性较广。 Homework 1、3、6、8