若当标准形计算.pdfVIP

第22 卷第5 期 大 学 数 学 V ol. 22, . 5 2006 年10 月 COLLEGE MAT HEMAT ICS Oct. 2006 若当标准形的计算 李尚志 ( 北京航空航天大学理学院, 北京 100083) 1 在线性代数与矩阵论中, 一个重要 题是如何将复数域C 上的任意n 阶方阵A 相似到尽可能简单 的方阵J . 也就是说: 寻找适当的可逆方阵P, 使P- 1AP = J. 在复数域上求出A 的特征多项式A ( ) = det( I - A) 的全部不同的根 1 , , t , 也就是A 的特征 值. 对每个特征值i , 解齐次线性方程组(A - iI ) X = 0, 可求得一个基础解系Li . 所有这些Li 的并集L = L1 ! !Lt 线性无关, 是由最多的特征向量组成的线性无关集合. 如果L 包含n 个向量, 它就是复数 域C 上n 维列向量空间C n∀ 1 一组基. 以L 中的向量为各列排成的矩阵P 是n 阶可逆方阵, 它将A 相似 到对角阵P - 1AP. 但是, 也可能按以上步骤求得的L 所含向量个数少于n, 此时A 不能相似于对角阵. 但可以对L 再 补充适当的向量组成C n ∀ 1 的一组基M , 以这组基的各向量为各列排成可逆矩阵P, 它将A 相似到如下 形式的准对角阵 i 1 J 1 - 1 P AP= J = , 其中每个块具有形式Ji = 1 J s i 这样的块J i 称为若当块, 这样的矩阵J 称为A 的若当标准形, P 称为过渡矩阵. 很自然提出的 题是: 怎样由A 计算它的若当标准形J , 并计算过渡矩阵P? 很多教科书都叙述了利用 矩阵I - A 的相抵标准形来计算A 的若当标准形J 的方法. 至于求过 渡矩阵P, 大多数高等代数、线性代数教材中都认为太难而不能加以讨论, 一些矩阵论的教材也只是给 出少量特殊例子而没有给出一般的算法. 有少数教材( 例如[ 2] ) 给出了利用 矩阵的相抵来求过渡矩阵 P 的算法: 求可逆 矩阵P ( ) , Q( ) 使P( ) ( I - A) Q( ) = I - J , 并求出 - 1 m m- 1 R ( ) = Q( ) = Rm Rm- 1 R 1 R0. 再将其中的替换成A 就得到过渡矩阵 P = R(A) = RmA m Rm- 1Am - 1 R 1A R0 . 但是, 对以多项式为元素的矩阵进行相抵变换, 与复系数矩阵的相抵变换比起来, 计算量要大得多, 要 实际用来对矩阵A 进行计算仍然过于繁琐. 本文作者在所编写的教材[ 1] 中给出了由任意复方阵A 计算若当标准形J 和过渡矩阵P 的一个新 算法, 在求得A 的特征根之后只要解齐次线性方程组和求极大线性无关组就可以求出J 和P . 本文通过