[论文解读] An Explicit Construction of Gauss-Jordan Elimination Matrix

Yi Li|ArXiv.org|Jul 29, 2009

Spectroscopy and Chemometric Analyses参考文献 6被引用 3

一句话总结

本文提出了一种基于行列式的显式构造方法，用于高斯-约旦消去法中的元素，表明在经过 $ k $ 步后，矩阵中每个非对角线元素 $ a^k_{ij} $（其中 $ i \neq j, j > k $）可以表示为原始矩阵 $ A $ 的两个行列式的比值。该方法提供了一种适用于任意线性系统的广义克拉默法则，包括含参数系数的系统，并可直接用于构造 $ AX = b $ 的解。

ABSTRACT

A constructive approach to get the reduced row echelon form of a given matrix $A$ is presented. It has been shown that after the $k$th step of the Gauss-Jordan procedure, each entry $a^k_{ij}(i<>j; j > k)$ in the new matrix $A^k$ can always be expressed as a ratio of two determinants whose entries are from the original matrix $A$. The new method also gives a more general generalization of Cramer's rule than existing methods.

研究动机与目标

• 开发一种利用行列式计算矩阵行最简形的构造性方法。
• 将克拉默法则推广至现有公式之外的范围，特别是适用于含参数系数的系统。
• 提供一种直接通过行列式表达式求解 $ AX = b $ 的计算路径。
• 建立一个理论框架，使得高斯-约旦消去法的每一步都能通过原始矩阵的行列式显式表示。

提出的方法

• 该方法将经过 $ k $ 次高斯-约旦消去步骤后的矩阵中每个元素 $ a^k_{ij} $ 表示为从原始矩阵 $ A $ 衍生出的两个行列式的比值。
• 利用基于巴雷伊斯（Bareiss）的无分式高斯消去法的递归行列式恒等式，推导出每个消去阶段的元素表达式。
• 该方法根据行索引 $ i $ 的不同分为三种情况：$ i > k $、$ i = k $ 和 $ i < k $，每种情况具有不同的行列式公式。
• 对于 $ i < k $ 的情况，表达式中引入符号因子 $ (-1)^{k-i+1} $，以考虑行交换和行列式展开的影响。
• 推导基于恒等式 $ \left| \begin{smallmatrix} M & U \\ R & a \end{smallmatrix} \right| = \text{det}(M) \cdot \text{det}(\text{子矩阵}) $，实现行列式的递归分解。
• 该方法确保 $ A^k $ 中的所有元素均表示为原始矩阵元素的有理函数，分母为顺序主子式。

实验结果

研究问题

• RQ1 在高斯-约旦消去过程中，矩阵的元素能否通过原始矩阵的行列式显式表达？
• RQ2 克拉默法则如何推广以适用于含参数系数或非方阵的系统？
• RQ3 是否存在一个统一的基于行列式的公式，适用于所有满足 $ j > k $ 的非对角线元素 $ a^k_{ij} $，而不受行位置影响？
• RQ4 该基于行列式的思路能否直接得出 $ AX = b $ 的解，而无需回代？

主要发现

• 经过 $ k $ 次高斯-约旦消去步骤后，每个满足 $ i \neq j $ 且 $ j > k $ 的元素 $ a^k_{ij} $ 均可表示为原始矩阵 $ A $ 的两个行列式的比值。
• 元素 $ a^k_{ij} $ 的表达式依赖于相对行索引 $ i $：当 $ i > k $ 时，其表达式为 $ \frac{a^{(k)}_{ij}}{a^{(k-1)}_{kk}} $；当 $ i = k $ 时，为 $ \frac{a^{(k-1)}_{ij}}{a^{(k-1)}_{kk}} $；当 $ i < k $ 时，包含符号因子 $ (-1)^{k-i+1} $。
• 该方法以区别于以往方法的方式推广了克拉默法则，使得 $ AX = b $ 的解可直接构造。
• 该方法在任意交换环上均成立，只要主子式非零，从而将适用范围扩展至实数或复数矩阵之外。
• 行列式表达式的递归结构与巴雷伊斯的无分式消去法一致，确保精确算术运算，避免浮点误差。
• 该方法支持参数化线性系统和整数规划的应用，通过基于行列式的解进行多项式插值。