広義固有空間の構造とジョルダン標準形 のバックアップ(No.9)

概要 †

現在推敲中で未完成です

ジョルダン標準形とは、対角化できない行列を「準対角化」した形である。 この標準形は行列の固有空間の概念を拡張した、「広義固有空間」が持つ構造を反映した形となる。

ここは発展項目なので、線形代数IIの内容を先取りして使う。
実のところ、線形代数IIでも扱わない内容なので、線形代数IIを学んでから戻ってきても良い。

固有空間とその次元 †

次正方行列 の、固有値 に対する固有空間 （固有ベクトルの集合が作る線形空間）とは、

　

の解が作る空間である。したがって、

　

と書ける。

すべての固有値について、その重複度と同数の一次独立な固有ベクトルを見つけられることが 行列が対角化可能な条件だった。

これは固有空間の次元が重複度と等しいことと同値である。 （重複度を超えることはない）

　

（ の重複度を とする）

広義固有空間 †

対角化できない場合には

　

となる固有値があるが、そのような場合にも必ず、

　

となることを以下のようにして示せる。すなわち、固有空間

　

に対して、「広義固有空間」を

　

と定義すれば、その次元は必ず重複度と等しくなる。

さらに、異なる固有値の広義固有空間はゼロ以外で重ならないため、 全空間を広義固有空間の直和に分解できることになる。

 & \ddots & * &\dots & \dots&* \\
 &  & \lambda & & & \vdots\\
 &  &  & \lambda'& &\vdots\\
 &  &  &  & \lambda''\\
 &  &  &  &  & \ddots\\

 & \ddots & * &\dots & \dots&* \\
 &  & 0 & & & \vdots\\
 &  &  & \lambda'-\lambda& &\vdots\\
 &  &  &  & \lambda''-\lambda\\
 &  &  &  &  & \ddots\\

 & \ddots & 0 &\dots & \dots&* \\
 &  & 0 & & & \vdots\\
 &  &  & (\lambda'-\lambda)^r& &\vdots\\
 &  &  &  & (\lambda''-\lambda)^r\\
 &  &  &  &  & \ddots\\

 & \ddots & * &\dots & \dots&* \\
 &  & \lambda-\lambda' & & & \vdots\\
 &  &  & 0& &\vdots\\
 &  &  &  & \lambda''-\lambda\\
 &  &  &  &  & \ddots\\

 & \ddots & * &\dots & \dots&* \\
 &  & (\lambda-\lambda')^{r'} & & & \vdots\\
 &  &  & 0& &\vdots\\
 &  &  &  & (\lambda''-\lambda')^{r'}\\
 &  &  &  &  & \ddots\\

広義固有空間の構造 †

定義の通り、広義固有空間に含まれるベクトルは、

　

を満たす。

そこでこの空間を、 「何回 をかけたらゼロになるか」で分類する。 これらの集合は部分空間となるから、

• ０回でゼロになる空間は
• １回でゼロになるのが通常の固有空間
• ２回かけたらゼロになる空間
• ３回かけたらゼロになる空間
• …
• 回かけたらゼロになる空間

とすると、以下が成り立つ。

回かけたらゼロになるベクトルは 回かけたらゼロになるベクトルでもあるので、

　&math( \underbrace{V^0(\lambda)}_{=\,\set{\bm 0}}\subset \underbrace{V^1(\lambda)}_{=\,V(\lambda)}\subset V^2(\lambda)\dots\subset \underbrace{V^r(\lambda)}_{=\,W(\lambda)} );

から へと、 増えた分の空間を と書こう。つまり、

　

回ではゼロにならないが、 回で初めてゼロになるベクトル（とゼロベクトル）の集合が である。 （ の場合もある）

すると、 は の直和で書けて、

　&math( W(\lambda)= \underbrace{V^0(\lambda)}_{=\,\set{\bm 0}} \,\dot+\,\underbrace{\Delta V(\lambda)}_{=\,V(\lambda)} \,\dot+\,\Delta V^2(\lambda) \,\dot+\, \dots \,\dot+\, \Delta V^m(\lambda) \,\dot+\, \underbrace{\Delta V^{m+1}(\lambda) \,\dot+\, \dots \,\dot+\, \Delta V^r(\lambda)}_{\set{\bm 0}} );

一般に、 に対して の はゼロ次元 {0} になる。 これは m 回かけたらすべてゼロになってしまうため、 ということ。

一方、このとき はゼロ次元にならない。

なぜなら一般に、 に含まれるベクトルに をかけたベクトルは に属すから、 とすれば、

　

しかも であれば には１回かけてゼロになるベクトルは含まれないから、

　

つまり、

　

このとき、 の基底ベクトルに をかけたものに、必要に応じていくつかベクトルを加えて の基底を作れる。

つまり、全体空間は広義固有空間の直和に分解でき、

　

それぞれの広義固有空間は次のような構造をしていることがわかった。

　&math( W(\lambda)=&\Delta V^m(\lambda)\dot + \underbrace{T\big(\Delta V^m(\lambda)\big)\dot +{\Delta V^{m-1}}^*(\lambda)} _{=\,\Delta V^{m-1}(\lambda)}\dot + \underbrace{T\big(\Delta V^{m-1}(\lambda)\big)\dot +{\Delta V^{m-2}}^*(\lambda)}_{=\,\Delta V^{m-2}(\lambda)}\dot + \dots\dot + \underbrace{T\big(\Delta V^2(\lambda)\big)\dot +{\Delta V^{1}}^*(\lambda)}_{=\,V(\lambda)}\\ =&\underbrace{{\Delta V^m}^*(\lambda)}_{=\,\Delta V^m(\lambda)}\dot +\\ &\underbrace{T\big({\Delta V^m}^*(\lambda)\big)\dot+{\Delta V^{m-1}}^*(\lambda)}_{=\,\Delta V^{m-1}(\lambda)}\dot +\\ &\underbrace{T^2\big({\Delta V^m}^*(\lambda)\big)\dot+T\big({\Delta V^{m-1}}^*(\lambda)\big)\dot+{\Delta V^{m-2}}^*(\lambda)}_{=\,\Delta V^{m-2}(\lambda)}\dot +\\ &\cdots\\ &\underbrace{T^{m-1}\big({\Delta V^m}^*(\lambda)\big)\dot+T^{m-2}\big({\Delta V^{m-1}}^*(\lambda)\big)\dot+\dots\dot+{\Delta V^{1}}^*(\lambda)}_{=\,\Delta V^{1}(\lambda)\,=\,V(\lambda)}\\ );

例１ ３次元の広義固有空間、固有空間が１次元 †

　&math(A=\begin{pmatrix}

4 & -2 & 0 \\
-3 & 4 & -2 \\
-11 & 9 & -2 \\

\end{pmatrix}

);

まずは固有値を求める。

　&math(|A-\lambda I|= \begin{vmatrix}

4-\lambda & -2 & 0 \\
-3 & 4-\lambda & -2 \\
-11 & 9 & -2-\lambda \\

\end{vmatrix}= 8-12\lambda+6\lambda^2-\lambda^3=(2-\lambda)^3 );

　 （３重解）

したがって、広義固有空間 は３次元。

次に固有空間 を求める。

　

　&math( \begin{pmatrix}

2 & -2 & 0 \\
-3 & 2 & -2 \\
-11 & 9 & -4 \\

\end{pmatrix}\sim \begin{pmatrix}

1 & -1 & 0 \\
0 & -1 & -2 \\
0 & -2 & -4 \\

\end{pmatrix}\sim \begin{pmatrix}

1 & 0 & 2 \\
0 & 1 & 2 \\
0 & 0 & 0 \\

\end{pmatrix} );

掃出せなかった をパラメータとして、

　&math(\bm x=z\begin{pmatrix}

• 2\\-2\\1 \end{pmatrix});

したがって、固有空間 は１次元。 基底として &math(\bm b_1=\begin{pmatrix}

• 2\\-2\\1 \end{pmatrix}); を取れる。

広義固有空間が３次元で固有空間が１次元ということは、

　&math( \dim V(2)=1\ge\dim \Delta V^2(2)\ge\dim \Delta V^3(2)\ge\dim \Delta V^4(2)\ge\dots);

かつ、

　&math( \dim V(2)+\dim \Delta V^2(2)+\dim \Delta V^3(2)+\dim \Delta V^4(2)+\dots=\dim W(2)=3 );

であるから、

　&math(&\dim V(2)=\dim \Delta V^2(2)=\dim \Delta V^3(2)=1\\ &\dim\Delta V^4(2)=\dim\Delta V^5(2)=\dots=0);

となるはず。

を求めるため、 「 をかけると に入るベクトル」を探す。

　

の拡大係数行列を変形して、

　&math( \begin{pmatrix}

2 & -2 & 0 & -2\\
-3 & 2 & -2 & -2\\
-11 & 9 & -4 & 1\\

\end{pmatrix}\sim \begin{pmatrix}

1 & -1 & 0 & -1\\
0 & -1 & -2 & -5\\
0 & -2 & -4 & -10\\

\end{pmatrix}\sim \begin{pmatrix}

1 & 0 & 2 & 4\\
0 & 1 & 2 & 5\\
0 & 0 & 0 & 0\\

\end{pmatrix} );

より、

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

掃出せなかった をパラメータとすれば、

　&math( \bm x=z\begin{pmatrix}

• 2\\-2\\1 \end{pmatrix}+s\begin{pmatrix} 4\\5\\0 \end{pmatrix} );

第１項は のベクトルだから、 は &math(\bm b_2=\begin{pmatrix} 4\\5\\0 \end{pmatrix}); を基底に持つ１次元空間。

次に を求めるために「 をかけると に入るベクトル」を探す。

　

の拡大係数行列を変形して、

　&math( \begin{pmatrix}

2 & -2 & 0 & 4\\
-3 & 2 & -2 & 5\\
-11 & 9 & -4 & 0\\

\end{pmatrix}\sim \begin{pmatrix}

1 & -1 & 0 & 2\\
0 & -1 & -2 & 11\\
0 & -2 & -4 & -22\\

\end{pmatrix}\sim \begin{pmatrix}

1 & 0 & 2 & -9\\
0 & 1 & 2 & -11\\
0 & 0 & 0 & 0\\

\end{pmatrix} );

掃出せなかった をパラメータとすれば、

　&math( \bm x=z\begin{pmatrix}

• 2\\-2\\1 \end{pmatrix}+s\begin{pmatrix}
• 9\\-11\\0 \end{pmatrix} );

第１項は のベクトルだから、 は &math(\bm b_3=\begin{pmatrix}

• 9\\-11\\0 \end{pmatrix}); を基底に持つ１次元空間。

念のため、 があるか調べるために 「 をかけると に入るベクトル」を探してみる。

　&math( \begin{pmatrix}

2 & -2 & 0 & -9\\
-3 & 2 & -2 & -11\\
-11 & 9 & -4 & 0\\

\end{pmatrix}\sim \begin{pmatrix}

2 & -2 & 0 & -9\\
-6 & 4 & -4 & -22\\
-22 & 18 & -8 & 0\\

\end{pmatrix}\sim \begin{pmatrix}

2 & -2 & 0 & -9 \\
0 & -2 & -4 & -49\\
0 & -4 & -8 & -99\\

\end{pmatrix}\sim \begin{pmatrix}

2 & -2 & 0 & -9 \\
0 & -2 & -4 & -49\\
0 & 0 & 0 & -1\\

\end{pmatrix} );

最後の行は を表し、これを満たすベクトルは存在しないから、 である。

したがって、

　&math( W(\lambda)= \underbrace{ \underbrace{ \underbrace{V(\lambda)}_{[\bm b_1]}\,\dot+\, \underbrace{\Delta V^2(\lambda)}_{[\bm b_2]} }_{V^2(\lambda)} \,\dot+\, \underbrace{\Delta V^3(\lambda)}_{[\bm b_3]} }_{V^3(\lambda)} );

の構造を持っている。

例２ 固有空間が２次元で非対称 †

　&math(A=\begin{pmatrix}

• 9&8&-4\\
• 15&13&-6\\
• 5&4&-1 \end{pmatrix} );

とすると、

　

は、 を導き、 掃出せなかった をパラメータで と置けば、

　&math(\bm x=s\begin{pmatrix} 4\\5\\0 \end{pmatrix}+t\begin{pmatrix}

• 2\\0\\5 \end{pmatrix});

は２次元で、 &math(\bm b_1=\begin{pmatrix} 4\\5\\0 \end{pmatrix},\bm b_2=\begin{pmatrix}

• 2\\0\\5 \end{pmatrix}); を基底に取れる。

を求めるために、

&math((A-\lambda I)\bm x=s\begin{pmatrix} 4\\5\\0 \end{pmatrix}+t\begin{pmatrix}

• 2\\0\\5 \end{pmatrix});

と置けば、

　&math( \begin{pmatrix}

• 10&8&-4&4&-2\\
• 15&12&-6&5&0\\
• 5&4&-2&0&5 \end{pmatrix}\sim \begin{pmatrix}
• 5&4&-2&0&5\\ 0&0&0&4&-12\\ 0&0&0&5&-15\\ \end{pmatrix}\sim \begin{pmatrix}
• 5&4&-2&0&5\\ 0&0&0&1&-3\\ 0&0&0&0&0\\ \end{pmatrix}\sim );

となって、&math( \begin{cases}

• 5x+4y-2z=5t\\ 0=s-3t\\ \end{cases} ); を得る。

すなわち、

　&math( \bm x=s'\begin{pmatrix} 4\\5\\0 \end{pmatrix}+t'\begin{pmatrix}

• 2\\0\\5 \end{pmatrix}+t\begin{pmatrix}
• 1\\0\\0 \end{pmatrix} );

に対して

　&math( (A-\lambda)\bm x= 3t\begin{pmatrix} 4\\5\\0 \end{pmatrix}+t\begin{pmatrix}

• 2\\0\\5 \end{pmatrix} =t\begin{pmatrix} 10\\15\\5 \end{pmatrix} =5t\begin{pmatrix} 2\\3\\1 \end{pmatrix} );

つまり、

　&math(\bm b_2'=\begin{pmatrix} 2\\3\\1 \end{pmatrix});

　&math(\bm b_3=\begin{pmatrix}

• 1\\0\\0 \end{pmatrix});

とすると、

も を張るから、 広義固有空間の構造は、

　&math( W(1)= \underbrace{ \underbrace{ \underbrace{V_1(1)}_{\bm b_1} \,\dot+\, \underbrace{V_2(1)}_{\bm b_2'} }_{V(1)} \,\dot+\, \underbrace{V_2^2(1)}_{\bm b_2} }_{V^2(1)} );

例３ 固有空間が２次元で対称 †

　&math( A=\begin{pmatrix}

3 & 8 & 4 & -12 \\
6 & 5 & 4 & -12 \\
-6 & -1 & -5 & 6 \\
5 & 9 & 4 & -15 \\

\end{pmatrix} );

　&math( \begin{pmatrix} \lambda&1\\ &\lambda&\\ &&\lambda&1\\ &&&\lambda\\ \end{pmatrix} );

一般のジョルダン標準形 †

ジョルダン標準形の応用 †

ジョルダン標準形を見れば、その行列の広義固有空間がどのような構造になっているか一目瞭然となる。 その意味でジョルダン標準形は非常に重要であるが、応用的にも重要な意味を持つ。