このページでは,2次〜3次の正方行列に対して,対角化,ジョルダン標準形を利用して行列のn乗を求める方法を調べる. 【ジョルダン標準形】 線形代数の教科書では,著者によって,[A] 対角行列を含めてジョルダン標準形と呼ぶ場合と,[B] 用語として対角行列とジョルダン標準形を分けている場合があるので,文脈を見てどちらの立場で書かれているかを見分ける必要がある. [A] ジョルダン標準形
[B]
対角行列
ジョルダン標準形 [A]はすべてのジョルダン細胞が1次正方行列から成る場合が正方行列であると考える.(言葉の違いだけ) |
3次正方行列の場合を例にとって,以下のこのページの教材に書かれていることの要約を示すと次の通り.
【要約】
はじめに与えられた行列に対する固有方程式を解いて,固有値を求める. (1) 固有値に重複がない場合(固有値が虚数であっても) となる固有ベクトルを求めると,これらは互いに1次独立になるので,これらの列ベクトルを束にしてできる変換行列を とおくと,この変換行列は正則になる(逆行列が存在する). 固有値を対角成分にした対角行列を とおくと …(1.1) もしくは …(1.2) が成り立つ. このとき,を(正則な)変換行列,を対角行列といい,は対角化可能であるという.「行列を対角化せよ」という問題に対しては,(1.1)または(1.2)を答えるとよい.
この教材に示した具体例
対角行列は行列の積としての累乗が容易に計算できるので,これを利用して行列の累乗を計算することができる.【例1.1】 【例1.2.1】 【例1.2.2】 【例1.3.1】 【例1.3.2】 |
(2) 固有方程式が重解をもつ場合,
@) 元の行列自体が対角行列であるとき これらの行列は,変換するまでもなく対角行列になっているから,n乗などの計算は容易にできる. A) 上記の@)以外で固有方程式が重複解をもつとき,次のようにジョルダン標準形と呼ばれる形にできる A) 重複度1の解と二重解が固有値であるとき
a) 任意のベクトル(ただし,後で求まるベクトルとは1次独立でなければならない)を選び
となる列ベクトルが求まるときは で定まる変換行列を用いて と書くことができる.
この教材に示した具体例
b) ≪2次正方行列≫ 【例2.1】(1) 【例2.2.1】【例2.2.2】【例2.2.3】【例2.2.4】 ≪3次正方行列≫ 【例2.1】(2) 【例2.3.1】 【例2.3.2】 となる列ベクトルが求まるときは で定まる変換行列を用いて対角化できる.すなわち と書くことができる.
この教材に示した具体例
【例2.3.3】 【例2.3.4】 【例2.3.5】 |
B) 三重解が固有値であるとき
なお,2次正方行列で固有値が重解となる場合において,1次独立な2つのベクトルについて
a) 任意のベクトル(ただし,後で求まるベクトルとは1次独立でなければならない)を選び
となるベクトルが定まるときは で定まる変換行列を用いて と書くことができる.
この教材に示した具体例
【例2.4.3】 【例2.4.4】
b) 任意のベクトル(ただし,後で求まるベクトルとは1次独立でなければならない)を選び
となるベクトルが定まるときは で定まる変換行列を用いて と書くことができる.
この教材に示した具体例
【例2.4.1】 【例2.4.2】 が成り立てば,平面上の任意のベクトルはと書けるから, となる.したがって となり,このようなことが起こるのは自体が単位行列の定数倍となっている場合に限られる. 同様にして,3次正方行列で固有値が三重解となる場合において,1次独立な3つのベクトルについて が成り立てば,空間内の任意のベクトルはと書けるから, となる.したがって となり,このようなことが起こるのは自体が単位行列の定数倍となっている場合に限られる. これらが(2)@)に述べたものである. |
1.1 対角化可能な行列の場合
与えられた行列から行列の累乗を求める計算は一般には難しい.しかし,次のような対角行列では容易にn乗を求めることができる.…(*1.1) ここで,だから,中央の掛け算が簡単になり 同様にして,一般に次の式が成り立つ. 両辺に左からを右からを掛けると …(*1.2) このように,が対角行列となるように変形できる行列は,対角化可能な行列と呼ばれ上記の(*1.1)を(*1.2)の形に変形することによって,を求めることができる. |
【例1.1】
(1)(1) (2) に対して,,とおくと すなわち が成り立つから (2) に対して, ,とおくと が成り立つ.すなわち ※上記の正則な変換行列および対角行列は固有ベクトルを束にしたものと固有値を対角成分に並べたものであるが,その求め方は後で解説する. |
1.2 対角化できる場合の対角行列の求め方(実際の計算)
2次の正方行列が,固有値,固有ベクトルをもつとは一次変換の結果がベクトルの定数倍になること,すなわち …(1) となることをいう. 同様にして,固有値,固有ベクトルをもつとは 一次変換の結果がベクトルの定数倍になること,すなわち …(2) となることをいう. (1)(2)をまとめると次のように書ける. |
固有値が相異なり重複解を持たないとき,すなわちのとき,固有ベクトルとは互いに1次独立に選ぶことができ,固有ベクトルを束にして作った変換行列は正則行列(逆行列が存在する行列)になる. そこで,を対角行列として すなわち の形で対角化できることになり,対角行列は累乗を容易に計算できるので によりが求められる. |
【例1.2.1】
(解答)(1)を対角化してください. 固有方程式を解く 固有ベクトルを求める ア) のとき すなわち より 1つの固有ベクトルとして,が得られる. イ) のとき すなわち より 1つの固有ベクトルとして,が得られる. ア)イ)より まとめて書くと …(答) |
【例1.2.2】
(解答)(2)を対角化してください. 固有方程式を解く 固有ベクトルを求める ア) のとき すなわち より1つの固有ベクトルとして,が得られる. 同様にして イ)のとき1つの固有ベクトルとして,が得られる. ウ)のとき1つの固有ベクトルとして,が得られる. 以上の結果をまとめると …(答) |
1.3 固有値が虚数の場合
正方行列に異なる固有値のみがあって,固有値に重複がない場合には,対角化できる.元の行列が実係数の行列であるとき,実数の固有値であっても虚数の固有値であっても重複がなければ対角化できる. 元の行列が実係数の行列であって,虚数の固有値が登場する場合でも行列のn乗の成分は実数になる---虚数の固有値と言っても共役複素数の対から成り,それらの和や積で表される行列のn乗は,実数で書ける.
【例題1.3.1】
(解答)次の行列が対角化可能かどうかを調べ,を求めてください. 固有方程式を解く 固有ベクトルを求める ア) のとき より イ) のとき より |
ア)イ)より ゆえに,行列は対角化可能…(答) は正の整数として,次の早見表を作っておくと後が楽
この表を使ってまとめると
1)n=4kのとき 2)n=4k+1のとき 3)n=4k+2のとき 4)n=4k+3のとき
原点の回りに角θだけ回転する1次変換
に当てはめると,となるから で左の計算と一致する |
【例題1.3.2】
(解答)次の行列が対角化可能かどうかを調べ,を求めてください. 固有方程式を解く 固有ベクトルを求める ア) のとき より イ) のとき より ア)イ)より ゆえに,行列は対角化可能…(答) |
ここで複素数の極表示を考えると ここで, だから 結局 以下 (nは正の整数,kは上記の1〜8乗) このように,元の行列の成分が実数であれば,その固有値や固有ベクトルが虚数であっても,(予想通りに)n乗は実数になることが示せる. (別解) 原点の回りに角θだけ回転して,次に原点からの距離をr倍することを表す1次変換の行列は であり,与えられた行列は と書けるから ※回転を表す行列になるものばかりではないから,前述のように虚数の固有値,固有ベクトルで実演してみる意義はある. |
2.1 対角化はできないがそれに近い形にできる場合
行列の固有値が重解になる場合などにおいて,対角化できない場合でも,次のように対角成分の1つ上の成分を1にした形を利用すると累乗の計算ができる.
【例2.1】
(1)(1) (2) に対して,,とおくと |
すなわち が成り立つから (2) に対して, ,とおくと が成り立つ.すなわち |
2.2 ジョルダン標準形の求め方(実際の計算)
【例題2.2.1】
(解答)(1) 次の行列のジョルダン標準形を求めてください. 固有方程式を解いて固有値を求める (重解) 固有ベクトルを求める のとき すなわち より1つの固有ベクトルとして,が得られる.
[以下の解き方@] となると1次独立なベクトルを求める.
いきなり,そんな話がなぜ言えるのか疑問に思うかもしれない.
実は,この段階ではとなる行列があるとは証明できていないが「求まったらいいのにな!」と考えて,その条件を調べている--方程式として解いているだけ.「もしこのような行列があれば右辺がジョルダン標準形になるから」対角化できなくてもn乗が計算できるから嬉しいのである.(実際には,必ず求まる!) |
両辺の成分を比較すると だから,…(*A)が必要十分条件 1つの固有ベクトルとして,が得られる. これにより すなわち …(答) (参考) この後,次のように変形すれば問題の行列Aのn乗が計算できる. |
[以下の解き方A] と1次独立な(が1次独立ならば行列は正則になり,逆行列が求まるが,そうでなければ逆行列は求まらない)ベクトル条件(*A)を満たせばよいから,必ずしもでなくてもよい.ここでは,他のベクトルでも同じ結果が得られることを示してみる. 1つの固有ベクトルとして,を使うと …(答) |
(参考) この後,次のように変形すれば問題の行列Aのn乗が計算できる. この結果は@の結果と一致する |
[以下の解き方B] 線形代数の教科書,参考書には,次のように書かれていることがある. 行列の固有値が(重解)で,これに対応する固有ベクトルがのとき, …(1) と1次独立なベクトルは,次の計算によって求められる. …(2) これらの式の意味は次のようになっている
(1)は固有値がで,これに対応する固有ベクトルがであることから
を移項すれば として(1)得られる. これに対して,(2)は次のように分けて考えると を表していることが分かる. を列ベクトルに分けると が(1)を表しており が(2)を表している. (2)はであるから と書ける.要するに(1)を満たす固有ベクトルを求めてそれをとして,次に を満たすを求めるという流れになる. 以上のことは行列とベクトルで書かれているので,必ずしも分かり易いとは言えないが,解き方@において となると1次独立なベクトルを求める. ・・・そのようながあったらいいのにな〜[対角成分の1つ上の成分が1になっている行列でもn乗ができるから]〜という「願いのレベル」で未知数を求めていることと同じになる. |
【解き方Bのまとめ】
[以下の解き方B]…(1) …(2) となるベクトルを2つの列ベクトルとして,それらを束にして行列にしたもの は,元の行列をジョルダン標準形に変換する正則な変換行列になる.すなわち が成り立つ. 実際に解いてみると・・・ 行列の固有値を求めると(重解) そこで,次の方程式を解いて,を求める. …(1) …(2) (1)より したがって,を満たすベクトル(ただし,零ベクトルでないもの)は固有ベクトル. そこで,とする. 次に(2)により したがって,を満たすベクトル(ただし,零ベクトルでないもの)は解のベクトル. そこで,とする. |
[解き方Bの2]・・・別の解説 線形代数の教科書,参考書によっては,次のように解説される場合がある. はじめに,零ベクトルでない(かつ固有ベクトルと平行でない)「任意のベクトル」を選ぶ.次に(2)式によってを求めたら,「は必ず(1)を満たす」ので,これらの組を解とするのである. …(1’) …(2’) 前の解説と(1’)(2’)の式は同じであるが,「は任意のベクトルでよい」「(2’)で求めた「は必ず(1’)を満たす」という所が,前の解説と違うように聞こえるが・・・実際に任意のベクトルを代入してみると,次のようになる. |
とおくとはAの固有ベクトルになっており,(1)を満たす. この場合,任意のベクトルは固有ベクトルの倍率を決めることだけに使われている. 例えば,任意のベクトルをとすると, となって が得られる. 初め慣れるまでは,考え方が難しいが,慣れたら単純作業で求められるようになる. |
【例題2.2.2】
(解答)次の行列のジョルダン標準形を求めて,を計算してください. (重解) のとき,固有ベクトルは より よって,1つの固有ベクトルは
(解き方@)
,正則行列,,ジョルダン標準形に対してこのベクトルと1次独立なベクトルを適当に選び とおくと となれば,対角化はできなくても,それに準ずる上三角化ができる. ゆえに,・・・(**) 例えば1つの解として とすると, すなわち となるから 結局 …(答) |
前述において,(解き方@)で示した答案は,(**)を満たす他のベクトルを使っても,同じ結果が得られる. (解き方A) 例えば1つの解として とすると, ,正則行列,,ジョルダン標準形に対して すなわち となるから 結局 …(答) となって,結果は等しくなる. (解き方B) 行列の固有値を求めると(重解) そこで,次の方程式を解いて,を求める. …(1) …(2) (1)より したがって,を満たすベクトル(ただし,零ベクトルでないもの)は固有ベクトル. そこで,とする. 次に(2)により そこで,とする. 以下は(解き方@)(解き方A)と同様になる. (解き方Bの2) 例えばとおくと, となり これを気長に計算すると,上記(解き方@)(解き方A)の結果と一致する. |
【例題2.2.3】
(解き方@1)次の行列のジョルダン標準形を求めて,を計算してください. 固有方程式を解く (重解) のとき,固有ベクトルは より よって,1つの固有ベクトルは そこで となるを求める ・・・(**) より ,正則行列,,ジョルダン標準形に対して すなわち となるから …(答) |
(解き方A) (**)においてを選んだ場合 ,正則行列,,ジョルダン標準形に対して すなわち となるから …(答) (解き方B) 行列の固有値を求めると(重解) そこで,次の方程式を解いて,を求める. …(1) …(2) (1)より したがって,を満たすベクトル(ただし,零ベクトルでないもの)は固有ベクトル. そこで,とする. 次に(2)により そこで,とする. 以下は(解き方@)と同様になる. (解き方Bの2) 固有ベクトルと1次独立な任意の(零ベクトルでない)ベクトルとしてを選び,によって定まるベクトルにより正則行列 を定めると ,正則行列,,ジョルダン標準形に対して すなわち となるから …(答) |
【例題2.2.4】
(解き方@)次の行列のジョルダン標準形を求めて,を計算してください. 固有方程式を解く (重解) のとき,固有ベクトルは より よって,1つの固有ベクトルは そこで となるを求める ・・・(**) より ,正則行列,,ジョルダン標準形に対して すなわち となるから …(答) |
(解き方A) (**)においてを選んだ場合 ,正則行列,,ジョルダン標準形に対して すなわち となるから …(答) (解き方Bの2) 固有ベクトルと1次独立な任意の(零ベクトルでない)ベクトルとしてを選び,によって定まるベクトルにより正則行列 を定めると ,正則行列,,ジョルダン標準形に対して すなわち となるから …(答) |
2.3 3次正方行列で固有値が二重解になる場合
3次正方行列をジョルダン標準形にすると,行列のn乗が次のように計算できる
【例題2.3.1】
(解き方@)次の行列のジョルダン標準形を求めてください. 固有方程式を解く (重複度1),(重複度2) 固有ベクトルを求める ア)(重複度1)のとき より イ)(重複度2)のとき より |
これら2つのベクトルと1次独立なベクトルをもう1つ求める必要があるから となるベクトルを求めるとよい. 以上により ,正則行列,ジョルダン標準形に対して となる (解き方Bの2) (重複度1),(重複度2)に対して, と1次独立になるように気を付けながら,任意のベクトルを用いて次の式から定まるを用いて,正則な変換行列 を定める. …(2) …(1) たとえば,,とおくと, に対しては,が定まるから,解き方@と同じ結果を得る. |
【例題2.3.2】
(解き方@)次の行列のジョルダン標準形を求めてください. 固有方程式を解く (重複度1),(重複度2) 固有ベクトルを求める ア)(重複度1)のとき より イ)(重複度2)のとき より |
これら2つのベクトルと1次独立なベクトルをもう1つ求める必要があるから となるベクトルを求めるとよい. 以上により ,正則行列,ジョルダン標準形に対して となる (解き方Bの2) (重複度1),(重複度2)に対して, と1次独立になるように気を付けながら,任意のベクトルを用いて次の式から定まるを用いて,正則な変換行列 を定める. …(2) …(1) たとえば,,とおくと, に対しては,が定まるから,解き方@と同じ結果を得る. |
2次正方行列が二重解をもつとき,元の行列自体が単位行列の定数倍である場合を除けば,対角化できることはなくジョルダン標準形 になる. これに対して,3次正方行列が1つの解と二重解をもつ場合,二重解に対応する側の固有ベクトルが1つしか定まらない場合は上記の【2.3.1】, 【2.3.2】のようにジョルダン標準形になるが,二重解に対応する側の固有ベクトルが独立に2個求まる場合には,この行列は対角化可能である.すなわち, となる列ベクトルが求まるときは で定まる変換行列を用いて対角化できる.すなわち と書くことができる.
【例題2.3.3】
固有方程式を解く次の行列が対角化可能かどうか調べてください. |
固有ベクトルを求める ア) のとき より イ) のとき より これを満たすベクトルは独立に2個できる 以上により 変換行列,対角行列により すなわち |
【例題2.3.4】
(略解)次の行列が対角化可能かどうか調べてください. 固有値に対する固有ベクトルは 固有値(二重解)に対する固有ベクトルは 対角化可能 |
【例題2.3.5】
(略解)次の行列が対角化可能かどうか調べてください. 固有値に対する固有ベクトルは 固有値(二重解)に対する固有ベクトルは 対角化可能 |
2.4 3次正方行列で固有値が三重解になる場合
3次正方行列をジョルダン標準形にすると,行列のn乗が次のように計算できる
三重解の場合,次の形が使えることがある. 次の形ではかなり複雑になる
【例題2.4.1】
(解き方@)次の行列のジョルダン標準形を求めてて,n乗を計算してください. 固有方程式を解く (重複度3) 固有ベクトルを求める (は任意) これを満たすベクトルは1次独立に2つ作れる |
正則な変換行列を作るには,もう1つ1次独立なベクトルが必要だから次の形でジョルダン標準形を求める (は任意) より 以上により n乗を計算するには,次の公式を利用する |
(解き方Bの3) 1次独立なベクトルの束から作った行列 が次の形でジョルダン標準形 となるようにベクトルを求める. 両辺を列ベクトルに分けると …(1) …(2) …(3) すなわち …(1’) …(2’) …(3’) そこで,任意の(ただし,後で求まるベクトルとは1次独立でなければならない)ベクトルを選び,(3)で定まるを求めると固有ベクトルになって(2)を満たしているので,これと独立にもう1つ固有ベクトルを定めるとよい. |
例えば,とおくと,となる. (1’)は次の形に書ける (は任意) と1次独立となるようにを選ぶと, このとき, について,だからは正則になる. 変換行列は解き方@と同じではないが,n乗の計算を同様に行うと,結果は同じになる |
【例題2.4.2】
(略解:解き方B)次の行列のジョルダン標準形を求めください. 固有方程式は三重解をもつ これに対応する固有ベクトルを求める より (は任意) これを満たすベクトルは独立に2つ選べる |
これらと独立にもう1つベクトルを定めるために となるベクトルを求める. (は任意) 以上により 正則な変換行列 ジョルダン標準形 として が成り立つ. |
【例題2.4.3】
(解き方Bの2)次の行列のジョルダン標準形を求めて,n乗を計算してくださいください. 固有方程式を解く (三重解) 次の形でジョルダン標準形を求める 正則な変換行列は3つの1次独立なベクトルを束にしたものとする 両辺を列ベクトルに分けると …(1) …(2) …(3) すなわち …(3’) …(2’) …(1’) 次の順に決める:任意の(ただし,後で求まるベクトルとは1次独立でなければならない)ベクトルを選び,(3’)で定まるを求める.さらに(2’)でを定める:(1’)は成り立つ. |
例えば とおくと となる. 以上がジョルダン標準形である n乗は次の公式を使って求める |
【例題2.4.4】
(解き方Bの2)次の行列のジョルダン標準形を求めて,n乗を計算してくださいください. 固有方程式を解く (三重解) 変換行列を求める. 任意のベクトル(ただし,後で求まるベクトルとは1次独立でなければならない)を選び となるを求めて,この作業を繰り返す. 例えば,次のように定まる. だから とおくと …(#1) により さらに …(#2) により なお …(#3) が成り立つ. |
(#1)は
以上によりすなわち …(#1’) を表している. (#2)は すなわち …(#2’) を表している. (#3)は すなわち …(#3’) を表している. (#1’)(#2’)(#3’)より変換行列をによって作ると (右辺のジョルダン標準形において,1列目のは単独,2列目,3列目のの上には1が付く) に対して,変換行列 ジョルダン標準形 とおくと が成り立つ.すなわち だから |