DSP数理综合.ppt

第5章 IIR数字滤波器的设计 5.1 数字滤波器的基本概念 5.2 模拟滤波器的设计 5.3 用脉冲响应不变法设计IIR数字低通滤波器 5.4 用双线性变换法设计IIR数字低通滤波器 5.5 数字高通、带通和带阻滤波器的设计 5.1 基本概念 1. 数字滤波器的分类 数字滤波器从实现的网络结构或者从单位脉冲响应分类，可以分成无限脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器。它们的系统函数分别为： 图5.1 理想低通、高通、带通、带阻滤波器幅度特性 2数字滤波器的技术要求 图5.2 低通滤波器的技术要求 5.3 用脉冲响应不变法设计IIR数字低通滤波器 为了保证转换后的H(z)稳定且满足技术要求，对转换关系提出两点要求： (1) 因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。 (2)数字滤波器的频率响应模仿模拟滤波器的频响，s平面的虚轴映射z平面的单位圆，相应的频率之间成线性关系。 设模拟滤波器的传输函数为Ha(s),相应的单位冲激响应是ha(t) 设模拟滤波器Ha(s)只有单阶极点，且分母多项式的阶次高于分子多项式的阶次，将Ha(s)用部分分式表示： (5.3.1) 式中si为Ha(s)的单阶极点。将Ha(s)进行逆拉氏变换得到ha(t)： (5.3.2) 式中u(t)是单位阶跃函数。对ha(t)等间隔采样，间隔为T，得到： (5.3.3) 对上式进行Z变换，得到数字滤波器的系统函数H(z)： (5.3.4) 设ha(t)的采样信号用 表示， 对 进行拉氏变换，得到: 式中ha(nT)是ha(t)在采样点t=nT时的幅度值，与序列h(n)的幅度值相等，即h(n)=ha(nT)，因此得到: (5.3.5) 上式表示采样信号的拉氏变换与相应的序列的Z变换之间的映射关系，可用下式表示： (5.3.6) 已知模拟信号ha(t)的傅里叶变换Ha(jΩ)和其采样信号 的傅 里叶变换 之间的关系满足： 将 s=σ+jΩ=0+jΩ=jΩ 代入上式，得 由(5.3.5)式和(5.3.8)式得到: (5.3.7) (5.3.8) (5.3.9) 上式表明：将模拟信号ha(t)的拉氏变换在s平面上沿虚轴按照周期Ωs=2π/T延拓后，再按照(5.3.6)式映射关系，映射到z平面上，就得到H(z)。(5.3.6)式可称为标准映射关系。 下面进一步分析这种 映射关系。设 按照(5.3.6)式: 最后得到： 于是有：σ=0, r=1 σ0, r1 σ0, r1 为整数 注意z=esT为周期函数，可写成： (5.3.10) 图5.3.1 z=esT, s平面与z平面之间的映射关系 假设 没有频率混叠现象，即满足 按照(5.3.9)式，将关系式s=jΩ、ω=ΩT代入得到： 图5.3.2 脉冲响应不变法的频率混叠现象 一般Ha(s)的极点si是一个复数，且以共轭成对的形式出现，(5.3.1)中将一对复共轭极点放在一起，即可构成一个二阶基本节。 极点 (5.3.11) 形式A对应的数字滤波器二阶基本节： (5.3.12) 模拟滤波器的二阶基本节形式A： 极点 (5.3.13) 形式B对应的数字滤波器二阶基本节： (5.3.14) 模拟滤波器的二阶基本节形式B： 例6.3.1 已知模拟滤波器的传输函数Ha(s)为 用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z)。 解：首先将Ha(s)写成部分分式： 极点为： 则H(z)的极点为： 按照(5.3.4)式，并经过整理，得到： T=1s T=0.1s 也可以直接按照(5.3.13),(5.3.14)式进行转换。首先将Ha(s)写成(5.3.13)式的形式，如极点s1,2=σ1±jΩ1,则 再按照(5.3.14)式，H(z)为 图5.3.3 例5.3.1的幅度特性 5.4 用双线性变换法设计IIR数字低通滤波器 正切变换实现频率压缩： (5.4.1