状态转移矩阵计算2讲课文档.pptVIP

塞尔维斯特内插法计算矩阵指数函数(1/4)3.塞尔维斯特内插法计算矩阵指数函数基于最小多项式(或特征多项式),塞尔维斯特内插法可以非常简洁、快速地计算出矩阵指数函数,其计算思想与过程可描述如下。若?(?)=?m+?1?m-1+…+?m-1?+?m为矩阵A的最小多项式,则由?(A)=0有Am=-?1Am-1-…-?m-1A-?mI即Am可用有限项Am-1,…,A,I的线性组合来表示。第27页，共41页。塞尔维斯特内插法计算矩阵指数函数(2/4)将上式两边乘以矩阵A,则有即Am+1可用有限项Am-1,…,A,I的线性组合来表示。第28页，共41页。塞尔维斯特内插法计算矩阵指数函数(3/4)依次类推,则可知,Ai(im)可用有限项Am-1,…,A,I的线性组合来表示。因此,我们有其中?i(t)(i=0,1,…,m-1)为待定的关于时间t的函数。即,矩阵指数函数eAt亦可以用有限项Am-1,…,A,I的线性函数组合表示。第29页，共41页。状态转移矩阵计算第1页，共41页。目录(1/1)目录概述3.1线性定常连续系统状态方程的解3.2状态转移矩阵及其计算3.3线性时变连续系统状态方程的解3.4线性定常连续系统的离散化3.5线性定常离散系统状态方程的解3.6Matlab问题本章小结第2页，共41页。状态转移矩阵计算(1/1)3.2状态转移矩阵计算在状态方程求解中,关键是状态转移矩阵?(t)的计算。对于线性定常连续系统,该问题又归结为矩阵指数函数eAt的计算。上一节已经介绍了基于拉氏反变换技术的矩阵指数函数eAt的计算方法,下面讲述计算矩阵指数函数的下述其他3种常用方法。级数求和法约旦规范形法化eAt为A的有限多项式矩阵函数法重点推荐第3页，共41页。级数求和法(1/3)级数求和法由上一节对矩阵指数函数的定义过程中可知:矩阵指数函数eAt的计算可由上述定义式直接计算。由于上述定义式是一个无穷级数,故在用此方法计算eAt时必须考虑级数收敛性条件和计算收敛速度问题。类似于标量指数函数eat,对所有有限的常数矩阵A和有限的时间t来说,矩阵指数函数eAt这个无穷级数表示收敛。第4页，共41页。级数求和法(2/3)显然,用此方法计算eAt一般不能写成封闭的、简洁的解析形式,只能得到数值计算的近似计算结果。其计算精度取决于矩阵级数的收敛性与计算时所取的项数的多少。如果级数收敛较慢,则需计算的级数项数多,人工计算是非常麻烦的,一般只适用于计算机计算。因此,该方法的缺点:计算量大精度低非解析方法,难以得到计算结果的简洁的解析表达式。第5页，共41页。级数求和法(3/3)—例3-4例3-4用直接计算法求下述矩阵的矩阵指数函数:解按矩阵指数函数的展开式计算如下:第6页，共41页。约旦规范形法(1/8)约旦规范形法上节给出了对角线矩阵、块对角矩阵和约旦块三种特殊形式矩阵的矩阵指数函数。由于任何矩阵都可经线性变换成为对角线矩阵或约旦矩阵,因此可通过线性变换将一般形式的矩阵变换成对角线矩阵或约旦矩阵,再利用上述特殊形式矩阵的矩阵指数函数来快速计算矩阵矩阵指数函数。下面讨论之。第7页，共41页。约旦规范形法(2/8)下面首先讨论矩阵指数函数的一条性质:对矩阵A,经变换矩阵P作线性变换后,有则相应地有如下矩阵指数函数的变换关系第8页，共41页。约旦规范形法(3/8)根据上述性质,对矩阵A,可通过线性变换方法得到对角线矩阵或约旦矩阵,然后利用该类特殊矩阵的矩阵指数函数,由矩阵指数函数的变换关系来求原矩阵A的矩阵指数函数。该结论可简单证明如下:第9页，共41页。约旦规范形法(4/8)—例3-5例3-5试求如下系统矩阵的矩阵指数函数解1.先求A的特征值。由特征方程可求得特征值为?1=-1?2=-2?3=-32.求特征值所对应的特征向量。由前述的方法可求得特征值?1,?2和?3所对应的特征向量分别为p1=[101]?p2=[124]?p3=[169]?第10页，共41页。约旦规范形法—例3-5故将A变换成对角线矩阵的变换矩阵P及其逆阵P-1为3.由系统矩阵和矩阵指数函数的变换关系,分别有第11页，共41页。约旦规范形法—例3-6例3-6试求如下系统矩阵的矩阵指数函数第12页，共41页。约旦规范形法(7/8)—例3-6解1.先求A的特征值。由特征方程可求得特征值为?1=2?2=?3=-12.由于矩阵A为友矩阵,故将A变换成约旦矩阵的变换矩阵P