第三章线性系统状态方程的解

  • 格式:doc
  • 大小:601.00 KB
  • 文档页数:19

3-1 第三章 线性系统的运动分析 §3-1线性连续定常齐次方程求解 一、齐次方程和状态转移矩阵的定义 1、齐次方程 状态方程的齐次方程部分反映系统自由运动的状况(即没有输入作用的状况),设系统

的状态方程的齐次部分为:)()(tAxtx 线性定常连续系统:Axx 2、状态转移矩阵的定义 齐次状态方程Axx有两种常见解法:(1)幂级数法;(2)拉氏变换法。其解为

)0()(xetxAt。其中Ate称为状态转移矩阵(或矩阵指数函数、矩阵指数),记为:Atet)(。

若初始条件为)(0tx,则状态转移矩阵记为:)(00)(ttAett 对于线性时变系统,状态转移矩阵写为),(0tt,它是时刻t,t0的函数。但它一般不能写成指数形式。

(1)幂级数法 设Axx的解是t的向量幂级数

kktbtbtbbtx2210)(

式中,,,,,kbbbb210都是n维向量,则

1232132)(kktkbtbtbbtx )(2210kktbtbtbbA 故而有:





00323021201!1!31312121bAkbbAAbbbAAbbAbb

KK

 3-2

且有0)0(bx。 故 kktbtbtbbtx2210)(

kktbAktbAtAbb020200!1!2

1

)0()!1!21(22xtAktAAtIkk

定义:022!1!1!21KkkkkAttAktAktAAtIe 则)0()(xetxAt。 (2)拉氏变换解法 将Axx两端取拉氏变换,有

)()0()(sAxxssx

)0()()(xsxAsI )0()()(1xAsIsx 拉氏反变换,有 )0(])[()(11xAsILtx 则 ])[()(11AsILetAt

【例3.1.1】 已知系统的状态方程为xx

00

10,初始条件为)0(x,试求状态转移矩阵

和状态方程的解。 解:(1)求状态转移矩阵

kkAttAktAAtIet!1!21)(22

此题中:

0010A, 

000032nAAA

所以 3-3



1010001001)(ttAtIetAt

(2)状态方程的解 )0(101)0()(xtxetxAt

【例3.1.2】 已知系统状态方程为xx

32

10,初始条件为)0(x,试求状态方程的

解。 解:)0()(xetxAt

321321000ssssAsI

2211221221112112213)2)(1(1)(1ssssssssssssAsI ttttttttAteeeeeeeeAsILet2222112222])[()( 故而 )0(2222)0()(2222xeeeeeeeexetxttttttttAt

二、状态转移矩阵Ate的性质 kkAttAktAAtIet!1!21)(22 (1)I)0( (2)AttAt)()()( A)0( (3))()()()()(122121tttttt 证明:)()()()()(1221)()()(212121tttteeetttAtAttA

%Example 3.1.2: %MATLAB syms s t x; A=sym('[0,1;-2,-3]'); I=eye(2); L=inv(s*I-A) lap=ilaplace(L) x=lap*x 3-4

(4))()(1tt,)()(1tt 证明:)()()()()()0(1ttItttt (5))()()(00txtttx 证明:)0()()(xttx )()()0()0()()(00100txtxxttx,代入上式 ∴)()()()()()(00001txtttxtttx 证毕。

(6))()()(011202tttttt 证明:)()()(0022txtttx………………………. …………………(1) )()()(0011txtttx……………………………………………(2) )()()()()()(001121122txtttttxtttx…………….(3) 比较(1)、(3)式,有)()()(011202tttttt成立。证毕。

(7))()(kttk 证明:)(][)()(kteeetktAkAtkAtk (8)若BAAB,则AtBtBtAttBAeeeee)( 若BAAB,则AtBtBtAttBAeeeee)( (9)设)(t为Axx的状态转移矩阵,引入非奇异变换xPx后的状态转移矩阵为: PePtAt1)( 证明:将xPx代入Axx中,有 xAPPx1 APtPet1)(

kkAPtPtAPPktAPPAPtPIe)(!1)(!21122111 3-5

kktAPPktAPPAPtPPP)(!1)(!21122111 PtAktAAtIPkk)!1!21(221 PePAt1 ∴PePtAt1)(。证毕。

(10)两种常见的状态转移矩阵 ①设],,,[21ndiagA,即A为对角阵,且具有互异元素。则

ttneet00)(1 ②设A为mm约当阵

mmA



11

,则ttmtttmttteetmteeetmetteet000)!2(10)!1(1!21)(212

【例3.1.3】 已知状态转移矩阵为 ttttttttAteeeeeeeee22222222

试求)(1t和A。 解:(1)根据状态转移矩阵的性质4,可知

tttttttteeeeeeeett222212222)()( (2)根据状态转移矩阵的性质2,可知 32100442222)0(2222teeeeeeeeAtttttttt 3-6

【例3.1.4】 已知 4401101A

试求状态转移矩阵Ate。 解:根据状态转移矩阵的性质10,可知

ttttttttttAteteeetteeetetteeet000002106121)(232

【例3.1.5】 验证如下矩阵是否为状态转移矩阵。 ttttsincos0cossin0001

解:利用性质(1)I)0( Ittttt0101000010sincos0cossin0001,所以该矩阵不是状态转移矩阵。 【例3.1.6】 已知系统状态方程为Axx, 当11)0(x时,tteetx22)(

当12)0(x时,tteetx2)( 试求系统矩阵A和状态转移矩阵Ate。 解:由性质(2)可知:)0(A 由已知,有

)0()(xetxAt 1121222Attttteeeee