北京邮电大学《数字信号处理》门爱东-DSP05-FIR 数字滤波器设计和实现.pptVIP

门爱东教授menad@bupt.edu.cn 数字信号处理 Digital Signal Processing 5.1 概述：IIR 和 FIR 比较 IIR与FIR性能特性比较 IIR数字滤波器： 幅频特性较好；但相频特性较差； 有稳定性问题； FIR数字滤波器： 可以严格线性相位，又可任意幅度特性 因果稳定系统 可用 FFT 计算 (计算两个有限长序列的线性卷积) 但阶次比 IIR 滤波器要高得多 IIR 与 FIR 设计方法比较 IIR DF： 无限冲激响应，H(Z) 是 z-1 的有理分式，借助于模拟滤波器设计方法，阶数低（同样性能要求）。其优异的幅频特性是以非线性相位为代价的。 缺点：模拟原型法只能设计特定类型的滤波器，不能逼近任意的频响。 信号通过线性滤波器时，其幅度和相位可能会发生改变，滤波器幅频特性 |H(ω)|和相频特性 θ(ω) 可能会随频率的变化而改变。 如：输入正弦信号 Acos(nω0) 则：输出为 |H(ω0)| Acos(nω0＋θ)，其中相移θ＝θ(ω0) 输出频率和输入频率相同，但幅度和相位都发生了变化 输出信号比输入信号滞后的样点数 n (相位延迟) 可由下式求得: 设：nω0＋θ＝0 5.2 线性相移FIR DF 约束条件和频率响应 三个内容： 约束条件 恒延时滤波 h(n) 偶对称：恒相延时和恒群延时同时成立 h(n) 奇对称：仅恒群延时成立 频率响应 Type I：h(n) 偶对称、N 为奇数 Type II：h(n) 偶对称、N 为偶数 Type III：h(n) 奇对称、N 为奇数 Type IV：h(n) 奇对称、N 为偶数 FIR DF 零极点分布 5.2.1 线性相移FIR DF 约束条件：恒延时滤波 5.2.1 线性相移FIR DF 约束条件：恒延时 5.2.1 线性相移FIR DF 约束条件：恒延时 5.2.1 线性相移FIR DF 约束条件：恒群延时 5.2.1 线性相移 FIR DF 约束条件 5.2.2 线性相移 FIR DF 频率响应：Type I 5.2.2 线性相移 FIR DF 频率响应：Type I 5.2.2 线性相移 FIR DF 频率响应：Type I 由此可以看出其线性相位特性。由于 cos(nω) 对于 ω=0、π、2π都是偶对称，所以振度响应 Hr(ω) 对ω=0、π、2π也是偶对称。 5.2.2 线性相移 FIR DF 频率响应：Type II h(n) 偶对称，N 为偶数（恒相时延、恒群时延〕 由于h(n) 序列的长度为偶数，因此滤波器的频率响应函数可拆分成如下两部分（前后对称部分，中心点处无值）： 5.2.2 线性相移 FIR DF 频率响应：Type II 5.2.2 线性相移 FIR DF 频率响应：Type II 与所设计的 b(n) 或 h(n) 无关，恒为 0。这种类型（即 h(n) 偶对称，N为偶数）不能用于高通或带阻滤波器。 2）由于 cos[(n-1/2)ω] 对于 ω=π是奇对称，所以，Hr(w) 对 ω=π也是奇对称；以ω =0、2π为偶对称。 线性相移 FIR DF 的零极点： 5.3 FIR DF 窗口法（傅里叶级数法） 5.3.1 窗口法：基本原理 设所要求的 DF 的频率响应是 Hd(ejω)，需要注意：它可能是低通、高通、带通和带阻 FIR DF，没有特指某种类型的数字滤波器。 不管是何种 FIR DF，它的频率响应是频域中的周期函数，周期为 2π，所以它可以展开为傅氏级数形式： 5.3.1 窗口法：基本原理 因此，所要求的 DF 的系统函数便可求得： 显然，Hd(z) 是非因果的，且 hd(n) 的持续时间为 -∞～ +∞，物理上不可实现。 可以采用逼近 Hd(ejω) 的方法 首先把 hd(n) 截短为有限项，为 2M+1 项，得： 5.3.1 窗口法：基本原理 然后把截短后的 hd(n) 右移，使之变成因果性的序列。 令 H(z) 等于 H1(z) 乘以 z-M 得： 令 h(n)= hd (n-M), n=0, 1, 2, ......, 2M，则 5.3.2 窗口法：性能分析 |H(ejω)| 对 |Hd(ejω)| 的逼近 5.3.2 窗口法：性能分析 5.3.2 窗口法：性能分析 5.3.2 窗口法：Gibbs 效应 5.3.2 窗口法：Gibbs 效应 5.3.2 窗口法：Gibbs 效应 5.3.2 窗口法：Gibbs 效应 5.3.2 窗口法：Gibbs 效应 现在分析几个特殊频率点的滤波器性能：