第三章线性系统状态方程的解
- 格式: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 )(2210kktbtbtbbA 故而有:
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 则 ])[()(11AsILetAt
【例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为对角阵,且具有互异元素。则
ttneet00)(1 ②设A为mm约当阵
mmA
11
,则ttmtttmttteetmteeetmetteet000)!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】 已知 4401101A
试求状态转移矩阵Ate。 解:根据状态转移矩阵的性质10,可知
ttttttttttAteteeetteeetetteeet000002106121)(232
【例3.1.5】 验证如下矩阵是否为状态转移矩阵。 ttttsincos0cossin0001
解:利用性质(1)I)0( Ittttt0101000010sincos0cossin0001,所以该矩阵不是状态转移矩阵。 【例3.1.6】 已知系统状态方程为Axx, 当11)0(x时,tteetx22)(
当12)0(x时,tteetx2)( 试求系统矩阵A和状态转移矩阵Ate。 解:由性质(2)可知:)0(A 由已知,有
)0()(xetxAt 1121222Attttteeeee