信号分析与处理 教学课件 作者 杨西侠 柯晶 5 1IIR.pptVIP

第5章 数字滤波器 5.1 基本概念 5.2 IIR数字滤波器设计 5.3 FIR数字滤波器设计 5.4 数字滤波器的结构 5.1 基本概念 5.2 IIR数字滤波器设计 5.2.2 冲激响应不变法 5.2.3 双线性变换法 将已知的?c = 400π，T=0.001代入，得 1. 基本原理 双线性变换的基本思路：让模拟滤波器与数字滤波器的输入、输出互相模仿，从而达到频率响应的相互模仿，也就是使数字滤波器的差分方程是模拟滤波器微分方程的近似解。 xa(t) o t ya(t) t x(n) o n digital filter 差分方程 x(n) y(n) y(n) n 模拟filter xa(t) 微分方程 ya(t) 数值近似的方法很多，例如，用差分直接代替微分，但结果使时域、频域上产生的误差大； 另一种近似方法是采用将微分方程作积分，再对积分进行数值近似（梯形近似积分）。 令数值抽样步长为T，t = nT，t0 = (n?1)T，则 用梯形法逼近积分项 双线性变换 2．稳定性与逼近程度 模拟滤波器是稳定 数字滤波器一定是稳定 设s = ? + j?，? 0，则 说明s平面的虚轴映射到z平面单位圆上，而且? = ??为? = ??，?与?呈单值映射，可以避免混叠。 对s平面的虚轴s = j?，? = 0，有? z ? = 1，且 非 线 性 ? π -π ? 单值映射 ω→0，tanω/2→ω/2 ，其频率转换关系近似成直线。 |H(e jΩ )| o Ωc Ωs Ω ?c ?s π |H(e j? )| o ? ? o Ω π ?=2arctan(ΩT)/2 双线性变换是非线性变换，频率频率特性会引起失真。双线性变换频率轴的这一非线性畸变问题，对于大量常见的一定频段内为常数幅频特性的滤波器来说，问题并不严重，如一般低通、高通、带通、带阻等，它们都要求在通带内逼近一个衰减为零的常数特性，阻带逼近衰减趋向无穷大的特性。通过双线性变换后，虽然频率发生了非线性畸变，但频率响应在通带和阻带的频率响应仍接近常数幅频特性，只是通带边界频率、阻带始点频率发生了非线性的变化。 解决这一问题可以采用所谓的“预畸”的方法，在指标转换过程中，不按线性关系求，而按非线性关系来求， 这样就保证了通过双线性变换后，使所设计的模拟滤波器截止频率正好映射在所要求的数字截止频率上。 * * 信号处理最广泛的应用是滤波。数字滤波，是指输入、输出均为离散时间信号，利用离散系统特性对输入信号进行加工和变换，改变输入序列的频谱或信号波形，让有用频率的信号分量通过，抑制无用的信号分量输入。 或者说，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的算法。因此，从概念上说，数字滤波与模拟滤波相同，只是信号的类型和实现方法不同。 数字滤波器相对模拟滤波器而言，在体积、重量、精度、稳定性、可靠性、存储功能、灵活性以及性能价格比等方面都显示明显的优点，而且，数字滤波器除利用硬件电路实现之外还可借助计算机以软件编程方式实现。 A/D xa(t) x(n) 数字滤波 H(z) y(n) D/A ya(n) 数字滤波器是一类重要的离散时间系统。可以从时域、复数域和频域的角度来描述。 一般离散时间系统可表示为N阶差分方程，即 系统函数可表示为 h(n) H(z) H(e j?) y(n) Y(z) Y(e j?) x(n) X(z) X(e j?) 如果x(n)，y(n)和h(n)满足绝对可积条件，则相应地X(z)、Y(z)和H(z)均在单位圆上收敛，离散时间系统频率响应为 ① 若1 ? r ? M，br = 0，即只有b0 ? 0时 系统只含有N个极点，无有限零点，记作AR模型。这种系统的单位抽样响应为无限长序列，习惯上称为无限冲激响应（IIR）离散系统。在系统辨识中，通常也称为自回归模型。 ② 若1 ? k ? M，ak = 0，有 系统只含有M个零点，无有限极点，记作MA模型。这种系统的单位抽样响应为有限长序列，习惯上称为有限冲激响应（FIR）离散系统。在系统辨识中，通常也称为滑动平均模型。 ③ 系统同时具有零点和极点，记作ARMA模型也称为自回归滑动平均模型。 【例5-1】一数字滤波器的差分方程为 试求其单位抽样响应h(n)及频率响应。 解：由差分方程可得 h(n) = z ?1 [H(z)] = 0.5?(n) ? 0.5?(n?1) 显然h