无限长单位冲激响应(IIR)数字滤波器的设计方法资料

1、数字(shz)信号处理 共五十七页 第六章无限长单位(dnwi)冲激响应（IIR）数字滤波器的设计方法 6.1引言 6.2最小与最大相位延时系统，最小与最大相位超前系统 6.3全通系统 6.4用模拟滤波器设计IIR数字(shz)滤波器 6.5冲激响应不变法 6.6阶跃响应不变法 6.7双线性变换法共五十七页本章主要讨论逼迫性性能要求或系统(xtng)函数的设计问题。共五十七页数字(shz)滤波器按频率特性划分为为低通、高通、带通、带阻、全通等。其理想幅度频率响应如图6.1。共五十七页数字滤波器在复频率响应下三个参量分析：1.幅度平方响应 (6.3)由于的极点既是共轭的，又是以单位圆对称的，故只

2、取单位圆内的极点作为(zuwi)的极点。如果选在z平面单位圆的零点作为零点，则得到是最小相位延迟滤波器。共五十七页2.相位(xingwi)响应由于(yuy)（6.4）所以（6.5）由于所以又有（6.6）共五十七页3.群延迟(ynch)响应它是滤波器平均延迟的一个度量（6.7）可以(ky)化为（6.8）由于所以因而又有（6.9）共五十七页同样可化为（6.10）当滤波器为线性相位(xingwi)响应特性时，则通带内延迟特性为常数。共五十七页IIR滤波器逼迫(bp)问题IIR系统函数（6.11）就是去求出滤波器的各系数，使得在规定意义上，逼近所要求特性。在z平面上逼近，得到数字滤波器。设计方法：（1

3、）先设计一个合适(hsh)的模拟滤波器，然后变换成满足指标的数字滤波器。（2）计算机辅助设计方法。共五十七页6.2最小与最大相位(xingwi)延时系统，最小与最大相位(xingwi)超前系统任一个线性时不变系统，其系统函数为（6.12）系统频率响应表达式（6.13）我们对进行(jnxng)研究。共五十七页的模为（6.14）其相角为（6.15）若分别表示(biosh)单位圆内外的零极点数则注：零（极）矢是指零（极）点指向z平面单位圆上频率点的矢量。当从0到时，只有单位圆内的零极点对相角有影响。共五十七页零点、极点的分布对系统相角影响的讨论（1）对因果稳定(wndng)系统，此时当从0变到时，则

4、辐角变化量它称为相位延时系统。a 当全部零点在单位圆内，即有其相位变化最小，称为最小相位系统。共五十七页b 当全部(qunb)零点在单位圆外，即有这时相位变化最大，又是负数，即最大相位延时系统，是因果稳定系统。（2）对逆因果移动系统，此时当从0变到时，则辐角变化量为称为相位超前系统。共五十七页a 当全部零点(ln din)在单位圆内时，即有这时相位变化最大，称最大相位超前系统，是逆因果稳定系统。b 当全部零点在单位圆外时即有这时相位超前最小，称最小相位超前系统，是逆因果移动系统。共五十七页表6.1四种系统及其因果性、稳定性、零点(ln din)、极点的关系。共五十七页最小相位(xingwi)系

5、统重要性质（1）在傅里叶变换(binhun)相同的所有系统中，它的负相位最小。（2）最小相位系统能量集中在n=0附近。（3）相同傅里叶变换幅度的各序列，最小相位序列的最大。（4）在幅度响应相同的系统中，只有唯一的一个最小相位延时系统。（5）利用级联全通的方法，可将最小相位系统的零点反射到单位圆外，构成幅度响应相同的非最小相位延时系统。共五十七页6.3全通系统(xtng)定义：是指系统(xtng)频率响应的幅度在所有频率下均为1或某一常数的系统(xtng)。满足简单一阶全通系统函数A为实数（6.16）零极点如图6.3共五十七页高阶全通系统(xtng)包含实零点实极点系统(xtng)；还包括复数零

6、点极点系统(xtng)；复数零点极点的全通节的系统(xtng)函数a为实数（6.17）如图6.4所示。共五十七页当h(n)是实函数，因而其系统函数的复数极点零点(ln din)必须共轭出现。实系数有理二阶全能系统函数（6.18）如图6.5。共五十七页一般来说，N阶数字全通系统(xtng)的系统(xtng)函数频率响应的模都为1。证明：N阶全通系统函数为式中当时，满足（6.21）所以有（6.22）共五十七页全通系统(xtng)应用（1）任何一个因果稳定的（非最小相位延时）系统的H(z)都可以表示(biosh)为全通系统和最小相位延时系统的时延。（6.24）它们频率响应的相位相同，相位不同。即（6

7、.25）共五十七页证明：设一个因果(yngu)稳定的非最小相位延时系统H(z)为（6.24）将6.24式表示为（6.25）由于，所以是最小相位延时，是全通级联。所以可表示为6.23式。共五十七页（2）如果设计出的滤波器是非稳定的，则可用级联全通函数的办法(bnf)将它变成一个稳定系统。例：原滤波器有一对极点在单位圆外级联一个全通系统则可将单位圆外极点抵消，但不改变系统幅度特性。共五十七页（3）可以作为相位均衡器（群时延均衡器）用，来得到线性相位，但不改变(gibin)幅度特性。设全通滤波器为，系统为级联后H(z)即相位关系共五十七页按得当通带中满足是常数则逼近误差(wch)的平方值共五十七页6

8、.4用模拟(mn)滤波器设计IIR数字滤波器方法：利用模拟滤波器设计数字滤波器就是要把s平面映射到z平面，使模拟系统函数变换成所需的数字滤波器的系统函数。基本(jbn)要求：（1）的频响特性应能模仿的频率响应（2）因果稳定的就能映射成因果稳定几种映射方法：冲激响应不变法、阶跃响应不变法、双线性变换法。共五十七页6.5冲激响应不变法(bin f)一、变换原理冲激响应不变法是使数字滤波器的单位冲激响应序列(xli)h(n)模仿模拟滤波器的单位冲激响应。即抽样满足（6.27）利用2.53式，得（6.28）它与2.5节从s平面到z平面2.51式一致。如图6.7共五十七页二、混叠失真(sh zhn)数字

9、滤波器和模拟滤波器频率响应之间关系（6.29）只有当（6.30）才能在折叠频率内重现模拟滤波器响应(xingyng)而不产生混叠失真，即（6.31）共五十七页注意：（1）模拟滤波器频率响应(pn l xin yn)在折叠频率以上衰减快、大，变换后混叠失真越小。（2）当增加，可减小混叠效应。（3）当滤波器指标用数字域给定时，若不变，用减小T的方法就不能解决混叠问题。共五十七页三、模拟(mn)滤波器的数字化方法冲激响应不变法所造成的s平面和z平面的对应关系设模拟滤波器的系统函数只有单阶极点，且假定分母的阶次大于分子的阶次。展开(zhn ki)部分分式（6.32）其相应的冲激响应是的拉普拉斯变换即（

10、6.33）对h(n)求z变换，得数字滤波器系统函数（6.34）共五十七页6.32式的和6.34式的H(z)比较结果分析：（1）s平面的单极点(jdin)变换到z平面上处的单极点(jdin)。（2）两者部分分式的系数是相同的，都是（3）如果模拟滤波器是稳定的，则数字滤波器也是稳定的。（4）能保证两者极点代数对应关系，但不保证整个s与z的代数对应关系。零点就没有这种关系。而是随的极点及系数两者而变化。共五十七页6.31式看出，数字滤波器频率响应还与抽样间隔(jin g)T成反比。当抽样频率高时，滤波器增益会太高。修正方法：令（6.35）则有（6.36）及（6.37）因共轭极点，其变成h(z)关系共

11、五十七页例6.1设模拟滤波器的系统函数为试利用冲激响应不变法(bin f)，设计IIR数字滤波器。解：直接利用6.36式可得数字滤波器的系统函数为设T1，则有共五十七页模拟滤波器与数字滤波器的频率响应(pn l xin yn)分别为画在图6.8上。共五十七页冲激不变响应法优缺点：数字滤波器的冲激响应完全(wnqun)模仿模拟滤波器的冲激响应，也就是时域逼近良好，模拟频率和数字频率之间呈线性关系。只适应限带的模拟滤波器。共五十七页6.6阶跃响应(xingyng)不变法该方法是使数字滤波器的阶跃响应g(n)模仿模拟(mn)滤波器的阶跃响应。使g(n)满足（6.38）因而满足将此式取z变换可得G (

12、 n ) = g ( n ) =z H ( z )/ ( z -1) （6.39）所以（6.40）共五十七页对于模拟(mn)滤波器其系统函数为，其阶跃响应如图6.10。因而满足共五十七页对上式两端取拉普拉斯变换可得（6.41）即（6.42）要满足(mnz)阶跃响应不变，则应有 （6.43）将6.43式取z变换，并代入6.40得数字滤波器系统函数（6.44）这就是阶跃响应不变法公式。共五十七页阶跃响应不变法频率响应周期延拓及混叠现象探讨(tnto)：理想抽样信号拉氏变换为（6.45）而阶跃响应序列的z变换为（6.46）比较6.45及6.46式，可知（6.47）共五十七页同样模拟(mn)阶跃信号与

13、理想抽样信号两者拉普拉氏变换关系（6.48）由6.47,48式，代入6.39,41式，得（6.49）可以看出，不是严格限带情况下仍有周期延拓及混叠失真现象，但比冲激响应不变法要小。共五十七页例6.2已知模拟二阶巴特沃思低通滤波器的归一化（3dB截止频率为1）的模拟系统函数为而3dB截止频率为50Hz的模拟滤波器，需将归一化中的s变量用来代替，即设系统抽样频率(pnl)为，要求从这一低通模拟滤波器设计一个低通滤波器，采用阶跃响应不变法。共五十七页解：按6.41式，模拟(mn)滤波器阶跃响应的拉普拉斯变换为由于共五十七页故则利用(lyng)以下z变换关系共五十七页且代入得阶跃响应的z变换(binh

14、un)再由6.40式可得数字低通滤波器的系统函数为共五十七页6.7双线性变换法目的：主要用来(yn li)克服混叠现象。共五十七页一、变换(binhun)原理双线性变换法是使数字滤波器的频率响应与模拟滤波器的频率响应相应的一种变换方法。首先把整个s平面压缩变换到某一中介(zhngji)的平面的一条横带里，其次再通过上面讨论过的标准变换关系将此横带变换到整个z平面上去，这样就使s平面与z平面一一对应来消除混叠。共五十七页将s平面整个 轴压缩变换(binhun)到上的到一段变换(binhun)关系：（6.50）这样，变到变到，可将6.50式变换成解析延拓到s平面和平面，令则（6.51）共五十七页再

15、将平面通过以下标准变换关系映射到z平面（6.52）从而得到s平面和z变换的单值映射关系（6.53）（6.54）为使模拟滤波器的某一频率(pnl)与数字滤波器的任一频率(pnl)有对应关系引入待定常数c，使50,51式变成（6.55）（6.56）共五十七页仍将代入6.56式，可得（6.57）（6.58）上两式是s平面与z平面之间的单值映射关系，这种变换(binhun)就称为双线性变换(binhun)。共五十七页二、变换(binhun)常数c的选择不同c可以调节(tioji)频带间的对应关系。有两种方法（1）使两种滤波器在低频处有较明确的对应关系。即低频时，当较小有由6.55式及可得因而得到（6.

16、59）此时两者低频特性近似。共五十七页（2）采用数字滤波器的某一特定频率与模拟原型滤波器的一个特定频率相对应，即则有（6.60）这种方法的主要优点是在特定的模拟和数字频率处，频率响应严格(yng)相等。可准确的控制截止频率的位置。共五十七页三、逼近(bjn)的情况6.57,58就满足要求：（1）首先把代入6.41式，可得（6.61）即s平面的确实(qush)与z平面的单位圆相对应。（2）其次将代入6.42式中，得，因此从上看出，s左半平面对应z平面的单位圆内，s右半对应z单位圆外，虚轴对应单位圆上。共五十七页四、优缺点由6.61式重写如下（6.62）其关系表见6.12所示。从12所示，当零频附

17、近呈线性关系，当增加时，就存在非线性关系。其产生(chnshng)两个问题：（1）经双线性变换不能保持原有相位。（2）这种非线性关系要求模拟滤波器的帧频响应必须是分段常数型的。注：图6.13是理想数字微分器变换后产生的畸变。共五十七页解决办法：将频率加以预畸，即利用将数字频率变换(binhun)成模拟频率，利用这组模拟频率来设计模拟带通滤波器，这是我们要求的模拟原型。这一预畸过程如图6.14。共五十七页五、模拟(mn)滤波器数字化方法利用双线性变换将模拟滤波器数字化的办法：（1）用6.57式将模拟系统函数(hnsh)变换为数字滤波器函数(hnsh)：（6.63）（2）先将模拟系统分解为并联子系统，然后再对子系统进行双线性变换。如子系统（6.64）变换后得离散系统函数（6.65）共五十七页其中(qzhng)：如果模拟系统分解为并联子系统，得（6.67）经双线性变换后，离散系统函数为：（6.68）其中（6.69）共五十七页内容摘要数字信号处理。任一个线性时不变系统(xtng)，其系统(xtng)函数为。当h(n)是实函数，因而