数字信号处理课件7.3.pptVIP

第七章 FIR数字滤波器的设计 7.1 引言 7.2 线性相位FIR数字滤波器的特性 7.3 窗口设计法 7.4 频率取样法 7.5 FIR数字滤波器的最优化设计 7.6 IIR与FIR数字滤器的比较 7.2窗函数设计法 一、窗函数设计法原理 一般先给出为理想频率响应 ，现要求设计一个N点( ) 的 FIR滤波器 去逼近 。 逼近方法有三种： 1、窗函数设计法（时域逼近） 2、频率采样法（频域逼近） 3、最优化设计（等波纹逼近） 7.2窗函数设计法 理想频率响应hd(n) 通过傅立叶反变换获得 7.2窗函数设计法 窗函数设计法是从单位冲激响应序列着手，使h(n)逼近理想的单位冲激响应序列hd(n)。我们知道hd(n)可以从理想频率响应 通过傅立叶反变换获得 7.2窗函数设计法 怎样用一个有限长序列h(n)来逼近无限长的hd(n)？ 最简单的办法是直接截取一段 hd(n) 代替 h(n) 。这种截取可以形象地想象为h(n)是通过一个“窗口” 看到一段hd(n) 。 因此 ， h(n)也可表示为hd(n)和一个“窗函数”的乘积，即 7.2窗函数设计法 窗函数设计法是从单位冲激响应着手，使h (n)逼近理想的单位冲激响应序列hd(n) 。 7.2窗函数设计法 二、窗函数设计法的步骤 已知理想频率响应 ，求hd(n) 。 7.2窗函数设计法 三、矩形窗函数设计方法 下面以一个截止频率为 ωc的线性相位理想低通滤波器为例，讨论FIR的设计问题。 给定的理想低通滤波器为 7.2窗函数设计法 下面是这个滤波器的理想频率特性 7.2窗函数设计法 第一步：，计算hd(n) 7.2窗函数设计法 clear;clc;close all; N_half=20; tao=10; n=-N_half:N_half+tao; wc=0.2*pi; h_d=(sin(wc*(n-tao))+eps)./(pi*(n-tao)+eps); h_d(N_half+1+tao)=wc/pi; stem(n,h_d,.) 7.2窗函数设计法 理想低通的hd(n) 7.2窗函数设计法 这是一个以为τ中心的偶对称的无限长非因果序列，如果截取一段n=0～N-1的hd(n)作为h(n) ，则为保证所得到的是线性相位FIR滤波器，延时τ应怎样选择？ 7.2窗函数设计法 第三步，检验求hd(n)的频率特性 7.2窗函数设计法 用幅度函数和相位函数来表示，则有 其线性相位部分 则是表示延时一半长度 7.2窗函数设计法 对频率响应起作用的是它的幅度函数 7.2窗函数设计法 加窗处理对矩形频率响应的影响 （1）窗函数的频率特性使矩形频率响应产生带内和带外波动。在 (在过渡带的两侧)波动最大，约为8.95%，与N无关，称为吉布斯效应。 （2）使理性频率特性不连续点处边沿加宽，形成过渡带，过渡带宽等于窗函数频率响应的主瓣宽度。 （3）增加N，只能改变窗谱的主瓣宽度、ω坐标的比例、 的绝对值大小等内容，但不能改变主瓣与旁瓣的相对比例。 7.2窗函数设计法 矩形窗函数自身频谱性能指标 （1）主瓣宽度2× （2）旁瓣峰值－13dB 加矩形窗后滤波器的性能指标 （1）过渡带宽1.8× （2）阻带最小衰减 7.2窗函数设计法 四、各种窗函数 （1） 矩形窗 7.2窗函数设计法 （2）巴特列特（Bartlett）窗（又称三角窗） 时域表达式和窗谱见课本P337 MATLAB指令：w = bartlett(n) 7.2窗函数设计法 时域和频域波形： 7.2窗函数设计法 （3）汉宁（Hanning）窗（又称升余弦窗） 时域表达式和窗谱见课本P338 MATLAB指令： 7.2窗函数设计法 汉宁（Hanning）窗时域和频域波形： 7.2窗函数设计法 （4）海明（Hamming）窗 时域表达式和窗谱见课本P338 MATLAB指令： 7.2窗函数设计法 海明（Hamming）窗时域和频域波形： 7.2窗函数设计法 （5）布拉克曼（Blackman）窗（又称二阶升余弦窗）