3差分方程Z变换解读

1、第 3 章 线性离散时间系统的描述及分析3.1 差分方程及其时域分析3.1.1 差分方程3.1.2 差分方程的解A 递推解B 古典解C Z 变换求解3.2 Z 变换3.2.1 Z 变换的定义3.2.2 Z 变换的性质3.2.3 Z 反变换A 长除法B 留数法C 部分分式法3.3 离散时间系统的 Z 域分析3.3.1 零输入响应3.3.2 零状态响应3.3.3 完全响应3.4 Z 传递函数及其求法3.4.1 Z 传递函数的定义3.4.2 离散系统的运算343由G(s)求G(z)连续时间系统的离散化A对G(s)的讨论B 对离散化方法的评价C 留数法D 直接代换法E系统等效法I 冲击响应不变法；F系

2、统等效法II阶跃响应不变法G 部分分式法3.4.4 离散化方法小结3.5 线性离散时间系统的稳定性分析3.5.1 闭环极点与输出特性之间的关系3.5.2 稳定判据3.6 线性离散时间系统的频率特性分析法3.6.1 线性离散时间系统的频率特性线性离散时间系统的频率特性分析法第3章 线性离散系统的描述及分析3.1差分方程及其时域分析3.1.1差分方程在线性离散时间动态系统中，输入激励序列u(k)与输出响应序列y(k)之间的动态关系在时域中用差分方程来描述，差分方程一般写成升序方式y(k+ n) +a(k + n-1)+| + a_iy(k+ 1)+ any(k(2 1)二 b°u(k m

3、) Qu(k m-1) | bm_iu(k 1) bmu(k)有始性：k - 0初始条件：y(0) = y。, y(1) = %，y(n-1) = y“-i 时间因果律：m岂n或写成mny(k n)二為 biu(k m - i) -' ajy(k n - j)i =0j =1上式表明某一离散时间点上输出值可能与当前时间点上的输入值(当以及此前若干个输入和输出值有关推论开来，当前的输出值是 此前”全部激励和内部状态共同作用的积累 效应。考虑实时控制系统的时间因果律，必须有 m命。当m=n时，表明当前时刻的输入会直接影响当前时刻的输出，可称为 直 传”当m<n时，表明当前时刻的输入不

4、会直接影响当前时刻的输出；当前时刻的输入对输出的影响会延时 n- m”拍。差分方程也可以写成降序方式式(2.1)中各项序号均减ny(k) a(k -1) a?y(k - 2)|- n 1) a(k - n)二(22)二 b°u(k) be(k -1) III bm_e(k - m 1) bmU(k - m)在降序方式中的n和m与升序方式中的n和m的含义不完全相同，因而对 n和m并无限制。在降序方式中，当bo工0时，相当于升序方式中 m=n的情况。此时 当前时 刻的响应与当前时刻的输入有关”。升序意味着超前，与连续时间系统中的微分相对应；当用 Z变换法求解差 分方程时，升序方式便于考虑

5、初始条件。降序意味着滞后，与连续时间系统中的积分相对应；当用 Z变换法求解差 分方程时，降序方式无法考虑初始条件。3.1.2差分方程的解5 1例：已知差分方程 x(k 2) x(k 1) x(k)二 r(k+1)+0.5r(k),其 中6 6r(k)=1, k>0 x(0)=1, x(1)=2试由迭代法求其全解的前5项；分别由古典法求其零输入解yzi(k)、零状态解yzs(k)，以及全解y(k)。给定一个差分方程，根据特定的输入时间序列u(k)和初始条件，来求得其输出序列y(k), 一般有三种方法。A. 递推解(迭代解)对式(2.1)差分方程可以写成mny(k n)八 bjU(k m-i

6、)-' ajy(k n-j)i =0j =1显然给定初始条件后，就可依次求出各点值。但是，式(2.1)差分方程中的n个初始条件x(0),x(10),x(n-1)仅仅是指“零 输入初始条件”，进行递推求解时的初始条件应该是“全解初始条件”；因而应 该先求出其“零状态初始条件”，“全解初始条件”是“零输入初始条件”与“零 状态初始条件”之和。上例已知零状态初始条件，由此可递推求得零输入解yzi(k)；可求零输入初始条件，由此可递推求得零状态解yzs(k)；以上初始条件之和为全解初始条件，由此递推即可直接求得全解 y(k)=yzi(k)+yzs(k)oB. 古典解法1) 零输入解在式(2.1

7、)中令输入为零，即u(k)=0, k为，则得齐次方程y(k n) a1 y(k n -1)any(k 1) any(k) = 0(2.3)类似于在解线性常微分方程时定义的微分算子P,对差分方程定义一个移序(增序)算子d,即dny(k)二 y(k n)n(2.4)dj(k)二 y(k-n)于是式(2.3)可以表示成(dn+ adnT 寺an)y(k)二 A(d)y(k)二 0以多项式A(d)存在n个单根为例，即nA d m (d - di) , di = 0, i = 1,2,., ni=1则有零输入解yzi(k)的“通解”式为n(2.5)y/k)二 CQk C2d2kCndnk 八 Gd, ,

8、 k0i丄其中Ci, C2,Cn是由n个(另输入)初始条件决定的n个待定常数。设给定初始条件为y(i)=yi, i=0, 1，n-1,分别代人上式可得yo1%d1d12.一ynd_d1n_1d2d22n -1d21C1dnC2dn2 C3n I I 3 I " I dnn_1Cn(2.6)可简记为矩阵方式Y)二 D*C以n个单根为例，矩阵D一定可逆。于是可得待定常数为C 二 D丫。当A(d)存在重根时，亦可得相应结果，不再赘述上例求得零输入解yzi (k)。2) 零状态解当“零输入初始状态”为零时，为求得式(2.1)在任意输入u(k)激励下的“零状态响应” yzs(k)，首先考虑单位

9、脉冲激励u(k)= (k)的特殊情况，此时的系统响应为单位脉冲响应，记为h(k),式(2.1)成为h(k n) a1h(k n T)an_(k 1) anh(k)二 b0 (k m) ty (k m T)b< (k 1) bm (k)可写成如下形式mnh(k + n) =bi§(k + m_i)_E aj h(k + n _ j), men(2.7)i=0j =0上式中依次令k=-n，-n +1,，-2, -1, 0,可求得前面n+1个点的结果,h(0) = ho人(1)=几，当 m<n 时，h(0)=h°=0h(n -1)= hn_!h(n)二 hn当k>

10、;0时，在式(2.7)中恒有k+m-i>0,即恒有、(k+m-i)=0,此时式(2.7)又成为一个齐次方程，等价为(2.8)h(k n) a1h(k n -1) an_(k 1) anh(k) = 0 人北何2)十2，h(n) = h n上式按差分方程的零输入解法求解，并考虑h(0)=0,即可得到式(2.1)的单位 脉冲响应序列h(k), k初。对于一个一般的输入序列u(k)= u(0), u(1), u(2),，可以写成cOu(k)二' u(i) (k - i)二 u(0) (k)u(1) (k -1)i=0按照线性系统的迭加原理，(k-1)所激励的响应为h(k-i)1(k-i

11、), i=0, 1,于是可得u(k)激励下的响应为y(k)二 u(0)h(k)1(k)u(1)h(k - 1)1(k T)u(k)h(0)kk(2.9)二 6 u(i)h(k _ i)二嘉 u(k - i)h(i)i =0i =0= u(k) h(k) , k- 0 称为u(k)和h(k)的卷和显然，卷和的定义与连续时间函数的卷积具有类似的形式。卷和计算例上例求得零状态解yzsi(k)。3) 全解1)和2)二者之和。上例y(k)=yzi(k)+ yzs(k)。C. Z变换解法一一后面再讲3.2 Z变换3.2.1 Z变换的定义Z变换是对离散序列定义的，设有y(k)八 y(0),y(1),，二 y

12、(o)(k)y(1) (k-1)y(2) (k - 2), k-O则y(k)的z变换定义为(单边)罗朗级数i (2.10) 丫(z)二 y(0)y(1)z 八 y(i)zi=0zZ变换域变量d增序算子两者在数字上具有完全相同的表现形式，但意义却不同，不能混淆。就像s S变换域(拉氏变换)变量p微分算子二者表现形式相同，但意义截然不同为什么要定义Z变换？Z变换把离散(等距时间点上)数值序列变换成有理分式；L变换把连续时间信号变换成有理分式；便于利用代数学的某些结论进行简单处理。Z变换的另一种“定义”对于时域信号y(t)=f(t),采样得离散信号y*(t)记得第1章中讨论过y*(t) 和y*(k)

13、的(冲量的)等价性，f*(t)二 f(0) (t) f(T) (t-T) f(2T) (T-2T)二' f(kT) (t-kT)kdD取其拉氏变换，得0 _kTs(2.11)(2.12)F*(s)二 Lf *(t)八 f (kT)ek =0再令Tsz 二 e即得,QOF(z)二Zf*(t)八 f(kT)z*k=Q二者的结果是一致的。但是，二者有两点区别, 前者是对y(k)定义的，后者是对y*(t)定义的。在离散时间系统中使用前者更符合工程实际。但是，对于首先熟悉了Laplace变换的工程技术人员而言,后者更容易理解 前者在数学上是严格的；而后者中的式(2.11)容易使得误解z和s之间的

14、关系。实时上z和s之间并没有式(2.11)所示的关系，仅仅是有时同一个被控对象的Z变换传递函数和L变换传递函数的特征根具有那个关系3.2.2 Z变换的性质A. 在简单的情况下，可直接按定义求得y(k)的Z变换Y(z)。CO(k)丨八、(i)z=1i =0oOoOZ1(k) 1 八1(i)z-i 八 zi =0i =01-z1Z e"八 e'iTz_ii =0八 eTz)i =0zeT11-eTz(2.13)(2.14)(2.15)做为线性离散系统的Z变换，它有许多与L变换类似的性质，不同的是按 照Z变换的定义，这些性质更容易被证明一些。B. 线性迭加性质:已知 Z f'

15、;k) = F1(z),Z f2(k) = F?(z),a,b R，下同。按定义可得,Zafdk) bf2(k)二 Zat(k) Zbf2(k)二 aZf1(k) bZf2(k) = aF1(z) bF2(z)C. 增序性质：(对应于L变换的微分性质)设 g(k)=f(k+n), k>0为什么？0Zf(k n)=Zg(k)八 g(k)zk =0二' f(jn)zjf(j n)z£ n)znj £j=0(2.17)7n -An -A八 f (i)z'zn C f (i)znf (i)izn)i hii =0i =0con 二 zn、f(i)zT-、f(i

16、)zi -0i -0=znF(z) - zn f (0) - zn f (1)- z2 f (n - 2) - zf (n - 1)(令 i=j+n )注意两点:是为什么要减去前面几项？因为按照定义g(k)中没有这几项!二是与L变换的微分性质相比，形式上多了一个“ z”D. 减序性质：(对应于L变换的积分性质)设 g(k)=f(k-n), k>0为什么?QOQOZf (k -n)八 f(i- n)z八 zf(i-n)z2i=0i 曲(2.18)=f(j)zf(j)zj(令 i -n =j)jj =0二 zF(z)为什么第一项没啦？因为按照定义f(k)中的这几项为零!E.卷和性质：(对应于

17、L变换的卷积性质)(2.19)F.初值性质:f(0) =1忸 f(k) pm：F (z)(2.20)Zfi(h)* f2(h)二 Fi(z)F2(z)证明：一一按照Z变换的定义。G.终值性质:(2.21)f L) Jim. f (k)=典(1 - zF(z) =1阿(z-1)F(z)当f(k)不收敛(F(z)中有单位圆外极点)时，终值性质不能使用!证明:'Z【f(k+1)-f (k)】 = zF(z) f (0) - F(z)= (z-1)F(z) zf(0)aZf(k+1)-f (k) I 八f (i+1)-f (i)zi =0同令ZT1得，I叫 z-1)F (Z)z j=f (0)

18、+ f (1H (0) f (2)-f .f (k)-f (k 一 1)二 f(：)其它略3.2.3 Z反变换(2.22)已知F(z)有理分式，求f(k)使得Z f(k)H F(z),记为 f(k)二 ZpF(z)A. 长除法一一罗朗级数展开如果F(z)是有理分式，必可展开为罗朗级数，,即有如果F(z)是真有理分式，必可展开为(单边)罗朗级数(有始函数)f(k), k>0如果F(z)是严格真有理分式，则一定有f(0)=0。例，B. 留数法在实时离散控制系统中有f(k), kQ则一定有Q0F二 f (0)zf z"+八 f(k)z*k=0按照复变函数的留数理论，考虑如下围线(逆时

19、针包围含全部极点)积分,CF(z)zkdz 二 -f(i)zzkdzCC y=f(0) f (1)zJ fCk-Jz"f(k)zf(k1)1+.z1dzC二cf(0)zkJ f(1)zkf (k -1)f (k)z-1f (k1)z'+.dz二 cf (k)zdz = 2二 jf (k)留数是如何定义的?f(k)、j1k 1cF(z)zdz称为F(z)zk的留数于是有n(2.23)f (k)二 Z-1F(z)二 ResF(z)zk八 ResF(z)zk即f(k)为T(z)zk在其所有极点Zi, i = 1, 2,，n,处的留数之和。按照留数计算规则，若zc是F(z)的单重极点

20、则有Res F(z)zk二二 lim( z _ zo)F(z)zkJzo 一zr.zo若zo是F(z)的m重极点,则有m X輕F(z)zk耐処萨(灯厂k -1C. 部分分式法一一留数法的特例一一一般都是直接查表设F(z)有n部分分式法是应用留数法得到的一些易于实际应用的特例情况,个单重根乙，zn，则可以写成部分分式形式n(2.24)F(z)=5： Ai # Z 一 Zi按照迭加原理，我们可以求得其中每一项的Z反变换，即f(k)= ZF(z)八 ZAi JzZ-ZiAz/按式(2.23)有,nnk 1Z k 1f(k)=E ResF(z)z = E A Res(z )i Aid z Z 乙(2.

21、25)z(z- Zi) z)z Zi正是所希望的结果3.3离散时间系统的Z域分析利用Z变换求解差分方程。3.3.1零输入响应对式(2.1)所示差分方程，当输入u(k)=O, k0时，成为齐次方程yzi(k n) yyzi(k n -1)anyzi(k 1) anyzi(k) = 0y(o)=y。, y(1)=y1,，y(n-1)=yn-1应用Z变换的增序性质，并注意给定的零输入初始条件，得整理可得znYzi(z)-zny° - zny1 -z%-+ a1znYzi (z) - zy。Dzyz+ anzY；i (z) - zy° + anYzi (z) = 0Y (z)B(y

22、o,y1,., y2,y1)z + a1z+ . + an /+ an于是可得式(2.1)的零输入响应为yzi(k)=Z1YZi(z)3.3.2零状态响应设式(2.1)所示系统在没有输入激励时，其内部初始能量积累为零，即所谓零状态，此时不考虑初始条件对式2.1的两边同时进行Z变换，可得Yzs(z)二b°zmb1Zmbm_1Z bmnn-1za1z . an/Z anU(z)定义G(z哼吩" bm z + a1z + + an十 an(2.26)称为离散动态系统式(2.1)的Z传递函数，则上式可写成Yzs(z)弋(z)U(z)则有yzs(k Z-1Yzs(z) = Z-1G(

23、z)U(z)按照卷和定理kyzs(k) = g(k)* u(k)八 g(k-i)u(i), k-O其中g(k)=zTG(z)g(k)是什么，以及如何求得g(k)?设u(k)=迭)是一个单位脉冲函数，已知，U(z)=Z迭)=1，即可得系统对u(k)= <k)的零状态响应，称为单位脉 冲响应，并记为h(k), k>0并有h(kz，G(z)g(k)现在，如欲解析求解式(2.1)所示的差分方程的零状态响应，主要有两种 方法。Z 域法：yzs(k) = ZG(z)U(z)时域法：yzs(k)二 h(k)* u(k)333完全响应对式(2.1)求Z变换时，同时考虑初始条件，即可得系统的完全响应

24、，与 分别求出yzi(k)和yzs(k)再相加是一致的。即:Y(z) =B(yo,y,yVni)zn af an_Z anG(z)U (z)=Yzi(z) Yzs(z) y(k) =yzi(k)yzs(k)(2.27)几点说明:在求零状态响应时，显然零状态解 yzs(k)的初始n个值并不一定为零，零状 态仅仅是说当输入为零时，系统初值为零。求零状态响应时，对式(2.1)两边求Z变换时，此时的yzs(k)与u(k)都是有初 值的，因此亦应考虑增序性质时的初值，但是在整理时两边的初值正好相互抵 消，因此在求零状态响应时的 Z变换时，可以不考虑初值。在求完全响应时，由u(k)引起的yzs(k)中的那

25、一部分初值效应必然由u(k)的初值效应所抵消，因此只考虑系统的零输入初值。5 1例：已知差分方程 x(k 2) x(k 1) x(k) = r(k+1)+0.5r(k)，其 中6 6r(k)=1, k>0 x(0)=1, x(1)=2。试由Z变换法求其全解。3.4 Z传递函数及其求法3.4.1 Z传递函数的定义定义一个离散时间被控对象的动态特性，或连续时间对象的离散控制动态特性。由输入-输出序列Z变换之比来定义。 传递函数描述一个动态系统的输入输出稳态传递特性(稳态的含义是 不包含初始条件的影响)。A对于离散时间系统u(k)y(k)离散时间系统G(z)U(z)G(z) =Y(z)U(z)

26、Y(z)图2.1离散时间被控对象传递函数比如这个离散时间系统原来是由差分方程描述的。对于式(2.1)描述的差分方程,(2.1)y(k n) a$(k n -1)any(k 1) 3ny(k)二，y(1)=,y(n-1)= %根据Z变换的性质，两边求Z变换(不考虑初始条件)，并化简可得m _1bz babmnn -1za1z a* an(2.28)如果差分方程是由式2.2描述的,y(k) a(k -1) a2y(k - 2)azy(k - n 1) any(k-二 b0u(k) b1u(k -1)bm1u(k - m 1) bmu(k - m)则同理可得G二n -1n-m<)n-mz bm

27、_1ZbmZzn 事' . an? an(0Zmbm 1Zbm)an-1zanY(z)bozn bU(z)n -mznn =za1z(2.2)(2.29)当n= m时,与式(2.28)相同=b0u(k m) b1u(k m -1)bm-1u(k 1) bmu(k) k = 0,1,2,，m 乞 n y(0) = y°,注意:2) 为什么上二式求Z变换时不考虑初始条件？传递函数只描述稳态特性，与初始条件无关！3) 式(2.28)和 (2.29)称为有理分式；n<m时称为(假)有理分式，反时间因果律，离散时间系统中不存在;n= m时称为真有理分式，输入-输出有直通分量；n&

28、gt;m时称为严格真有理分式，输入-输出至少延时一拍B对于一个连续时间的采样控制系统对于一个连续时间系统，对其进行离散时间控制时前面必须加一个零阶保 持器(ZOH)。只有对其输入和输出采样得到响应的输入-输出离散时间序列时, 才能对其定义Z传递函数。G(s)二 丫(s)/ U Uu(k)丫离散时间系统G(z)y(k)G(z)二丫U(z)图2采样控制的连续时间系统的离散时间传递函数3.4.2离散系统的运算流图化简，与连续时间系统完全相同。A串联U(z)* Gi(z) G2f Gn(Z)Y(z)B并联4 G(z)> Y(z)nG (z)二 ii Gi (z)i =1图3 离散时间系统的串联C

29、反馈系统鸣G(z)七nG(z)八 Gi(z)i =1图4 离散时间系统的并联U(z)G(z)丫(z)G(z)1 Gdz)G2(z)G,z)图5 离散时间反馈系统对于任意的复杂系统，可由梅森公式求得343 由G(s)求G(z)连续时间系统(或信号)的离散化A 对 G(s) 的讨论一般来说，G(s)的含义可能有以下三种情况：1) G(s)为时域信号g(t)的Laplace变换此时，应该由G(s)求的g(t)，对g(t)离散化得g(k)，最后再求G(z)。2) G(s)为控制器的传递函数一一它只是一个数字模型G(s)既可以由连续时间系统(模拟)实现，输入输出为连续时间变量；G(s)也可以由离散时间系

30、统(数字)实现、输入输出为离散时间变量； 此时，对G(s)直接离散化即可，不需要 ZOH。3) G(s)是一个(连续时间)被控对象离散化后的输入时离散时间的，但是 G(s)只能接受连续时间激励信号， 因此必须在输入端需增设一个保持器(例如零阶保持器ZOH)，将离散序列转化为连续时间函数。G(s)的输出一定是连续时间函数，需对其进行采样y(k)> Gi 匚)G(»zG()J对连续时间被控对象的离散化B 对离散化方法的评价离散化方法不是唯一的，它们各有其特点和适用范围。因而需要对离散化 方法建立评价指标体系。对信号的离散化结果应该是唯一的，严格的。就是说在采样点上的取值严 格等于原

31、函数。对调节器传递函数G(s)的离散化结果G(z),应与G(s)的频率特相一致。这 时会因所用方法的不同而有差异。对被控对象传递函数G(s)的离散化结果G(z),在不同情况下有不同的要求， 后面会详细讨论。这时也会因方法的不同而有差异。评价一个离散化方法，大概有如下5项指标。但是在不同的应用场合有不同的要求。1) 易操作性。2) 从S平面到Z平面的映射关系。包括映射的单值性和 稳定性的遗传 性。3) 频率特性畸变。指G(z)的频率特性与G(s)的频率特的一致性。4) 稳态增益畸变。指G(z)的稳态增益与G(s)的稳态增益的一致性。5) 时域(采样点)响应的一致性。指在采样点上G(z)和G(s)

32、取值的一致性。C 留数法适用于G(s)为时域信号g(t)的Laplace变换的情况。这时，G(z)和G(s)在采样点上的取值是完全一致的。=F g(t)k=0G(z) =' g(k)z*kRkzt=kT-jskTk一产 G(s)e ds z按定义带入g(t)1- jskT _k、云j “G(sh(e z )ds交换求和求积分的顺序级数和的闭式按留数定理即可得,m1G( Z = ResG(s) r-r30)匸1 i1-ezD 直接代换法操作简单，但却有误差。直接代换法既适用于对控制器的离散化，亦适用于对被控对象的离散化。但是不适用于对信号的离散化(在采样点上取值不严格)。使用直接代换法对

33、被控对象离散化时，一方面物理上需要引入ZOH，两一方面代换是并不包括ZOH。直接代换法有很多种，下面介绍常用的几种。1)后向差分法设连续时间描述为：ddtG(s)二X (s)U (s)用差分代替微分，采样周期取为 T,x(k 1)- x(k)Tu(k 1),G(z)X(z)UTzz-1(为什么叫“后向”差分？ ？)1z 二1 - sT(2.31)比较G(s)和 G(z)，可得代换式,zT s 二TzS-； z映射关系：单值一一对应S平面上左半平面稳定域二= Z平面上单位圆内正实轴上小圆G(s)稳定=G(z)稳定cc：)S平面ZZac0b7z平面-1cc：)图7 后向差分法的稳定性遗传显然稳定性

34、的遗传不是可逆的，但 S稳定” “z稳定”因此常被采用（S平面上除了 aef小圆外，所有的s映射到Z平面都是稳定的）频轴畸变较大。稳态增益无畸变。即:G(s)»訂 G(z)|z.cc：)cc：)不能保证时域（采样点）响应的一致性。2）前向差分法连续时间系统描述为X = U, dtG( As用差分代替微分x(k 1-)x(ku( k),cc：)cc：)（为什么叫“后向”差分？ ？）比较G（s）和 G（z），可得代换式,z-1(2.32)s到z映射关系：单值一一对应。事实上就是一个平移图8前向差分法的稳定性遗传G(s)稳定 G稳定显然，G(s)稳定很难保证G是稳定的，固很少采用 频轴畸变

35、较大。稳态增益无畸变，即：g(s)| 二g (z)|不能保证时域(采样点)响应的一致性。3) 双线性变换法(Tustin法)连续时间系统描述为G(s)二sdx = u, dt用差分代替微分，x(k 1)- x(k) u(k 1) u(k)T2比较的代换式，G(zT z 12 z- 1(2.33)(为什么叫“双线性变换？ ？)图9双线性变换法的稳定性遗传s到z的映射关系：单值对应；S平面上左半平面稳定域二-Z平面上单位圆内稳定域G(s)稳定 =G(z)稳定当T足够小时(即当 匕足够大时)频轴畸变很小；稳态增益无畸变；事实上在程序化处理的G(s)显然，在直接代换法中，双线性变换是最好的 到G(z)

36、变换中都采用双线性变换法，应用最为广泛。E 系统等效法I冲激响应不变法提法：设有(被控对象)G(s)和G(z),若G(s)在 吐)的激励下的响应g(t)在 kT处的采样值g(kT)与G(z)在S (k)的激励下所得之响应相等，即称 G(z)和G(s) 是冲击响应不变(等价)的。但是，事实上 迭)和 我)并不等价。原因是，吐)的冲量为1,而(加上零阶 保持器之后)Sk)的冲量为“ T”二者差一个系数“ T”使得G(z)的稳态增益随 着T大幅变化，这是不允许的。为什么还要讲这种方法？按定义，在St)激励下，有冲激响应g(t)1g(t)二 L'G(s)-：G(s)estds2兀j按采样周期T

37、采样即得g(k)二 g(kT )-j： G(s)eskTds2花j按照输入输出等效原则，在单位脉冲输入Sk)的激励下，应有输出g(k)如上 式所示。根据Z变换的定义，即有对上式求Z变换G(z)= ' g(k)z"k=0乂1<j+j 閃skT_k交换和积顺序求级数和的闭式八：G(s)e dszkk=0 2 二 j1；j "skT -k、.j：G(s)' (e z )ds2 二 jk =0壷二Gwds按留数定理G( z)ResG(s) i=1 i1-sT -1e z(2.34)因此，冲击响应等效法也是留数计算法。显然，此式与式(2.30)的留数法相同。此式

38、用来对信号的G(s)求其G(z)时是严格正确的，但是，用来对被控对象的G(s)求其G(z)时却是不对的。此代换不易操作，特别是不易计算机实现。S到z的映射关系分析如下。若G(s)有一个极点gS =则G(z)定有一个极点2兀iT j( J -门.)Te ij diT其中 r = eiT 1,di(T显然，S平面=z平面,单值映射Z平面=S平面,多值映射主域图10冲击响应等效法的稳定性遗传3T 31 灼 1如果只考虑S平面的主值域，即i ' T"T ，则有一一对应的关系。在主值域内有di 一 '，因此，频轴无畸变。求式(2.34)的稳态增益lim G (z) Lt =旳T

39、 ： 0z -1可见G(z)的稳态增益受采样周期T响应很大。因此，稳态增益畸变严重一使得本法很少使用。当T足够小时，一定可使所有 S域极点均落在主域之内，此时的映射可相当于一一对应的。主域一一整个Z平面；左半平面一一单位圆内；右半平面一一单位圆外；虚轴单位圆；容易理解，如果在§(k)的激励下也引入零阶保持器时，&k)和8(t)就成为等 价的了(为什么？)，于是式(2.34)成为,m1- e"Ts1(2.35)G( z)=' ResG(s) st ti=1 i s1- e z.m 1 1=(1 - z_ 比 Res G( s)srTr1 i s 1- e z由

40、下式可以证明稳态增益无畸变G(z) |z=_i、J1 1-(1 -G(s)s+ |(V z庄 R丿(1 - z )迟 Res-G(s)、i s1 - Z1 - z二 G(s)|esTzzG(s)(2.36)s Z0f系统等效法n阶跃响应不变法提法：设有G(s)和G(z),若G(s)在1(t)的激励下的响应e(t)在kT处的采样值 e(kT)与G(z)在1(k)的激励下所得之响应相等，即称G(z)和G(s)是阶跃响应不变(等价)的。在阶跃输入的特殊情况下，在1(t)的后面有无零阶保持器是无区别的(?)。有.1 1.，11Z G(z)=L皿)两边求z变换，得G匚1 1ZL G(s) I sRes-

41、G(s)- Si s 1 - e z可得,G(z)= (1 -z1)'i =11 .Re于(s)E(2-37)sTS到z的映射关系与冲激响应不变法相同;从变换关系式可知，无频轴畸变。由下式可知，无增益畸变m11、G1C(1-z花 Res1G(s)PL-11 - z1_zG(s)+s=0(1- Z)' Ress=0G(s) 1s(VesTr1h1二 G(s)s -0对比式(2.36)和式(2.37)可知，引入零阶保持器时的冲激响应等效法 式(2.36) 与不引入零阶保持器时的阶跃响应等效法式(2.37)二者是等价的。G 部分分式法事实上，部分分式法是留数计算法的一个变形，也是留数

42、法的一种使用形 式。一般教科书中都给出相应的表格以供查照。3.4.4离散化方法小结1) 对于表示信号的G(s)的离散化必须直接使用留数法(部分分式法)2) 在物理上，表示调节器的G(s)不需要ZOH，表示被控对象的G(s)必需要加ZOH3) 无论对于表示调节器还是表示被控对象的 G(s)的离散化，都可以使用直 接代换法，也可以使用留数法(部分分式法)。但是在数学上，使用直 接代换法时不需要ZOH，使用留数法(部分分式法)时需要先加上ZOH3.5线性离散时间系统的稳定性分析3.5.1离散系统的闭环极点(特征值)与系统输出特性的关系设线性离散时间系统G(z)，G(z)二= k n Am(z)Bn(

43、z)n (z- pji 二其中Am(z)为m阶首一多项式，并设pi为单实根或单共轭复根的情况，且设G(z)中没有z=1的极点，即有pi工1当存在复根或z=1的极点时，如下各项分析 结论仍然成立。当存在一对共轭复根时，有j-j 打Pi 二 re , Pi i 二 re当输入为单位阶跃序列即魁)z ，此时输出为Z-1Y (z)二 G (z) R(z)二 kAm(z)nH (z - Pi)i 2由上一节讨论可知，求上式的Z反变换，可得y(k)=人 ' kr(pjk ' ks(e%kejks /%隹八)prps二 kkr(p)、ksrsTerprps.e（k s飞）二心、K(Pr)k

44、' 2ksskcos(k spprs上式中，ko为与阶跃输入相对应的稳态响应项Pr为单重实根极点，kr为与Pr相对应的输出项系数ps为单重共轭复极点，其中rs为其幅值，申s为其幅角ks为与极点Ps相对应的输出项的系数幅值，二为其相位角s由上式可知，如果|>1，贝U随着kT=O ，Pi的对应输出项发散，不稳定如果|p|=1，贝吓随着 匕珀,Pi的对应输出项为恒值（实根）或等幅振荡 （共轭复根），临界稳定。如果Opp|<1，贝U随着kToO , Pi的对应项收敛，稳定。再考察共轭复根对应输出项的相角特性（周期振荡），令k$ = 2二，则一个振荡周期对应的周期数为kd2, o.:

45、(考虑共轭复数)s9s显然，質越接近零，kd越大，即振荡周期越长，当兀时，kdT 2，输出正负交替。震荡周期为两个采样周期。图例：Pi对应输出，发散， ® s = 22 5 ° kd = 16 P8对应输出，稳定，申=15kd=24稳定，/s= 180 ,kd = 2稳定，.:s=0 K八临界稳定，s =135 ,kd =2.7临界稳定，s = 45 , kd = 8临界稳定，90 *厂4P6对应输出,P7对应输出,R对应输出,P2对应输出,P3对应输出,试分别画出与上述各特征根对应的输出模态波形的示意图。由上分析可知：线性离散系统的稳定域是其传递函数的全部特征根均落在Z平面上的单位圆之内3.5.2线性离散时间系统的稳定判据不需要解出特征方程的根，只需根据特征方程的系数就可判定特征根是否 全部位于单位圆内。1.舒