平面上の 4 点 に、直線 を最小二乗法であてはめる。すなわち、実数 を動かして
を最小にする。ここで は上の 4 点である。 を最小にする を、組 の形で求めよ。
直交射影と最小二乗法
部分空間への直交射影を求め、正規方程式 で最小二乗解を求められる
前提
- 内積空間とグラム・シュミットの直交化 (大学数学 線形代数)
この技能までの道筋(37 技能)
読み込み中…
右端(太枠)がこの技能です。左へドラッグすると、さらに前の前提をさかのぼれます。技能を押すと、そこまでの道筋が光ります。
一覧で見る(前提が先)
- 三角形の角の大きさ (算数 平面図形)
- 対頂角と平行線の同位角・錯角 (中学数学 平面図形)
- 三角形の外角と角の組み合わせ (中学数学 平面図形)
- 小数の四則計算 (算数 計算と数のきまり)
- 分数のたし算・ひき算 (算数 計算と数のきまり)
- 分数のかけ算・わり算 (算数 計算と数のきまり)
- 小数と分数のまじった計算 (算数 計算と数のきまり)
- 逆算(□ を求める) (算数 計算と数のきまり)
- 正負の数の加法・減法 (中学数学 数と式)
- 四則の混じった計算の順序 (算数 計算と数のきまり)
- 正負の数の乗除・累乗と四則混合 (中学数学 数と式)
- 分配法則を使った計算のくふう (算数 計算と数のきまり)
- 文字式の計算(同類項・分配法則) (中学数学 数と式)
- 1 次方程式の解き方 (中学数学 方程式)
- 二等辺三角形の角 (中学数学 平面図形)
- 平行四辺形と特別な四角形の性質 (中学数学 平面図形)
- 三角形・四角形の面積 (算数 平面図形)
- 約数の見つけ方 (算数 整数の性質)
- 素因数分解と約数の個数 (算数 整数の性質)
- 平方根の意味と根号の簡単化 (中学数学 数と式)
- 三平方の定理 (中学数学 相似と三平方の定理)
- 乗法公式による展開 (中学数学 数と式)
- 折り返した図形の角と長さ (中学数学 平面図形)
- 三平方の定理と平面図形 (中学数学 相似と三平方の定理)
- ベクトルの演算と成分 (数学C ベクトル)
- 行列の和・積と 2 次の逆行列 (大学数学 線形代数)
- 掃き出し法と連立 1 次方程式 (大学数学 線形代数)
- 行列の階数と解の存在 (大学数学 線形代数)
- ベクトル空間・一次独立・基底・次元 (大学数学 線形代数)
- 特別な直角三角形の辺の比 (中学数学 相似と三平方の定理)
- 最大公約数と最小公倍数 (算数 整数の性質)
- 比の性質と簡単な比 (算数 割合と比)
- 相似な三角形と長さの比 (算数 平面図形)
- 三角形の相似条件と対応する辺の長さ (中学数学 相似と三平方の定理)
- 三角比の定義と値 (数学I 図形と計量)
- 内積となす角 (数学C ベクトル)
- 内積空間とグラム・シュミットの直交化 (大学数学 線形代数)
説明
部分空間 の外にあるベクトル に対して、 の中で に最も近いベクトル(直交射影)を求められるようになります。これを応用すると、方程式の数が未知数より多くて厳密には解けない連立方程式 に対し、誤差の二乗和 を最小にする (最小二乗解)が求まります。測定データに直線や曲線を当てはめる場面で使います。
なぜそうしてよいのか
この項では の標準内積 を使います。ここで は の行と列を入れ替えた行列(転置行列)で、縦ベクトル の転置 は横ベクトルです。 を 行列とし、その列ベクトル が張る部分空間を とします。 なので、 を動かしたとき が動く範囲はちょうど です。よって「 を最小にする」ことは「 の中で に最も近い点を探す」ことと同じです。いま で、 が のすべてのベクトルと直交しているとします。 の任意の について なので と は直交し、三平方の定理から です。等号は のときだけなので、最も近い点は ただ 1 つです。この を の への直交射影といいます。
このような は必ず存在します。グラム・シュミットの直交化で の正規直交基底 を作り、 とおけば、 となり、 は基底すべて、したがって 全体と直交するからです。一方、正規直交基底を作らずに済ませるには、直交条件を の列で書きます。 が と直交することは、各列と直交すること、つまり と同値です。 の第 行は なので、これらをまとめると 、すなわち正規方程式 になります。 だから となる があり、それが正規方程式の解なので、正規方程式はいつも解をもちます。逆に正規方程式のどの解 も残差が と直交するので となり、最小二乗解です。
さらに の列が一次独立なら解はただ 1 つです。 とすると より 、列が一次独立なので です。つまり 行列 は逆行列をもち(正則)、 と書けます。列が一次従属のときは最小二乗解が無数にありますが、射影 は 1 つに決まります。
手順
問題を の形に書く。求めたい係数を にまとめ、 の列が射影先の部分空間を張るようにする(直線 の当てはめなら、1 列目はすべて 1、2 列目は各点の 座標)
( 行列)と ( 次元ベクトル)を計算する。成分はそれぞれ と
正規方程式 を掃き出し法などで解き、最小二乗解 を求める( の列が一次独立なら解はただ 1 つ)
必要なら直交射影 を計算し、残差 と の各列との内積が 0 になることで検算する
例題
3 点 に直線 を当てはめる。誤差の二乗和 を最小にする を求めよ。
例題の解答
手順 1:各点で としたいので、、、 とおくと、誤差の二乗和は です。 の 2 つの列は比例していないので一次独立です。
手順 2:、。
手順 3:正規方程式は 。2 式目から 1 式目を引くと より 。1 式目に代入して より 。
手順 4:、残差は 。1 列目との内積は 、2 列目との内積は で、確かに直交しています(二乗和の最小値は )。
よって 。
答え:
よくある誤り
が正方行列でないのに と逆行列で解こうとする誤り。 行列に逆行列はない。正しくは両辺に左から を掛けた形の正規方程式 を解く。
の代わりに を使ってしまう誤り。例題では は 行列で階数 2 なので正則でなく、未知数 の個数とも合わない。正しくは 行列 を使う。
直交していない基底 に対して を使ってしまう誤り。例えば の基底 と では となり、本来の射影である 自身にならない。正しくは、先にグラム・シュミットで直交化してからこの公式を使うか、正規方程式を解く。
直線 の当てはめで定数項に対応する「すべて 1 の列」を忘れ、 を を並べた 1 列だけにしてしまう誤り。これでは原点を通る直線 への当てはめになる。正しくは の 1 列目を 、2 列目を とする。
この説明は AI が書き、別の AI が例題を独立に解いて答えを確かめ、説明を一歩ずつ査読しました。
入口の問題
平面上の 4 点 に、直線 ( は実数)を最小二乗法であてはめる。すなわち、 を最小にする実数 を求める。ここで は上の 4 点である。
とおくと、( は の標準内積によるノルム)である。正規方程式 を解いて、 を最小にする を組 の形で答えよ。
ルベーグ測度に関する実ヒルベルト空間 を、内積 で考える。部分空間 をとり、定数関数 とする。 と との距離 を求めよ。
この技能を使う問題
に標準内積を考える。ベクトル
に対し、 で張られる部分空間を とする。 と の距離、すなわち が 全体を動くときの の最小値を求めよ。に標準内積を考え、
とする。 の張る部分空間を とし、 の への直交射影を とするとき、 を求めよ。の列ベクトル , に対して
と定めると、これは の内積である。この内積に関するノルムを とする。, で張られる部分空間を とし、 とする。この内積に関する と の距離、すなわち から への直交射影を としたときの を求めよ。
ルベーグ測度に関する実ヒルベルト空間 において、 次以下の多項式全体のなす部分空間を とする。関数 と との距離の 乗
を求めよ。ルベーグ測度に関する実ヒルベルト空間 (内積 )を考える。実数 を動かすとき、
の最小値を求めよ。実ヒルベルト空間 を、内積 とそれから定まるノルム で考える。部分空間 を で定め、 とする。 と との距離 について、 の値を求めよ。
に標準内積 を考え、ノルムを とする。
とし、 を で張られる の部分空間とする。 が 全体を動くときの の最小値( と の距離)を求めよ。に標準内積 を考え、そのノルムを とする。
とし、 を で張られる の部分空間とする。 と の距離
を求めよ。を正の整数とし、 には標準内積 を考える。 次実正方行列 が次の 2 つの条件を満たすとする。
このとき、 は対称行列である(すなわち である)ことを証明せよ。ただし、 の部分空間の正規直交基底を 全体の正規直交基底に延長できること(グラム・シュミットの直交化法)、 が像 の次元に等しいこと、トレースの性質 は用いてよい。
実ヒルベルト空間
を、内積 とそれから定まるノルム で考える。部分空間 を
で定め、 を とする。このとき
の値を求めよ。とし、 を 上のルベーグ測度とする。 上の正則関数 で
をみたすもの全体を とする(ベルグマン空間)。これは内積 により複素ヒルベルト空間になる。
部分空間
を考え、 とする。 と との距離
を求めよ。実数係数の の多項式で次数が 以下のもの全体のなす実線形空間に、内積
を考える。実数 を動かすとき
の最小値を求めよ。
この技能の次に学べること
- ヒルベルト空間と射影定理 (大学院数学 関数解析)