控制系统的运动分析.ppt
《控制系统的运动分析.ppt》由会员分享,可在线阅读,更多相关《控制系统的运动分析.ppt(41页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、Ch.3 Ch.3 控制系统的运动分析控制系统的运动分析状态转移矩阵计算状态转移矩阵计算(1/1)3.2 状态转移矩阵计算状态转移矩阵计算 q在状态方程求解中,关键是状态转移矩阵(t)的计算。对于线性定常连续系统,该问题又归结为矩阵指数函数eAt的计算。上一节已经介绍了基于拉氏反变换技术的矩阵指数函数eAt的计算方法,下面讲述计算矩阵指数函数的下述其他3种常用方法。级数求和法级数求和法约旦规范形法约旦规范形法 化化eAt为为A的有限多项式矩阵函数法的有限多项式矩阵函数法重点推荐级数求和法级数求和法(1/3)3.2.1 级数求和法级数求和法 q由上一节对矩阵指数函数的定义过程中可知:矩阵指数函数
2、eAt的计算可由上述定义式直接计算。q由于上述定义式是一个无穷级数,故在用此方法计算eAt时必须考虑级数收敛性条件和计算收敛速度问题。类似于标量指数函数eat,对所有有限的常数矩阵A和有限的时间t来说,矩阵指数函数eAt这个无穷级数表示收敛。级数求和法级数求和法(2/3)q显然,用此方法计算eAt一般不能写成封闭的、简洁的解析形式,只能得到数值计算的近似计算结果。其计算精度取决于矩阵级数的收敛性与计算时所取的项数的多少。如果级数收敛较慢,则需计算的级数项数多,人工计算是非常麻烦的,一般只适用于计算机计算。因此,该方法的缺点:计算量大精度低非解析方法,难以得到计算结果的简洁的解析表达式。级数求和
3、法级数求和法(3/3)例例3-4q例例3-4 用直接计算法求下述矩阵的矩阵指数函数:q解 按矩阵指数函数的展开式计算如下:约旦规范形法约旦规范形法(1/8)3.2.2 约旦规范形法约旦规范形法 q上节给出了对角线矩阵、块对角矩阵和约旦块三种特殊形式矩阵的矩阵指数函数。由于任何矩阵都可经线性变换成为对角线矩阵或约旦矩阵,因此可通过线性变换将一般形式的矩阵变换成对角线矩阵或约旦矩阵,再利用上述特殊形式矩阵的矩阵指数函数来快速计算矩阵矩阵指数函数。下面讨论之。约旦规范形法约旦规范形法(2/8)q下面首先讨论矩阵指数函数的一条性质:对矩阵A,经变换矩阵P作线性变换后,有则相应地有如下矩阵指数函数的变换
4、关系约旦规范形法约旦规范形法(3/8)q该结论可简单证明如下:q根据上述性质,对矩阵A,可通过线性变换方法得到对角线矩阵或约旦矩阵,然后利用该类特殊矩阵的矩阵指数函数,由矩阵指数函数的变换关系来求原矩阵A的矩阵指数函数。约旦规范形法约旦规范形法(4/8)例例3-5q例3-5 试求如下系统矩阵的矩阵指数函数q解解 1.先求先求A的特征值。的特征值。由特征方程可求得特征值为1=-1 2=-2 3=-32.求特征值所对应的特征向量。求特征值所对应的特征向量。由前述的方法可求得特征值1,2和3所对应的特征向量分别为p1=1 0 1 p2=1 2 4 p3=1 6 9约旦规范形法约旦规范形法例例3-5q
5、故将A变换成对角线矩阵的变换矩阵P及其逆阵P-1为3.由系统矩阵和矩阵指数函数的变换关系,分别有约旦规范形法约旦规范形法例例3-6q例3-6 试求如下系统矩阵的矩阵指数函数约旦规范形法约旦规范形法(7/8)例例3-6q解解 1.先求先求A的特征值。的特征值。由特征方程可求得特征值为1=2 2=3=-12.由于矩阵A为友矩阵,故将A变换成约旦矩阵的变换矩阵P和其逆阵P-1分别为3.由系统矩阵和矩阵指数函数的变换关系,分别有约旦规范形法约旦规范形法(8/8)-例例3-6塞尔维斯特内插法塞尔维斯特内插法(1/1)3.2.3 塞尔维斯特内插法塞尔维斯特内插法q在讨论塞尔维斯特(Sylvester)内插
6、法计算矩阵指数函数eAt时,需要用到关于矩阵特征多项式的凯莱-哈密顿(Cayley-Hamilton)定理以及最小多项式的概念。因此,首先给出凯莱-哈密顿定理及最小多项式的概念,再讨论塞尔维斯特内插法。下面依次介绍:凯莱凯莱-哈密顿定理哈密顿定理最小多项式最小多项式塞尔维斯特内插法计算矩阵指数函数塞尔维斯特内插法计算矩阵指数函数凯莱-哈密顿定理(1/4)1.凯莱-哈密顿定理q凯莱-哈密顿定理是矩阵方程分析和求解中非常重要的定理,其表述和证明如下。q定理3-1(凯莱-哈密顿定理)设nn矩阵A的特征多项式为f()=|I-A|=n+a1n-1+an-1+an则矩阵A必使由上述特征多项式决定的矩阵多项
7、式函数f(A)=An+a1An-1+an-1A+anI=0上述特征多项式亦称为矩阵A的零化特征多项式。凯莱-哈密顿定理(2/4)q证明 因为I=(I-A)-1(I-A)=adj(I-A)/|I-A|(I-A)故|I-A|I=adj(I-A)(I-A)由伴随矩阵的定义可知,伴随矩阵adj(I-A)可表示为如下多项式矩阵函数:adj(I-A)=n-1I+n-2B2+Bn-1+Bn其中矩阵B2,B3,Bn为nn维的常数矩阵。凯莱-哈密顿定理(3/4)因此由前面两式,有(n+a1n-1+an-1+an)I=(n-1I+n-2B2+Bn-1+Bn)(I-A)整理得 (n+a1n-1+an-1+an)I
8、=nI+(B2-A)n-1+(Bn-Bn-1A)-BnA凯莱-哈密顿定理(4/4)上式中,令等号两边的同幂次项的系数相等,则有a1I-B2+A=0a2I-B3+AB2=0 an-1I-Bn+ABn-1=0anI+ABn=0因此,将上述各等式从上至下依次右乘以An-1,A,I,然后将各等式相加,即得An+a1An-1+an-1A+anI=0故矩阵A满足其本身的零化特征多项式。最小多项式最小多项式(1/3)2.最小多项式最小多项式 q根据凯莱-哈密尔顿定理,任一nn维矩阵A满足其自身的特征方程,即特征多项式为A的一个零化多项式。然而特征多项式不一定是A的最小阶次的零化多项式。将矩阵A满足的最小阶次
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 控制系统 运动 分析
限制150内