n次正方行列の逆行列を部分ピボットLU分解で計算

使い方

1. ページ先頭の『行列A』へ、n次行列の逆行列で計算したい実数行列を貼り付けます。1行を行列の1行として、成分は空白または半角コンマで区切り、各行の列数をそろえます。
2. 行数・列数、負号、小数点、指数表記を確認します。正方行列が必要なページでは行数と列数を一致させ、矩形対応ページでは不要なゼロ行やゼロ列を追加しません。
3. 許容誤差は通常1e-12のまま計算し、特異・低ランク・収束境界を調べる場合だけ理由を記録して変更します。値を変えてエラーを隠す目的では使いません。
4. 計算後は結果表の行列名または指標、行、列、値を確認し、元の式へ代入した残差、既知の小さな例、単位行列やゼロ行列などの境界例と照合します。
5. 必要なら見出し付きTSVをコピーし、元の入力、計算条件、計算日、使用目的と一緒に保存します。重大な用途では独立した数値ライブラリでも再計算してください。

計算式・考え方

正方行列Aを部分ピボット付きで P A=L U と分解し、各標準基底e_jについて L y=P e_j、U x_j=y を解いて列x_jを並べます。A A⁻¹とA⁻¹ Aが単位行列へ戻る残差も検算します。

注意事項

• 逆行列A⁻¹は、A A⁻¹=A⁻¹ A=Iを満たす正方行列です。正方であっても行や列が一次従属なら逆行列は存在しません。本ページは余因子を展開するのではなく、部分ピボット付きLU分解を一度作り、単位行列の各列を右辺として三角方程式を解きます。次数が増えたときも式の意味を保ち、連立方程式の解法と共通の数値基盤で検算できます。
• 例としてA=[[4,7],[2,6]]ではdet(A)=10で、A⁻¹=[[0.6,-0.7],[-0.2,0.4]]です。結果表は行列名A^-1、行番号、列番号、値の順で4成分を返します。小数表示を手作業で丸めてから積を取ると単位行列からずれるため、検算には内部値を使い、表示された有限桁だけを再入力しないでください。
• 特異行列[[1,2],[2,4]]は2行目が1行目の2倍で、異なる右辺に対する一意な対応を作れません。ピボットが許容誤差以下になった時点で逆行列を返さず、入力行列と許容誤差を見直すエラーにします。極めて小さなピボットを無理に割ると巨大な成分が現れ、見かけ上の積だけが不安定になるためです。
• 逆行列を使ってAx=bを解くことは数学的にはx=A⁻¹bですが、数値計算では逆行列を明示的に作らずLU分解から直接解く方が計算量と誤差の面で適切です。連立方程式だけが目的なら専用ページを使ってください。逆行列ページは変換そのもの、複数右辺への再利用、教材上の性質確認に向いています。
• 単位行列の逆行列は同じ単位行列、対角行列の逆行列はゼロでない対角成分の逆数を並べた対角行列です。これらは境界検算に適しています。三角行列でも対角に0があれば特異です。行列式が小さいという事実だけで誤りとは決められず、入力スケール、残差、必要精度を一緒に判断してください。
• FAQでよくある『逆行列が非常に大きい』という結果は、入力ミスだけでなく行列が特異に近い可能性を示します。桁を増やせば必ず解消する問題ではありません。行や列の尺度をそろえ、条件数を扱える専門環境でも確認し、元データのわずかな変化で結果が大きく変わらないか感度を調べてください。
• 入力は実数だけを含む行列です。1行が行列の1行に対応し、成分は空白または半角コンマで区切ります。空行は前後だけを無視し、途中の空行、列数の不一致、数値以外の文字、NaN、Infinity、指数表記として解釈できない断片は暗黙に補正しません。表計算ソフトから貼り付ける場合も、貼り付け後に行数・列数と負号、小数点、桁区切りの有無を確認してください。
• このページは可変サイズを扱いますが、ブラウザーを停止させないよう成分数と文字数を有界にします。大きな疎行列を密行列として貼り付ける用途や、数千行規模の科学計算を置き換えるものではありません。上限を超える行列は途中まで計算せず、入力エラーとして停止します。大規模計算ではLAPACKを利用する科学計算環境など、目的に合う専門ソフトを使用してください。
• 計算はIEEE 754の倍精度相当で行います。おおむね15〜16桁の有効数字を持ちますが、入力した全桁が正しい結果として保存される意味ではありません。参照ページにある最大50桁表示をそのまま再現したとは表示せず、本ツールでは最大有効桁と丸め誤差を明記します。非常に大きい値と小さい値を混在させると有効桁が失われるため、必要なら単位や尺度をそろえてください。
• 部分ピボットは各消去段階で候補列の絶対値が大きい行を選び、ゼロまたは極端に小さい割り算を避ける手順です。多くの一般行列で安定性を改善しますが、すべての悪条件を解消する保証ではありません。行交換が起きた場合は置換情報を計算へ一貫して反映し、行列式の符号、右辺の並べ替え、Lの既計算部分にも同じ交換規約を適用します。
• 許容誤差は、計算上の小さな値をゼロと判定するための尺度です。既定値1e-12は万能な物理定数ではなく、入力行列の最大成分や次数と組み合わせて使います。誤差を大きくすると特異・低ランクと判定しやすくなり、小さくすると丸め雑音を有効な成分と誤認しやすくなります。結果を比較するときは、行列だけでなく許容誤差も同じにしてください。
• 検算では結果の見た目だけでなく、元の式へ戻した残差を確認します。逆行列ならA A⁻¹−I、LU分解ならP A−L U、連立方程式ならA x−bが対象です。行列の大きさが変わると同じ絶対残差でも意味が変わるため、必要に応じて入力と結果のノルムに対する相対的な大きさも見ます。残差が小さいことは有力な確認ですが、悪条件問題の真の解に近いことを単独で保証しません。
• 結果表の行番号と列番号は1から始めます。行列名、行、列、値の各列を持つ表は、見出しを含むTSVとしてコピーできる設計です。画面幅が狭い場合は横スクロールしても列の意味が失われないよう固定した見出しを使います。コピー後に並べ替えると行列の形が変わるため、元の行・列番号を残し、入力、許容誤差、計算日も一緒に記録してください。
• 正方行列が必要な計算へ矩形行列を入れたり、係数行列と右辺の行数が違ったりする場合は実行しません。足りない成分を0で埋める、余分な列を捨てる、右辺を繰り返すといった推測は結果を別問題へ変えてしまうためです。エラーが出たら、各行の成分数、末尾の区切り文字、全角記号、貼り付け時の改行を見直してください。
• 一次資料としてNetlib LAPACKのDGETRF、DGETRI、DGESV、DLANGE、DGESVDとLAPACK Users' Guideの一般連立方程式の説明を2026年7月14日に照合しました。これらは実装方針と検算規約の根拠であり、ブラウザーから資料へ通信して計算する仕組みではありません。資料名と確認日は本文に残し、URLは配布定義ではなく別の監査記録で管理します。
• 入力検査、数値計算、表の描画、コピー用TSVの生成はclient-onlyでブラウザー内だけで完結します。行列や右辺をサーバー、外部API、外部CDNへ送信せず、数値ライブラリを実行時に取得しません。そのため最新データを自動取得する機能もありません。共有端末、画面共有、ブラウザー拡張、クリップボード履歴、貼り付け先の保存設定は利用者側で管理してください。
• 教育目的では、まず2×2や3×3の手計算可能な例で操作し、次に行交換が起きる例、ゼロ行を含む例、ほぼ従属な例を試すと違いを確認できます。通常例だけでなく、1×1、ゼロ行列、単位行列、対角行列、特異行列という境界を分けて扱ってください。エラーを回避するため値を少し変えるのではなく、なぜ計算が定義できないかを式から確認することが重要です。
• 本ツールは入力した数学問題を計算するもので、測定値の品質、モデルの妥当性、単位系、座標系、境界条件を判定しません。構造設計、制御、電気、医療、金融など重大な判断へ使う場合は、独立した実装、専門家のレビュー、既知解、残差、感度解析を組み合わせてください。表示値を根拠に安全率や判定基準を自動で決める機能はありません。
• 倍精度浮動小数点による補助計算です。特異または悪条件の行列、厳密演算、高精度保証、重大な設計判断には単独で使用できません。

よくある質問