IIR滤波器的理论与设计培训讲义.pptVIP

§6.1 IIR数字滤波器的设计概述 一、设计步骤 1、按照实际任务确定滤波器的性能要求； 2、用一个因果稳定的系统函数去逼近这一要求 （IIR、FIR两类系统）； 3、利用有限精度算法来实现这个系统函数 （选择运算结构，选择合适字长，选择有效的 数字处理方法）； 4、实际的技术实现（软件、硬件或二者结合）。 二、设计方法 H(z)的设计就是要确定零、极点ci、di，以使滤波器满足给定的性能要求。 三、设计原理 利用模拟滤波器来设计数字滤波器要先根据给出的性能指标设计出相应的模拟滤波器系统函数Ha(s)，然后由Ha(s)经变换得到所需的数字滤波器的系统函数H(z)。 时域变换法有：冲激响应不变法、阶跃响应不变法、匹配z变换法； 频域变换法有：双线性变换法、微分映照法 §6.2 模拟滤波器的数字仿真 数字仿真就是要设计出数字滤波器Ld，当其输入为模拟滤波器La输入x(t)取样时，输出也为La输出y(t)的取样。 下面从时域与频域两方面分析实现数字仿真的条件 一、时域的数字仿真 设有一模拟滤波器（线性非移变的因果系统），x(t)与y(t)关系为 现对输出y(t)抽样，令t=nT代入上式 二、频域的数字仿真 取样信号的频谱X^a(W) 是原模拟信号的频谱Xa(W)的周期延拓。 ha(t)也是一模拟信号，对其取样即有 如上图，如果 Ha(W) 在－Ws/2 ～Ws/2之外为零，则在此区间内H(ejω)与Ha(W)完全一致，否则H(ejω) 将产生混叠失真。 综上所述： 若有一个脉冲响应为ha(t)，频率响应为Ha(W)的模拟滤波器， 当|W|Wm时，Ha(W)=0，则可得到一个与之仿真的数字滤波器； 其脉冲响应为h(n)=Tha(nT)，抽样周期满足 Ws=2p/T≥2Wm ； 频率响应为H(ejw) ，是 Ha(W) 的周期延拓，且在－Ws/2 ～Ws/2 区间内与Ha(W)完全一致。 §6.3 脉冲响应不变法（时域变换） 有了刚才的基础（数字仿真），如何在给定模拟滤波器系统函数Ha(s)的情况下，求出数字滤波器系统函数H(z)。Ha(s)是模拟滤波器冲激响应的拉氏变换 一、设计方法 模拟滤波器系统函数 化为部分分式合形式 交换求和顺序 例6-3-1：设模拟滤波器的系统函数为 特性图 令z=gejw（极坐标）s=s+jW（复坐标），代入式 z=esT 故有 g=esT w=WT； 即可知z平面的模和幅角与s平面实部与虚部的关系。 2、w与W的关系 w=WT ⑴ W=0（s 平面实轴）对应于 w=0（ z平面正实轴） 将二者关系合并 三、特点 ⑴优点：由于是根据h(n)=Tha(nT)进行设计，保持了模拟滤波器的时域瞬态特性； ⑵数字滤波器频率w与模拟滤波器频率W之间呈线性关系，w=WT，不会发生线性失真； ⑶缺点：如果Ha(W)不是在 -p/T～p/T 严格带限，则会产生混叠失真。实际上，任何一个模拟滤波器其频响不可能是真正限带的，但只要在p/T以上锐减，保证其混叠失真足够小，仍能满足工程实际需要； ⑷此法不能设计高通和带阻滤波器。 §6.4 双线性变换法（频域变换） 它是为克服混叠失真即多值映射这一缺点而提出的 一、变换原理 首先把整个s平面压缩到某一中介平面（s1平面的）一条横带里（宽度为2π/T，从-π/T～π/T ）再由标准变换关系z=es1T，将此横带变换到整个z平面上，这样使s平面与z平面是一一对应关系，可消除混叠失真。 首先将s平面上jΩ轴压缩到s1上jΩ1轴一段[-π/T～π/T ] 压缩关系为 将z=es1T代入 系数c的确定 在低频处，为使Ω≈Ω1（使模拟滤波器与数字滤波器在低频处有较好的确切的对应关系） 即当Ω1较小时，有 可求得 例6-4-1：已知模拟滤波器Ha(W)=1/(1+s)，利用双线性变换法将其转换成数字滤波器H(z)。（设T=1） 解： Ha(jΩ)与H(ejw)幅频特性如图 三、特点 1、优点：模拟滤波器经过双线性变换后，不存在频率特性的混叠失真。因而对Ha(Ω)无限带要求，且能直接用于设计低通、高通、带通、带阻等各种类型滤波器 2、缺点：由于映射关系为w= 2arctan(WT/2) ，当W很小时，近似呈线性关系，当W增加，变换关系就是非线性的了。会产生相频特性失真。 怎样预畸变？ 书例6.1(P124)用双线性变换法设计一个三阶Butterworth数字低通滤波器，采样频率为 fs=1.2KHZ，截止频率为fc=400HZ； 数字滤波器截止频率 wc=WcT=2pfc/fs=2p×400/1200=2p/3 用双线性变换法，则相应的模拟滤波器截止频率为 说明：各种频率间关系 1、在考虑数字滤波器