数字信号处理—基于计算机的方法2014cph10FIR的设计.ppt

* (P381) 以上窗函数以0为中心。需要对其进行右移。得到形式略为不同的函数形式。 * * M=25 * 矩形窗的幅度谱 W N p 2 N p 2 - N p 4 N p 4 - N p p - 主瓣 旁瓣 ) ( W W 主瓣宽度 一个窗函数的性能主要取决于它的两个参数，即主瓣宽度和相对旁瓣级（最大旁瓣与主瓣以dB为单位的幅度差值）; * Rectangular(矩形)窗: * * 0 0.25 0.5 0.75 1 -40 -30 -21 -10 0 M=14 M=30 Gain db 用矩形窗设计的Wc=p/2 FIR滤波器的幅度响应 * Bartlett (巴特利特)窗 (三角形窗): * * Hanning（汉宁）窗（升余弦窗）： 倒余弦 * * (M=38) 0 0.25 0.5 0.75 1 -80 -60 -44 -20 0 Square Hanning Gain dB 用Hanning窗设计的Wc=p/2 FIR滤波器的幅度响应 * Hamming（海明）窗（改进的升余弦窗）: * * 0 0.25 0.5 0.75 1 -80 -60 -52 -20 0 Square Hamming Gain db (M=38) 用Hamming窗设计的Wc=p/2 FIR滤波器的幅度响应 * Blackman（布莱克曼）窗: * * 0 0.25 0.5 0.75 1 -100 -75 -60 -40 -20 0 Square Blackman Gain db (M=38) 用Blackman窗设计的Wc=p/2 FIR滤波器的幅度响应 * Hanning Hamming Blackman * * 一个窗函数的性能主要取决于它的两个参数，即主瓣宽度和相对旁瓣级（最大旁瓣与主瓣以dB为单位的幅度差值）; 最大通带偏移和最小阻带值之间的距离近似等于窗的主瓣宽度（见下页图）; 为了保证从通带快速过渡到阻带，窗函数应该有一个非常小的主瓣，另一方面，为了减小通带和阻带波纹，旁瓣下的面积也要求非常小。遗憾的是，这两个要求是相互矛盾的。 * * 矩形窗 三角形窗 汉宁窗 海明窗 布莱克曼窗 凯泽窗 * 几种窗函数是各以一定主瓣加宽为代价， 来换取某种程度的旁瓣抑制。 * * 固定窗函数的特性 阻带最小衰减只由窗形状决定 过渡带宽则与窗形状和窗宽 N 都有关 * An important step in the development of a digital filter is the determination of a realizable transfer function G(z) approximating the given frequency response specifications. 73 —— 数字信号处理 —— 李 永 杰 Liyj@uestc.edu.cn 电子科技大学生命科学与技术学院 第10章 FIR数字滤波器的设计 Design of FIR Digital Filter * FIR的系统函数： 或 FIR滤波器总是稳定的！ IIR： * IIR数字滤波器： 可以利用模拟滤波器设计 但相位非线性 FIR数字滤波器： 可以严格线性相位，又可任意幅度特性 因果稳定系统 但阶次比IIR滤波器要高得多 * 设FIR滤波器单位冲激响应 h(n)长度为N，其系统函数H(z)为： 按照一定的准则选取 N 和系数 h(n). FIR滤波器的设计任务 * 10.1 预备知识 FIR数字滤波器设计的基本方法 FIR数字滤波器阶数的估计 10.2 基于加窗傅立叶级数的FIR滤波器设计 10.5 用Matlab设计FIR滤波器 主要内容 * 10.1 预备知识 FIR滤波器的设计是基于对指定幅度响应的直接逼近，通常还加上线性相位的条件限制； 一些研究者提出了一些先进的方法对给定性能指标的数字滤波器进行阶数估计。 两种直接的FIR滤波器设计方法是加窗傅立叶级数法（时域）和频率采样法（频域）； * 基于加窗傅立叶级数的FIR滤波器设计——基于对指定频率响应的傅立叶级数进行截短来设计。 为了使滤波器具有线性相位特性，h[n]需满足： * FIR滤波器阶数的估计 低通滤波器：归一化通带截止频率ωp，归一化阻带截止频率ωs，最大通带波纹δp，最大阻带波纹δs Kaiser方程： Bellanger方程： Hermann方程：(略)（P425~426） * 对通带宽度适中的FIR滤波器： 对通带宽度较窄的FIR滤波器： 对通带宽度较宽的FIR滤波器： 这样设计出来的FIR滤波器的频率响应特性可能满足也可能不满足给定的指标。如果不满足指标，建议逐渐增加滤波器阶数直到满足指标为止。 FIR滤