CH3无限长单位脉冲响应滤波器[IIR]的设计方法.pptVIP

第３章　 无限长单位脉冲响应(IIR) 数字滤波器的设计方法 内容 §3.2 常用模拟低通滤波器（LPF）特性 §3.1 根据模拟滤波器设计IIR滤波器 §3.3 从模拟LPF原型到各种滤波器的频率变换 §3.4 从数字LPF原型到各种数字滤波器的频率变换 引 言 许多信息处理过程，如 信号的过滤，检测、预测等 都要用到滤波器，数字滤波器是数字信号处理中使用得最广泛的一种线性系统环节，是数字信号处理的重要基础。 数字滤波器的功能是将一组输入的数字序列通过一定的运算后转变为另一组输出的数字序列。 数字滤波器——线性时不变系统。 实现方法主要有两种： 数字信号处理机 计算机软件 设计数字滤波器的步骤： 一般包括以下三步： (1)按照任务的要求，确定滤波器的性能指标任务包括： 需要滤除哪些频率分量 保留哪些频率分量 保留的部分允许有多大的幅度或相位失真 (2)用一个因果稳定的离散线性时不变系统的系统函数去逼近这一性能要求 系统函数可以分为IIR和FIR两类系统 DF的分类（补充） 系统函数 递归系统 非递归系统 系统响应 IIR FIR 频率响应 高通、低通 带通、带阻 (3)数字滤波器的实现 选择运算结构 确定运算和系数存储的字长 选用 通用计算机及相应的软件 专用数字滤波器硬件实现这一系统。 数字滤波器的数学描述 1)差分方程 2)系统函数 系统的组成 一般，M≤N,这类系统称为N阶系统 当M＞N时，H(z)可看成是一个N阶IIR子系统与一个(M-N)阶的FIR子系统的级联。 IIR (N阶) FIR (M-N阶) X(n) y(n) 数字低通滤波器频率响应幅度特性的容限图 |ω|≤ωc ωr≤|ω|≤π IIR滤波器的逼近问题 寻找滤波器的各系数ai和bi，使其逼近一个所要求的特性。 通常有以下两种方法： 1.模拟滤波器=预定指标的DF。 优点：经典方法，成熟，方便，准确 要求：掌握 2.最优化设计方法 (1)确定一种最优准则 例: 设计出的实际频率响应的幅度特性|H(ejω)|与所要求的理想频率响应|Hd(ejω)|的均方误差最小准则，或最大误差最小准则 (2)求此准则下的滤波器系数ai和bi。 特点： (1)直接法: 不需要模拟滤波器作为之间环节 (2)现代方法，需要大量的迭代运算 3.1 根据模拟滤波器来设计IIR滤波器 教材p100 从已知的模拟滤波器传递函数Ha(s)设计DF传递函数H(z)。 由s平面到z平面的变换，满足两条基本要求： (1)H(z)的频响要能模仿Ha(s)的频响 s平面的虚轴映射到z平面的单位圆ejω上。 (2)因果稳定的Ha(s)能映射成因果稳定的H(z) S平面的左半平面Re[s]＜0=z平面的单位圆内|z|＜1 模拟滤波器的设计方法 1.巴特沃兹滤波器（了解） 2.切比雪夫滤波器（了解） 3.椭圆滤波器（*） 两种映射方法 由DF模仿AF的特性，也即从AF映射成DF的问题。 1.脉冲响应不变法（重点） 2.双线性变换法（重点） 各有优势 3.1.1 脉冲响应不变法 教材P101 使DF的h(n)模仿模拟滤波器的ha(t)， T为采样周期。即 h(n)=ha(nT) (1) 假设 Ha(s)=L[ha(t)]，H(z)=Z[h(n)] (2) 从Ha(s)到H(z) 结论 通过N阶模拟滤波器的Ha(s)的Ai和si，可以求出DF的H(z)。 (3) (4) (5) 教材P101 ∴理想采样的拉氏变换与采样序列h(n)的z变换H(z)之间存在着s平面与z平面的映射关系 (3.7) z的模r仅对应于s的实部σ z的幅角ω仅对应于s的虚部Ω。 (1)σ=0时，r=1，s平面虚轴映射为z平面的单位圆。 (2)σ＜0时，r＜1；当σ＞0时，r＞1。s左半平面映射为z平面的单位圆内部，而s右半平面则映射为z平面单位圆外部 (3)ω=ΩT，当ω自0至±π变化时，Ω的对应值为0至±π/T。 s平面上每一条宽为2π/T的横带部分，都将重叠地映射到z平面的整个平面上。 每一横带的左半部分映射在z平面单位圆内；横带的右半部分映射在单位圆以外 jΩ轴映射到单位圆上，但jΩ轴上的每一段2π/T都对应于绕单位圆一周。 图3.1 脉冲响应不变法的映射关系 混叠失真不可避免 任何一个实际的模拟滤波器的频响都不可能真正是带限的，这就不可避免地会产生混叠失真 因而AF的频响在折叠频率以上处衰减越大，这个失真就越小。 (1)稳定性 (1)如果Ha(s)是稳定的，即其极点全部都在s左半平面内 由映射关系可知，对应的H(z)的极点也全部都在z平面