メインコンテンツへスキップ

不変部分空間と不変部分空間に関する線形変換の表現行列、不変部分空間の視覚化

··4937 文字·10 分·
数学 線形代数 不変部分空間 表現行列 線形変換 線形写像 直和
著者
gf
目次

不変部分空間(invariant subspace)とは
#

ベクトル空間VVと線形変換T:VVT: V \to Vについて、VVの部分空間WWTTに関する不変部分空間である(TT-不変, TT-invariant)とは、 {T(w) wW}W \{T(w)\thinspace|\space \forall w \in W\} \subset W を満たすことをいう。

日本語で言うと、部分空間WWのどのベクトルw\bm{w}を線形変換TTで写しても、行き先はWW内に留まるということだ。

定義から分かるとおり、ベクトル空間VVと線形変換T:VVT: V \to Vについて、VV自身と{0}\{\bm{0}\}は常にTTに関する不変部分空間となる。

そして、Ker T, Im T\mathrm{Ker} \space T, \space \mathrm{Im} \space Tも不変部分空間になる。

補足(クリックで展開)

Ker T\mathrm{Ker} \space T のベクトルはTTで写すと全て0Ker T\bm{0} \in \mathrm{Ker} \space Tに写る。

Im T\mathrm{Im} \space Tに関しては、いま、TTの定義域と値域はVVであることに注意。

Im T\mathrm{Im} \space TのベクトルはTTで写すと、当たり前だが、Im T\mathrm{Im} \space T内に写る。

不変部分空間に関する線形変換の表現行列
#

V, W, TV, \space W, \space Tなどの記号は上の項のものをそのまま引用する。
VVの次元はnnWWの次元はmmとする。

WWの基底をw1, ,wm\bm{w}_1, \space \dots, \bm{w}_mとする。

この基底を基底の延長定理で延長して、VVの基底w1, , wm, v1, ,  vnm\bm{w}_1, \space \dots, \space \bm{w}_m, \space \bm{v}_{1} , \space \dots, \space \space \bm{v}_{n-m}を得る。

いま、WWTT-不変なので、w1, ,wm\bm{w}_1, \space \dots, \bm{w}_mをそれぞれTTで写したものは、WWの基底の一次結合で表される。
よって、 T(wi)=a1iw1+a2iw2++amiwm(p)(axyは適当なスカラー, i=1  m) T(\bm{w}_i) = a_{1i} \bm{w}_1 + a_{2i}\bm{w}_2 + \dots + a_{mi} \bm{w}_m \quad(p) \\ (a_{xy}は適当なスカラー, \space i = 1 \space \dots \space m) である。

他方、TTの定義より、VVのベクトルw1, , wm, v1, ,  vnm\bm{w}_1, \space \dots, \space \bm{w}_m, \space \bm{v}_{1} , \space \dots, \space \space \bm{v}_{n-m} をそれぞれTTで写したものは、VVの基底の一次結合で表される。

よって、先程WWの基底を延長して得たVVの基底を用いて、 T(wi)=a1iw1+a2iw2++amiwm+am+1 iv1++anm ivnm(q)T(vj)=a1jw1+a2jw2++amjwm+am+1 jv1++anm jvnm(i=1  m, j=m+1  nm) \small T(\bm{w}_{i}) = a_{1i} \bm{w}_1 + a_{2i}\bm{w}_2 + \dots + a_{mi} \bm{w}_m + a_{m+1 \space i}\bm{v}_{1} + \dots + a_{n-m \space i}\bm{v}_{n-m} \quad (q) \\ T(\bm{v}_{j}) = a_{1j} \bm{w}_1 + a_{2j}\bm{w}_2 + \dots + a_{mj} \bm{w}_m + a_{m+1 \space j}\bm{v}_{1} + \dots + a_{n-m \space j}\bm{v}_{n-m} \\ (i = 1 \space \dots \space m, \space j = m +1 \space \dots \space n-m) と表せる。

この関係は行列を用いてまとめて以下のように表すことができる。
(T(w1)T(w2)T(wm)T(v1)T(vnm))=(w1wmv1vnm)(a11a12a1ma1 m+1a1 nma21a22a2ma2 m+1a2 nmam1am2ammam m+1am nmam+1 1am+1 2am+1 mam+1 m+1am+1 nmanm 1anm 2anm manm m+1anm nm) \tiny \begin{pmatrix} T(\bm{w}_{1}) & T(\bm{w}_{2}) & \dots & T(\bm{w}_{m}) & T(\bm{v}_{1}) & \dots & T(\bm{v}_{n-m}) \end{pmatrix} \\ = \begin{pmatrix} \bm{w}_1 & \dots & \bm{w}_m & \bm{v}_{1} & \dots & \bm{v}_{n-m} \end{pmatrix} \begin{pmatrix} a_{11} & a_{12} & \dots & a_{1m} & a_{1 \space m+1} & \dots & a_{1 \space n-m} \\ a_{21} & a_{22} & \dots & a_{2m} & a_{2 \space m+1} & \dots & a_{2 \space n-m} \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ a_{m1} & a_{m2} & \dots & a_{mm} & a_{m \space m+1} & \dots & a_{m \space n-m} \\ a_{m+1 \space 1} & a_{m+1 \space 2} & \dots & a_{m+1 \space m} & a_{m+1 \space m+1} & \dots & a_{m+1 \space n-m} \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ a_{n-m \space 1} & a_{n-m \space 2} & \dots & a_{n-m \space m} & a_{n-m \space m+1} & \dots & a_{n-m \space n-m} \end{pmatrix}

この行列は線形変換(線形写像)TTの表現行列である。

🔗関連記事: 抽象ベクトル空間の間の線形写像とその表現行列

ここで、(q)(q)と基底(一次独立なベクトルの組)の一次結合の一意性と、(p)(p)に関して、 (p)=a1iw1+a2iw2++amiwm+0 (p) = a_{1i} \bm{w}_1 + a_{2i}\bm{w}_2 + \dots + a_{mi} \bm{w}_m + \bm{0} であることより、 T(wi)=(p)=a1iw1+a2iw2++amiwm+0=a1iw1+a2iw2++amiwm+am+1 iv1++anm ivnm=(q)=T(wi) T(\bm{w}_{i}) = (p) = a_{1i} \bm{w}_1 + a_{2i}\bm{w}_2 + \dots + a_{mi} \bm{w}_m + \bm{0} \\ = a_{1i} \bm{w}_1 + a_{2i}\bm{w}_2 + \dots + a_{mi} \bm{w}_m + a_{m+1 \space i}\bm{v}_{1} + \dots + a_{n-m \space i}\bm{v}_{n-m} \\ =(q)=T(\bm{w}_{i}) となり、ここで登場するベクトルは一次独立なので、ゼロベクトルは存在しない。

つまり、am+1 ianm i(i=1  m)a_{m+1 \space i} \dots a_{n-m \space i} \quad (i = 1 \space \dots \space m)は全て0であることが判明する。

よって、TTの表現行列は基底をうまく取ると以下のようになる。 (a11a12a1ma1 m+1a1 nma21a22a2ma2 m+1a2 nmam1am2ammam m+1am nmam+1 m+1am+1 nmOanm m+1anm nm) \begin{pmatrix} a_{11} & a_{12} & \dots & a_{1m} & a_{1 \space m+1} & \dots & a_{1 \space n-m} \\ a_{21} & a_{22} & \dots & a_{2m} & a_{2 \space m+1} & \dots & a_{2 \space n-m} \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ a_{m1} & a_{m2} & \dots & a_{mm} & a_{m \space m+1} & \dots & a_{m \space n-m} \\ & & & & a_{m+1 \space m+1} & \dots & a_{m+1 \space n-m} \\ & & \huge{O} & & \vdots & \vdots & \vdots \\ & & & & a_{n-m \space m+1} & \dots & a_{n-m \space n-m} \end{pmatrix}

線形変換の定義域を不変部分空間に制限した際の表現行列
#

先程の線形変換T:VVT: V \to Vの定義域を不変部分空間WWに制限(restriction)すると、

restrictionについて
🔗https://en.wikipedia.org/wiki/Restriction_(mathematics)
日本語版: 🔗https://ja.wikipedia.org/wiki/%E5%88%B6%E9%99%90_(%E6%95%B0%E5%AD%A6)

TW:WVT|_{W}: W \to Vを得る。

いま、WWTT不変なので、不変部分空間の定義より、TW:WWT|_{W}: W \to Wとなることがわかる。

TWT|_{W}の表現行列はどのような形になっているのだろうか?

実は、先程の議論でもうほとんど答えが分かっている。 (T(w1)T(w2)T(wm)T(v1)T(vnm))=(w1wmv1vnm)(a11a12a1ma1 m+1a1 nma21a22a2ma2 m+1a2 nmam1am2ammam m+1am nmam+1 m+1am+1 nmOanm m+1anm nm) \tiny \begin{pmatrix} T(\bm{w}_{1}) & T(\bm{w}_{2}) & \dots & T(\bm{w}_{m}) & T(\bm{v}_{1}) & \dots & T(\bm{v}_{n-m}) \end{pmatrix} \\ = \begin{pmatrix} \bm{w}_1 & \dots & \bm{w}_m & \bm{v}_{1} & \dots & \bm{v}_{n-m} \end{pmatrix} \begin{pmatrix} a_{11} & a_{12} & \dots & a_{1m} & a_{1 \space m+1} & \dots & a_{1 \space n-m} \\ a_{21} & a_{22} & \dots & a_{2m} & a_{2 \space m+1} & \dots & a_{2 \space n-m} \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ a_{m1} & a_{m2} & \dots & a_{mm} & a_{m \space m+1} & \dots & a_{m \space n-m} \\ & & & & a_{m+1 \space m+1} & \dots & a_{m+1 \space n-m} \\ & & \huge{O} & & \vdots & \vdots & \vdots \\ & & & & a_{n-m \space m+1} & \dots & a_{n-m \space n-m} \end{pmatrix} の右辺の積を行列のmm列目までだけ計算すると、 T(w1)=a11w1+a21w2++am1wm+0v1++0vnmT(w2)=a12w1+a22w2++am2wm+0v1++0vnmT(wm)=a1mw1+a2mw2++ammwm+0v1++0vnm T(\bm{w}_1) = a_{11} \bm{w}_1 + a_{21} \bm{w}_2 + \dots + a_{m1} \bm{w}_{m} + 0 \bm{v}_{1} + \dots + 0 \bm{v}_{n-m} \\ T(\bm{w}_2) = a_{12} \bm{w}_1 + a_{22} \bm{w}_2 + \dots + a_{m2} \bm{w}_{m} + 0 \bm{v}_{1} + \dots + 0 \bm{v}_{n-m} \\ \vdots \\ T(\bm{w}_m) = a_{1m} \bm{w}_1 + a_{2m} \bm{w}_2 + \dots + a_{mm} \bm{w}_{m} + 0 \bm{v}_{1} + \dots + 0 \bm{v}_{n-m} \\ を得る。

この対応は、TTの表現行列内の左上のm×mm \times mの行列はTWT|_{W}の表現行列であることを示している。

先程の議論内でも暗に判明していたが、同時に、不変部分空間WWの元はVVの基底の一次結合で表すと、 k1w1+k2w2++kmwm+0v1++0vnm k_1 \bm{w}_1 + k_2 \bm{w}_2 + \dots + k_m \bm{w}_{m} + 0 \bm{v}_{1} + \dots + 0 \bm{v}_{n-m} の形になることも分かる。

まとめると、結局、TTの表現行列ATA_{T}基底をうまく取ると AT=(ATWBm×nmOnm×mCnm×nm) A_{T}= \begin{pmatrix} A_{T|_{W}} & B_{m \times n-m} \\ O_{n-m \times m} & C_{n-m \times n-m} \end{pmatrix} の形になる事がわかる。

不変部分空間の直和と線形変換の表現行列
#

更に、線形変換T:VVT: V \to V に関する不変部分空間X, YX, \space Yに関して VVXXYYの直和で表されるとき、つまり、 V=XY V = X \oplus Y のとき、TTの表現行列は基底をうまく取ると AT=(ATXOOATY) A_{T}= \begin{pmatrix} A_{T|_X} & O \\ O & A_{T|_Y} \end{pmatrix} の形になる。

なぜなら直和であるということはVVの任意のベクトルv\bm{v}v=x+y(xX, yY) \bm{v} = \bm{x} + \bm{y} \quad (\bm{x} \in X, \space \bm{y} \in Y) の形で一意的に表されるということであり、XX, YYはそれぞれ不変部分空間なので、 先程の項で出てきたVVの基底をXXの基底にYYの基底を付け加えた x1, , xdimX, y1, , ydimY \bm{x}_1 , \space \dots, \space \bm{x}_{\mathrm{dim} X}, \space \bm{y}_{1} , \space \dots , \space \bm{y}_{\mathrm{dim} Y} で考えることによって、上の基底のうちのYYのベクトルをそれぞれTTで写したものであるT(yj)T(\bm{y}_{j})YYのベクトルであると同時にVVの基底の一次結合で表されることを考えると、 T(yj)=a1 jx1++adimX jxdimX+adimX+1 jy1++adimVdimX jydimY=0x1++0xdimX+adimX+1 jy1++adimVdimX jydimYj=(dimX+1, , dimY) \small T(\bm{y}_{j}) = a_{1 \space j} \bm{x}_1 + \dots + a_{\mathrm{dim}X \space j} \bm{x}_{\mathrm{dim} X} + a_{\mathrm{dim} X + 1 \space j} \bm{y}_{1} + \dots + a_{\mathrm{dim} V - \mathrm{dim} X \space j} \bm{y}_{\mathrm{dim} Y} \\ = 0 \bm{x}_1 + \dots + 0 \bm{x}_{\mathrm{dim} X} + a_{\mathrm{dim} X + 1 \space j} \bm{y}_{1} + \dots + a_{\mathrm{dim} V - \mathrm{dim} X \space j} \bm{y}_{\mathrm{dim} Y} \\ j = ( \mathrm{dim} X +1, \space \dots, \space \mathrm{dim} Y) と表されるので、結果的にTTの表現行列の右上側もOOになることが分かる。

添字がごちゃごちゃしていて分かりにくいので分からなければ他のリソースも参照することをおすすめする。

[おまけ] 不変部分空間の具体例の一例とその視覚化
#

ある平面でベクトルを反射させるような線形写像に関して、不変部分空間を観察してみる。 このような線形写像の行列はcomputer graphicsではreflection行列と呼ばれる。

ベクトル空間VVR3\Reals^3とする。

A=(100010001) A = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & -1 \end{pmatrix} としておく。

Vの基底 (100),(010),(001) \begin{pmatrix} 1 \\ 0 \\ 0 \end{pmatrix}, \begin{pmatrix} 0 \\ 1 \\ 0 \end{pmatrix}, \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix} VVのベクトルは

((abc)V)=a(100)+b(010)+c(001) (\begin{pmatrix} a \\ b \\ c \end{pmatrix} \in V) = a \begin{pmatrix} 1 \\ 0 \\ 0 \end{pmatrix} + b \begin{pmatrix} 0 \\ 1 \\ 0 \end{pmatrix} + c \begin{pmatrix} 0 \\ 0 \\ 1 \end{pmatrix} と基底の一次結合の形で一意的に表すことができる。

このとき、成分ベクトルは(abc)\begin{pmatrix} a \\ b \\ c \end{pmatrix}である。

一般のベクトル空間上の線形写像は基底を固定した際の、定義域のベクトルの成分と値域のベクトルの成分との間の関係であり、その関係は表現行列で表すことができた。

Aで成分ベクトルがどのように写るか見てみよう。
A (abc)=(100010001)(abc)=(abc) A \begin{pmatrix} a \\ b \\ c \end{pmatrix} = \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & -1 \end{pmatrix} \begin{pmatrix} a \\ b \\ c \end{pmatrix} = \begin{pmatrix} a \\ b \\ -c \end{pmatrix} と写ることが分かる。

この計算結果の成分ベクトルと先程の基底の一次結合でVVのベクトルを表すことができるので、結局、 VVのベクトル(abc)\begin{pmatrix} a \\ b \\ c \end{pmatrix}VVのベクトル(abc)\begin{pmatrix} a \\ b \\ -c \end{pmatrix}に写ることがわかる。

VVのベクトル(xyz)\begin{pmatrix} x \\ y \\ z \end{pmatrix}を3次元空間の点x, y, zx,\space y, \space zに対応させて考えてみる。

一応、理屈を説明しておこう。

VVのベクトル(xyz)\begin{pmatrix} x \\ y \\ z \end{pmatrix}に対して、3次元空間上の点(x, y, z)(x, \space y, \space z)を対応させる。
続いて、全てのスカラー倍k(xyz)k \begin{pmatrix} x \\ y \\ z \end{pmatrix}も写すと、点の集まりは3次元空間上の原点と(x, y, z)(x, \space y, \space z)の2点を通る直線になることがわかる。

以上を踏まえて、いくつかのVVの元を視覚化すると以下の動画のようになる。

VVのあるベクトルv1のスカラー倍を全て3次元空間上の点に対応させたものが緑色の直線で、v1スカラー倍をそれぞれ行列AAで写した後に3次元空間上の点に対応させたものが黄色の直線v1 reflectedだ。

v1v1 reflectedではz座標の符号が逆なので対称的になっていることに注目。

青(v2)と赤(v2 reflected)の直線も同様の関係にある。

ここで、AAで定まる線形写像fA:VVf_A: V \to Vに関する2つの不変部分空間が視覚化された様子が以下の動画だ。

動画の概要の方にテイストの違う説明を書いたので分からなければそちらも合わせて参照してみてほしい。

黒い平面に含まれる任意の点に着目すると、z座標の値が00になっていることがわかる。
つまり、その点に対応するVVの元は(xy0)\begin{pmatrix} x \\ y \\ 0 \end{pmatrix}の形になっていて、 (100010001)(xy0)=(xy0) \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & -1 \end{pmatrix} \begin{pmatrix} x \\ y \\ 0 \end{pmatrix} = \begin{pmatrix} x \\ y \\ 0 \end{pmatrix} と線形変換(reflection)の影響を受けないことがわかる。

つまり、(xy0)\begin{pmatrix} x \\ y \\ 0 \end{pmatrix}の形のベクトルの全体は不変部分空間であり、その視覚化が動画内の黒い平面である。

続いて、先程の黒い平面に垂直な、原点を通るオレンジ色の直線に注目すると、対応するVVの元は(00z)\begin{pmatrix} 0 \\ 0 \\ z \end{pmatrix}の形になっており、 (100010001)(00z)=(00z) \begin{pmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & -1 \end{pmatrix} \begin{pmatrix} 0 \\ 0 \\ z \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \\ -z \end{pmatrix} となることがわかる。

reflectionの影響を受けているが、VVの部分空間、(00z)\begin{pmatrix} 0 \\ 0 \\ z \end{pmatrix}全体(オレンジ色の直線)を考えると、線形変換しても行き先は部分空間をはみ出ないことがわかる。

よって、黒い平面とオレンジ色の直線はVVの不変部分空間の一例が視覚化されたものであることが分かる。

編集後記
#

TTの表現行列の右下の行列Cnm×nmC_{n-m \times n-m}TV/W:V/WV/WT|_{V/W}: V/W \to V/Wの表現行列らしい。

T:VVT: V \to VからV/WV/WV/W \to V/W🔗universal propertyを利用してinduceすることはできたが、具体的にその写像の表現行列がCnm×nmC_{n-m \times n-m}と一致するかどうか分からなかった、示せなかったので割愛した。

もっとレベルが上がったら再チャレンジしようと思う。

関連記事

抽象ベクトル空間の間の線形写像とその表現行列
··3039 文字·7 分
数学 線形代数 ベクトル空間 線形写像 表現行列
一般のベクトル空間の間の線形写像もスカラーから成る行列で表現することができる。 すげえよ。
抽象ベクトル空間と部分空間
··5342 文字·11 分
数学 線形代数 ベクトル空間 部分空間 線形写像
俺自身が部分空間になることだ。
基底の取り替えと線形写像の表現行列の関係
··1535 文字·4 分
数学 線形代数 基底 表現行列 ベクトル空間 行列