数字信号处理第七章.pptVIP

第七章 FIR数字滤波器的设计方法 IIR数字滤波器： 一、线性相位FIR滤波器的特点 FIR滤波器的单位冲激响应： 一、FIR滤波器具有线性相位的条件 对于长度为N的h(n),传输函数为： 1、H(ejw)线性相位 2、FIR滤波器具有线性相位的条件 二、线性相位条件的证明 我们将根据单位冲激响应h(n)的奇对称和偶对称进行讨论。 由下面讨论可知： 1、线性相位条件 h(n)为实序列时，其频率响应： 2、线性相位FIR滤波器频率响应的特点 频率响应： 频率响应： 3、幅度函数的特点 1）h(n)偶对称，N为奇数 2）h(n)偶对称，N为偶数 3）h(n)奇对称，N为奇数 4）h(n)奇对称，N为偶数 4、零点位置 二、窗函数设计法 1、设计方法 以低通滤波器为例讨论： 线性相位理想低通滤波器的频率响应： 取矩形窗： 加窗处理后对频率响应的影响： 时域乘积相当于频域卷积 理想滤波器的频率响应： 加窗函数的影响： 不连续点处边沿加宽形成过渡带，其宽度（两肩峰之间的宽度）等于窗函数频率响应的主瓣宽度。 2、各种窗函数 窗函数的要求： 三角形（Bartlett）窗 汉宁（Hanning）窗 （升余弦窗） 海明（Hamming）窗 （改进的升余弦窗） 布莱克曼（Blackman）窗 （二阶升余弦窗） 凯泽（Kaiser）窗 3、窗函数法的设计步骤 4、线性相位FIR低通滤波器的设计 2）求hd(n) 5）确定FIR滤波器的h(n) 5、线性相位FIR高通滤波器的设计 6、线性相位FIR带通滤波器的设计 7、线性相位FIR带阻滤波器的设计 五、IIR和FIR数字滤波器的比较 IIR滤波器 FIR滤波器 4）确定N 值 3）选择窗函数：由 确定海明窗（-53dB） 6）求 ，验证 若不满足，则改变N或窗形状重新设计 其单位抽样响应： 理想高通的频响： 其单位抽样响应： 理想带通的频响： 其单位抽样响应： 理想带阻的频响： 其中： h(n)为奇对称时，有900相移，适用于微分器和900移相器，而选频滤波器采用h(n)为偶对称时 得： 由 1）若 z = zi 是H(z)的零点，则 z = zi-1 也是零点 2）h(n)为实数，则零点共轭成对 线性相位滤波器的零点是互为倒数的共轭对 即共轭成对且镜像成对 1） 零点： 2） ，即零点在单位圆上 零点： 3） ，即零点在实轴上 零点： 4） 即零点既在实轴上，又在单位圆上 零点： w(n)：窗函数序列 要选择合适的形状和长度 其理想单位抽样响应： 中心点为 α 的偶对称无限长非因果序列 则FIR滤波器的单位抽样响应： 按第一类线性相位条件，得 而矩形窗的频率响应： 其幅度函数： 则FIR滤波器的频率响应： 幅度函数： 在 处出现肩峰值，两侧形成起伏振 荡，振荡的幅度和多少取决于旁瓣的幅度和多少 改变N只能改变窗谱的主瓣宽度，但不能改变主瓣与旁瓣的相对比例。其相对比例由窗函数形状决定，称为Gibbs效应 窗谱主瓣尽可能窄以获得较陡的过渡带 尽量减少窗谱最大旁瓣的相对幅度 以减小肩峰和波纹 矩形窗 主瓣宽度最窄： 旁瓣幅度大 窗谱： 幅度函数： 主瓣宽度宽： 旁瓣幅度较小 窗谱： 幅度函数： 主瓣宽度宽： 旁瓣幅度小 幅度函数： 主瓣宽度宽： 旁瓣幅度更小 幅度函数： 主瓣宽度最宽： 旁瓣幅度最小 幅度函数： ：第一类变形零阶 贝塞尔函数 -21 -25 -44 -53 -74 -80 0.9 2.1 3.1 3.3 5.5 5 2 4 4 4 6 -13 -25 -31 -41 -57 -57 矩形窗 三角形窗 汉宁窗 海明窗 布拉克曼窗 凯泽窗 阻带最小衰减 /dB 过渡带宽 主瓣宽度 旁瓣峰值 /dB 加窗后滤波器性能指标 窗谱性能指标 窗函数 阻带最小衰减只由窗形状决定 过渡带宽则与窗形状和窗宽N都有关 给定理想的频率响应函数 及技术指标 求出理想的单位抽样响应 根据阻带衰减选择窗函数 计算频率响应 ，验算指标是否满足要求 根据过渡带宽度确定N值 求所设计的FIR滤波器的单位抽样响应 公式法： IFFT法： 计算其IFFT，得： 对 M点等间隔抽样： 解： 1）求数字频率 例：设计一个线性相位FIR低通滤波器，