现代控制理论(第2版)
上QQ阅读APP看本书,新人免费读10天
设备和账号都新为新人

2.5 线性时变系统状态方程的解

线性时变系统的结构参数随时间而变化,其一般形式的状态方程为时变非齐次状态方程,即

式中,At)、Bt)分别为n×nn×r维时变实值矩阵。若输入控制u=0,式(2-47)则变为时变齐次状态方程,即

若矩阵A仅为一阶,即At)=at),则向量-矩阵时变齐次状态方程式(2-48)变为式(2-49)所示的标量时变齐次微分方程,即

可应用分离变量法求解式(2-49),即

式(2-50)表明,xt)也可视为初值xt0)的转移,但时变系统与定常系统转移特性不同之处在于,其不仅与系统特性at)及t有关,而且与初始时刻t0有关,但与tt0之差无关,即对于标量时变系统,其状态转移函数为

则标量时变齐次微分方程式(2-49)的解可表示为

仿照标量时变齐次微分方程解的表达式(2-52),时变齐次状态方程式(2-48)的解为

式中,Φtt0)为式(2-48)所描述时变系统状态转移矩阵。将式(2-53)代入式(2-48)得

由式(2-53)及式(2-54)可推知状态转移矩阵Φtt0)满足如下矩阵微分方程和初始条件

应用经典控制理论分析时变系统较为困难,而采用状态空间分析法的优点之一在于可将线性定常系统的求解方法推广到线性时变系统,且应用状态转移矩阵的概念和性质,可使时变系统的解在形式上与定常系统统一,即自由运动均可视为初始状态的转移。应该指出,时变系统状态转移矩阵用Φt,t0)表示,反映其为tt0的函数;但定常系统状态转移矩阵用Φt-t0)表示,反映其为t-t0的函数。

2.5.1 线性时变系统状态转移矩阵的求解

由式(2-11)可知,线性定常系统状态转移矩阵可用矩阵指数表示,即

但时变系统状态转移矩阵一般不能用矩阵指数给出,只有当At)与

满足矩阵相乘可交换条件,即成立时,Φtt0)才可用如下矩阵指数及其幂级数展开式表示,即

式(2-56)两边对t求导数,得

则由式(2-57)有

以上推导证明,若At)与

满足矩阵相乘可交换条件式(2-58),状态转移矩阵Φtt0)可用式(2-56)所示的矩阵指数表示,此时可得式(2-48)闭合形式的解为

下面进一步分析At)与

满足矩阵相乘可交换条件对时变系统状态矩阵At)的要求。由式()得

显然,若对任意的t1t2,下式

成立,则At)与

满足矩阵相乘可交换条件。

应该指出,时变系统的系统矩阵At)一般并不满足式(2-62),这时Φtt0)就不能采用简便方法求解,通常也得不到闭合形式的Φtt0),但可表示成递推形式,采用数值计算近似求解。由式(2-55)得

dΦtt0)=AtΦtt0)dt

t0t对上式两边取积分,得

反复应用式(2-63),可将Φtt0)展成无穷级数,即

式(2-64)所示级数称为Peano-Baker级数,若At)的元素在积分区间有界,则该级数收敛,但难以表示成封闭形式的解析式,可根据精度要求采用数值计算方法近似求解。

2.5.2 线性时变系统状态转移矩阵的性质

1.传递性

证明 由时变齐次状态方程解的表达式(2-53),有

故有

Φt2t0)=Φt2t1Φt1t0

2.可逆性

证明 由式(2-65)及式(2-55),则有

Φt0tΦtt0)=Φt0t0)=IΦtt0Φt0t)=Φtt)=I

故有

Φ-1tt0)=Φt0t

2.5.3 线性时变非齐次状态方程的解

设线性时变非齐次状态方程式(2-47)的解为

将式(2-67)代入式(2-47),并根据式(2-55)得

则有

上式中的ξt0)可据式(2-67)、式(2-55)求得,即

ξt0)=Φ-1t0t0xt0)=xt0

则线性时变非齐次状态方程式(2-47)的解为

式(2-68)表明,由于线性系统满足叠加原理,线性时变系统状态的全响应xt)由源于系统初始状态xt0)的零输入响应Φtt0xt0)和源于系统输入ut)控制作用的零状态响应两部分构成。应该指出,由于通常得不到闭合形式的Φtt0),故式(2-68)右边一般得不到闭合形式,需在数字计算机上根据精度要求采用数值计算方法近似计算。

【例2-7】 已知线性时变齐次状态方程为,求当t0=1

时状态方程的解。

At)与

满足矩阵相乘可交换条件,系统状态转移矩阵Φtt0)可由式(2-56)所示的。

【例2-8】 求线性时变系统t0=0的状态转移矩阵Φt,0)。

可见,At1At2)≠At2At1),At)与

不可交换,应采用式(2-64)所示Peano-Baker级数表示Φtt0),即

【例2-9】 已知线性时变系统状态空间表达式为

试求初始时刻t0=0,初始状态xt0)=

时,输入为单位阶跃信号ut)=1(t)系统的输出响应。

可见,At1At2)=At2At1),即At)与

可交换,则可由式(2-56)的矩阵指数求系统状态转移矩阵Φtt0),即

由以上计算可以看出,若t0≠0,则

则系统的输出响应为