線形代数 · 第13回 · 動画3 · インタラクティブ教材

解けない連立方程式に Aᵀ を掛ける――正規方程式と AᵀA が可逆になる条件

I重要な公式

公式 名称 使う場面
ATAx^=ATbA^T A \hat{x} = A^T b 正規方程式 解けない Ax=bAx = b を、解ける連立方程式に変える
x^=(ATA)−1ATb\hat{x} = (A^T A)^{-1} A^T b 最小二乗解 AA の列が一次独立のとき
AA が m×n⇒ATAm \times n \Rightarrow A^T A は n×nn \times n サイズの規則 長方形の連立方程式 → 正方形の連立方程式
ATAA^T A が可逆   ⟺  \iff AA の列が一次独立 可逆性の判定 x^\hat{x} がただ1つに決まるかを判定する

II可視化1 — 列空間への射影

Col(A)\text{Col}(A) を直線として描くと、最も短い残差 b−pb - p は、Col(A)\text{Col}(A) と直角に交わるものです。

💡 bb がはじめから Col(A)\text{Col}(A) の上にあるときは、この最短の残差は 00 です。そのとき Ax=bAx = b には厳密な解があり、x^\hat{x} がそのままその解になります。

III可視化2 — ATAA^T A が特異になるとき

A(t)A(t) の列が t=1t = 1 で一次従属になるのに向かって動くにつれ、det⁡(ATA)=26(1−t)2\det(A^T A) = 26(1 - t)^2 は 00 まで下がり、ATAA^T A は逆行列を失います。

💡 ATAA^T A が可逆になるのは、ちょうど AA の列が一次独立のときです。ですから x^\hat{x} がただ1つに決まるかどうかは、その一次独立性しだいです。

IV可視化3 — 最小二乗法は誤差を最小にする

正規方程式の解は、残差の二乗和 ∑ri2\sum r_i^2(RSS)をできるだけ小さくする、ただ1本の直線です。

💡 GPS もこの仕組みで動いています。スマートフォンは、未知数を決めるのに必要な数より多くの衛星の方程式に自分の位置を当てはめ、正規方程式によって各信号の雑音をならしているのです。

V確認クイズ

問1 · 正規方程式のサイズ

AA が 7×37 \times 3 行列(方程式7つ、未知数3つ)のとき、ATAA^T A のサイズはいくつですか?

✅ 正解です! (3×7)(7×3)=3×3(3 \times 7)(7 \times 3) = 3 \times 3 です。正規方程式からは、いつでも n×nn \times n の正方形の連立方程式が得られます!

❌ 惜しい! AA が m×nm \times n なら、ATA^T は n×mn \times m です。内側のサイズがそろうように掛けると、(n×m)(m×n)=n×n(n \times m)(m \times n) = n \times n になります。

解答を見る

解答:

AA は 7×37 \times 3 なので、ATA^T は 3×73 \times 7 です。したがって:

ATA=(3×7)⏟AT⋅(7×3)⏟A=(3×3)A^T A = \underbrace{(3 \times 7)}_{A^T} \cdot \underbrace{(7 \times 3)}_{A} = (3 \times 3)

規則:AA が m×nm \times n なら、ATAA^T A はいつでも n×nn \times n の正方行列です。サイズは方程式の数ではなく、未知数の数に一致します。

問2 · ATAA^T A の対称性

正しいか誤りか:行列 ATAA^T A は、AA の形や成分にかかわらず、いつでも対称である。

✅ 正解です! (ATA)T=AT(AT)T=ATA(A^T A)^T = A^T(A^T)^T = A^T A です。どんな行列 AA についても、いつでも対称です!

❌ 惜しい! (ATA)T(A^T A)^T を、(XY)T=YTXT(XY)^T = Y^T X^T を使って一歩ずつ計算してみましょう。

解答を見る

解答:正しい

ATAA^T A の転置をとります:

(ATA)T=AT(AT)T=ATA(A^T A)^T = A^T (A^T)^T = A^T A

ここでは2つの規則を使っています:

  1. (XY)T=YTXT(XY)^T = Y^T X^T(積の転置は順序が逆になる)
  2. (AT)T=A(A^T)^T = A(2回転置するともとに戻る)

ATAA^T A の転置は、それ自身に等しくなります。これがまさに対称の定義です! この性質は、どんな形のどんな行列 AA についても成り立ちます。

動画の例では ATA=[38830]A^T A = \begin{bmatrix}3 & 8\\8 & 30\end{bmatrix} で、非対角成分はどちらも 8 です。✓

問3 · det(AᵀA) を求める

A=[111215]A = \begin{bmatrix}1 & 1\\1 & 2\\1 & 5\end{bmatrix} のとき、det⁡(ATA)\det(A^T A) はいくつですか?

(ヒント:まず ATAA^T A を計算し、次に det⁡ ⁣[abcd]=ad−bc\det\!\begin{bmatrix}a & b\\c & d\end{bmatrix} = ad - bc を使います。)

✅ 正解です! 3×30−8×8=90−64=26≠03 \times 30 - 8 \times 8 = 90 - 64 = 26 \neq 0 です。可逆なので、x^\hat{x} はただ1つに決まります!

❌ 惜しい! まず ATA=[38830]A^T A = \begin{bmatrix}3 & 8\\8 & 30\end{bmatrix} を求め、それから det⁡=3(30)−8(8)\det = 3(30) - 8(8) を計算しましょう。

解答を見る

解答:

ステップ1:ATAA^T A を計算します。

AT=[111125]A^T = \begin{bmatrix}1 & 1 & 1\\1 & 2 & 5\end{bmatrix}

ATA=[1 ⁣⋅ ⁣1+1 ⁣⋅ ⁣1+1 ⁣⋅ ⁣11 ⁣⋅ ⁣1+1 ⁣⋅ ⁣2+1 ⁣⋅ ⁣51 ⁣⋅ ⁣1+2 ⁣⋅ ⁣1+5 ⁣⋅ ⁣11 ⁣⋅ ⁣1+2 ⁣⋅ ⁣2+5 ⁣⋅ ⁣5]=[38830]A^T A = \begin{bmatrix}1\!\cdot\!1+1\!\cdot\!1+1\!\cdot\!1 & 1\!\cdot\!1+1\!\cdot\!2+1\!\cdot\!5\\1\!\cdot\!1+2\!\cdot\!1+5\!\cdot\!1 & 1\!\cdot\!1+2\!\cdot\!2+5\!\cdot\!5\end{bmatrix} = \begin{bmatrix}3 & 8\\8 & 30\end{bmatrix}

ステップ2:行列式を求めます。

det⁡(ATA)=3×30−8×8=90−64=26\det(A^T A) = 3 \times 30 - 8 \times 8 = 90 - 64 = 26

26≠026 \neq 0 なので、ATAA^T A は可逆です → x^\hat{x} がただ1つ存在します!

誤答を選びやすい理由:

  • 90 = adad だけ(bcbc を引き忘れた)
  • 64 = bcbc だけ(adad の項を忘れた)
  • 0 = 列が一次従属であることを意味する(ここでは列は一次独立)

問4 · ただ1つに決まらなくなるのはいつ?

A=[131313]A = \begin{bmatrix}1 & 3\\1 & 3\\1 & 3\end{bmatrix} のとき、正規方程式 ATAx^=ATbA^T A \hat{x} = A^T b を作るとどうなりますか?

✅ 正解です! 列が一次従属 → det⁡(ATA)=81−81=0\det(A^T A) = 81 - 81 = 0 → 特異 → x^\hat{x} はただ1つには決まりません。

❌ 惜しい! 2列目 = 3×3 \times 1列目 になっていることに注目しましょう。列は一次従属です。一次従属な列と ATAA^T A について、重要な定理は何と言っていたでしょうか?

解答を見る

解答:

AA の2列目は [333]=3⋅[111]\begin{bmatrix}3\\3\\3\end{bmatrix} = 3 \cdot \begin{bmatrix}1\\1\\1\end{bmatrix} = 3 × 1列目です。列は一次従属です!

ATA=[111333][131313]=[39927]A^T A = \begin{bmatrix}1&1&1\\3&3&3\end{bmatrix}\begin{bmatrix}1&3\\1&3\\1&3\end{bmatrix} = \begin{bmatrix}3 & 9\\9 & 27\end{bmatrix}

det⁡=3×27−9×9=81−81=0\det = 3 \times 27 - 9 \times 9 = 81 - 81 = 0

ATAA^T A は特異(可逆でない)です。正規方程式に解はありますが、ただ1つには決まりません。

これで、動画の重要な定理が確かめられます:ATAA^T A が可逆になるのは、AA の列が一次独立のとき、かつそのときに限ります。

ほかの選択肢が誤りである理由:

  • ATAA^T A は、AA の形にかかわらず、いつでも作れます。サイズは n×nn \times n です。
  • ここでの ATAA^T A は 2×22 \times 2 です(3×33 \times 3 ではありません)。AA の列は2本だからです。
  • 可逆性は保証されません。列が一次独立かどうかで決まります。

正解数:0 / 4