dsp7无限长单位冲激响应数字滤波器的设计方法研讨.ppt

* 第7章 无限长单位冲激响应数字滤波器的设计方法 7.1 线性相位FIR数字滤波器的特点 7.2 FIR数字滤波器的窗函数设计法 7.3 频率抽样法设计FIR数字滤波器 7.4 IIR和FIR数字滤波器的比较 7.1 线性相位FIR数字滤波器的特点 7.1.1 FIR数字滤波器具有线性相位的条件 对于离散线性时不变系统，系统频率响应 可表示为 其中 称为滤波器的幅度响应， 称 为相位响应。 幅度响应反映了滤波器对信号通过该滤波器后各频率成分 衰减的强弱，而相位响应则反映了各频率成分通过滤波器 时间上延时的多少。 滤波器具有线性相位是指相位响应 与频率 呈线性关系， 满足： 或者 其中 都是常数。 当 也标志滤波器具有线性相位。 对于FIR滤波器，设其单位冲激响应为h(n)，长度为N，则 对应的系统函数为： 下面讨论有限长单位冲激响应h(n)为实序列，并关于 偶对称或者奇对称两种情况的相位特性。 1. 有限长单位冲激响应h(n)为偶对称 令 则有 因此，系统的频率响应 幅度响应： 相位响应： 2. 有限长单位冲激响应h(n)为奇对称 幅度响应： 相位响应： 结论：当FIR滤波器的单位冲激响应h(n)为实序列，具有 偶对称性或者奇对称，那么该FIR滤波器具有线性相位 特性。 7.1.2 线性相位FIR数字滤波器的幅度特点 1. h(n)偶对称，N为奇数 由于N为奇数,中间项为 , 其余项偶对称 , 关于 偶对称，可实现任意形式滤波器。 2. h(n)偶对称，N为奇数 FIR滤波器的幅度响应如右 图所示。 幅度响应具有以下特点： (1)当 时， 对 呈奇对称， 对 呈偶对称。 因此，具有h(n)偶对称，N 为偶数的FIR滤波器不能用 于高通滤波器或者带阻滤 波器。 3. h(n)奇对称，N为奇数 幅度响应 具有以下特点： (1)当 时， 对 呈奇对称。 因此，具有h(n)奇对称，N为奇数的FIR滤波器只能实现带通滤波器。 4. h(n)奇对称，N为偶数 幅度响应 具有以下特点： (1)当 时， 对 呈偶对称。 因此，具有h(n)奇对称，N为偶数的FIR滤波器不能实现低通、带阻滤波器。 对 呈奇对称 7.1.3 线性相位FIR数字滤波器零点分布特点 线性相位FIR滤波器的系统函数： 线性相位FIR数字滤波器零点分布特点：互为共轭对称 线性相位FIR滤波器零点分布图 7.2 FIR滤波器的窗函数设计法 7.2.1 设计原理 1.根据滤波器技术要求直接设计一个理想数字滤波器的频率响应。 2.求具有上述频率响应的滤波器对应的冲激响应序列理想频率响应的滤波器具有矩形频率特性，因而其时域响应为无限长的。 3.将求得的冲激响应序列截短，并使之具有对称性。 下面以低通滤波器为例说明设计过程 理想低通数字滤波器的频率响应为 式中 表示群延时，其单位冲激响应为： 从上式可以看出，理想低通滤波器的单位冲激响应 是 无限长非因果序列，且具有偶对称性，其中心点在 如下图所示。 理想单位冲激响应 用一个有限长的窗函数序列（比如说矩形窗）来截取 截取后的实际单位冲激响应h(n)如下图所示，它是一个以 为对称点的对称图形。 实际单位冲激响应 这时FIR滤波器的频率响应为： 从上式中可以看出，FIR滤波器的频率响应实际上是 的傅里叶级数的有限项之和，所以窗函数法也称为傅里叶 级数法，显然，N越大， 和 的差别就越小，但所 需的计算量就越大。 7.2.2 窗函数设计法的截断效应 根据卷积定理可知，时域相乘则频域相卷，实际FIR滤波器 的频率特性应为： 式中 是窗函数w(n)的傅里叶变换 对于理想低通滤波器： 为了保证截取的h(n)中心对称，需保证其相移 满足 实际滤波器的频率响应为： 说明：1.加窗后的滤波器仍然具有线性相位； 2.加窗后滤波器的幅度特性等于理想低通滤波器的 幅度特性 与矩形窗的幅度特性 的卷积， 即 理想