数字信号处理Digital Signal Processing Lecture 18-FIR数字滤波器设计.ppt

4、线性相位FIR低通滤波器的设计 解： 1）求数字频率 例：设计一个线性相位FIR低通滤波器， 给定抽样频率为 ， 通带截止频率为 ， 阻带起始频率为 ， 阻带衰减不小于-50dB，幅度特性如图所示 2）求hd(n) 4）确定N 值 3）选择窗函数：由 , 确定海明窗（-53dB） 5）确定FIR滤波器的h(n) 6）求 ，验证 若不满足，则改变N或窗形状重新设计 5、线性相位FIR高通滤波器的设计 其单位抽样响应： 理想高通的频响： 6、线性相位FIR带通滤波器的设计 其单位抽样响应： 理想带通的频响： 7、线性相位FIR带阻滤波器的设计 其单位抽样响应： 理想带阻的频响： 五、IIR和FIR数字滤波器的比较 IIR 滤波器 FIR 滤波器 h(n)无限长 h(n)有限长 极点位于z平面任意位置 滤波器阶次低 非线性相位 递归结构 不能用FFT计算 可用模拟滤波器设计 用于设计规格化的选频滤波器 极点固定在原点 滤波器阶次高得多 可严格的线性相位 一般采用非递归结构 可用FFT计算 设计借助于计算机 可设计成各种幅频特性和相频特性的滤波器 离散时间信号—序列 数字信号处理DIGITAL SIGNAL PROCESSING 学习目标 掌握FIR数字滤波器窗函数设计法 理解IIR与FIR数字滤波器各自的优缺点 数字IIR与FIR滤波器设计方法的比较 IIR数字滤波器： 可以利用模拟滤波器的设计结果 但一般相位具有非线性性 FIR数字滤波器： 严格设计成线性相位， 可设计成任意幅度响应特性 是因果稳定系统 可用FFT计算 但阶次比IIR滤波器要高得多 FIR滤波器具有严格线性相位的条件 FIR滤波器单位抽样响应h(n)为实数， 且满足： 偶对称： 或奇对称： 即对称中心在 (N-1) / 2处 则FIR滤波器的频率响应具有严格线性相位。 窗函数设计法 1、设计方法——已知 ，求逼近此频率响应的 w(n)：窗函数序列 需要选择合适的形状和长度 以低通滤波器为例讨论： 线性相位理想低通滤波器的频率响应： 其理想单位抽样响应： 是中心点为 α 的偶对称、无限长、非因果序列 取矩形窗： 则 FIR 滤波器的单位抽样响应： 按第一类线性相位条件，得 加窗处理对频率响应的影响： 时域乘积相当于频域卷积 而矩形窗的频率响应： 理想低通滤波器的频率响应： 其幅度函数： 则FIR滤波器的频率响应： 幅度函数： 加窗函数对频率响应的影响： 不连续点处边沿加宽, 形成过渡带，过渡带宽度（两肩峰之间的宽度）等于窗函数频率响应的主瓣宽度。 在 处出现肩峰值，两侧形成起伏振荡，振荡的幅度和大小取决于旁瓣的幅度和大小。称为 Gibbs效应 改变 N 只能改变窗谱的主瓣宽度，但不能改变主瓣与旁瓣的相对比例。主瓣与旁瓣的相对比例由窗函数形状决定。 2、各种窗函数的选择 对窗函数的要求： 窗谱主瓣尽可能窄,以获得较陡的过渡带 尽量减少窗谱最大旁瓣的相对幅度以减小肩峰和波纹 矩形窗函数 主瓣宽度最窄： 旁瓣幅度大 窗谱： 幅度函数： 三角形（Bartlett）窗 主瓣宽度宽： 旁瓣幅度较小 窗谱： 幅度函数： 汉宁（Hanning）窗 （升余弦窗） 主瓣宽度宽： 旁瓣幅度小 幅度函数： 海明（Hamming）窗 （改进的升余弦窗） 主瓣宽度宽： 旁瓣幅度更小 幅度函数： 布莱克曼（Blackman）窗 （二阶升余弦窗） 主瓣宽度最宽： 旁瓣幅度最小 幅度函数： 凯泽（Kaiser）窗 ：第一类变形零阶 贝塞尔函数 窗函数 窗的谱性能指标 加窗后滤波器性能指标 旁瓣峰值 /dB 主瓣宽度 过渡带宽 阻带最小衰减 /dB 矩形窗 -13 2 0.9 -21 三角形窗 -25 4 2.1 -25 汉宁窗 -31 4 3.1 -44 海明窗 -41 4 3.3 -53 布拉克曼窗 -57 6 5.5 -74 凯泽窗 -57 5 -80 1）阻带最小衰减只由窗形状决定 2）过渡带宽则与窗形状和窗宽 N 都有关 3、窗函数法的设计步骤 给定理想滤波器的频率响应函数 及技术指标 求出理想的单位抽样响应 根据阻带衰减选择窗函数 计算频