3.2g 状转移矩阵计算.pptVIP

状态转移矩阵计算(1/1) 3.2 状态转移矩阵计算 在状态方程求解中,关键是状态转移矩阵?(t)的计算。 对于线性定常连续系统,该问题又归结为矩阵指数函数eAt的计算。 上一节已经介绍了基于拉氏反变换技术的矩阵指数函数eAt的计算方法,下面讲述计算矩阵指数函数的下述其他3种常用方法。 级数求和法 约旦规范形法 化eAt为A的有限多项式矩阵函数法 级数求和法(1/3) 3.2.1 级数求和法 由上一节对矩阵指数函数的定义过程中可知: 级数求和法(2/3) 显然,用此方法计算eAt一般不能写成封闭的、简洁的解析形式,只能得到数值计算的近似计算结果。 其计算精度取决于矩阵级数的收敛性与计算时所取的项数的多少。 如果级数收敛较慢,则需计算的级数项数多,人工计算是非常麻烦的,一般只适用于计算机计算。 因此,该方法的缺点: 计算量大 精度低 非解析方法,难以得到计算结果的简洁的解析表达式 。 级数求和法(3/3)—例3-4 例3-4 用直接计算法求下述矩阵的矩阵指数函数: 约旦规范形法 (1/8) 3.2.2 约旦规范形法 上节给出了对角线矩阵、块对角矩阵和约旦块三种特殊形式矩阵的矩阵指数函数。 由于任何矩阵都可经线性变换成为对角线矩阵或约旦矩阵,因此 可通过线性变换将一般形式的矩阵变换成对角线矩阵或约旦矩阵, 再利用上述特殊形式矩阵的矩阵指数函数来快速计算矩阵矩阵指数函数。 下面讨论之。 约旦规范形法(2/8) 下面首先讨论矩阵指数函数的一条性质: 对矩阵A,经变换矩阵P作线性变换后,有 则相应地有如下矩阵指数函数的变换关系 约旦规范形法(4/8)—例3-5 例3-5 试求如下系统矩阵的矩阵指数函数 约旦规范形法—例3-5 故将A变换成对角线矩阵的变换矩阵P及其逆阵P-1为 约旦规范形法—例3-6 例3-6 试求如下系统矩阵的矩阵指数函数 约旦规范形法(7/8)—例3-6 解 1. 先求A的特征值。由特征方程可求得特征值为 ?1=2 ?2=?3=-1 2. 由于矩阵A为友矩阵,故将A变换成约旦矩阵的变换矩阵P和其逆阵P-1分别为 约旦规范形法(8/8)--例3-6 塞尔维斯特内插法(1/1) 3.2.3 塞尔维斯特内插法 在讨论塞尔维斯特(Sylvester)内插法计算矩阵指数函数eAt时,需要用到关于矩阵特征多项式的凯莱-哈密顿(Cayley-Hamilton)定理以及最小多项式的概念。 因此,首先给出凯莱-哈密顿定理及最小多项式的概念,再讨论塞尔维斯特内插法。 下面依次介绍: 凯莱-哈密顿定理 最小多项式 塞尔维斯特内插法计算矩阵指数函数 凯莱-哈密顿定理(1/4) 1. 凯莱-哈密顿定理 凯莱-哈密顿定理是矩阵方程分析和求解中非常重要的定理,其表述和证明如下。 定理3-1(凯莱-哈密顿定理) 设n?n矩阵A的特征多项式为 f(?)=|?I-A|=?n+a1?n-1+…+an-1?+an 则矩阵A必使由上述特征多项式决定的矩阵多项式函数 f(A)=An+a1An-1+…+an-1A+anI=0 上述特征多项式亦称为矩阵A的零化特征多项式。 □ 凯莱-哈密顿定理(2/4) 证明 因为 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为n?n维的常数矩阵。 凯莱-哈密顿定理(3/4) 因此由前面两式,有 (?n+a1?n-1+…+an-1?+an)I=(?n-1I+?n-2B2+…+?Bn-1+Bn)(?I-A) 整理得 (?n+a1?n-1+…+an-1?+an)I =?nI+(B2-A)?n-1+…+(Bn-Bn-1A)?-BnA 凯莱-哈密顿定理(4/4) 上式中,令等号两边?的同幂次项的系数相等,则有 a1I-B2+A=0 a2I-B3+AB2=0 … an-1I-Bn+ABn-1=0 anI+ABn=0 因此,将上述各等式从上至下依次右乘以An-1,…,A,I,然后将各等式相加,即得 An+a1An-1+…+an-1A+anI=0 故矩阵A满足其本身的零化特征多项式。  最小多项式 (1/3) 2. 最小多项式 根据凯莱-哈密尔顿定理,任一n×n维矩阵A满足其自身的特征方程,即特征多项式为A的一个零化多项式。 然而特征多项式不一定是A