数字信号处理及MATLAB实现数字滤波器设计方法.ppt

数字信号处理 第六章　无限长单位冲激响应（IIR）数字滤波器的设计方法 6.1　引言 6.2　最小与最大相位延时系统，最小与最大相位超前系统 6.3　全通系统 6.4　用模拟滤波器设计IIR数字滤波器 6.5　冲激响应不变法 6.6　阶跃响应不变法 6.7　双线性变换法 本章主要讨论逼迫性性能要求或系统函数的设计问题。 数字滤波器按频率特性划分为为低通、高通、带通、带阻、全通等。其理想幅度频率响应如图6.1。 　 数字滤波器在复频率响应　　下三个参量分析： 1.　幅度平方响应 (6.3) 由于　　　　　的极点既是共轭的，又是以单位圆对称的，故只取单位圆内的极点作为　　的极点。 如果选　　　　　在z平面单位圆的零点作为　　零点，则得到是最小相位延迟滤波器。 2.相位响应 由于 　　　　　　　　　　　　　　　　　　（6.4） 所以　　　　　　　　　　　　　　　　（6.5） 由于 所以又有 　　　　　　　　　　　　　　　　　　（6.6）　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　 3.群延迟响应 它是滤波器平均延迟的一个度量 　　　　　　　　　　　　　　　　　（6.7） 可以化为 　　　　　　　　　　　　　　　　　（6.8） 由于 所以 因而又有 　　　　　　　　　　　　　　　　　（6.9） 同样可化为 　　　　　　　　　　　　　　　（6.10） 当滤波器为线性相位响应特性时，则通带内延迟特性为常数。 IIR滤波器逼迫问题 IIR系统函数 　　　　　　　　　　　　　　　（6.11） 就是去求出滤波器的各系数　　，使得在规定意义上，逼近所要求特性。在z平面上逼近，得到数字滤波器。 设计方法： （1）先设计一个合适的模拟滤波器，然后变换成满足指标的数字滤波器。 （2）计算机辅助设计方法。 6.2　最小与最大　相位延时系统，最小与最大相位超前系统 任一个线性时不变系统，其系统函数为 　　　　　　　　　　　　　　　　（6.12） 系统频率响应表达式 　　　　　　　　　　　　　　　　（6.13） 我们对　　进行研究。 　　　的模为 　　　　　　　　　　　　　　　　（6.14） 其相角为 　　　　　　　　　　　　　　　　（6.15） 若　　　　　分别表示单位圆内外的零极点数则 注：零（极）矢是指零（极）点指向z平面单位圆上频率点的矢量。当　从0到　　时，只有单位圆内的零极点对相角有影响。 零点、极点的分布对系统相角影响的讨论 （1）对因果稳定系统，此时 当　从0变到　　时，　　　　　则辐角变化量 它称为相位延时系统。 a 当全部零点在单位圆内，即　　　　　　有 其相位变化最小，称为最小相位系统。 b 当全部零点在单位圆外，即　　　　　　　有 这时相位变化最大，又是负数，即最大相位延时系统，是因果稳定系统。 （2）对逆因果移动系统，此时 当　从0变到　　时，　　　　，则辐角变化量为 称为相位超前系统。 a 当全部零点在单位圆内时，即　　　　　有 这时相位变化最大，称最大相位超前系统，是逆因果稳定系统。 b 当全部零点在单位圆外时即　　　　　　有 这时相位超前最小，称最小相位超前系统，是逆因果移动系统。 表6.1　四种系统及其因果性、稳定性、零点、极点的关系。 最小相位系统重要性质 （1）在傅里叶变换　　相同的所有系统中，它的负相位最小。 （2）最小相位系统能量集中在n=0附近。 （3）相同傅里叶变换幅度的各序列，最小相位序列的　　　最大。 （4）在幅度响应　　　相同的系统中，只有唯一的一个最小相位延时系统。 （5）利用级联全通的方法，可将最小相位系统的零点反射到单位圆外，构成幅度响应相同的非最小相位延时系统。 6.3　全通系统 定义：是指系统频率响应的幅度在所有频率下均为1或某一常数的系统。满足 简单一阶全通系统函数 A为实数　　（6.16） 零极点如图6.3 高阶全通系统包含实零点－实极点系统；还包括复数零点－极点系统；复数零点－极点的全通节的系统函数 　　　　　　　　　　　a为实数　（6.17） 如图6.4所示。 当h(n)是实函数，因而其系统函数的复数极点零点必须共轭出现。实系数有理二阶全能系统函数 　　　　　　　　　　　　　　　（6.18） 如图6.5。 一般来说，N阶数字全通系统的系统函数频率响应的模都为1。 证明：N阶全通系统函数为 式中 当　　　时，满足 　　　　　　　　　　　　　　　　（6.21） 所以有　　　　　　　　　　　　　（6.22） 全通系统应用 （1）任何一个因果稳定的（非最小相位延时）系统的H(z)都可以表示