擬似逆行列計算ツールの使い方
擬似逆行列計算ツールは、任意のフルの実数行列に対するムーア・ペンローズ擬似逆行列 A⁺ をオンラインで即座に計算します。行列のサイズ（行数・列数）と各成分を入力すると、本ツールが最適な計算公式を自動選択し、計算手順やランクとともに結果を表示します。
- サイズ（次元）を設定 – 行数と列数（それぞれ 1〜5）を指定します。
- 行列 A の要素を入力 – 各マスに値を入力します。負の数や小数も入力可能です。
- 擬似逆行列 A⁺ を確認 – 擬似逆行列計算ツールが算出された擬似逆行列と採用された計算手法を表示します。
- 計算手順を確認 – 中間行列（AᵀA や (AᵀA)⁻¹ など）の段階的な計算プロセスを確認できます。
- 方程式系への適用 – 線形方程式系 Ax = b の最小二乗解 x = A⁺b の算出に活用できます。
擬似逆行列の計算公式と理論
擬似逆行列計算ツールでは、ムーア・ペンローズ擬似逆行列の算出において、行列の性質に応じた以下の閉形式（Closed-form）公式を実装しています。
ケース 1: 正則な正方行列（フルランク）
A⁺ = A⁻¹
ケース 2: フル列ランク（rank = n ≤ m）
A⁺ = (AᵀA)⁻¹ Aᵀ
ケース 3: フル行ランク（rank = m ≤ n）
A⁺ = Aᵀ (AAᵀ)⁻¹
ケース 4: ランク落ち・階数不足（フルランクでない場合）
A⁺ ≈ (AᵀA + εI)⁻¹ Aᵀ （ティホノフ正則化）| 記号 | 意味・定義 |
|---|---|
| A⁺ | ムーア・ペンローズ擬似逆行列 |
| AᵀA | グラム行列（列数 × 列数） |
| AAᵀ | 双対グラム行列（行数 × 行数） |
| ε | 正則化用の微小定数（10⁻⁹） |
| r | rank(A) （行列 A のランク・階数） |
擬似逆行列は以下の ペンローズの4条件 をすべて満たします：
- AA⁺A = A
- A⁺AA⁺ = A⁺
- (AA⁺)ᵀ = AA⁺
- (A⁺A)ᵀ = A⁺A
方程式 Ax = b において、最小二乗誤差 ||Ax − b||₂ を最小化し、かつその中でノルム ||x||₂ を最小にする解は x = A⁺b として一義的に求められます。
前提条件と計算精度
本擬似逆行列計算ツールは最大 5×5 までの行列をサポートしています。ランクの判定は、閾値 10⁻¹⁰ を用いたガウスの消去法（掃き出し法）によって行われます。ランク落ちのケースでは、AᵀA が可逆になるよう微小な正則化項（ε = 10⁻⁹）を加算しています。これは近似計算であり、特異値分解（SVD）による厳密な擬似逆行列の計算結果と僅かに異なる場合があります。
擬似逆行列計算ツールの活用事例
擬似逆行列は、応用数学、データサイエンス、工学分野で幅広く利用されています。
- 最小二乗回帰（データフィッティング） – 方程式の数が未知数より多い過剰決定系 Ax = b の最適解を x = A⁺b で導出。
- 欠損・不足決定系の方程式解法 – 未知数が方程式より多い場合に、ノルムが最小となる唯一の解を取得。
- 機械学習・統計モデル – 条件数が悪く正則行列でない計画行列の重みパラメータ計算。
- 現代制御理論 – 可制御性行列の擬似逆行列を用いた最小エネルギー制御入力の設計。
- 画像処理・信号処理 – 逆問題やデコンボリューション（逆畳み込み）演算における復元。
本ツールで得られた擬似逆行列 A⁺ を方程式 x = A⁺b に代入することで、様々な線形システムの計算をスムーズに行うことができます。