第二章控制系统状态空间分析概要.doc

该【第二章控制系统状态空间分析概要 】是由【泰山小桥流水】上传分享，文档一共【29】页，该文档可以免费在线阅读，需要了解更多关于【第二章控制系统状态空间分析概要 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。精品文档精品文档2精品文档《自动控制原理》 (下部)第二章第二章 控制系统的状态空间分析在讨论了状态方程的描述、标准形和模型转换后,本章将讨论线性多变量系统的运动分析,包括线性状态方程的求解和系统的能控与能观性分析。 线性系统的解给定线性定常系统非齐次状态方程为Σ:x(t)Ax(t)Bu(t)()其中,x(t)Rn,u(t)Rr,ARnn,BRnr,且初始条件为x(t)t0x(0)。将方程()写为x(t)Ax(t)Bu(t)在上式两边左乘e-At,可得eAt[x(t)Ax(t)]d[eAtx(t)]eAtBu(t)dt将上式由O积分到t,得eAtx(t)x(0)teABu()do故可求出其解为x(t)teA(t)Bu()deAtx(0)()o或x(t)(t)x(0)t(t)Bu()d()o式中(t)eAt为系统的状态转移矩阵。对于线性时变系统非齐次状态方程,x(t)A(t)x(t)B(t)u(t)()类似可求出其解为x(t)(t,0)x(0)t精品文档精品文档2精品文档(t,)B()u()d()o一般说来,线性时变系统的状态转移矩阵 (t,t0)只能表示成一个无穷项之和,只有在特殊情况下,才能写成矩阵指数函数的形式。1精品文档精品文档29精品文档《自动控制原理》 (下部)(t,t0)是满足如下矩阵微分方程和初始条件(t,t0)A(t)(t,t0)()(t0,t0) I的解。下面不加证明地给出线性时变系统状态转移矩阵的几个重要性质:1、 (t,t) I;艙浏獲诂鳎页剀鹞貰镐锓电腸繞龔。2、(t2,t1)(t1,t0)(t2,t0);3、1(t,t0)(t0,t);4、当A给定后,(t,t0)唯一;5、计算时变系统状态转移矩阵的公式(t,t0)tA()dtA(1)12)d2d1()It0A(t0t0上式一般不能写成封闭形式,可按精度要求,用数值计算的方法取有限项近似。特别地,只有当满足A(t)tt)dA(t)A()dA(t0t0即在矩阵乘法可交换的条件下, (t,t0)才可表示为如下矩阵指数函数形式(t,t0)texpA()d()t0显然,定常系统的状态转移矩阵(tt0)不依赖于初始时刻t0,其性质仅是上述时变系统的特例。[]试求如下线性定常系统x101x1x223x2的状态转移矩阵Ф(t)和状态转移矩阵的逆Ф-1(t)。[解]对于该系统,0 1A2 3其状态转移矩阵由下式确定(t) eAt L1[(sI A)1]由于精品文档精品文档4精品文档2精品文档精品文档29精品文档《自动控制原理》 (下部)第二章sIAs001s10s232s3其逆矩阵为(sIA)11s31(s2)2s1)(ss31(s1)(s2)(s1)(s2)2s(s1)(s2)(s1)(s2)因此(t)eAtL1[(sIA)1]=2ete2tete2t2et2e2tet2e2t由于Ф-1(t)=Ф(-t),故可求得状态转移矩阵的逆为1(t)eAt2ete2tete2t2et2e2tet2e2t[]求下列系统的时间响应:x101x10x223x2u1式中,u(t)为t=0时作用于系统的单位阶跃函数,即u(t)=1(t)。[解]对该系统A01B0231状态转移矩阵(t),即(t)eAt2ete2tete2t2et2e2tet2e2t因此,系统对单位阶跃输入的响应为:x(t)eAtt2e(t)e2(t)e(t)e2(t)01(t)dx(0)(t)2(t)(t)2(t)1o2e2e2ee或x1(t)2ete2tete2tx1(0)1et1e2tx2(t)2et2e2tet2e2tx2(0)2t2t2ee精品文档精品文档7精品文档精品文档精品文档7精品文档3精品文档精品文档29精品文档《自动控制原理》 (下部)第二章如果初始状态为零,即 X(0)=0,可将X(t)简化为x1(t)1et1e2tx2(t),状态方程的解实质上可归结为计算状态转移矩阵,即矩阵指数函数Ate。如果给定矩阵A中所有元素的值,MATLAB将提供一种计算eAT的简便方法,其中T为常数。除了上述方法外,对eAt的计算还有几种分析方法可供使用。这里我们将介绍其中的四种计算方法。方法一:直接计算法 (矩阵指数函数 )讞癮絷驿靈鄭员跷糧鮮顎鐘顰坟恹。eAtIAtA2t2A3t31Aktk()2!3!k0k!可以证明,对所有常数矩阵A和有限的t值来说,这个无穷级数都是收敛的。方法二:对角线标准形与Jordan标准形法若可将矩阵A变换为对角线标准形,那么Ate可由下式给出e1t0e2teAtPetP1PP1()0ent式中,P是将A对角线化的非奇异线性变换矩阵。类似地,若矩阵A可变换为Jordan标准形,则At可由下式确定出eeAt=SeJtS–1()[]考虑如下矩阵A010A001133[解]该矩阵的特征方程为|IA|33231(1)30因此,矩阵A有三个相重特征值λ=1。可以证明,矩阵A也将具有三重特征向量(即有两个广义特征向量)。易知,将矩阵A变换为Jordan标准形的变换矩阵为1 0 0S 1 1 01 2 1袭櫛瘋诼泺励虾厩遞滎銳羈殒嫔塏。精品文档精品文档9精品文档精品文档精品文档29精品文档4精品文档精品文档11精品文档《自动控制原理》 (下部)第二章矩阵S的逆为1 0 0S1 1 1 01 2 1于是100010100S1AS110001110121133121110011J001注意到ettet1t2et2eJt0ettet00et可得eAt=SeJtS–1即100ettet1t2et100et21100tet11012100et121et tet 1t2ettett2et 1t2et2 21t2et ettett2et tet1t2et2 2tet 1t2et 3tet t2et et 2tet 1t2et2 2方法三:拉氏变换法eAt L1[(sI A)1] ()为了求出 eAt,关键是必须首先求出( sI-A)的逆。一般来说,当系统矩阵 A的阶次较高时,可采用递推算法。[] 考虑如下矩阵 A亲笕鴯觴闡镞玀组鴆盏谅钼幣顥骟。精品文档精品文档12精品文档5精品文档精品文档29精品文档《自动控制原理》 (下部)第二章A0102试用前面介绍的两种方法计算eAt。[解]方法一由于A的特征值为0和-2(λ1=0,λ2=-2),故可求得所需的变换矩阵P为11P=20因此,由式)可得11eo01111(1e2t)At2e2020e2t100e2t2方法二由于sIAs001s10s020s2可得11(sIA)1ss(s2)01s2因此AtL1[(sIA)1]11(1e2t)e20e2t方法四:化eAt为A的有限项法(Caley-Hamilton定理法)第四种是利用凯莱-哈密尔顿定理,化eAt为A的有限项,然后通过求待定时间函数获得eAt的方法。凯莱-哈密尔顿(Caley-Hamilton)定理在证明有关矩阵方程的定理或解决有关矩阵方程的问题时,凯莱-哈密尔顿定理是非常有用的。考虑n×n维矩阵A及其特征方程|IA|na1n1an1an0凯莱-哈密尔顿定理指出,矩阵A满足其自身的特征方程,即Ana1An1an1AanI0()为了证明此定理,注意到(λI-A)的伴随矩阵adj(λI-A)是λ的n-1次多项式,即精品文档精品文档14精品文档adj(IA)B1n1B2n2Bn1Bn精品文档精品文档15精品文档6精品文档精品文档29精品文档