自动控制原理(下)
上QQ阅读APP看书,第一时间看更新

8.3z变换

在线性连续控制系统中,以拉普拉斯变换作为数学工具,将系统的微分方程化为代数方程,建立系统的传递函数,可以非常方便地对系统进行分析和设计。与此相似,在线性离散控制系统中,应用z变换将描述系统的线性差分方程转化为代数方程,建立系统的脉冲传递函数,从而对系统进行分析和设计。因此,z变换是研究线性离散控制系统的重要数学工具。

8.3.1z变换的定义

设连续信号xt)存在拉普拉斯变换,其象函数为Xs)。xt)经过等速采样后,得到离散信号xt),由式(8-2)可知其表达式为

978-7-111-66176-4-Part01-40.jpg

式中,T为采样周期。对上式表示的脉冲序列进行拉普拉斯变换,可得采样函数

978-7-111-66176-4-Part01-41.jpg

式(8-13)最后一个等式中每一项均含有复变量s的指数函数eTs,直接运算不方便。若引入复变量

978-7-111-66176-4-Part01-42.jpg

则式(8-13)可写成以z为自变量的函数

978-7-111-66176-4-Part01-43.jpg

式(8-15)为离散信号xt)的z变换的定义式。Z[xt)]表示离散信号xt)的z变换,记为Xz)。通过引入复变量z,将采样函数拉普拉斯变换的表达式转换为z的幂级数形式。

z变换,又称为采样拉普拉斯变换,它是从拉普拉斯变换直接引申出来的变换方法,它实际上是采样函数拉普拉斯变换的一种变形。对一连续函数xt)取z变换,只考虑这个函数在采样时刻的采样值,即时间序列x(0),xT),x(2T),…。也就是说

Xz)=Z[xt)]=Z[xt)]=Z[xkT)]=Z[Xs)]=Z[Xs)]

都表示对离散信号xt)的z变换。

8.3.2z变换的方法

常用的z变换方法有级数求和法、部分分式法和留数计算法。

1.级数求和法

根据z变换的定义,将式(8-15)写成级数展开式

978-7-111-66176-4-Part01-44.jpg

显然,只要知道连续函数xt)在各采样时刻t=kTk=0,1,2,…,∞)的采样值xkT),便可求得离散函数xt)的z变换的级数展开形式。这种级数展开式是开放式的,如果不能写成闭合式形式,实际应用就不方便。一些常用函数z变换的级数形式,都可以写成闭合式形式。

例8-1 求单位阶跃函数1(t)的z变换。

解:单位阶跃函数在各采样时刻的采样值均为1,即1(kT)=1 (k=0,1,2,…,∞),由式(8-16)可得

978-7-111-66176-4-Part01-45.jpg

若|z|>1,则上式的无穷级数是收敛的,利用等比级数求和公式可得对应的闭合式形式为

978-7-111-66176-4-Part01-46.jpg

例8-2 求指数函数e-ata>0,t≥0)的z变换。

解:根据式(8-16)可得

978-7-111-66176-4-Part01-47.jpg

若|z|>a,则无穷级数是收敛的,利用等比级数求和公式可得对应的闭合式形式为

978-7-111-66176-4-Part01-48.jpg

2.部分分式法

利用部分分式法求z变换时,先求出已知连续函数xt)的拉普拉斯变换Xs)。Xs)通常是s的有理分式,然后将其展开成部分分式之和的形式,使每一部分分式对应简单的时间函数,再分别求出或查表得到每一项的z变换。最后作通分化简运算,于是可方便地求出Xs)对应的z变换Xz)。

Xs)展开成部分分式和的形式,即

978-7-111-66176-4-Part01-49.jpg

式中,nXs)的极点个数;Ai为常系数;siXs)的极点。由拉普拉斯反变换可知Ai/(s-si)对应的原函数为978-7-111-66176-4-Part01-50.jpg,利用在例8-2中求得的指数函数的z变换,可得

978-7-111-66176-4-Part01-51.jpg

因此,函数xt)的z变换由象函数Xs)求得为

978-7-111-66176-4-Part01-52.jpg

例8-3 已知连续函数xt)的拉普拉斯变换978-7-111-66176-4-Part01-53.jpg,求对应的z变换Xz)。

解:将Xs)展开成如下部分分式

978-7-111-66176-4-Part01-54.jpg

于是,与式(8-17)对照可知,A1=1,A2=-1,s1=0,s2=-a。根据式(8-18),可得连续函数xt)的z变换

978-7-111-66176-4-Part01-55.jpg

例8-4 已知连续函数xt)=sinωt,求对应的z变换Xz)。

解:对xt)取拉普拉斯变换,得

978-7-111-66176-4-Part01-56.jpg

Xs)展开成部分分式

978-7-111-66176-4-Part01-57.jpg

根据式(8-18),可得

978-7-111-66176-4-Part01-58.jpg

3.留数计算法

已知连续函数xt)的拉普拉斯变换Xs)及其全部极点sii=1,2,3,…,n),则其z变换可通过下列留数计算式求得,即

978-7-111-66176-4-Part01-59.jpg

式中,siXs)的彼此不相等的极点;ri为重极点si的阶数;n为彼此不相等的极点个数。

例8-5 已知连续函数xt)=t2,求对应的z变换Xz)。

解:xt)的拉普拉斯变换为

978-7-111-66176-4-Part01-60.jpg

由上式可知,n=1,s1=0,r1=3。根据式(8-19)可得

978-7-111-66176-4-Part01-61.jpg

例8-6 已知拉普拉斯变换978-7-111-66176-4-Part01-62.jpg,求对应的。变换X(z)。

解:由拉普拉斯变换式Xs)可知,n=2,s1=-1,r1=2,s2=-2,r2=1。根据式(8-19)可得

978-7-111-66176-4-Part01-63.jpg

常用时间函数的z变换见表8-2。由表可知,这些函数的z变换都是z的有理分式,且分母多项式的次数大于或等于分子多项式的次数。值得指出,表中所列常用时间函数z变换式的分母z多项式的阶次与相应拉普拉斯变换式的分母s多项式的阶次相等。

表8-2z变换表

978-7-111-66176-4-Part01-64.jpg

(续)

978-7-111-66176-4-Part01-65.jpg

8.3.3z变换的基本定理

z变换中,有一些与拉普拉斯变换类似的基本定理。熟悉了这些定理,可以更加简便地应用z变换。

1.线性定理

X1z)=Z[x1t)],X2z)=Z[x2t)],a1a2为常数,则

978-7-111-66176-4-Part01-66.jpg

证明:由z变换定义

978-7-111-66176-4-Part01-67.jpg

线性定理表明,时域函数的线性组合的z变换等于各函数z变换的线性组合,z变换是一种线性变换。

2.实数位移定理

实数位移定理又称平移定理。实数位移是指整个采样序列在时间轴上左右平移若干采样周期,其中向左平移为超前、向右平移为滞后。实数位移定理如下:

若连续函数xt)的z变换为Xz),则有滞后定理

978-7-111-66176-4-Part01-68.jpg

以及超前定理

978-7-111-66176-4-Part01-69.jpg

证明:1)式(8-21)的证明。由z变换定义

978-7-111-66176-4-Part01-70.jpg

m=n-k,则有

978-7-111-66176-4-Part01-71.jpg

由于z变换的单边性,当m<0时,有xmT)<0,所以上式可写为

978-7-111-66176-4-Part01-72.jpg

式(8-21)得证。

2)式(8-22)的证明。由z变换定义

978-7-111-66176-4-Part01-73.jpg

m=n+k,则有

978-7-111-66176-4-Part01-74.jpg

式(8-22)得证。

算子z有明确的物理意义:z-k代表时域中的滞后算子,它将采样信号滞后k个采样周期;zk代表时域中的超前算子,它将采样信号超前k个采样周期。但是,zk仅用于运算,在实际物理系统中并不存在。实数位移定理是一个重要的定理,可将描述离散系统的差分方程转换为z域的代数方程。

3.复数位移定理

若连续函数xt)的z变换为Xz),a为常数,则有

978-7-111-66176-4-Part01-75.jpg

证明:由z变换定义

978-7-111-66176-4-Part01-76.jpg

z1=ze±aT,则有

978-7-111-66176-4-Part01-77.jpg

式(8-23)得证。

复数位移定理的含意是函数xt)乘以指数函数eatz变换,就等于在xt)的z变换式Xz)中以ze±aT取代原算子z

4.初值定理

若连续函数xt)的z变换为Xz),且极限978-7-111-66176-4-Part01-78.jpg存在,则有

978-7-111-66176-4-Part01-79.jpg

证明:由z变换定义

978-7-111-66176-4-Part01-80.jpg

z→∞时,等式右边只有x(0),所以

978-7-111-66176-4-Part01-81.jpg

5.终值定理

若连续函数xt)的z变换为Xz),Xz)不含有z=1的二重及其以上的极点且在z平面的单位圆外无极点,则有

978-7-111-66176-4-Part01-82.jpg

证明:由z变换的线性定理,有

978-7-111-66176-4-Part01-83.jpg

由实数位移定理

Z[xt+T)]=zXz)-zx(0)

于是

978-7-111-66176-4-Part01-84.jpg

上式两边取z→1的极限,得

978-7-111-66176-4-Part01-85.jpg

978-7-111-66176-4-Part01-86.jpg

所以有

978-7-111-66176-4-Part01-87.jpg

如果已知xt)的z变换Xz),不需要求其z反变换,利用初值定理和终值定理可以方便地求出xt)的初值和终值。

8.3.4z反变换

所谓z反变换,就是已知z变换表达式Xz),求取相应的离散函数xt)或离散时间序列xkT)的过程。记为

978-7-111-66176-4-Part01-88.jpg

通过z反变换只能求出连续信号在采样时刻的数值,而不能给出连续信号在采样时刻之间的有关信息。常用的z反变换方法有幂级数法、部分分式法和留数计算法。

1.幂级数法

幂级数法就是利用长除法将z变换表达式Xz)展开成z-1的幂级数,求取离散函数xt)或离散时间序列xkT)的数值。

z变换表达式Xz)是z的有理分式函数,即

978-7-111-66176-4-Part01-89.jpg

式中,aii=0,1,2,…,n)和bjj=0,1,2,…,m)均为常数。将Xz)的分子、分母多项式表示为z-1的升幂形式,则可以直接用分母去除分子,得到无穷幂级数的展开式

978-7-111-66176-4-Part01-90.jpg

根据实数位移定理,离散函数xt)为

978-7-111-66176-4-Part01-91.jpg

考虑到式(8-27),由z变换定义式(8-15)可知,式(8-28)中的系数ckk=0,1,2,…,∞)即为xt)在采样时刻t=kT的值xkT)。

在实际应用中,通常计算有限的几项就够了,因此用幂级数法得到xt)较简便,但不容易求出xt)的通项表达式。

例8-7 已知z变换表达式978-7-111-66176-4-Part01-92.jpg,试求其z反变换。

解:将Xz)表示为

978-7-111-66176-4-Part01-93.jpg

长除法算式如下

978-7-111-66176-4-Part01-94.jpg

利用长除法得

Xz)=1+3.5z-1+4.75z-2+6.375z-3+…

则离散时间函数为

xt)=δt)+3.5δt-T)+4.75δt-2T)+6.375δt-3T)+…

2.部分分式法

已知z变换表达式Xz),考虑到在z变换表中z变换表达式Xz)的分子都含有因子z,所以应将Xz/z展开为部分分式,逐项查z变换表,就可以得到离散时间函数xt)或离散时间序列xkT)。

z变换表达式Xz)无重极点,将Xz/z展开为如下部分分式

978-7-111-66176-4-Part01-95.jpg

式中,nXz)的极点个数;Ai为常系数且978-7-111-66176-4-Part01-96.jpgziXz)的极点。

式(8-29)的两端同乘以z,得到Xz)的部分分式展开式为

978-7-111-66176-4-Part01-97.jpg

逐项查z变换表,得到

978-7-111-66176-4-Part01-98.jpg

最后写出Xz)对应的采样函数

978-7-111-66176-4-Part01-99.jpg

例8-8 已知z变换表达式978-7-111-66176-4-Part01-100.jpg,试求其z反变换。

解:将Xz/z展开成如下部分分式

978-7-111-66176-4-Part01-101.jpg

上式两端同乘以z,得

978-7-111-66176-4-Part01-102.jpg

z变换表得离散时间序列

xkT)=10(2k-1) (k=0,1,2,3,…,∞)

于是,离散时间函数为

978-7-111-66176-4-Part01-103.jpg

3.留数计算法

在实际问题中遇到的z变换表达式Xz),除了有理分式外,也可能是超越函数,此时无法应用部分分式法及幂级数法来求z反变换,而采用留数计算法却比较方便。已知z变换表达式Xz),求取离散时间序列xkT)的留数计算式如下:

978-7-111-66176-4-Part01-104.jpg

式中,ziXz)的彼此不相等的极点;ri为重极点zi的阶数;n为彼此不相等的极点个数。则z变换函数Xz)对应的离散函数xt)为

978-7-111-66176-4-Part01-105.jpg

例8-9 已知z变换表达式978-7-111-66176-4-Part01-106.jpg,试求其z反变换。

解:由Xz)可知,n=2,z1=1,r1=1,z2=5,r2=2。根据式(8-32),可得

978-7-111-66176-4-Part01-107.jpg

于是,由式(8-33)得到Xz)对应的离散函数xt)为

978-7-111-66176-4-Part01-108.jpg