数字信号处理第七章-2(精品·公开课件).ppt

数字信号处理Digital Signal Processing    窗函数法的主要工作是计算hd(n)和w(n)的选取。 　用窗函数进行FIR设计的步骤： (1) 确定待求滤波器的理想单位脉冲取样响应hd(n)； (2) 根据过渡带及阻带衰减的要求，选择窗函数的形状； (3) 估计窗长N (4)确定 (5) 计算 ，验证是否满足指标。 过度带宽 A取决于窗的形状，如矩形窗 ，哈明窗 等。 主瓣带宽 * * §7.2 利用窗函数法设计FIR滤波器（时域方法） 如果希望得到的滤波器的理想频率响应为 ，那么 FIR滤波器的设计就在于寻找一个频率响应 去逼近 ，逼近方法有三种： 窗口设计法（时域逼近） 频率采样法（频域逼近） 最优化设计（等波纹逼近） 时间窗函数设计法是从单位脉冲响应序列着手，使h(n)逼近理想的单位脉冲响应序列hd(n)。我们知道hd(n)可以从理想频响通过傅氏反变换获得： 但理想频响 是分段恒定，在边界频率处有突变点，这样得到的理想单位脉冲响应hd(n)是无限长序列，而且是非因果的。 FIR的h(n)是有限长的，问题是怎样用一个有限长的序列去近似无限长的hd(n)。最简单的办法是直接截取一段hd(n) 代替 h(n) 。这种截取可以形象地想象为h(n)是通过一个“窗口”所看到的一段hd(n)，因此，h(n)也可表达为hd(n)和一个“窗函数”的乘积，即 h(n) = w(n) hd(n) 在这里窗口函数就是矩形脉冲函数RN(n)，为了改善所设计滤波器的特性，窗函数还可以有其它的形式。 设计步骤： 设计FIR滤波器就是根据要求，找到有限个傅里叶级数系数，以有限项傅立叶级数去近似代替无限项傅里叶级数。窗函数法也称傅里叶级数法。 为理想低通滤波器的传输函数。 则滤波器的幅度特性等于理想低通滤波器的幅度特性与窗函数的幅度特性的卷积。 如果对截取后的信号进行傅里叶变换，假设采用矩形窗截取，对截取后信号进行傅里叶变换得： 线性相位 后面有证明 频域卷积定理 一. 矩形窗口法 ：低通滤波器的延时 以一个截止频率为 的线性相位理想低通滤波器为例,讨论FIR的设计问题。 1. 对于给定的理想低通滤波器 ，计算 理想特性的hd(n)和Hd(ω) （a） （b） hd(n)是一个以(N-1)/2为中心的偶对称的无限长非因果序列，如果截取一段n=0～N-1的hd(n)作为h(n)，则为保证所得到的是线性相位FIR滤波器，延时 应为h(n)长度N的一半，即 hd(n)的最大值是多少？ 其中 2. 计算 3. 计算  窗口函数的频谱 4. 加窗后的截断效应 矩形窗函数及其幅度函数（见P94图4.4） 幅度函数为 7.2.1 主瓣 旁瓣 w(n)=RN(n) 把理想频响写成幅度函数和相位函数的表示形式 其中幅度函数为 两个信号时域的乘积对应于频域卷积，所以有 则实际FIR滤波器的幅度函数H (ω) 为 正好是理想滤波器幅度函数与窗函数幅度函数的卷积。 取样函数 矩形窗 矩形窗的卷积过程（P95的图4.5来说明） 7.2.3 课本P205图7.2.3 H(0) 0.5H(0) H(ω)max H(ω)min 4 个特殊频率点看卷积结果： ，故 H(0) 近似为 在 （1）ω=0时，H(0)等于 在[-ωc, ωc]内的积分面积 因一般 [-π, π]内的积分面积； （2）ω=ωc时，仅一半波形重叠，H(ωc)= 0.5 H(0)； （3）ω=ωc –2π/N时，主瓣完全在[-ωc, ωc]内，最大 的第一旁瓣（负数）在通带外，在该点 H(ω) 出 现正肩峰； （4）ω=ωc +2π/N 时，主瓣完全在 [-ωc, ωc] 外，最