擬似逆行列計算機
フルランクの長方形行列または正方行列のムーア・ペンローズ擬似逆行列を計算します。
行列の擬似逆
行をセミコロンで、値をカンマまたは空白で区切って入力してください。
例:1, 2; 3, 4 は二行二列の行列を表します。
擬似逆行列計算機について
ムーア・ペンローズ擬似逆行列は、逆行列の概念を長方形行列や通常の一意な逆行列を持たない系へ拡張したものです。A+ と表し、A、A+、それらの転置に関する四つの整合条件を満たす唯一の行列です。可逆な正方行列では、擬似逆行列は通常の逆行列と一致します。縦長または横長のフルランク行列では、最小二乗計算や最小ノルム計算に必要な、逆に相当する自然な変換を与えます。
この計算機は行ごとに記述した行列を受け付けます。成分はカンマまたは空白で、行はセミコロンで区切ってください。行数が列数以上で列が独立な場合は、(A の転置 A) の逆行列に A の転置を掛けます。列数が行数より多く行が独立な場合は、A の転置に (A A の転置) の逆行列を掛けます。いずれもムーア・ペンローズ逆行列の標準的なフルランク公式です。必要な正方行列の逆行列には、ピボット選択付きガウス・ジョルダン消去法を使います。
代表的な用途は、未知数より方程式が多く、すべての式を厳密に満たすベクトルがない過剰決定系 A x = b です。b に左から A+ を掛けると、残差の二乗和を最小にする最小二乗解が得られます。行フルランクの劣決定系では、A+ はユークリッドノルムが最小の厳密解を与えます。この性質により、回帰、信号処理、制御系、ロボット工学、数値最適化、データ近似、機械学習で役立ちます。
ランクは重要です。便利な正規方程式の公式には、列または行の独立性が必要です。ランク不足では A の転置 A または A A の転置が特異となるため、この実装は不安定な値を返さず、制限を報告します。一般的な特異値分解なら、数値許容値を超える特異値のみ逆数にすることでランク不足を扱えます。大規模、悪条件、ランク不足の問題には、SVD ベースの専用数値ライブラリを使い、許容値の設定を確認してください。
浮動小数点演算も結果に影響します。表示成分は小数点以下六桁に丸めますが、内部計算では JavaScript の倍精度を保持します。正規方程式では条件数が二乗になるため、桁の大きさの差が大きいと丸め誤差が増幅されます。列のスケール調整が役立つこともありますが、専門的な数値計算では QR 法や SVD 法を優先してください。授業の例や小さなフルランク行列なら、この計算機で何度もの行列積や消去を手計算せず、擬似逆行列を分かりやすく求めて確認できます。
擬似逆行列の計算例
| 行列 A | 擬似逆行列 A+ | 意味 |
|---|---|---|
| [1, 0]; [0, 1] | [1, 0]; [0, 1] | 単位行列の逆行列と擬似逆行列は、それ自身です。 |
| [1]; [2] | [0.2, 0.4] | 列ベクトルを、その長さの二乗である 5 で割ります。 |
| [1, 0, 0]; [0, 1, 0] | [1, 0]; [0, 1]; [0, 0] | 横長の行列から最小ノルムの右逆行列が得られます。 |
擬似逆行列の計算方法
- 各行の成分をカンマまたは空白で区切って記述します。
- 行同士をセミコロンで区切り、各行の成分数が等しいことを確認します。
- 「擬似逆行列を計算」を選び、適切なフルランク公式を適用します。
- 出力の各行を読み取り、後続の計算で使う場合は十分な精度を保持します。
擬似逆行列計算機のよくある質問
擬似逆行列は逆行列とどう違いますか?
通常の逆行列は非特異な正方行列にのみ存在します。ムーア・ペンローズ擬似逆行列は長方形行列にも定義され、逆行列が存在するときはそれと一致します。
擬似逆行列で何が解けますか?
整合しない過剰決定系の最小二乗解を求められます。劣決定系ではユークリッドノルムが最小の解を選びます。
この計算機は特異行列に対応していますか?
行または列がフルランクの長方形行列と正方行列に対応します。ランク不足の行列には SVD ベースの方法が必要なため、不確かな近似値ではなく制限を報告します。
行列はどの形式で入力しますか?
同じ行の値はカンマまたは空白で、行はセミコロンで区切ります。すべての行に同じ数の数値成分が必要です。
なぜ結果を丸めるのですか?
読みやすさのため小数点以下六桁で表示します。内部では倍精度浮動小数点数を使いますが、悪条件行列では精度が失われることがあります。