形式陈述
不知道特征值,仍能完整判断相似吗
设 A ∈ F n × n ,n ≥ 1 。存在唯一一列首一、次数为正的多项式
f 1 ∣ f 2 ∣ ⋯ ∣ f r , ∑ i deg f i = n , 使 A 经换基 理路 换基与坐标变换 Change of basis · Coordinate transformation 用可逆过渡矩阵在不同基之间转换向量坐标与算子矩阵。 相似于
R = diag ( C ( f 1 ) , … , C ( f r ) ) , 其中 C ( f ) 是列约定的伴随矩阵 理路 循环子空间与伴随矩阵 Cyclic subspace and companion matrix · Cyclic vector · Polynomial companion matrix 把一个向量的精确幂递推表示为多项式商模和伴随矩阵,给出列向量约定下的换基证书。 。f i 称为不变因子,R 称为有理标准形。“有理”表示构造只需底域运算,并不要求 F = Q 。
两个同阶矩阵在 F 上相似,当且仅当这列不变因子相同。若 P 的列是上述循环块基,则
A P = P R , det P ≠ 0 是可直接复核的完整证书。只列出 R 而没有给出可逆换基,并不等于已经验证某个具体输入的分类。
模结构怎样变成矩阵分类
在 F n 上令 x ⋅ v = A v ,得到 F [ x ] -模 V A 。标准向量空间基也生成这个模;极小多项式消去全部向量,所以它是有限生成挠模。应用PID 模结构定理 理路 PID 上有限生成模结构定理 Structure theorem for finitely generated modules over a PID PID 上每个有限生成模唯一分解为有限秩自由部分与满足整除链的循环挠模。 ,得到
V A ≅ ⨁ i F [ x ] / ( f i ) . 每个商模的幂基成为一个伴随块,组合各块基便给出 P 。没有自由 F [ x ] 部分,因为一个非零自由副本作为 F -向量空间是无限维的。
一个 F -线性同构 P : V R → V A 是模同构,恰好要求它与 x 的作用相容,即 A P = P R 。因此模同构与矩阵相似完全对应,模结构定理的唯一性才真正成为相似类的唯一性。
各块的特征多项式 理路 特征多项式 Characteristic polynomial 由 det(λI−T) 定义并编码线性算子特征值的多项式。 为 f i ,而消去全部块的最低首一多项式是它们的最小公倍式,故
χ A = ∏ i f i , m A = f r . 首尾信息不能一般恢复中间分配,这正是必须保留整列不变因子的原因。
从输入矩阵怎样计算不变因子
矩阵 x I − A 是 V A 的一个关系矩阵。具体地,把多项式列 g ( x ) = ∑ k x k v k 送到 ∑ k A k v k ;映射满射,并把 ( x I − A ) h 送到零。反过来,逐项使用 x v ≡ A v ( mod ( x I − A ) F [ x ] n ) ,每个 g 都等价于它的求值常向量,核因此恰为这组关系。
对 x I − A 计算Smith 形 理路 PID 上的 Smith 正规形 Smith normal form over a PID · Smith normal form PID 上的矩阵可经可逆行列变换化为满足整除链的对角形,且对角因子在相伴意义下唯一。 ,删除单位对角元后,剩余首一因子就是 f i 。这里的左右变换是多项式矩阵等价 ,不能把某一个 Smith 左因子直接当作常矩阵 P 。恢复 P 要把商模生成元映回 V A ,再展开各自的幂轨道;也可以直接提出循环向量后检查 A P = P R 。
直觉
Jordan 形按根和链长编码;有理标准形按底域上的多项式关系编码。两者都保留完整相似类,但前者需要相应多项式在底域分裂,后者不用。
每个循环块内部是一条幂轨道。不变因子的整除顺序不是装饰:任意拼起的伴随块虽可描述同一算子,却未必已成为唯一标准形。例如互素的两块可以通过中国剩余结构合并为一块。
图片加载失败 两条轨道对应二阶与三阶伴随块。图中每支箭头表示再作用一次 A,灰色框是已有向量的线性组合,不增加基维数。
例子与边界
五阶有理输入的完整证书
取
A = ( 0 − 1 0 1 1 1 0 0 0 − 1 0 0 0 1 3 0 0 1 1 − 1 0 0 0 1 1 ) . 令 u = e 1 、v = e 1 + e 3 。直接计算
A u = e 2 , A 2 u = − u ; A v = e 2 + e 4 , A 2 v = e 3 + e 4 + e 5 , 以及
A 3 v = 2 v − A v + 2 A 2 v . 因此两条关系分别是 p = x 2 + 1 和 q = x 3 − 2 x 2 + x − 2 = ( x − 2 ) p 。把轨道排成列,得
P = [ u , A u , v , A v , A 2 v ] = ( 1 0 1 0 0 0 1 0 1 0 0 0 1 0 1 0 0 0 1 1 0 0 0 0 1 ) . P 上三角且 det P = 1 ,所以两条轨道独立并覆盖全空间。前述五个作用关系逐列给出
A P = P diag ( ( 0 − 1 1 0 ) , ( 0 0 2 1 0 − 1 0 1 2 ) ) . 由于 p ∣ q ,这已经是有理标准形,而非任意块分解。不变因子为 ( p , q ) ,χ A = p 2 ( x − 2 ) ,m A = p ( x − 2 ) 。x I − A 的 Smith 对角元则为 ( 1 , 1 , 1 , p , q ) 。
可再独立核对 A 2 + I 的秩为 1 ,且 ( A − 2 I ) ( A 2 + I ) = 0 ;这与二次主分量四维、特征值 2 的主分量一维一致。但这项低阶校验不能替代已经给出的换基证书。
相同的两个多项式为什么仍可能不够
C ( x 2 ) ⊕ C ( x 2 ) 与 C ( x ) ⊕ C ( x ) ⊕ C ( x 2 ) 都有 χ = x 4 、m = x 2 ,不变因子却分别为 ( x 2 , x 2 ) 与 ( x , x , x 2 ) ,所以不相似。这与Jordan 页 理路 Jordan 标准形 Jordan canonical form · Jordan normal form 在特征多项式分裂时,把有限维算子表示成 Jordan 块直和。 的核维数反例是同一结构,换了一套不必分裂的语言。
推论与应用
A 有循环向量,当且仅当不变因子只有一个,也等价于 m A = χ A :整除链中若有其他正次数因子,乘积的次数必大于最后一个。此判别是整个算子的性质,不能由某个未充分探索的起始向量判断。
有理标准形适合精确分类,浮点矩阵的小扰动会改变精确因式关系。需要数值全谱时复用Schur 分解 理路 Schur 分解 Schur decomposition · Schur triangularization · Real Schur form 用酉或正交相似变换把一般矩阵化为上三角或实准上三角形,作为浮点全谱计算的稳定结构目标。 ;需要分类证明时则保留底域、首一约定、块次序及换基方向,不把小浮点残差当成精确相似关系。
参考资料