数字信号处理第七章 有限长数字滤波器的设计.ppt

凯塞窗的幅度谱函数为： 凯塞窗参数对滤波器的性能影响 六种窗函数的基本参数 三、窗函数法的设计步骤 给定理想的频率响应函数 及技术指标 求出理想的单位抽样响应 根据阻带衰减选择窗函数 计算频率响应 ，验算指标是否满足要求 根据过渡带宽度确定N值 求所设计的FIR滤波器的单位抽样响应 奇对称单位冲击响应 h(n)=-h(N-1-n) 三、线性相位FIR滤波器零点分布特点 1）若 z = zi 是H(z)的零点，则 z = zi-1 也是零点 2）h(n)为实数，则零点共轭成对 线性相位滤波器的零点是互为倒数的共轭对 这样H(z)的零点有4种情况 (1) 既不在实轴上，也不在单位圆上,则零点是互为倒数的两组共轭对,如图z1 (2) 不在实轴上,但在单位圆上,共轭对的倒数就是它们本身,如图z3 (3) 在实轴上,不在单位圆上,只有倒数。如图z2 (4) 既在实轴上也在单位圆上。此时，只有一个零点，且有两种可能，或位于z=1，或位于z=-1。如图z4 四、线性相位FIR滤波器网络结构 1）N为偶数 2)N为奇数，则将中间项单独列出 先给定理想filter的频响 ,设计一个FIR filter的频响为 ,使 逼近 设计是在时域进行的,先用傅氏反变换求出理 想filter的单位抽样响应 ,然后加时间窗 对 截断,以求得FIR filter的单位抽样响应h(n)。 §7.2 窗函数设计法 2、设计过程 1、设计思想 一、设计方法 因而窗的形状和长度的选择很关键。 1 0 0 以低通滤波器为例讨论： 这是一个以为 中心的偶对称的无限长非因果序列， 如果截取一段n=0～N-1的hd(n)作为h(n)，则为 保证所得到的是线性相位FIR滤波器，必须使  ，即使关于 对称，故必须使 即： 其中 现在考虑，对理想的 用矩形窗简单地截断后，对其频响的影响，即 去逼近 的质量如何？ 矩形窗 其频谱 具有线性相位 其中： 其在 内有一主瓣 理想滤波器频率响应： 其中： 实际FIR滤波器的幅度函数H(ω)为 (1)ω=0时, H(0)等于 在[-ωc, ωc]内的积分面积,因一般 ,故H(0)近似为 在[-π, π]内的积分面积 (2)ω=ωc时，一半重叠，H(ωc)=0.5 H(0); (3) ω=ωc –2π/N时，第一旁瓣(负数)在通带外，出现正肩峰； (4) ω=ωc +2π/N 时,第一旁瓣(负数)在通带内,出现负肩峰。 加矩形窗处理后，对理想频率响应产生了两点影响： 1）使理想频率特性不连续点ω=ωc 处，形成一个过渡带, 过渡带宽等于矩形窗的WR(ω)的主瓣宽度△ω=4π/N ； 2) 在截止频率ωc的两边ω=ωc ±2π/N 处，H(ω)出现最大的肩峰值，肩峰的两侧形成起伏振荡，其振荡幅度取决于旁瓣的相对幅度，而振荡快慢则取决于WR(ω)波动的快慢。 在主瓣附近 该函数的性质：随着x加大（即N加大时），函数曲线波动的频率加快，主瓣幅度加高，旁瓣幅度也同样加高，主瓣与旁瓣的相对比例保持不变。 其中： 在矩形窗情况下,最大相对肩峰值为8.95%，N增加时,过渡带宽4π/N减小,起伏振荡变密,但最大肩峰则总是8.95%,这就是吉布斯(Gibbs)效应。 由于肩峰的存在，影响到H(ω)通带的平坦和阻带的衰减，使阻带最小衰减只有21dB左右，因此在实际中，矩形窗很少采用。 怎样才能改善阻带的衰减特性呢？只能从改善窗函数的形状上找出路。要求： ①主瓣宽度要窄，以获得较陡的过渡带； ②相对于主瓣幅度，旁瓣要尽可能小，使能量尽量集中在主瓣中，这样就 可以减小肩峰和余振，以提高阻带衰减和通带平稳性。 但实际上这两点不能兼得，一般总是通过增加主瓣宽度来换取对旁瓣的抑制。 二、几种常用的窗函数： 1. 矩形窗 时域表达式： 频域表达式（频谱）： 幅度函数： 主瓣宽度最窄： 旁瓣幅度大 2.三角窗（Bartlett）窗 时域表达式： 频谱： 幅度函数： 主瓣宽度 ,比矩形宽一倍。旁瓣幅度较小 3、汉宁(Hanning)窗(又称升余弦窗) 频谱： 时域： 幅度函数： 由于汉宁窗的