[論文レビュー] Condensation of Determinants

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

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

ひとこと要約

本稿では、n×n行列の行列式を2×2小行列式を要素とする(n−1)×(n−1)行列の行列式に縮約する、新しい行列式縮約法を提示する。主な貢献は、元の行列式を(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行列の行列式を、2×2小行列式を要素とする(n−1)×(n−1)行列の行列式に縮約する手法を開発すること。
• ドジソンの縮約法を一般化し、完全な余因子展開を回避する再帰的アルゴリズムを構築すること。
• 元の行列式と2×2小行列式から構成された縮約行列の行列式を結びつける数学的に厳密な公式を提供すること。ただし、ピボット要素がゼロの場合でも正しく保たれるようにすること。
• 任意の正方行列の行列式を、反復的縮約によって簡素化されたアルゴリズム的手法で計算すること。

提案手法

• 本手法は、元の行列Aの最初の行と他の行から形成された2×2行列式を要素とする(n−1)×(n−1)行列Bを構築する。
• 最初の行における最初の非ゼロ要素の位置に基づいてBの要素を定義し、数値的安定性を確保するとともに、ゼロ除算を回避する。
• 元の行列式は、det(A) = det(B) / (a₁ₗ)ⁿ⁻²として回復され、ここでa₁ₗは最初の行における最初の非ゼロ要素である。
• 行列式の等価性を維持するために、ブロック行列の分解と行・列インデックスに基づく符号調整を用いる。
• 縮約プロセスを再帰的に適用し、2×2行列に達するまで繰り返す。2×2行列は直接計算される。
• ピボット要素がゼロの場合を扱うための補題を証明し、すべての状況下でも手法の有効性を保証する。

実験結果

リサーチクエスチョン

• RQ1 n×n行列の行列式を、2×2小行列式を要素とする(n−1)×(n−1)行列の行列式に、2×2小行列式のみを用いて体系的かつ一貫して縮約できるか？
• RQ2 元の行列式と2×2小行列式から構成された縮約行列の行列式を結びつける正確な代数的公式は何か？
• RQ3 ピボット要素（例：a₁₁）がゼロの場合、どのようにして縮約プロセスを強靭に保つか？
• RQ4 再帰的アルゴリズム実装に適した、縮約プロセスの再帰的構造は何か？
• RQ5 この手法は、完全な余因子展開を回避しつつも行列式不変性を保つように一般化可能か？

主な発見

• 本稿では、(a₁ₗ)ⁿ⁻² det(A) = det(B) という恒等式を確立した。ここでBは、最初の行と他の行から形成された2×2行列式を要素とする(n−1)×(n−1)行列である。
• ピボット要素がゼロの場合でも、a₁₁ = 0 ならば、その小行列式の行列式がゼロになることを示す補題により、一貫性が保たれる。
• n = 7の場合、元の行列式の32倍が、2×2小行列式から構成された6×6行列の行列式に等しいことを示した。
• アルゴリズム的実装により、完全な余因子展開を回避し、再帰的縮約が可能になり、計算量が低減された。
• 適切に符号変更を加えることで、行および列の置換に対しても、行列式不変性が保証されることを証明した。
• 構造的にドジソンの方法に類似しているが、最初の行からピボットを選択可能に一般化されており、数値的安定性が向上している。