将矩阵化为约当标准型.docxVIP

增补： sIA 1 ~ 1 1 PsI ~1 P 1 sIPAP A adjsI ~ = A P 1 P ~ detsI A = =  (s 1)(s 2) (s 1)(s2) P (s 1)2(s (s1)2 P1 2) s 2 s 2 P (s1)(s (s 1) P1 2) 矩阵 A的特点多项式det(IA) ( 1)2(2)有双重根 12 1和单特点值3 2，凡是在adjsI A 中的公因子则必定和 detsI A能够相消。 经过线性变换后，系统矩阵成为对角线矩阵形式的状态空间表达式， 特别指出，假如nn维矩阵A由下式给出 0 1 0 0 1 A 0 0 0 1 0 1 2 n1 而且其特点值 1,2 , n互异，作非奇怪线性变换x ~ ，则化A为对 Px 角线标准型矩阵 1 0 0 0 2 0 P1AP 0 0 0 n 此中，P为范德蒙德（Vandermond）矩阵。即 1 1 1 1 P 1 2 3 n n1 n1 n1 n1 1 2 3 n 增补： 设约当块数为q和q个mi（约当块的阶数）。A矩阵唯一决定的约当型矩 阵式 J1 J2 J Jq 设变换矩阵p与J拥有相同阶数组的分块矩阵型 令P [p1 p2 pq] 即，pi是n mi 阶矩阵。则AP PJ J1 A[p1 p2 pq][p1p2 J2 pq] Jq 依据分块矩阵的乘法例则，有 [Ap1  Ap2  Apq]  [p1J1  p2J2  pqJq] 上式其实是  q个等式，即  Api  piJi,  i  1,2,  ,q 将n  mi  阶矩阵  pi  写成列向量形式，于是有 Pi  [pi1  pi2  pimq] i  1 i  1 A[pi1pi2pimi][pi1pi2pimi] 1 i 即 Api1 ipi1 Api2 pi1ipi2 Apimipimi1ipimi 也可写成 (iI A)pi1 0 (iI A)pi2 pi1 (iI A)pimi pimi1 次序解以上方程组就能够确立 pi的mi个列向量。这些列向量中只有 第一个pi1是对应于 i的特点向量，而其他的mi1个向量pi2,,pimi,称 之为对应特点值i的广义特点向量，可由上式递推解出。 设矩阵A的重特点值为 1，代入式(iI A)pi10中，即由 (1IA)p11 0 可求出A的对应于1的特点向量。有上式解出的线性独立特点向量的个数，就 是该特点值对应的约当块数，或表示为 11nrank1IA 降秩数11就是对应1的线性没关特点向量个数，或许是对应1的约当块块数。 换句话说，矩阵A的特点值分组1,2,q中，有12a11。 将式 上当算p12 的式子， (1IA)p12 p11 两头同时乘以(1I A)，(1I A) 2p12 (1I A)p110 该方程线性没关的解的个数是 n rank1I A 2 ，但这个数量中包含 11 的个数， p 即11。因此，解出线性没关的列向量p12的个数， 12rank1IArank1IA2 也就是对应1的大于或等于2阶约当块的块数。 例将已知矩阵A化为约当型。 2i10 021 002 A 1 3 2 21 2 解：先求A的特点多项式，由于A矩阵是对角分块矩阵，因此特点多 项式是每个对角分块矩阵特点多项式的乘积，即 (sIA)(sI3A1)(sI2A2)(sI1A3)(s2)6 将特点值 2代入式11n rank 1IA，求约当型中的约当块数。 0 1 0 0 0 1 0 0 0 11n rank2IA6 rank 1 1 =6-3=3 2 2 1 0 由此，由A矩阵化成的约当型共有 3个约当块。 而后，将 2代入12 rank1I Arank 1I A2求出约当型中大于等于2 阶的块数 0 1 0 0 0 0 123rank 2I A2 3 rank 0 0 0 0 0 3 12 0 0 0 因此，由A矩阵化成的约当型将有一个 1阶块，两个大于或等于2阶的块。 再将将 2代入13 rank 1I A 2 rank 1IA 3，求出约当型中大于等 于3阶的块数 13 1 rank 2 A 3 1 0 1 I 因此，由A矩阵化成的约当型共有 3个约当块，此中一个1阶块，一个2阶块， 一个3阶块，即 2 1 0 0 2 1 J P1AP 0 0 2 1 2 0 2 2 设i是系统的一个特点值，若存在一个 n维非零向量pi，知足 Api i pi 或 (iI A)pi 0 则称pi为系统相关于特点值 i的特点向量。 比如：系统矩阵为 A 0 1 2 3 试求其特点值和一组特点向量。 解：由系统的特点方程 1 2 3 20 IA 3 2 系统的特点值为1 1,2 2 设关于特点值1,2的特点向量p