3.2状态转移矩阵计算资料.pptVIP

塞尔维斯特内插法计算矩阵指数函数(4/4) 利用上式去计算矩阵指数函数eAt的关键是如何计算待定函数?i(t)。 下面分 A的特征值互异 A有重特征值 两种情况来讨论如何计算?i(t)以及eAt。 A的特征值互异(1/4) (1) A的特征值互异 设矩阵A的n个互异特征值为?1,?2,…,?n,则矩阵A的最小多项式?(?)等于特征多项式f(?)=|?I-A|=?n+a1?n-1+…+an-1?+an。 因系统的所有特征值?i使特征多项式f(?i)=0,故与前面证明过程类似,我们亦有 其中待定函数?i(t)(i=0,1,…,n-1)与矩阵指数函数eAt的表达式中的?i(t)一致。 A的特征值互异(2/4) 因此,可得如下待定函数?i(t)(i=0,1,…,n-1)的线性方程组: 求解上述方程得函数?i(t)后,由式(3-49)可计算得矩阵指数函数eAt。 A的特征值互异(3/4)--例3-7 例3-7 试求如下系统矩阵的矩阵指数函数 解 由于矩阵A的3个特征值互异,并分别为-1,-2和-3,因此解方程组(3-52)可得 A的特征值互异(4/4) 则系统的状态转移矩阵为 A有重特征值(1/4) (2) A有重特征值 由于矩阵A与它的约旦矩阵 具有相同的最小多项式?(?),因此由前面的推导过程可知,约旦矩阵 也满足 设A与 的特征值?i的代数重数为mi,则由上