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

一つの解、二つの絵:最良近似直線は射影の影だった

I重要な公式

公式 名称 説明
ATAx^=ATbA^T A \hat{x} = A^T b 正規方程式 最小二乗問題を解く
p=Ax^\mathbf{p} = A\hat{x} 射影 Col(A)\text{Col}(A) の中で b\mathbf{b} に最も近いベクトル
e=b−p\mathbf{e} = \mathbf{b} - \mathbf{p} 誤差ベクトル e⊥Col(A)\mathbf{e} \perp \text{Col}(A) を満たさなければならない
ATe=0A^T \mathbf{e} = \mathbf{0} 直交条件 誤差は AA のすべての列に垂直

II可視化1 — データ空間での最良近似

3つの点すべてを通る直線はありません。最良近似とは、残差の二乗の和 ∑ei2\sum e_i^2 を最小にする直線のことです。

III可視化2 — 列空間への射影

p=Ca1+Da2\mathbf{p}=C\mathbf{a}_1+D\mathbf{a}_2 は、平面 Col(A)\text{Col}(A) に落ちた b\mathbf{b} の影です。誤差 e=b−p\mathbf{e}=\mathbf{b}-\mathbf{p} がいちばん短くなるのは、ちょうど平面と直角に交わるときです。

IV可視化3 — 誤差の地形

SSE(C,D)\text{SSE}(C,D) は下に凸なお椀の形をしていて、正規方程式はそのただ一つの最低点へまっすぐ飛び込みます。

V確認クイズ

問1

データ点 (1,1)(1,1)、(2,2)(2,2)、(3,2)(3,2) とモデル y=C+Dty = C + Dt から、連立方程式を Ax=bA\mathbf{x} = \mathbf{b} の形に書きます。右辺のベクトル b\mathbf{b} はどれですか。

✅ 正解です! ベクトル b\mathbf{b} は、各データ点の yy の値を集めたものです。

❌ 惜しい! 各方程式は C+Dti=yiC + Dt_i = y_i なので、b\mathbf{b} は yy の値を並べた列です。

解答を見る

解答:

各データ点 (ti,yi)(t_i, y_i) を y=C+Dty = C + Dt に代入します:

C+D(1)=1,C+D(2)=2,C+D(3)=2C + D(1) = 1, \quad C + D(2) = 2, \quad C + D(3) = 2

行列の形で書くと:

(111213)⏟A(CD)⏟x=(122)⏟b\underbrace{\begin{pmatrix} 1 & 1 \\ 1 & 2 \\ 1 & 3 \end{pmatrix}}_{A} \underbrace{\begin{pmatrix} C \\ D \end{pmatrix}}_{\mathbf{x}} = \underbrace{\begin{pmatrix} 1 \\ 2 \\ 2 \end{pmatrix}}_{\mathbf{b}}

右辺 b=(1,2,2)\mathbf{b} = (1, 2, 2) は yy の値を集めたものです。

よくある間違い:

  • (1,2,3)(1, 2, 3) = tt の値 = AA の2列目
  • (1,1,1)(1, 1, 1) = 1 が並んだ列 = AA の1列目
  • (2/3,1/2)(2/3, 1/2) = 解 x^\hat{\mathbf{x}} であって、b\mathbf{b} ではない

問2

正しいか誤りか: p⋅e=0\mathbf{p} \cdot \mathbf{e} = 0(射影が誤差に垂直)ならば、 最小二乗解が得られていることが保証される。

✅ 正解です! 条件 p⋅e=0\mathbf{p} \cdot \mathbf{e} = 0 は必要ですが、十分ではありません。必要なのは ATe=0A^T\mathbf{e} = \mathbf{0} です。

❌ 惜しい! p\mathbf{p} に垂直であるだけでは、列空間の平面全体に垂直であることは保証されません。

解答を見る

解答:誤り

p⋅e=0\mathbf{p} \cdot \mathbf{e} = 0 を確かめても、1つのベクトル(p\mathbf{p})に垂直であることしかわかりません。

最小二乗法の条件は次のとおりです:

ATe=0A^T\mathbf{e} = \mathbf{0}

つまり e\mathbf{e} は、列空間のすべてのベクトルに垂直でなければなりません。 a1\mathbf{a}_1 にも a2\mathbf{a}_2 にも、そしてそれらのあらゆる一次結合にもです。

R3\mathbb{R}^3 には、p\mathbf{p} には垂直でも平面全体には垂直でないベクトルがたくさんあります。 ∥e∥\|\mathbf{e}\| を最小にするのは、列空間全体に垂直なベクトルだけです。

問3

最適な誤差ベクトルは e=(−1/6,  1/3,  −1/6)\mathbf{e} = (-1/6, \; 1/3, \; -1/6) で、AA の2列目は a2=(1,2,3)\mathbf{a}_2 = (1, 2, 3) です。

e⋅a2\mathbf{e} \cdot \mathbf{a}_2 はいくつですか。

✅ 正解です! 内積は 0 です。ATe=0A^T\mathbf{e} = \mathbf{0} が要求するとおり、e⊥a2\mathbf{e} \perp \mathbf{a}_2 が確かめられました。

❌ 惜しい! (−1/6)(1)+(1/3)(2)+(−1/6)(3)(-1/6)(1) + (1/3)(2) + (-1/6)(3) を丁寧に計算しましょう。すべての項を6分のいくつかにそろえてみてください。

解答を見る

解答:

e⋅a2=(−16)(1)+(13)(2)+(−16)(3)\mathbf{e} \cdot \mathbf{a}_2 = \left(-\frac{1}{6}\right)(1) + \left(\frac{1}{3}\right)(2) + \left(-\frac{1}{6}\right)(3)

6分のいくつかにそろえると:

=−16+46−36=−1+4−36=06=0= -\frac{1}{6} + \frac{4}{6} - \frac{3}{6} = \frac{-1 + 4 - 3}{6} = \frac{0}{6} = 0

これで e⊥a2\mathbf{e} \perp \mathbf{a}_2 が確かめられました。これは ATe=0A^T\mathbf{e} = \mathbf{0} に含まれる2つの条件のうちの1つです。 (もう一方の条件 e⋅a1=0\mathbf{e} \cdot \mathbf{a}_1 = 0 も成り立ちます。)

問4

「二つの絵」の対応で、ベクトル空間の絵の列空間 Col(A)\text{Col}(A) は、 データ空間の絵では何に対応しますか。

✅ 正解です! Col(AA) の各点は、直線 y=C+Dty = C + Dt の当てはめ値に対応します。列空間全体が、考えられるすべての直線の当てはめを表しています。

❌ 惜しい! CC と DD を動かすと AxA\mathbf{x} は何を生み出すか考えましょう。選び方ごとに、違う直線と違う当てはめ値が得られます。

解答を見る

解答:

列空間 Col(A)\text{Col}(A) は、CC と DD を動かしたときのベクトル Ca1+Da2C\mathbf{a}_1 + D\mathbf{a}_2 全体からなります。そのようなベクトルはどれも、次の当てはめ値を与えます:

(C+D,  C+2D,  C+3D)(C + D, \; C + 2D, \; C + 3D)

これは、直線 y=C+Dty = C + Dt の t=1,2,3t = 1, 2, 3 での値です。

したがって、Col(A)\text{Col}(A) = 当てはめ可能なすべての直線 y=C+Dty = C + Dt の集合です。

二つの絵の対応をまとめると:

データ空間 ベクトル空間
データ点 (y1,y2,y3)(y_1, y_2, y_3) ベクトル b\mathbf{b}
考えられるすべての直線 y=C+Dty = C + Dt 列空間 Col(A)\text{Col}(A)
最良近似直線上の当てはめ値 射影 p\mathbf{p}
縦の残差 誤差ベクトル e\mathbf{e}
∑ei2\sum e_i^2 を最小化 ∥e∥2\|\mathbf{e}\|^2 を最小化

正解数:0 / 4