[论文解读] Jordan Normal and Rational Normal Form Algorithms

Bernard Parisse, Morgane Vaughan|arXiv (Cornell University)|Dec 2, 2004

Polynomial and algebraic computation参考文献 7被引用 3

一句话总结

本文提出了一种确定性算法，利用Fadeev公式和(λI − A)的余子式矩阵计算矩阵的若尔当代数标准形和有理标准形。该方法利用了B(λ₀)的非零列即为特征值λ₀对应的特征向量这一事实，通过多项式矩阵运算实现高效计算，当特征多项式已因式分解时，其复杂度为O(n⁴)域运算。

ABSTRACT

In this paper, we present a determinist Jordan normal form algorithms based on the Fadeev formula: \[(λ\cdot I-A) \cdot B(λ)=P(λ) \cdot I\] where $B(λ)$ is $(λ\cdot I-A)$'s comatrix and $P(λ)$ is $A$'s characteristic polynomial. This rational Jordan normal form algorithm differs from usual algorithms since it is not based on the Frobenius/Smith normal form but rather on the idea already remarked in Gantmacher that the non-zero column vectors of $B(λ_0)$ are eigenvectors of $A$ associated to $λ_0$ for any root $λ_0$ of the characteristical polynomial. The complexity of the algorithm is $O(n^4)$ field operations if we know the factorization of the characteristic polynomial (or $O(n^5 \ln(n))$ operations for a matrix of integers of fixed size). This algorithm has been implemented using the Maple and Giac/Xcas computer algebra systems.

研究动机与目标

• 开发一种在代数闭域上计算矩阵若尔当代数标准形的确定性算法。
• 将该方法扩展至在任意域上计算有理标准形，无需引入域扩张。
• 通过使用(λI − A)的余子式矩阵，为Frobenius/Smith标准形方法提供一种替代方案。
• 在已知特征多项式因式分解的前提下，实现O(n⁴)的时间复杂度，使其适合并行化处理。
• 将该算法实现在Maple和Giac/Xcas等计算机代数系统中，以支持实际应用。

提出的方法

• 利用Fadeev算法计算特征多项式P(λ)和矩阵多项式B(λ)，即(λI − A)的余子式矩阵。
• 利用恒等式(λI − A)B(λ) = P(λ)I，通过基于迹的公式推导出B(λ)系数的递推关系。
• 应用Horner方法高效计算在特征值λ₀处的B(λ₀)，以提取特征向量。
• 通过求解涉及Q(A)的线性系统来构造广义特征向量，其中Q是P(λ)的不可约因式。
• 利用二项式系数恒等式，通过递推式(9)从基向量重构高阶循环向量。
• 在Maple和Giac/Xcas中实现了一种混合方法，结合复数若尔当代数标准形的计算与对次数大于1的不可约因式的有理标准形转换。

实验结果

研究问题

• RQ1 能否仅通过(λI − A)的余子式矩阵和特征值求值，确定性地计算若尔当代数标准形？
• RQ2 如何仅通过矩阵多项式运算，在不引入域扩张的前提下构造有理标准形？
• RQ3 与基于Frobenius/Smith标准形的方法相比，基于余子式的算法的计算复杂度如何？
• RQ4 该算法是否能够高效并行化，或在计算机代数系统中以最小符号计算开销实现？
• RQ5 如何系统地通过多项式矩阵运算，从基特征向量重构广义特征向量？

主要发现

• 当特征多项式完全因式分解时，该算法以O(n⁴)次域运算完成若尔当代数标准形的计算。
• 对于固定大小的整数矩阵，由于多项式因式分解的代价，复杂度上升至O(n⁵ ln n)。
• 有理标准形通过先计算伪有理形式，再经由结构化变换转换得到。
• 该方法避免了对Frobenius或Smith标准形的依赖，为有理标准形的计算提供了新途径。
• 该算法已在Maple中成功实现，并原生集成于Giac/Xcas中，用户可直接调用jordan和rat_jordan等函数。
• 一个关键洞察是：B(λ₀)的非零列即为特征值λ₀对应的特征向量，从而可直接提取特征向量，无需求解(A − λ₀I)v = 0。