射影・直和・直交直和のバックアップ差分(No.36)

更新
追加された行はこの色です。
削除された行はこの色です。
[[前の単元 ＜＜＜>線形代数II/内積と計量空間]]
　　　　　　　　　　　　　　　[[線形代数II]]
　　　　　　　　　　　　　　　[[＞＞＞ 次の単元>線形代数II/固有値問題・固有空間・スペクトル分解]]

* 目次 [#xf03612b]

#contents
#mathjax
&katex();

* ベクトルの成分 [#x61c2758]

#ref(射影.png,right,around,40%);

規格化されたベクトル &math(\bm e); に対して、ベクトル &math(\bm x); を 
- &math(\bm e); に平行な成分 &math(\bm x_{\parallel}=x_\parallel \bm e); と、
- &math(\bm e); に垂直な成分 &math(\bm x_{\perp}); とに分け、

&math(\bm x=\bm x_{\parallel}+\bm x_{\perp}); としたい。

両辺に左から &math(\bm e); をかければ、

　&math((\bm e,\bm x)=x_{\parallel}(\bm e,\bm e)+(\bm e,\bm x_{\perp})=x_{\parallel});

が得られ、

　&math(\bm x_{\parallel}=(\bm e,\bm x)\bm e);~
　&math(\bm x_{\perp}=\bm x-\bm x_\parallel=\bm x-(\bm e,\bm x)\bm e);

としてこれらのベクトルを求められる。~
（同じことをグラム・シュミットの直交化で行った）

この &math(\bm x_\parallel); を &math(\bm x); の &math(\bm e); への射影と呼ぶ。

&math(\bm e); に垂直な光を &math(\bm x); に当てたとき、
&math(\bm e); 軸上にできる影が &math(\bm x_\parallel); 
であるという気持ちが込められている → 「射影」

#clear
*** 注意１ [#r5af3ba1]

規格化されていない &math(\bm v); 方向の成分を求めるなら、&math(\bm e=\bm v/\|\bm v\|); だから、

　&math(\bm x_{\parallel}=(\frac{\bm v}{\|\bm v\|},\bm x)\,\frac{\bm v}{\|\bm v\|}=\frac{(\bm v,\bm x)}{\|\bm v\|^2}\bm v);~

*** 注意２ [#s8228082]

複素ベクトルに対しては &math((\bm x,\bm e)\ne(\bm e,\bm x)); なので、
どちらから掛けるかが重要である。

&math((\bm e,\bm x)=(\bm e,x_\parallel \bm e)=x_{\parallel}); だが、~
&math((\bm x,\bm e)=(x_\parallel \bm e,\bm e)=\overline{x_{\parallel}}); となってしまう。

*** 注意３ [#p72a2888]

この授業では &math((\bm a,k\bm b)=k(\bm a,\bm b)); となる内積の公理を採用しているため
上記が正しいが、

多くの教科書では &math((k\bm a,\bm b)=k(\bm a,\bm b)); を採用しているため、
そのような公理系では左ではなく右から掛ける必要がある。

* 射影演算子 [#c9c77b82]

&math(\bm x); から &math(\bm x_\parallel); を求める演算、

　&math(P_{\bm e}:\bm x\mapsto\bm x_\parallel);

は線形変換であり、&math(P_{\bm e}); は射影変換あるいは射影演算子と呼ばれる。

正規直交基底 &math(A); の下での数ベクトル表現を考えれば、

　&math(
(\bm a_A,\bm b_A)&=\sum_k^n \overline{a_k}b_k=\begin{pmatrix}\overline a_1&\overline a_2&\dots&\overline a_n\end{pmatrix}
\begin{pmatrix}b_1\\b_2\\\vdots\\b_n\end{pmatrix}={}^t\overline {\bm a_A}\bm b_A=\bm a_A^\dagger\bm b_A
); 

となることを用いて、

　&math(
(\bm x_{\parallel})_A
&=(\bm e_A,\bm x_A)\bm e_A\\
&=\{\bm e_A^\dagger \bm x_A\}\bm e_A\\
&=\bm e_A\{\bm e_A^\dagger \bm x_A\}\hspace{1cm}\because \{\ \}内はスカラー\\
&=\{\bm e_A\bm e_A^\dagger\}\bm x_A\hspace{1cm}\because 結合法則\\
&=(P_{\bm e})_A\,\bm x_A\\
);

すなわち &math(P_{\bm e}); の表現は、任意の正規直交基底 &math(A); 
に対して次のように表せる。

　&math(
(P_{\bm e})_A&=\bm e_A\bm e^\dagger_A=
\begin{pmatrix}
e_1\\e_2\\\vdots\\e_n
\end{pmatrix}
\begin{pmatrix}
\overline{ e_1}&\overline{ e_2}&\dots&\overline{ e_n}
\end{pmatrix}\\
&=\begin{pmatrix}
e_1\overline{e_1}&e_1\overline{e_2}&\cdots&e_1\overline{e_n}\\
e_2\overline{e_1}&e_2\overline{e_2}&&\vdots\\
\vdots&&\ddots&\vdots\\
e_n\overline{e_1}&\cdots&\cdots&e_n\overline{e_n}
\end{pmatrix}
);


*** 例 [#p94d13e4]

&math(\bm x=\begin{pmatrix}x\\y\\z\end{pmatrix}); の &math(\bm v=\begin{pmatrix}1\\-1\\0\end{pmatrix}); 方向成分を求めよう。

&math(\bm e=\frac{1}{\|\bm v\|}\bm v=\frac{1}{\sqrt 2}\begin{pmatrix}1\\-1\\0\end{pmatrix}); を用いて、

&math(
P_{\bm v}=P_{\bm e}=\bm e\bm e^\dagger
=\frac{1}{2}\begin{pmatrix}1\\-1\\0\end{pmatrix}\begin{pmatrix}1&-1&0\end{pmatrix}
=\frac{1}{2}\begin{pmatrix}1&-1&0\\-1&1&0\\0&0&0\end{pmatrix}
);

&math(\bm x_\parallel&=P_{\bm v}\bm x\\
&=\frac{1}{2}\begin{pmatrix}1&-1&0\\-1&1&0\\0&0&0\end{pmatrix}\begin{pmatrix}x\\y\\z\end{pmatrix}
=\frac{1}{2}\begin{pmatrix}x-y\\-x+y\\0\end{pmatrix}
=\frac{x-y}{2}\begin{pmatrix}1\\-1\\0\end{pmatrix}\hspace{1cm}\parallel \bm v
);

&math(\bm x_\perp&=\bm x-\bm x_\parallel\\
&=\begin{pmatrix}x\\y\\z\end{pmatrix}-\frac{x-y}{2}\begin{pmatrix}1\\-1\\0\end{pmatrix}
=\begin{pmatrix}(x+y)/2\\(x+y)/2\\z\end{pmatrix}\hspace{1cm}\perp \bm v
);

** 射影演算子はエルミートになる。 [#tee56cfe]

　&math(
(P_{\bm e})_A^\dagger=\big(\bm e_A\bm e_A^\dagger\big)^\dagger=\big(\bm e_A^\dagger\big)^\dagger\bm e_A^\dagger=\bm e_A\bm e_A^\dagger=(P_{\bm e})_A
);

より、射影演算子の表現行列はエルミートである。

このとき、任意のベクトル &math(\bm x,\bm y); に対して 

　&math(
(\bm x,P_{\bm e}\bm y)=(P_{\bm e}^\dagger\bm x,\bm y)=(P_{\bm e}\bm x,\bm y)
);

が成り立ち、このような演算子はエルミート演算子と呼ばれる。

* $n$ 次元空間への射影を考える [#oec8bd2b]

ここまで、あるベクトルに平行な直線（一次元空間）への射影を考えたが、
以下では平面への射影や、もっと一般に &math(n); 次元空間への射影を考える。

そのためにまずはいくつか準備を行う。
** 復習１：線形空間 [#y64f3022]

&math(K); 上の線形空間とは、ベクトルの和とスカラー倍について閉じた集合のことだった。

- 任意の &math(\bm x,\bm y\in V); に対して、必ず &math(\bm x+\bm y\in V);~
　　↔ &math(\bm x,\bm y\in V); から &math(\bm x+\bm y\in V); を導けるということ
- 任意の &math(\bm x\in V,k\in K); に対して、必ず &math(k\bm x\in V);~
　　↔ &math(\bm x\in V,k\in K); から &math(k\bm x\in V); を導けるということ

** 復習２：部分空間 [#t397a1fd]

線形空間の部分集合 &math(W\subset V); がベクトルの和とスカラー倍について閉じている場合、
&math(W); も線形空間となり、&math(W); は &math(V); の部分空間であるという。

&math(\mathbb R^3); の部分空間：

- ０次元の部分空間は原点のみからなる集合 &math(\set{\bm 0});
- １次元の部分空間は原点を通る直線　　　 &math(\set{\bm p=s\bm a|s\in K});
- ２次元の部分空間は原点を通る平面　　　 &math(\set{\bm p=s\bm a+t\bm b|s,t\in K});
- ３次元の部分空間は &math(\mathbb R^3); そのもの

同じ直線的でも、原点を通らない &math(\set{\bm p=s\bm a+\bm b|s\in K}); 
は線形空間にならない。（和やスカラー倍が元の集合からはみ出す）

** 復習３：集合の積と和 [#sea08144]

集合 &math(A); と集合 &math(B); の積と和は、

#multicolumns
:積（交わり）|
&math(A\cap B=\set{x|x\in A\,\mathrm{かつ}\,x\in B});~
&math(A); および &math(B); の両方に含まれる要素の集合~
&math(A); キャップ &math(B); と読む。(&math(\cap); は帽子の形)
#multicolumns
:和（結び）|
&math(A\cup B=\set{x|x\in A\,\mathrm{または}\,x\in B});~
&math(A); あるいは &math(B); の少なくとも片方に含まれる要素の集合~
&math(A); カップ &math(B); と読む。(&math(\cup); はカップの形)
#multicolumns(end)

#ref(積集合和集合.png,center,33%);

~

記号の覚え方：
-「&math(x\in A\,\mathrm{かつ}\,x\in B);」は英語では「&math(x\in A\,\mathrm{and}\,x\in B);」
- And の A と &math(\cap); とは似ている（でしょ？）

* 以下、$U$ の部分空間 $V,W$ について考える [#i78b84b7]

&math(K); 上の線形空間 &math(U); の部分空間 &math(V,W); を考え、~
&math(\{\bm v_1,\bm v_2,\dots,\bm v_n\},);, &math(\{\bm w_1,\bm w_2,\dots,\bm w_m\}); をそれぞれの基底とする。~
(&math(\dim V=n,\ \dim W=m);)

* 交空間 $V\cap W$ [#oa2a08fa]

２つの線形空間の、集合としての交わり &math(V\cap W); は常に線形空間になり、交空間と呼ばれる。

>証明：~
&math(\bm x,\bm y\in V\cap W,\ k\in K); とする。~
&math(\bm x,\bm y\in V); かつ &math(\bm x,\bm y\in W); であるから、~
&math(\bm x+\bm y\in V); かつ &math(\bm x+\bm y\in W); また
&math(k\bm x\in V); かつ &math(k\bm x\in W); ~
すなわち、&math(\bm x+\bm y, k\bm x\in V\cap W); であり、~
&math(V\cap W); はベクトルの和とスカラー倍に対して閉じている。

交わり &math(V\cap W); が空集合になることはない。

線形空間は必ず &math(\bm 0); を含むから、常に &math(\bm 0\in V\cap W); である。

&math(V\cap W=\set{\bm 0}); のとき、&math(\dim(V\cap W)=0);

&attachref(1D_cap_2D.png,,25%);　
&attachref(2D_cap_2D.png,,25%);

* 和集合 $V\cup W$ はベクトル和に対して閉じていないことがある [#xd00301e]

#ref(1D_cup_1D.png,right,around,25%);

例えば図のように、２つの１次元空間 &math(V,W); 
の和集合 &math(V\cup W); は、
原点で交わる２本の直線の形をしている。

&math(V,W); 上から２つのベクトルを取り
&math(\bm v\in V, \bm w\in W); とすれば、~
&math(\bm v, \bm w\in V\cap W); でない限り、
&math(\bm v+\bm w\notin V\cup W); である。

すなわち、和集合は必ずしも線形空間にならない

#clear
* 和空間 $V+W$ [#i53e231d]

和集合をベクトル和について閉じるように拡大した線形空間が和空間 &math(V+W); である。

これは &math(V); の元と &math(W); の元の和で表せるすべてのベクトルからなる集合である。

　&math(V+W\equiv\set{\bm x=\bm x_V+\bm x_W|\bm x_V\in V,\bm x_W\in W});

&math(\bm x_V\in V,\bm x_W\in W); はそれぞれ、&math(V,W); の基底 
&math(\set{\bm v_k},\set{\bm w_k}); の線形結合として表せるから、
&math(\bm x\in V+W); は

　&math(
\bm x&=\bm x_V+\bm x_W\\
&=\underbrace{\sum_{k=1}^n c_k\bm v_k}_{\bm x_V}+
  \underbrace{\sum_{k=1}^m d_k\bm w_k}_{\bm x_W}
);

のように &math(V,W); の基底を合わせた線形結合として表せる。

すなわち、&math(V,W); の基底ベクトルすべてで「張られる」空間が和空間である。

ここから &math(\dim(V+W)\le\dim V+\dim W); が言える。

&math(V,W); の基底ベクトルを合わせたものが一次独立であるときに限り、
それがそのまま &math(V+W); の基底となるから、
&math(\dim(V+W)=\dim V+\dim W); となる。

** 和空間の次元 [#iadf0b68]

厳密な証明は省くが、

　&math(\dim (V+W)=\dim V+\dim W-\dim(V\cap W));

の関係がある。

これは、&math(V\cap W); の基底にいくつかベクトルを加えて &math(V); の基底を作成し、
同じ &math(V\cap W); の基底にいくつかベクトルを加えて &math(W); の基底を作成したならば、
それらすべてのベクトルを合わせると &math(V+W); の基底となる、という事実による。

#ref(2D_cap_2D.png,right,around,25%);
>例：
>
>右図の平面状の &math(V,W); の和空間は３次元空間全体となる。また２平面の交線が &math(V\cap W); に相当する。すなわち、
>
>&math(
\underbrace{\dim(V+W)}_3=\underbrace{\dim V}_2+\underbrace{\dim W}_2-\underbrace{\dim(V\cap W)}_1
);

#clear


* 直和 $V\dot +W$ [#u527a7fd]

上記より、&math(V\cap W=\{\bm 0\}); のとき、&math(\dim(V+W)=\dim V+\dim W); となる。

このとき「和空間 &math(V+W); は &math(V); と &math(W); の直和になっている」と言い、

　&math(V+W=V\dot +W); 

と書く。

- 直和は新たな演算ではない
- 「～～の場合に &math(V+W); は直和となる」「～～の場合には直和にならない」といった文脈で用いられる。
- &math(V); の基底と &math(W); の基底を合わせると、そのまま &math(V\dot +W); の基底になる~
↔ &math(\dim V+\dim W=\dim(V\dot +W));~

直和となるのは &math(V); と &math(W); （の基底）が一次独立なときである。

** 成分分解の一意性 [#o401c83d]

&math(V\dot +W); のときに限り &math(\bm x\mapsto\bm x_V); 
の射影演算を定義できる。

もし &math(\bm\delta\in V\dot+W); が &math(\bm\delta\ne\bm 0); であれば、

　&math(\bm x&=\bm x_V+\bm x_W\\
&=(\bm x_V+\bm \delta)+(\bm x_W-\bm \delta)\\
&=\bm x_V'+\bm x_W'
);

のように成分分解が一意に定まらないためだ。

** 線形独立な空間 [#h781440c]

直和は「線形独立な空間」の和空間のイメージになる。

　&math(\bm x_V+\bm x_W=\bm 0); から &math(\bm x_V=\bm x_W=\bm 0); を導ける。

** 成分の値はもう一方の空間に依存する [#e7949533]

成分分解のイメージは下図のようなものになる。

&attachref(1d_dplus_1d.png,,25%);

同じベクトル &math(\bm x); を~
　　&math(V); と &math(W); に分解したときの &math(\bm x_V); と、~
　　&math(V); と &math(W'); に分解したときの &math(\bm x'_V); とは~
一般には異なる値になる。

すなわち、ある部分空間の成分は、その部分空間だけでは決まらずに、他の部分空間の取り方にも依存する。

すなわち、上記の &math(P_{\bm e}); とは違って、&math(V); の情報のみから &math(\bm x_V); を求めることはできない。

&math(V); が２次元の時の成分分解のイメージは次の通り。

&attachref(2D-1D.png,,20%);

* 直交する空間 [#he5eb453]

&math(V); の任意の元が、
&math(W); の任意の元と直交するとき、
&math(V); と &math(W); とは直交すると言う。

&math(V); のすべての基底ベクトルが、~
&math(W); のすべての基底ベクトルと直交することと同義。

例えば $xy$ 平面からなる空間 $V$ と $yz$ 平面からなる空間 $W$ とは図形的には直交しているが、
$y$ 軸上のベクトル $\bm v=(0,1,0)$ は $\bm v\in V$ かつ $\bm v\in W$ であり、
当然 $\bm v\perp \bm v$ は成り立たないので、$V$ と $W$ は直交する空間とは呼ばない。

$V,W$ が直交する空間であれば、$V\cap W=\{\bm 0\}$ である。

* 直交直和 $V \oplus W$ [#d251a548]

２つの空間が直交する時、&math(V + W); を &math(V); と &math(W); の「直交直和」であるといい、

　&math(V+W=V \oplus W); 

と書く。

#ref(perpendicular.png,right,around,25%);

このとき、&math(V,W); の''正規直交''基底を合わせると &math(V \oplus W); の''正規直交''基底となる。

∵ &math(V); の正規直交基底が &math(W); の正規直交基底とも直交するから

当然、直交直和は直和でもある。

** 直交直和の成分分解 [#d2df2673]

直交直和の成分分解は簡単である。&math(\set{\bm v_k},\set{\bm w_k}); 
が正規直交基底であるとすると、

&math(
\bm x&=\underbrace{\sum_{k=1}^n c_{k}\bm v_{k}}_{\,\bm x_V}
      +\underbrace{\sum_{k=1}^m d_{k}\bm w_{k}}_{\,\bm x_W}\\
&=c_k\bm v_k+\underbrace{\sum_{k'\ne k}^n c_{k'}\bm v_{k'}+\sum_{k'=1}^m d_{k'}\bm w_{k'}}_{\perp\bm v_k}
);

のように、&math(\bm c_k\bm v_k); を取り出せば、残りの部分は &math(\bm v_k); と直交するから、
&math(\bm c_k\bm v_k); は &math(\bm x); の &math(\bm v_k); 方向成分である。

すると、上で見た任意の &math(\bm x); から &math(\bm v_k); 方向成分を取り出す１次元射影演算子
&math(P_{\bm v_k}); を使って &math(c_k\bm v_k=P_{\bm v_k}\bm x); と書けるから、

&math(
\bm x_V
&=\sum_{k=1}^n P_{\bm v_k}\bm x\\
&=\left(\sum_{k=1}^n P_{\bm v_k}\right)\bm x\\
&=P_V\bm x
);

すなわち、

　&math(P_V=\sum_{k=1}^n P_{\bm v_k});

が &math(V\oplus W); から &math(V); への射影演算子となる。

数ベクトルに対しては上で見たとおり

　&math(P_V=\sum_{k=1}^n \bm v_k\bm v_k^\dagger);

である。

射影演算子は &math(V); の情報だけから定まり、
&math(W); に依存しないことに注意せよ。

エルミート演算子の和はエルミート演算子になるから、
&math(P_{\bm v_k}); の和である &math(P_V); もエルミートである。

空間が直交しない一般の直和の場合にも「逆基底」を考えることにより、
直交直和の場合とほとんど同じように射影や成分分解が可能である。
→ 発展：[[線形代数II/非直交基底の成分分解]]

** 直交補空間 [#td0c1fcb]

#ref(complementary-space.svg,around,right);

全体空間 &math(U); が &math(U=V\oplus W); と表されるとき、
&math(W); を &math(V); の「直交補空間」と呼び、&math(W=V^\perp); と書く。

ある線形空間 &math(V); に対してその直交補空間は一意に定まる。

　&math(V^\perp=\set{\bm x\in U|\forall\bm y\in V,(\bm x,\bm y)=0});

つまり全体集合を、ある空間と、それに直交する補空間と、に分解することはいつも可能である。

あるベクトル &math(\bm x); を &math(\bm e); 
に平行な成分 &math(\bm x_\parallel); と垂直な成分 &math(\bm x_\perp); 
に分ける問題は、それぞれ線形空間 &math(V=\set{\bm p=t\bm e|t\in K}); 
とその直交補空間 &math(V^\perp); の成分への分解を表わしていたことになる。

一方、全体空間 &math(U); を &math(U=V\dot + W); と表せるとき、
&math(W); を &math(V); の（単なる）「補空間」と呼ぶ。
ある空間の直交補空間が一意に決まるのに対して、
補空間にはさまざまな取り方がある。


** 射影演算子の性質 [#ae98dccc]

- &math(\bm x\in V); のとき &math(P_V\bm x=\bm x);
- &math(\bm x\in V^\perp); のとき &math(P_V\bm x=\bm 0);
- &math(E=P_V+P_{V^\perp}); ← ∵&math(\bm x=P_V\bm x+P_{V^\perp}\bm x=\bm x_\parallel+\bm x_\perp);
- &math(P_V^2=P_V); あるいは &math(P_V(E-P_V)=O);~
∵ &math(P_V\bm x\in V); だから、&math(P_V^2\bm x=P_V\bm x);
- これは、&math(P_{V^\perp}=E-P_V); であり、&math(P_VP_{V^\perp}=O); であることからも理解できる

** 例 [#mfbb4345]

&math(\mathbb R^3); の部分空間として
&math(\bm a=\begin{pmatrix}1\\2\\3\end{pmatrix},\bm b=\begin{pmatrix}-1\\0\\1\end{pmatrix}); で張られる空間 &math(V=\big[\bm a,\bm b\big]\subset \mathbb R^3);を考える。

(1) &math(\mathbb R^3); から &math(V); への射影演算子を求めよ。

(2) 直交補空間 &math(V^\perp); に正規直交基底を定めよ。

*** 解答 (1) [#b1fdc171]

&math(\bm a,\bm b); から正規直交基底を作る。

&math(\bm b); と垂直なのは &math(\begin{pmatrix}s\\t\\s\end{pmatrix}); 
の形のベクトルであることに注意して、
&math(\bm c=\bm a-\bm b=\begin{pmatrix}2\\2\\2\end{pmatrix}); とすれば 
これは &math(V); 内にあり &math(\bm b); と垂直なベクトルである。

これらを正規化すれば、

　&math(\bm e_1=\frac{1}{\sqrt 2}\begin{pmatrix}-1\\0\\1\end{pmatrix}, 
\bm e_2=\frac{1}{\sqrt 3}\begin{pmatrix}1\\1\\1\end{pmatrix});

として正規直交基底が得られる。

したがって、求める射影演算子は

&math(
P_V&=\bm e_1\bm e_1^\dagger+\bm e_2\bm e_2^\dagger\\
&=\frac{1}{2}\begin{pmatrix}-1\\0\\1\end{pmatrix}\begin{pmatrix}-1&0&1\end{pmatrix}
+ \frac{1}{3}\begin{pmatrix}1\\1\\1\end{pmatrix}\begin{pmatrix}1&1&1\end{pmatrix}\\
&=\frac{1}{2}\begin{pmatrix}1&0&-1\\0&0&0\\-1&0&1\end{pmatrix}
+ \frac{1}{3}\begin{pmatrix}1&1&1\\1&1&1\\1&1&1\end{pmatrix}\\
&=\frac{1}{6}\begin{pmatrix}5&2&-1\\2&2&2\\-1&2&5\end{pmatrix}
);

各射影演算子がエルミート（実数行列では対称）になっていることにも注目せよ。

~

*** 解答 (1) 別解 [#r5bbcc18]

&math(\bm b,\bm a); からシュミットの直交化を用いて正規直交系を作る。

&math(\bm e_1=\frac{1}{\|\bm a\|}\bm a=\frac{1}{\sqrt 14}\begin{pmatrix}1\\2\\3\end{pmatrix});

&math(
\bm f_2
&=\bm b-(\bm e_1,\bm b)\bm e_1\\
&=\begin{pmatrix}-1\\0\\1\end{pmatrix}
-\frac{1}{14}\cdot 2\cdot\begin{pmatrix}1\\2\\3\end{pmatrix}\\
&=\frac{1}{7}\begin{pmatrix}-8\\-2\\4\end{pmatrix}
=\frac{2}{7}\begin{pmatrix}-4\\-1\\2\end{pmatrix}
);

&math(
\bm e_2=\frac{1}{\|\bm f_2\|}\bm f_2=\frac{1}{\sqrt{21}}\begin{pmatrix}-4\\-1\\2\end{pmatrix}
);

&math(
P_V&=
\frac{1}{14}\begin{pmatrix}1\\2\\3\end{pmatrix}\begin{pmatrix}1&2&3\end{pmatrix}+
\frac{1}{21}\begin{pmatrix}-4\\-1\\2\end{pmatrix}\begin{pmatrix}-4&-1&2\end{pmatrix}\\
&=\frac{1}{14}\begin{pmatrix}1&2&3\\&4&6\\&&9\end{pmatrix}+
\frac{1}{21}\begin{pmatrix}16&4&-8\\&1&-2\\&&4\end{pmatrix}\\
&=\frac{1}{42}\begin{pmatrix}3+32&6+8&9-16\\&12+2&18-4\\&&27+8\end{pmatrix}\\
&=\frac{1}{42}\begin{pmatrix}35&14&-7\\&14&14\\&&35\end{pmatrix}\\
&=\frac{1}{6}\begin{pmatrix}5&2&-1\\2&2&2\\-1&2&5\end{pmatrix}\\
);

射影演算子はエルミートになるため、左下部分の計算は省略した。

&math(P_V); の形は正規直交基底の取り方によらないことに注目せよ。

*** 解答 (2) [#n59bccf5]
 
&math(\mathbb R^3); が３次元、&math(V); が２次元なので、&math(V^\perp); は１次元となる。

&math(\bm e_1,\bm e_2); に垂直なベクトルを１つ挙げれば例えば、&math(\begin{pmatrix}1\\-2\\1\end{pmatrix}); 

したがって、

　&math(V^\perp=\Big[\begin{pmatrix}1\\-2\\1\end{pmatrix}\Big]);

である。正規直交基底はこれを正規化して、

　&math(\Big\{\frac{1}{\sqrt 6}\begin{pmatrix}1\\-2\\1\end{pmatrix}\Big\});

このとき、

　&math(
P_{V^\perp}&=\frac{1}{6}\begin{pmatrix}1\\-2\\1\end{pmatrix}\begin{pmatrix}1&-2&1\end{pmatrix}\\
&=\frac{1}{6}\begin{pmatrix}1&-2&1\\-2&4&-2\\1&-2&1\end{pmatrix}
);

であり、&math(P_V+P_{V^\perp}=E); となることが確かめられる。

** 演習 [#tac0cf13]

３次元空間に原点を通る平面 &math(x+y+z=0); を考える。
この平面への射影演算子を求めよ。
またその直交補空間を求めよ。

** 解答例 [#kee7e633]

*** まず平面内に基底を取る [#l8835428]

この「原点を通る平面」は２次元部分空間となるから、
平面内に２つの一次独立なベクトルを取れば、
それが平面に対応する線形空間の基底となる。

&math(x+y+z=0); を満たせば良いから、例えば、 
&math(\begin{pmatrix}
1\\-1\\0
\end{pmatrix},
\begin{pmatrix}
0\\-1\\1
\end{pmatrix}
); 
など、条件を満たすベクトルを「&ruby(めのこ){目の子};」で探しても良いが、
どんな場合にも通用する一般的なやり方を使うなら、

&math(x+y+z=0); を 
&math(\begin{pmatrix}
1&1&1
\end{pmatrix}\begin{pmatrix}
x\\y\\z
\end{pmatrix}=\bm 0
); の形に書いて、掃出し法により係数行列を階段化する。

今の場合は元の
&math(\begin{pmatrix}
1&1&1
\end{pmatrix}); がすでに階段行列であり、１列目は掃出しの完了した形になっているから、
掃出しの行えなかった列に対応する &math(y,z); をパラメータと見て、

　&math(\begin{cases}
x=-y-z\\
y=y\\
z=z
\end{cases});

すなわち、

　&math(\bm x=y\begin{pmatrix}
-1\\1\\0
\end{pmatrix}+z\begin{pmatrix}
-1\\0\\1
\end{pmatrix});

とすれば、２つのベクトル &math(\bm b_1=\begin{pmatrix}
-1\\1\\0
\end{pmatrix}, \bm b_2=\begin{pmatrix}
-1\\0\\1
\end{pmatrix}); がこの空間の基底となることが明らかである。

*** 基底を正規直交化する [#me88ef96]

これらを直交化するのも暗算で行っても良いが、
シュミットの直交化を使えばどんな場合にも必ず実行できて、

　&math(\bm f_1=b_1=\begin{pmatrix}
-1\\1\\0
\end{pmatrix});

　&math(\bm e_1=\frac{\bm f_1}{\|\bm f_1\|}=\frac{1}{\sqrt 2}\begin{pmatrix}
-1\\1\\0
\end{pmatrix});

　&math(
\bm f_2&=\bm b_2-\big(\bm e_1,\bm b_2\big)\bm e_1\\
&=\begin{pmatrix}
-1\\0\\1
\end{pmatrix}-\frac{1}{2}\underbrace{\overline{\begin{pmatrix}
-1&1&0
\end{pmatrix}}\begin{pmatrix}
-1\\0\\1
\end{pmatrix}}_{=\,1}\cdot\begin{pmatrix}
-1\\1\\0
\end{pmatrix}\\
&=\frac{1}{2}\begin{pmatrix}
-1\\-1\\2
\end{pmatrix}
);

　&math(
\bm e_2=\frac{\bm f_2}{\|\bm f_2\|}=\frac{1}{\sqrt 6}\begin{pmatrix}
-1\\-1\\2
\end{pmatrix}
);

途中で、転置されたベクトルの上に線が引いてあるのは
複素共役を取る演算であるが、ここでは実ベクトルなので値は変わらない。

したがって、&math(V); の正規直交基底は、

　&math(\set{\bm e_1,\bm e_2}=\Bigg\{\ 
\frac{1}{\sqrt 2}\begin{pmatrix}
-1\\1\\0
\end{pmatrix},\ \frac{1}{\sqrt 6}\begin{pmatrix}
-1\\-1\\2
\end{pmatrix}\ \Bigg\}
);

*** 正規直交基底から射影演算子を作る [#hc1e08df]

　&math(
P_V&=\bm e_1\bm e_1^\dagger+\bm e_2\bm e_2^\dagger=
\frac{1}{2}\begin{pmatrix}
-1\\1\\0
\end{pmatrix}\overline{\begin{pmatrix}
-1&1&0
\end{pmatrix}}+
\frac{1}{6}\begin{pmatrix}
-1\\-1\\2
\end{pmatrix}\overline{\begin{pmatrix}
-1&-1&2
\end{pmatrix}}\\
&=\frac{1}{2}\begin{pmatrix}
1&-1&0\\
-1&1&0\\
0&0&0\\
\end{pmatrix}+
\frac{1}{6}\begin{pmatrix}
1&1&-2\\
1&1&-2\\
-2&-2&4\\
\end{pmatrix}=
\frac{1}{6}\begin{pmatrix}
4&-2&-2\\
-2&4&-2\\
-2&-2&4\\
\end{pmatrix}\\
&=
\frac{1}{3}\begin{pmatrix}
2&-1&-1\\
-1&2&-1\\
-1&-1&2\\
\end{pmatrix}
);

この &math(P_V); に任意の &math(\bm x); をかければ、

　&math(
P_V\begin{pmatrix}
x\\y\\z
\end{pmatrix}=
\frac{x}{3}\begin{pmatrix}
2\\-1\\-1
\end{pmatrix}+
\frac{y}{3}\begin{pmatrix}
-1\\2\\-1
\end{pmatrix}+
\frac{z}{3}\begin{pmatrix}
-1\\-1\\2
\end{pmatrix}
);

となるが、右辺に現れる３つのベクトルはすべて &math(x+y+z=0); を満たしており、
確かに &math(P_V\bm x\in V); となることが確認できる。

*** 直交補空間を見つける [#r2a7e72f]

直交補空間 &math(V^\perp); は、&math(V); の任意の元と直交するベクトルを集めた集合である。

　&math(
V^\perp=\set{\bm x|\forall\bm y\in V,(\bm x,\bm y)=0}
);

&math(\bm y\in V); は &math(\bm y=y_1\bm e_1+y_2\bm e_2); と表せるから、

&math((\bm x,\bm y)=0); は &math(y_1(\bm x,\bm e_1)+y_2(\bm x,\bm e_2)=0); を表し、
任意の &math(\bm y); すなわち任意の &math(y_1,y_2); についてこれが成り立つには、

　&math((\bm x,\bm e_1)=(\bm x,\bm e_2)=0);

が必要十分条件となる。すなわち、

　&math(
V^\perp=\set{\bm x|\bm x\perp\bm e_1\ \mathrm{and}\ \bm x\perp\bm e_2}
);

&math(V); のすべての基底と直交するベクトルを集めた集合が &math(V^\perp); である。

　&math(\begin{cases}
-x+y=0\\
-x-y+2z=0
\end{cases});

の係数行列を同値変形して、

　&math(\begin{pmatrix}
-1&1&0\\
-1&-1&2
\end{pmatrix}\sim
\begin{pmatrix}
1&-1&0\\
0&-2&2
\end{pmatrix}\sim
\begin{pmatrix}
1&-1&0\\
0&1&-1
\end{pmatrix}\sim
\begin{pmatrix}
1&0&-1\\
0&1&-1
\end{pmatrix}
);

すなわち、&math(\begin{cases}
x-z=0\\
y-z=0
\end{cases});

掃出せなかった列に対応する &math(z); をパラメータとすれば、

　&math(\bm x=z\begin{pmatrix}
1\\1\\1
\end{pmatrix}
);

すなわち &math(V^\perp); の正規直交基底は 
&math(\Bigg\{\ \frac{1}{\sqrt 3}\begin{pmatrix}
1\\1\\1
\end{pmatrix}\ \Bigg\});

そもそも &math(V); を与える条件式
&math(x+y+z=0); は 
&math(
\Big(\ \bm x,\ \begin{pmatrix}
1\\1\\1
\end{pmatrix}\ 
\Big)=0
); という条件であるから、

&math(V^\perp); が
&math(\begin{pmatrix}
1\\1\\1
\end{pmatrix}); に平行な１次元空間となることは当然のことである。

* 一般化 [#jaf6d678]

以上の話は２つ以上の部分空間がある場合にも拡張できて、以下の通りである。

|交空間  |&math(V_1\cap V_2\cap \dots\cap V_r);      |全空間の共通部分|
|和空間  |&math(V_1+V_2+\dots+V_r);                  |一般には一次従属な空間たちを内包する空間|
|直和    |&math(V_1\dot +V_2\dot +\dots\dot +V_r);   |一次独立な空間たちの和空間|
|直交直和|&math(V_1\oplus V_2\oplus \dots\oplus V_r);|直交する空間たちの和空間  |

たとえば &math(V_1\cap V_2\cap V_3\cap V_4=(((V_1\cap V_2)\cap V_3)\cap V_4); などの意味であるが、これらの演算子には結合法則や交換法則が成り立ち、&math((V_1\oplus V_2)\oplus V_3=V_1\oplus (V_2\oplus V_3));, &math(V_1\oplus V_2=V_2\oplus V_1); などとなる。

&math(\cap); や &math(+); が複数の部分空間から新しい部分空間を作る演算子であるのに比べて、
&math(\dot +); や &math(\oplus); は 「線形空間同士の演算」 ではなく、
和空間を形成する空間が特殊な条件を満たすことを表現しているに過ぎない。

この違いに注意せよ。



* 質問・コメント [#m8a71c0e]

#article_kcaptcha
Counter: 135881 (from 2010/06/03), today: 6, yesterday: 0