[論文レビュー] A New Algorithm for Inverting General Cyclic Heptadiagonal Matrices Recursively

A. A. Karawia|arXiv (Cornell University)|Nov 10, 2010

Scientific Research and Discoveries参考文献 12被引用数 3

ひとこと要約

本稿では、行列要素に制限的な仮定を課さずに一般の循環型ヘプタダイアゴナル行列を逆行列化するための記号的で再帰的なアルゴリズム—CHINV—を提示する。DoolittleのLU分解、並列計算、再帰的関係式を用いて、逆行列を列ごとに計算し、最終的な逆行列は記号的簡略化とパラメータの代入によって得られる。分母がゼロになる可能性がある場合の処理のため、パラメータを導入している。

ABSTRACT

In this paper, we describe a reliable symbolic computational algorithm for inverting general cyclic heptadiagonal matrices by using parallel computing along with recursion. The algorithm is implementable to the Computer Algebra System(CAS) such as MAPLE, Matlab and Mathematica . An example is presented for the sake of illustration.

研究の動機と目的

• 行列要素に制限的な仮定を課さずに、一般の循環型ヘプタダイアゴナル行列の信頼性の高い記号的逆行列化アルゴリズムの開発。
• 既存のピアノダイアゴナルおよびヘプタダイアゴナル行列用のアルゴリズムを、完全に一般化された循環型ヘプタダイアゴナルケースに拡張すること。
• 記号的計算を用いることで除算ゼロエラーを回避し、数値的安定性を確保すること。
• 循環型ヘプタダイアゴナル行列を含む線形方程式系を効率的に解く再帰的メソッドの提供。
• ピアノダイアゴナル行列用の既存の再帰的アルゴリズムをヘプタダイアゴナルケースに一般化すること。

提案手法

• DoolittleのLU分解を用いて、循環型ヘプタダイアゴナル行列Hを、要素に記号的表現を用いたLおよびU行列に分解する。
• LU分解と記号的簡略化から導出された再帰的関係式を用いて、逆行列H⁻¹の最後の4列を計算する。
• 逆行列計算における分母がゼロになる可能性を処理するため、記号的パラメータtを導入し、後にt = 0に代入することで最終結果を得る。
• 再帰的列関係式（式2.39–2.40）を適用して、既知の最後の4列から残りの(n−5)列をH⁻¹に計算する。
• MAPLE、MATLAB、MATHEMATICAなどのコンピュータ代数システム（CAS）でアルゴリズムを実装し、記号的実行を可能にする。
• 逆行列計算の前に、非特異性を検証するDETCHERPTAサブアルゴリズムを適用する。

実験結果

リサーチクエスチョン

• RQ1 一般の循環型ヘプタダイアゴナル行列は、その要素に制限的な条件を課さずに逆行列化可能か？
• RQ2 記号的計算をどのように活用することで、行列逆行列化アルゴリズムの数値的信頼性を確保できるか？
• RQ3 標準的なピアノダイアゴナルケースを超えて、循環型ヘプタダイアゴナル行列の逆行列を効率的に計算する再帰的構造は何か？
• RQ4 LU分解と再帰的関係式を用いて、逆行列を列ごとにどのように計算できるか？
• RQ5 除算ゼロエラーを回避するための記号的パラメータ化戦略は何か？

主な発見

• CHINVアルゴリズムは、n ≥ 8の一般の循環型ヘプタダイアゴナル行列の逆行列を、行列要素に制限的な仮定を課さずに正しく計算できる。
• パラメータtを用いた記号的パラメータ化により、アルゴリズムは頑健で、失敗が生じない。最終的な逆行列を得るために、後でt = 0に代入する。
• 逆行列は列ごとに計算される：最後の4列は再帰的関係式（式2.15–2.38）により導出され、残りの列は式2.39–2.40により計算される。
• 本手法は、ピアノダイアゴナル行列用の既存の再帰的アルゴリズムを一般化し、ヘプタダイアゴナル系への適用可能性を拡張している。
• n=10の具体例において、与えられた線形方程式系に対して、解ベクトルX = [1,2,3,4,5,6,7,8,9,10] を正しく回復している。
• MAPLE、MATLAB、MATHEMATICAなどのCAS環境で本アルゴリズムを実装可能であり、正確な記号的計算が可能である。