矩阵论课程论文

西安理工大学研究生课程论文报告课程名称：矩阵论课程代号：任课教师：论文报告题目：矩阵函数在线性定常系统状态转移矩阵求解中的应用完成日期：2015年10月25日学科：电力电子与电力传动学号：姓名：成绩：矩阵函数在线性定常系统状态转移矩阵求解中的应用摘要控制系统的运动是系统性能定量分析的重要内容。“运动”是物理学上的一个概念，它是通过求系统方程的解)(tx、)(ty来分析研究的。由于状态方程是矩阵微分（差分）方程，输出方程式为矩阵代数方程，因此求系统方程的解主要是求状态方程的解。而求状态方程的解的关键是求状态转移矩阵。本文主要介绍了矩阵对角化标准型，约当标准型，凯莱-哈密顿定理及矩阵函数知识在线性定常系统的齐次状态方程的状态转移矩阵求解中的应用。关键词：状态转移矩阵，约当标准型，凯莱-哈密顿定理，矩阵函数.1.问题提出线性系统有线性定常系统和线性时变系统，最为基本的是线性定常系统。而线性定常系统根据有无初始输入，分为线性定常齐次方程，和线性定常非齐次方程。本文只给出线性定常系统的齐次状态方程的状态转移矩阵的求解。线性定常系统齐次方程的解亦即系统的自由解，是指系统输入为零时，由初始状态引起的自由运动。线性定常系统齐次状态方程为tAxtx1-1其中，x是n维状态向量；A为nn系数矩阵。设初始时刻00t，系统的初始状态00xtx。仿照标量微分方程求解的方法求方程1-1的解。设方程1-1的解为t的向量幂级数形式，即)(tx=kktbtbtbtbb3322102-1式中，,2,1,0ibi为n维向量。式2-1代入方程1-1得kkkktbtbtbbbAtkbtbtbb33221012321323-1既然式2-1是方程1-1的解，则式3-1对任意的t都成立。因此，式3-1的等式两边t的同次幂项的系数应相等，有0!11103!31231302!21121201bAAbbbAAbbbAAbbAbbkkkkk4-1Ⅰ.当0t时，由式2-1可得到00xb5-1将式4-1和式5-1代入式2-1，得到齐次状态方程的解0!122!21xtAtAAtItxkkk6-1上边右边括号内的级数是nn矩阵指数函数，记成Ate，即kkkAttAtAAtIe!122!217-1所以式6-1可写成0xetxAt8-1Ⅱ.如果初始时刻00t，初始状态为0tx，则齐次状态方程的解为00xetxttA9-1由上式可知，系统在状态空间的任一时刻t的状态tx，可视为系统的初始状态0tx通过矩阵指数函数0ttAe的转移而得到的。因此，矩阵指数函数0ttAe又称为状态转移矩阵。从上面的分析看，求状态方程的解tx，关键是求矩阵指数Ate。2.问题求解2.1矩阵指数的基本性质在介绍求矩阵指数Ate的方法之前，先介绍Ate的一些主要性质和几个特殊的指数函数：（1）0!kkkAtktAe，该无穷级数在有限时间时绝对收敛的（2）AtAtAeedtd（3）2121AtAtttAeee（4）AtAtee1（5）若BAAB，则tBABtAteee；若BAAB，则tBABtAteee（6）若P为非奇异矩阵，A通过非奇异变换成对角阵，即APPA1,则有11PPeeAPpAt(7)若A为对角阵nA0021，则tttAtneeee00212.2状态转移矩阵t的几种计算方法1.根据Ate的定义直接计算nnnAttAtAAtIet!122!212.拉普拉斯变换法对于线性定常系统的齐次状态方程tAxtx两边求拉普拉斯变换，得sAxxssX0，即0xsXAsI，有01xAsIsX因此，011xAsILetxAt若初始时刻00t，初始状态为0x，则对上式进行拉普拉斯变换，得11AsILetAtL3.非奇异线性变换（1）矩阵A经线性变换化为对角线矩阵求Ate当矩阵A的n个特征值互异或者虽有重根但是仍有n个独立的特征向量时，经过线性变换，将A化为对角形矩阵，即nPAP00211此时，系统的状态转移矩阵tttnnnnAtneeettttttttttIe00000000002122!21222!212221!2112221!212122!21111111由于PPA1所以矩阵A的状态转移矩阵PePPttIPPtPPtPPPtPPPtPPPtPPPtPIeettPtPAt122!21122!21111221!2111221!2111（2）矩阵A经线性变换化为约当形矩阵J求Ate当矩阵A的n个特征值均相同，且为1时，经过线性变换，可化为约当形矩阵J111111100JPAP则tnnJtettnttnte11!211!111210所以，系统的状态转移矩阵为PePetJtAt14.应用凯莱-哈密顿定理首先介绍一下凯莱-哈密顿定理：nn矩阵A满足自身的特征方程，即矩阵A的特征多项式是A的零化多项式。0det012211aaaaAInnnnn即012211-----aaaannnnn根据凯莱-哈密顿定理，有0012211IaAaAaAaAAnnnnn于是IaAaAaAaAnnnnn012211-----上式表明，nA是1-nA，2-nA，…A，I的线性组合。显然有AaAaAaAaAAAnnnnnn0211211-----则IaaAaaaAaaaAaaAnnnnnnnnnn0101122211-2121)()(-）（依次类推，可得2-1-,nnAA，…均是1-nA，2-nA，…A，I的线性组合。那么，Ate就化成一个A的最高幂次为1n的n项幂级数的形式，即1110!122!21nnnnnAtAtaAtaItatAtAAtIet）（(1)A的特征值),,2,1(nii互异应用凯莱-哈密顿定理，i和A均是特征多项式的零根。因此，nitatataeniniti,,2,11110那么，tatataeeennnnnnntttn11012122221121111121于是，tttnnnnnnnneeetatata211-121222211211110111(2)A的特征值均相同设A的特征值为1，待定系数tai的计算公式如下ttttnntnnnnnnnnneetetetetnnnntatatatata111111!112!212!211!111-112131211212113111123101!113210!22131001110000(3)A的n个特征值有重特征值和互异特征值当A的n个特征值有重特征值和互异特征值时，待定系数tai可以根据21综合得出，然后求出状态转移矩阵t。3.举例计算已知3-2-10A，分别用上述四种方法求解状态转移矩阵t。解：1定义法根据定义计算ttttttttkkkAteeeeeeeettttttttttttttAtAAtIet22223252273372367223322!122!212222313213210!21321010012拉普拉斯变换法22112212211121122121332111ssssssssssssssAsI那么，ttttttttAteeeeeeeeAsILe22221122223化矩阵为对角线标准型由021-AI得特征根2-1-21，211-1-222-1-111111-21APPPadjPPP所以，tttttttttttAteeeeeeeeeePPee2222212001222211120021114应用凯莱-哈密顿定理已知特征根2-1-21，，两两相异，则有ttttttttttttAttttttttttteeeeeeeeeeeeAtaItaeeeeeeeeeeetata22222210222211211022223210100122111221111121可见，四种算法的计算结果是一样的。4.应用小结本文在给出线性定常系统的齐次状态方程的状态转移矩阵的求解时，用到的矩阵论知识主要有矩阵的特征多项式，特征值，特征向量，矩阵逆，矩阵乘法，幂级数等一些基础知识。其次，在给出求解状态转移矩阵的四种方法中，其中利用线性变换将矩阵A化成对角线形，约当形；应用凯莱-哈密顿定理，再根据矩阵指数函数的性质对矩阵函数Ate进行求解，从而计算出状态转移矩阵t。通过这些矩阵论知识的应用，使得复杂的线性定常系统状态转移矩阵的求解变得方便了许多，对解决本专业问题有很大的意义。参考文献【1】王孝武主编.现代控制理论基础[M].2版.北京：机械工业出版社,2006.【2】郑大钟．线性系统理论[M].2版.北京：清华大学出版社，2002．【3】刘豹，唐丙生.现代控制理论[M].2版.北京：机械工业出版社，2005．【4】施颂椒，陈学中，杜秀华.现代控制理论基础[M].北京：高等教育出版社，2007．【5】王枞.控制系统理论及应用[M].北京：北京邮电大学出版社，2005.【6】王显正，