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

更新


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

* 目次 [#xf03612b]

#contents
#mathjax

* ベクトルの成分 [#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
*** 注意1 [#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);~

*** 注意2 [#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}}); となってしまう。

*** 注意3 [#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(
(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}
);

のような &math(n); 次正方行列になる。

** 射影演算子はエルミートになる。 [#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); 次元空間への射影を考える。

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

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

- 任意の &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);

** 復習2:部分空間 [#t397a1fd]

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

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

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

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

** 復習3:集合の積と和 [#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); と読む。
#multicolumns
:和(結び)|
&math(A\cup B=\set{x|x\in A\,\mathrm{または}\,x\in B});~
&math(A); あるいは &math(B); の少なくとも片方に含まれる要素の集合~
&math(A); カップ &math(B); と読む。
#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]

2つの線形空間の交わり &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%);

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

&math(V,W); 上から2つのベクトルを取り
&math(\bm v\in V, \bm w\in W); とすれば、~
&math(\bm v, \bm w\in V\cap 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(V,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); の基底ベクトルすべてで「張られる」空間が和空間である。

** 和空間の次元 [#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); の和空間は3次元空間全体となる。また2平面の交線が &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));
** 成分分解の一意性 [#o401c83d]

任意の &math(\bm x\in V\dot +W); に対して、

 &math(\bm x=\bm x_V+\bm x_W\ \ \ \bm x_V\in V, \bm x_W\in W);

の分解は一意に定まる。

これは、&math(\bm \delta\in V\cap W); を使えば

 &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'
);

となり、&math(\bm x_V'\in V, \bm x_W'\in W); も &math(V); と &math(W); 
への分解を与えるが、直和であれば &math(\bm delta=\bm 0); に限るから、
のように異なる分解 &math(\bm x_V'\in V, \bm x_W'\in W); を作れるが、
直和であれば &math(\bm \delta=\bm 0); に限るから、
常に &math(\bm x_V=\bm x_V', \bm x_W=\bm x_W'); となる、ということ。

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

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

** 成分の値はもう一方の空間に依存する [#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); が2次元の時の成分分解のイメージは次の通り。

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

* 直交する空間 [#he5eb453]

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

&math(V); のすべての基底ベクトルが、~
&math(W); のすべての基底ベクトルと直交することと同義。
* 直交直和 $V \oplus W$ [#d251a548]

2つの空間が直交する時、&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); の正規直交基底とも直交するから
∵ &math(V); の正規直交基底が &math(W); の正規直交基底とも直交するから

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

直交直和の成分分解は簡単である。

&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});

に対して、&math(c_k=(\bm v_k,\bm x)); より、
に対して、&math(\bm c_k\bm v_k); は &math(\bm x); の &math(\bm v_k); 方向成分であるから、
上で見た、任意の &math(\bm x); から &math(\bm v_k); 方向成分を取り出す1次元射影演算子
&math(P_{\bm v_k}); を使って、

&math(
\bm x_V&=\sum_{k=1}^n (\bm v_k,\bm x)\bm v_k\\
\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\bm x=\sum_{k=1}^n P_{\bm v_k});

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

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

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

である。

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

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

** 直交補空間 [#td0c1fcb]

全体空間 &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]

(1)

&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) &math(\mathbb R^3); が3次元、&math(V); が2次元なので、&math(V^\perp); は1次元となる。
*** 解答 (2) [#n59bccf5]
 
&math(\mathbb R^3); が3次元、&math(V); が2次元なので、&math(V^\perp); は1次元となる。

&math(\bm e_1,\bm e_2); に垂直なベクトルを1つ挙げれば例えば、&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]

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

** 解答例 [#kee7e633]

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

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

&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}); がすでに階段行列であり、1列目は掃出しの完了した形になっているから、
掃出しの行えなかった列に対応する &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});

とすれば、2つのベクトル &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}
);

となるが、右辺に現れる3つのベクトルはすべて &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}); に平行な1次元空間となることは当然のことである。

* 一般化 [#jaf6d678]

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

|交空間  |&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: 117127 (from 2010/06/03), today: 10, yesterday: 0