状态转移矩阵计算概述.ppt

最小多项式(2/3) 最小多项式在矩阵多项式的分析与计算中起着重要作用。 定理3-2给出了特征多项式与最小多项式的关系。 定理3-2 设首一多项式d(?)是?I-A的伴随矩阵adj(?I-A)的所有元素的最高公约式,则最小多项式为 ?(?)=?m+?1?m-1+…+?m-1?+?m 最小多项式(3/3) 证明 由假设知,矩阵adj(?I-A)的最高公约式为d(?),故 adj(?I-A)=d(?)B(?), 式中,B(?)的n2个元素(为?的函数)的最高公约式为1。 由于 (?I-A)adj(?I-A)=|?I-A|I 可得 d(?)(?I-A)B(?)=|?I-A|I 由上式可知,特征多项式|?I-A|可被整除d(?)。 因此设d(?)整除|?I-A|得到的因式记为?(?),故有 |?I-A|=d(?)?(?), 最小多项式(4/3) 由于首一多项式d(?)的最高阶次的系数为1,所以?(?)的最高阶次的系数也应为1。 因此,综合上两式,可得 (?I-A)B(?)=?(?)I 因而 ?(A)=0 即?(?)亦为A的零化多项式。 设?(?)为A的最小多项式,因此零化多项式?(?)可写为 ?(?)=g(?)?(?)+e(?) 其中g(?)和e(?)分别是多项式?(?)除以?(?)的商和余项,且e(?)的阶次低于?(?)。 最小多项式(5/3) 由于?(A)=0和?(A)=0,所以必然有e(A)=0。 考虑到?(?)为矩阵A的最小多项式,所以不存在比?(?)阶次还低的A的零化多项式,故e(?)必为零,即有 ?(?)=g(?)?(?) 又因为?(A)=0,所以?(?)可写为 ?(?)I=(?I-A)H(?) 式中,H(?)为?(?)的一个因子矩阵,故 ?(?)I=g(?)?(?)I=g(?)(?I-A)H(?) 将上式与(?I-A)B(?)=?(?)I比较,有 B(?)=g(?)H(?) 最小多项式(6/3) 又因为B(?)的n2个元素的最高公约式为1,因此 g(?)=1 于是 ?(?)=?(?) 因此,由前面证明的|?I-A|=d(?)?(?)而证明了最小多项式?(?)为 最小多项式(7/3) 根据上述定理3-2,n×n维矩阵A的最小多项式可按以下步骤求出。 1) 根据伴随矩阵adj(?I-A),写出作为?的因式分解多项式的adj(?I-A)的各元素; 2) 确定作为伴随矩阵adj(?I-A)各元素的最高公约式d(?)。 选取d(?)的最高阶次系数为1。 如果不存在公约式,则d(?)=1; 3) 最小多项式?(?)可由|?I-A|除以d(?)得到。 塞尔维斯特内插法计算矩阵指数函数(1/4) 3. 塞尔维斯特内插法计算矩阵指数函数 基于最小多项式(或特征多项式),塞尔维斯特内插法可以非常简洁、快速地计算出矩阵指数函数,其计算思想与过程可描述如下。 若?(?)=?m+?1?m-1+…+?m-1?+?m为矩阵A的最小多项式,则由?(A)=0有 Am=-?1Am-1-…-?m-1A-?mI 即Am可用有限项Am-1,…,A,I的线性组合来表示。 塞尔维斯特内插法计算矩阵指数函数(2/4) 将上式两边乘以矩阵A,则有 即Am+1可用有限项Am-1,…,A,I的线性组合来表示。 Understand? 塞尔维斯特内插法计算矩阵指数函数(3/4) 其中?i(t)(i=0,1,…,m-1)为待定的关于时间t的函数。 即,矩阵指数函数eAt亦可以用有限项Am-1,…,A,I的线性函数组合表示。 依次类推,则可知,Ai(im)可用有限项Am-1,…,A,I的线性组合来表示。 因此,我们有 关键喔 塞尔维斯特内插法计算矩阵指数函数(4/4) 利用上式去计算矩阵指数函数eAt的关键是如何计算待定函数?i(t)。 下面分 A的特征值互异 A有重特征值 两种情况来讨论如何计算?i(t)以及eAt。 (1) A的特征值互异 设矩阵A的n个互异特征值为?1,?2,…,?n,则矩阵A的最小多项式?(?)等于特征多项式f(?)=|?I-A|=?n+a1?n-1+…+an-1?+an。 因系统的所有特征值?i使特征多项式f(?i)=0,故与前面证明过程类似,我们亦有 A的特征值互异(1/4) 其中待定函数?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的特征值