第4章 无限脉冲响应数字滤波器设计.pptVIP

第4章 无限脉冲响应数字滤波器的设计 4.1滤波器的基本原理 4.2模拟滤波器的设计方法 4.3根据模拟滤波器设计IIR滤波器 4.4从模拟滤波器低通原型到各个数字滤波器的频率变换 4.5从低通数字滤波器到各种数字滤波器的频率变换 4.6IIR数字滤波器的最优化设计方法 离散系统可以用差分方程来描述 两边求Z变换得： 整理后得系统函数： 系统函数 M≤N,只要有一项bi ≠0,i＞o,对应的h(n）就无限长，那么该系统就是IIR系统。当bi =0,i＞o就是FIR系统。 4.1.2数字滤波器的特性 我们通常用的数字滤波器一般属于选频滤波器。假设数字滤波器的传输函数H(e jω)用下式表示： 通带内和阻带内允许的衰减一般用dB数表示，通带内允许的最大衰减用αp表示，阻带内允许的最小衰减用αs表示，αp和αs分别定义为： 3. 数字滤波器设计方法概述 IIR滤波器和FIR滤波器的设计方法是很不相同的。IIR滤波器设计方法有两类，经常用的一类设计方法是借助于模拟滤波器的设计方法进行的。其设计步骤是：先设计模拟滤波器得到传输函数Ha(s)，然后将Ha(s)按某种方法转换成数字滤波器的系统函数H(z)。还有一种是最优化设计，也叫直接设计。 4.3 根据模拟滤波器IIR滤波器 为了保证转换后的H(z)稳定且满足技术要求，对转换关系提出两点要求： (1) 因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。 (2)数字滤波器的频率响应模仿模拟滤波器的频响，s平面的虚轴映射z平面的单位圆。 4.3.1脉冲响应不变法 本节首先如何有模拟滤波器的得到特性相近的数字滤波器，从模拟滤波器映射成数字滤波器的问题，主要映射方法有两种：脉冲映射不变法和双线性变换法。 设模拟滤波器的传输函数为Ha(s),相应的单位冲激响应是ha(t) 式中u(t)是单位阶跃函数。对ha(t)进行等间隔采样，采样间隔为T，得到： 和h(z）变换的关系 对 进行拉氏变换，得到: 上式表示采样信号的拉氏变换与相应的序列的Z变换之间的映射关系可用下式表示： 我们知道模拟信号ha(t)的傅里叶变换Ha(jΩ)和其采样信号 的傅里叶变换 之间的关系满足(1.5.5)式，重写如下： 上式表明将模拟信号ha(t)的拉氏变换在s平面上沿虚轴按照周期Ωs=2π/T延拓后，再按照(6.3.6)式映射关系，映射到z平面上，就得到H(z)。(6.3.6)式可称为标准映射关系。下面进一步分析这种映射关系。设 那么 σ=0,r=1 σ0,r1 σ0,r1 另外，注意到z=esT是一个周期函数，可写成 4.3.2用双线性变换法 正切变换实现频率压缩： 再通过 转换到z平面上，得到： 下面分析模拟频率Ω和数字频率ω之间的关系。 令s=jΩ,z=e jω，并代入(6.4.3)式中，有 例6.4.1试分别用脉冲响应不变法和双线性不变法将图6.4.4所示的RC低通滤波器转换成数字滤波器。 解 首先按照图6.4.4写出该滤波器的传输函数Ha(s)为 利用双线性变换法转换，数字滤波器的系统函数H2(z)为 下面我们总结利用模拟滤波器设计IIR数字低通滤波器的步骤。 (1)确定数字低通滤波器的技术指标：通带截止频率ωp、通带衰减αp、阻带截止频率ωs、阻带衰减αs。 (2)将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。 (3)按照模拟低通滤波器的技术指标设计模拟低通滤波器。 (4)将模拟滤波器Ha(s)，从s平面转换到z平面，得到数字低通滤波器系统函数H(z)。 §4．4 从模拟滤波器低通原型到