[论文解读] Condensation of Determinants

Abdelmalek Salem, Saïd Kouachi|ArXiv.org|Dec 5, 2007

Matrix Theory and Algorithms参考文献 1被引用 77

一句话总结

本文提出了一种新颖的矩阵行列式约化方法，通过将 n×n 矩阵的行列式简化为 (n−1)×(n−1) 矩阵的行列式，其中新矩阵的元素为 2×2 子式。关键贡献在于推导出一个公式，将原始行列式乘以 (a₁ₗ)ⁿ⁻² 后等于由第一行与其他行构成的 2×2 子式作为元素的矩阵的行列式，从而实现复杂度更低的递归计算。

ABSTRACT

In this paper we tried to condense the determinant of n square matrix to the determinant of (n - 1) square matrix with the mathematical proof.

研究动机与目标

• 开发一种将 n×n 矩阵的行列式约化为 (n−1)×(n−1) 矩阵行列式的方法，仅使用 2×2 子式作为元素。
• 通过构建避免完整余子式展开的递归算法，推广 Dodgson 的约化方法。
• 提供一个数学上严谨的行列式约化公式，即使在主元元素为零时也能保持正确性。
• 通过迭代约化提供一种简化且算法化的任意方阵行列式计算方法。

提出的方法

• 该方法构建一个 (n−1)×(n−1) 矩阵 B，其元素 bᵢ,ⱼ 为原始矩阵 A 的第一行与其他行构成的 2×2 行列式。
• 基于第一行中首个非零元素的位置定义 B 的元素，以确保数值稳定性并避免除零错误。
• 原始矩阵的行列式通过公式 det(A) = det(B) / (a₁ₗ)ⁿ⁻² 恢复，其中 a₁ₗ 为第一行中首个非零元素。
• 该方法利用分块矩阵分解，并根据行和列索引调整符号，以保持行列式等价性。
• 递归应用约化过程，直至得到 2×2 矩阵，再直接计算其行列式。
• 证明了一个引理，以处理主元元素为零的情况，确保方法在所有条件下均有效。

实验结果

研究问题

• RQ1 n×n 矩阵的行列式能否通过仅使用 2×2 子式作为元素，系统地约化为 (n−1)×(n−1) 矩阵的行列式？
• RQ2 原始行列式与由 2×2 子式构成的约化矩阵行列式之间的精确代数公式是什么？
• RQ3 当主元元素（如 a₁₁）为零时，如何使约化过程保持鲁棒性？
• RQ4 约化过程的递归结构是什么，使其能够高效地实现算法化？
• RQ5 该方法能否被推广，以避免完整余子式展开，同时保持行列式不变性？

主要发现

• 本文建立了恒等式 (a₁ₗ)ⁿ⁻² det(A) = det(B)，其中 B 是一个 (n−1)×(n−1) 矩阵，其元素为第一行与其他行构成的 2×2 行列式。
• 该方法通过引理成功处理了主元元素为零的情况，证明若 a₁₁ = 0，则相应子式矩阵的行列式为零，从而保持一致性。
• 对于 n = 7，本文证明 32 倍的原始行列式等于由 2×2 子式构成的 6×6 矩阵的行列式。
• 算法实现通过避免完整余子式展开，实现了计算复杂度的降低，并支持递归约化。
• 该方法在正确考虑子式行列式符号变化的前提下，被证明在行和列置换下保持不变性。
• 约化过程在结构上等价于 Dodgson 的方法，但被推广为允许从第一行选择主元，从而提升了数值稳定性。