矩阵论课程论文

西安理工大学

研究生课程论文报告

课程名称：矩阵论

课程代号：

任课教师：

论文报告题目：矩阵函数在线性定常系统

状态转移矩阵求解中的应用完成日期：2015 年10 月25 日学科：电力电子与电力传动

学号：

姓名：

成绩：

矩阵函数在线性定常系统状态转移矩阵

求解中的应用

摘 要

控制系统的运动是系统性能定量分析的重要内容。“运动”是物理学上的一个概念，它是通过求系统方程的解)(t x 、)(t y 来分析研究的。由于状态方程是矩阵微分（差分）方程，输出方程式为矩阵代数方程，因此求系统方程的解主要是求状态方程的解。而求状态方程的解的关键是求状态转移矩阵。本文主要介绍了矩阵对角化标准型，约当标准型，凯莱-哈密顿定理及矩阵函数知识在线性定常系统的齐次状态方程的状态转移矩阵求解中的应用。

关键词：状态转移矩阵，约当标准型，凯莱-哈密顿定理，矩阵函数.

1.问题提出

线性系统有线性定常系统和线性时变系统，最为基本的是线性定常系统。而线性定常系统根据有无初始输入，分为线性定常齐次方程，和线性定常非齐次方程。本文只给出线性定常系统的齐次状态方程的状态转移矩阵的求解。

线性定常系统齐次方程的解亦即系统的自由解，是指系统输入为零时，由初始状态引起的自由运动。

线性定常系统齐次状态方程为

()()t Ax t x

= ()1-1

其中，x 是n 维状态向量；A 为n n ⨯系数矩阵。设初始时刻00=t ，系统的初始状态()()00x t x =。仿照标量微分方程求解的方法求方程()1-1的解。

设方程()1-1的解为t 的向量幂级数形式，即

)(t x = ++++++k k t b t b t b t b b 332210 ()2-1

式中，() ,2,1,0=i b i 为n 维向量。

式()2-1代入方程()1-1得

()

+++++=+++++-k k k k t b t b t b b b A t kb t b t b b 3322101232132 ()3-1

既然式()2-1是方程()1-1的解，则式()3-1对任意的t 都成立。因此，式()3-1的等式两边t 的同次幂项的系数应相等，有

⎪⎪⎪⎭

⎪⎪

⎪⎬⎫=======-0!1110

3!31

231302!21

12120

1b A Ab b b A Ab b b A Ab b Ab b k k k k k

()4-1 Ⅰ.当0=t 时，由式()2-1可得到

()00x b = ()5-1

将式()4-1和式()5-1代入式()2-1，得到齐次状态方程的解

()()

()0!1

22!21x t A t A At I t x k k k +++++= (

)6-1 上边右边括号内的级数是n n ⨯矩阵指数函数，记成At e ，即

+++++=k k k At t A t A At I e !

1

22!21 ()7-1 所以式()6-1可写成

()()0x e t x At = ()8-1

Ⅱ.如果初始时刻00≠t ，初始状态为()0t x ，则齐次状态方程的解为

()()()00x e t x t t A -= ()9-1

由上式可知，系统在状态空间的任一时刻t 的状态()t x ，可视为系统的初始状态()0t x 通过矩阵指数函数()0t t A e -的转移而得到的。因此，矩阵指数函数()0t t A e -又称为状态转移矩阵。

从上面的分析看，求状态方程的解()t x ，关键是求矩阵指数At e 。

2.问题求解

2.1 矩阵指数的基本性质

在介绍求矩阵指数At e 的方法之前，先介绍At e 的一些主要性质和几个特殊的指数函数：

（1）∑∞

==0!

k k

k At

k t A e ，该无穷级数在有限时间时绝对收敛的

（2）At

At Ae e dt

d =

（3）()2121At At t t A e e e ⋅=+ （4）[]

At At

e e --=1

（5）若BA AB =，则()t B A Bt At e e e +=⋅； 若BA AB ≠，则()t B A Bt At e e e +≠⋅

（6）若P 为非奇异矩阵，A 通过非奇异变换成对角阵，即AP P A 1-∧

=,则有 11

--=P Pe e AP p At (7) 若A 为对角阵

⎥

⎥⎥⎥⎦⎤⎢⎢⎢

⎢⎣⎡=n A λλλ0021

，则⎥

⎥⎥⎥⎥⎦⎤

⎢⎢⎢⎢⎢⎣⎡=t t

t At n e e e e λλλ0021

2.2状态转移矩阵()t φ的几种计算方法

1.根据At

e 的定义直接计算

() +++++==n n n At t A t A At I e t !1

22!21φ

2.拉普拉斯变换法

对于线性定常系统的齐次状态方程

()()t Ax t x

= 两边求拉普拉斯变换，得

()()()s Ax x s sX =-0，

即()()()0x s X A sI =-， 有()()()01

x A sI s X --=

因此，()()[]

()01

1x A sI L e t x At ---==

若初始时刻00=t ，初始状态为()0x ，则对上式进行拉普拉斯变换，得

()()[]1

1---==A sI L e t At L

φ 3.非奇异线性变换

（1）矩阵A 经线性变换化为对角线矩阵Λ求At

e

当矩阵A 的n 个特征值互异或者虽有重根但是仍有n 个独立的特征向量时，经过线性变换，将A 化为对角形矩阵Λ，即

⎥

⎥⎥⎥⎦⎤

⎢⎢⎢⎢⎣⎡=Λ=-n PAP λλλ00 211 此时，系统的状态转移矩阵