コンテンツにスキップ

対角化とは

対角行列をベクトルに掛けると、各成分を別々の倍率で伸び縮みさせるだけです。固有ベクトルを座標軸に選べれば、元の複雑な変換もこの単純な形で見られます。

AA が一次独立な固有ベクトル v1,,vn\boldsymbol{v}_1,\ldots,\boldsymbol{v}_n をもち、対応する固有値が λ1,,λn\lambda_1,\ldots,\lambda_n なら

P=(v1  vn),D=diag(λ1,,λn).P=(\boldsymbol{v}_1\ \cdots\ \boldsymbol{v}_n),\qquad D=\operatorname{diag}(\lambda_1,\ldots,\lambda_n).

PP の列と DD の対角成分は、必ず同じ順序で対応させます。

固有値が重複しても、独立な固有ベクトルが必要本数あれば対角化できます。たとえば I2I_2 は固有値1だけを2重にもっていますが、すべての非零ベクトルが固有ベクトルなので対角化可能です。

A=PDP1A=PDP^{-1} なら途中の P1PP^{-1}P が消えるため

Ak=PDkP1.A^k=PD^kP^{-1}.

上の例では

Dk=(3k001),D^k=\begin{pmatrix}3^k&0\\0&1\end{pmatrix},

したがって

Ak=12(3k+13k13k13k+1).A^k=\frac12\begin{pmatrix}3^k+1&3^k-1\\3^k-1&3^k+1\end{pmatrix}.

これは漸化式、連立微分方程式、長時間の変化を調べるときに役立ちます。

対角化できない行列も、複素数上ではジョルダン標準形というほぼ対角な形に整理できます。ここでは名前の紹介に留めます。

Q対角行列を組み立てる★☆☆

A=(4013)A=\begin{pmatrix}4&0\\1&3\end{pmatrix} を対角化してください。

HINTヒントを見る
固有値は 4,34,3。それぞれの固有ベクトルを同じ順序で PP に並べます。
ANS解答を見る

λ=4\lambda=4 では固有ベクトル (1,1)T(1,1)^{\mathsf T}λ=3\lambda=3 では (0,1)T(0,1)^{\mathsf T} を選べます。したがって

P=(1011),P1=(1011),D=(4003).P=\begin{pmatrix}1&0\\1&1\end{pmatrix},\quad P^{-1}=\begin{pmatrix}1&0\\-1&1\end{pmatrix},\quad D=\begin{pmatrix}4&0\\0&3\end{pmatrix}.

直接掛けると P1AP=DP^{-1}AP=D です。

Qべき乗を求める★★☆

上の AA について AkA^k を求めてください。

HINTヒントを見る
Ak=PDkP1A^k=PD^kP^{-1} を使います。
ANS解答を見る
Dk=(4k003k)D^k=\begin{pmatrix}4^k&0\\0&3^k\end{pmatrix}

なので

Ak=(1011)(4k003k)(1011)=(4k04k3k3k).A^k =\begin{pmatrix}1&0\\1&1\end{pmatrix} \begin{pmatrix}4^k&0\\0&3^k\end{pmatrix} \begin{pmatrix}1&0\\-1&1\end{pmatrix} =\begin{pmatrix}4^k&0\\4^k-3^k&3^k\end{pmatrix}.

k=1k=1 を代入すると元の AA に戻り、検算できます。

Q対角化できないことを示す★★★
A=(1101)A=\begin{pmatrix}1&1\\0&1\end{pmatrix}

が対角化できないことを示してください。

HINTヒントを見る
唯一の固有値 11 に対する固有空間の次元を調べます。
ANS解答を見る

特性多項式は (1λ)2(1-\lambda)^2 です。λ=1\lambda=1 に対し

AI=(0100)A-I=\begin{pmatrix}0&1\\0&0\end{pmatrix}

なので固有ベクトルは (x,0)(x,0) の形で、固有空間は1次元です。R2\mathbb R^2 の基底となる2本の固有ベクトルを取れないため、AA は対角化できません。