うさぎでもわかる線形代数 第22羽 ジョルダン標準形

スポンサードリンク

こんにちは、ももやまです。

今回は、線形代数の中でもかなりの難易度を誇り、期末試験や院試などで出題されるジョルダン標準形がどんなものなのかを簡単に説明し、3次ジョルダン標準形までの求め方を例題や練習問題を用意し、(たぶん)わかりやすくまとめています。

※注意

ジョルダン標準形の求め方が知りたい方は、項目3以降をご覧ください。

前回の記事「うさぎでもわかる線形代数 第21羽」はこちら↓

www.momoyama-usagi.com

スポンサードリンク

1.ジョルダン標準形とは

行列の中には対角化ができないものもありました。例えば、(2114)は、2次正方行列なのに固有ベクトルが1つしか求められないため、対角化ができませんね*1

しかし、対角化できない行列でもできる限り対角行列に近い形にすることはできます。この形のことをジョルダン標準形と呼ばれます。

ジョルダン標準形は、ジョルダンブロック(下で説明します)を対角線上に並べた行列を表します。

ジョルダンブロック(ジョルダン細胞)

ジョルダンブロックは下のように対角成分に同じ数(固有値)を並べ、対角成分の1つ上が1となる行列のことを表し、

f:id:momoyama1192:20190831142057g:plain

ジョルダンブロックは、対角成分に並んでいる数(固有値)を α、ジョルダンブロックのサイズ(行列の大きさ)を n とし、 J(α,n) と表すことができます。

ジョルダン標準形とジョルダンブロック

ジョルダン標準形は下のようにジョルダンブロックを対角線上に並べ、残りの成分が0であるような行列を表します。

f:id:momoyama1192:20190831142116g:plain

ジョルダン標準形を記号で表現する際には、どんな

例えば、図の一番左の行列は、J(3,2)J(4,1) のジョルダンブロックで構成されていますね。なので、J(3,2)J(4,1)と表現します。

行列が対角化できる場合、ジョルダン標準形はすべてサイズ1のジョルダンブロック J(x,1) で構成されます。

対角化はジョルダン標準形の特殊なパターンなのです。

(もし対角化がまだわからない人は早急にこちらの記事で復習しましょう!)

スポンサードリンク

2.2次・3次のジョルダン標準形とジョルダンブロック

2次平方行列、3次平方行列の場合のジョルダン標準形とジョルダンブロックの組み合わせを下に記しておきます。

2次正方行列の場合

f:id:momoyama1192:20190831142101g:plain

3次正方行列の場合

ジョルダン標準形が対角行列

f:id:momoyama1192:20190831142105g:plain

対角化できない場合

f:id:momoyama1192:20190831142110g:plain

スポンサードリンク

3.ジョルダン標準形の求め方(2次)

ではまずは2次正方行列のジョルダン標準形を求める方法を例題を使って説明していきましょう。

例題1

行列(2114)を正則行列 P を用いて P1AP をジョルダン標準形にしなさい。

解説1

(1)

Step1:固有値を求める

固有値を t とすると、固有方程式は、|AtE|=|2t114t|=(t4)(t2)+1=t26t+9=(t3)2=0より固有値は3(2重解)

Step2:固有ベクトルを求める

固有値が3のときの固有ベクトルは、(A3E)=(1111) (1100)となる。xy=0を解くと、任意定数 k を用いて(xy)=k(11)と表せるので、固有ベクトルは1本あり、固有ベクトル p1 は、p1=(11)となる。

ここまでは通常の対角化と全く同じ流れです。

Step3:広義固有ベクトルを求める

このままでは固有ベクトルが1本足りないため、ジョルダン標準形がつくれませんね。

そこで、(A3E)p2=p1を満たすような p2 を探します。1つ前に求めた固有ベクトルを借りてくるのです。これを広義固有ベクトルと呼びます。

通常の固有ベクトルも広義固有ベクトルの中に含まれるので「広義固有ベクトルを求めなさい」という問題が出た際には通常の固有ベクトルを答えに入れ忘れないように注意しましょう。

固有値3に対する広義固有ベクトル p2 は、(111111) (111000)となる。xy=1を解くと、任意定数 k を用いて(xy)=(01)+k(11)と表せますね。
(実はここは任意定数とおかなくても xy=1 を満たすような解を1つ見つければOK、例題2からは任意定数とおかずに解を1つ見つける書き方に変えます。)

よって広義固有ベクトル p2k=0 としてp2=(01)となります。

Step4:ジョルダン標準形を求める

よって、{(A3E)p1=0Ap1=3p1   {(A3E)p2=p1Ap2=p1+3p2が成立するので、P=(p1,p2)=(1011)とすると、AP=A(p1,p2)=(3p1,p1+3p2)=(p1,p2)(3103)=P(3103)となります。

ここで、J=(3103)とすると、AP=PJP1AP=J=(3103)となり、うまくジョルダン標準形が求められますできますね。
(ジョルダン標準形の検算も対角化と同じように AP=PJ で検算ができます。)

ジョルダンブロック:J(3,2)

広義固有ベクトルと広義固有空間

広義固有ベクトルも固有ベクトルと同じようにベクトル空間を用いて表すことができます。

例題1の広義固有空間は、V~(3)={(11),(01)}となります。

下のほうに一応定義は示しておきますが、参考程度にどうぞ。
固有値ごとに広義固有ベクトルをまとめたものが広義固有空間だと思ってもらってOKです。

広義固有空間を用いた表し方n 次正方行列 A の固有値 t に対し、固有空間 V(t)V(t)={xRn | 1 k (AtE)kx=0}と表せる。(固有空間は k=1 のときの特殊なパターン)

さきほどの例題で広義固有ベクトルを求めるときに(A3E)p2=p1を解きましたね。実は、(A3E)p1=(A3E)2p2=0と変形できるので、k=2 のときの広義固有ベクトル(広義固有空間)を求めていることと同じなのです!

4.ジョルダン標準形の求め方(3次)

3次のジョルダン標準形も2次と同じように計算をすることができます。

ただし、重解のパターンによって解き方が3パターンにわかれるので注意してください。

パターン1 2重解の固有値があり、固有ベクトルが1本足りない場合

例題2

行列(101020011)を正則行列 P を用いて P1AP をジョルダン標準形にしなさい。

解説2

固有方程式は、|AtE|=|1t0102t0011t|=(1t)|2t011t|=(1t)2(2t)=0より固有値は1(2重解)、2となる。

つぎに、広義固有ベクトル(および固有ベクトル)を求めていく。

(1) 固有値が1のとき

固有値1に属する固有ベクトルは、(AE)=(001010010) (010001000)となる。{y=0z=0を解くと(要注意)、(xyz)=k(100)となるので、固有ベクトルは1本あり、固有ベクトル p1 は、p1=(100)となる。

2重解に対し、固有ベクトルが1本しかないので、(AE)p2=p1を満たすような p2 を求める。(001101000100) (00110100)となる。y=0,z=1の解の1つに(xyz)=(001)があるので*2、広義固有ベクトル p2p2=(001)となります。

(ii)

固有値が2のときの固有ベクトル

固有値2に属する固有ベクトルは、(A2E)=(101000011) (101011000)となる。{x+zyz=0を解くと(要注意)、(xyz)=k(111)となるので、固有ベクトルは1本あり、固有ベクトル p3 は、p3=(111)となる。

よって、{(AE)p1=0Ap1=p1   {(AE)p2=p1Ap2=p1+p2   {(A2E)p3=0Ap3=2p3が成り立つ。ここでP=(p1,p2,p3)=(101001011)とすると、AP=A(p1,p2,p3)=(p1,p1+p2,p3)=(p1,p2,p3)(110010002)=P(110010002)となります。

よって、P=(101001011)   J=P1AP=(110010002)とジョルダン標準形にできます。

ジョルダンブロックの構成:J(1,2)J(2,1)

[検算] AP=PJ=(112002012)

※広義固有空間はV~(1)={(100),(001),(111)}となる。

パターン2 3重解の固有値があり、固有ベクトルが2本足りない場合

つぎに固有値が3重解となり、固有ベクトルが1本しか出てこないパターンをみていきましょう。

例題3

行列(711812611)を正則行列 P を用いて P1AP をジョルダン標準形にしなさい。

解説3

変形がめんどくさいの

固有多項式は、|AtE|=|7t1181t2611t|=|1t02t81t2611t|=(1t)|1t211t|+(2t)|81t61|=(1t)(t22t1)+(2t)(8+(66t))=(t39t2+27t27)=(t3)(t26t+9)=(t3)3と変形できる。

よって固有値は3(3重解)となる。

(1) 固有値が3のとき

固有値3に属する固有ベクトルは、(A3E)= (411822612) (201402612) (201011000) となる。{2x+z=0y+z=0を解くと(要注意)、(xyz)=k(122)となるので、固有ベクトルは1本あり、固有ベクトル p1 は、p1=(122)となる。

広義固有ベクトルが2本足りないので、まずは(A3E)p2=p1を満たすようなベクトルを求める。

(411182226122) (201140226122) (201100000111)となる。

{2x+z=1y+z=1の解の1つに(xyz)=(001)があるので、広義固有ベクトル p2p2=(001)となります。

しかし、まだ1つ広義固有ベクトルが足りませんね。

なので、今度は、(A3E)p3=p2を満たすベクトルを求めます*3

(411082206121) (201140226121) (201100000112)となる。

2x+z=1,y+z=2の解の1つに(xyz)=(011)があるので、広義固有ベクトル p3p3=(001)となります。

※解の1つを選ぶ際に p1, p2, p3 が1次独立かを確認すること!

よって、{(A3E)p1=0Ap1=3p1   {(A3E)p2=p1Ap2=p1+3p2   {(A3E)p3=p2Ap3=p2+3p3が成り立つ。ここでP=(p1,p2,p3)=(100201211)とすると、AP=A(p1,p2,p3)=(p1,p1+p2,p2+p3)=(p1,p2,p3)(310031003)=P(310031003)となります。

よって、P=(100201211)   J=P1AP=(310031003)とジョルダン標準形にできます。

ジョルダンブロックの構成:J(3,3)

[検算] AP=PJ=(310623612)

※広義固有空間はV~(3)={(122),(001)}となる。

パターン3 3重解の固有値があり、固有ベクトルが1本足りない場合

例題4

行列(311131226)を正則行列 P を用いて P1AP をジョルダン標準形にしなさい。

解説4

Step1:固有値を求める

固有値を t とする。

固有多項式は、|AtE|=|3t1113t1226t|=|4t4t013t1226t|=(4t)|11013t1226t|=(4t)|11002t1046t|=(4t)|2t146t|と変形できる。

ここで、|2t146t|=(t2)(t6)+4=t28t+16=(t4)2となるので、固有方程式は|AtE|=(4t)(t4)2=(t4)3を満たす t が固有値となり、固有値は4(3重解)となる。

Step2:広義固有ベクトル(広義固有空間)を求める

固有値4に属する固有ベクトルは、(A4E)= (111111222) (111000000)となる。xyz=0を解くと(要注意)、任意定数 s,t を用いて(xyz)=s(110)+t(101)となります。なので、固有ベクトルを2本求めることができます。
(ただしこの時点では解が出せないので注意)

広義固有ベクトルが1本足りないので先ほどの例題のように固有ベクトルを持ってくるのですが、p1, p2 と候補の固有ベクトルが2つあるため、どのように固有

なので、定義にしたがって、(A4E)2p3=0を満たすベクトルを探します。(A4E)2=(000000000)=Oとなるので、どのような広義固有ベクトル p3 を持ってきてもOKなことがわかりますね。

例えば、p3=(100)とします。つぎに、(A4E)p3=p2を満たすような固有ベクトル p2 を求めます。

すると、(A4E)p3=(211111222)(100)=(112)=p2となります。これは、s(110)+t(101)s=1, t=2 に相当します。

よって、(A4E)p3=p2を成立させるためには広義固有ベクトル p2, p3p2=(112),   p3=(100)とする必要があります。

残りの固有ベクトル p1 は先ほどの解の (s,t)=(1,0) もしくは (s,t)=(0,1) どちらを選んでもOKです*4

なので、p1=(110)とします。

よって、{(A4E)p1=0Ap1=4p1   {(A4E)p2=0Ap2=4p2   {(A4E)p3=p2Ap3=p2+4p3が成り立つ。ここでP=(p1,p2,p3)=(111110020)とすると、AP=A(p1,p2,p3)=(4p1,4p2,p2+4p3)=(p1,p2,p3)(400041004)=P(400041004)となります。

よって、P=(111110020)   J=P1AP=(400041004)とジョルダン標準形にできます。

ジョルダンブロックの構成:J(4,2)J(4,1)

※広義固有空間はV~(4)={(110),(112),(100)}となる。

このように同じ固有値から複数の固有値を求められ、かつ固有ベクトルが重解の数存在しない場合の広義固有ベクトルの求め方は少しめんどくさいです。

5.ジョルダン標準形の求め方のまとめ

ジョルダン標準形がどんなものなのか、そしてどうやって求めるかについてを簡単にですがまとめました。

確認用にご覧ください。

ジョルダン標準形行列 A の固有値、広義固有ベクトルの組が (t1,p1), (t2,p2), … (tn,pn) のとき、対角化ができなくても、正則行列 P を用いて、対角成分の上に0もしくは1が並んだジョルダン標準形 J=P1AP を下のように作成できる。2次正方行列の場合

P=(p1,p2)   J=(t110t2)

3次正方行列の場合

P=(p1,p2,p3)   J=(t1100t2000t3)

n 次正方行列の場合

P=(p1,p2,p3,,pn)   J=(t11000t21000t30000tn)

ただし、対角成分の上の成分は、(AtE)pi=0として固有ベクトル p を求めた場合は0が、(AtE)pi=pi1のように1つ前の固有ベクトルを利用して広義固有ベクトル pi を求めた場合は1が入る。

3次までのジョルダン標準形を求めるためのステップを簡単にまとめたいと思います。

  1. 行列の固有値を求める
  2. 行列の固有ベクトルを求め、重解分だけ固有ベクトルが出るかを確認
  3. 重解分だけ固有ベクトルが求められなかった固有値に対し、(AtE)pi=pi1のように1つ前の固有ベクトル(or 広義固有ベクトル)を用いて広義固有ベクトルを出す
  4. 求めた広義固有ベクトルを並べて正則行列 P を作成し、ジョルダン標準形 J=P1AP を作成する。

注意点としては、固有ベクトルが複数本求められたにも関わらず、広義固有ベクトルを用いないとジョルダン標準形が求められない場合は、(AtE)pi=pi1がうまく成立するように pi1 を決める必要があるため、(AtE)2pi=0を解き、pi を決めてから pi1 を決める必要があります。

6.練習問題

では、ジョルダン標準形を求める練習をしましょう。

練習1→練習2→練習3の順に難易度があがります。

練習1

行列(3411)についてつぎの問いに答えなさい。

(1) 行列の固有値を求めなさい。
(2) 固有ベクトル(もしくは固有空間)を求めなさい。
(3) 広義固有ベクトル(もしくは広義固有空間)を求めなさい。
(4) 正則行列 P を用いて P1AP をジョルダン標準形にしなさい。

練習2

行列A=(1583146613)についてつぎの問いに答えなさい。

(1) 行列 A のの固有値を求めなさい。
(2) 固有ベクトル(もしくは固有空間)を求めなさい。
(3) 広義固有ベクトル(もしくは広義固有空間)を求めなさい。
(4) 正則行列 P を用いて P1AP=J をジョルダン標準形にしなさい。

練習3

行列A=(331111264)についてつぎの問いに答えなさい。

(1) 行列 A のの固有値を求めなさい。
(2) 固有ベクトル(もしくは固有空間)を求めなさい。
(3) 広義固有ベクトル(もしくは広義固有空間)を求めなさい。
(4) 正則行列 P を用いて P1AP をジョルダン標準形にしなさい。

7.練習問題の答え

解答1

(1) 固有値を求める

固有値を t とすると、固有方程式は、|AtE|=|3t411t|=(t3)(t+1)+4=t22t+1=(t1)2=0より固有値は1(2重解)となる。

(2) 固有ベクトルを求める

固有値が1のときの固有ベクトルは、(AE)=(2412) (1200)となる。x2y=0を解くと、任意定数 k を用いて(xy)=k(21)と表せるので、固有ベクトルは1本あり、固有ベクトル p1 は、p1=(21)となる。

固有空間はV(1)={(21)}となる。

(3) 広義固有ベクトルを求める

固有ベクトル p1 はもちろん広義固有ベクトルである。(AE)p2=p1を満たすような p2 を求める。

(242121) (121000)となる。x2y=1の解の1つに(xy)=(10)があるので、広義固有ベクトル p2p2=(10)となります。
(広義固有ベクトルは p1, p2 の2つ存在)

広義固有空間はV~(1)={(21),(10)}となる。

(4) ジョルダン標準形を求める

よって、{(AE)p1=0Ap1=3p1   {(AE)p2=p1Ap2=p1+p2が成立するので、P=(p1,p2)=(2110)とすると、AP=A(p1,p2)=(p1,p1+p2)=(p1,p2)(1101)=P(1101)となります。

よって、P=(2110),   J=P1AP=(1101)とジョルダン標準形にできる。

[検算] AP=PJ=(2311)

解答2

(1) 固有値を求める

行列式では列基本変形も行えるというのを忘れないようにしましょう。

固有値を t とする。固有多項式は、|AtE|=|1t5831t46613t|=|1t4t834t46013t|=(4t)|1t183146013t|=(4t)|4t0123146013t|=(4t)|4t12613t|と変形できる。

ここで、|4t12613t|=(t13)(t+4)+72=t29t+20=(t4)(t5)=0となるので、固有方程式は|AtE|=(4t)(t4)(t5)=(t4)2(t5)を満たす t が固有値となり、固有値は4(2重解)、5となる。

(2) 固有ベクトルを求める

固有値が4のときの固有ベクトルは、(A4E)=(558334669) (558334001) (550330001) (110001000) となる。{xy=0z=0を解くと、任意定数 k を用いて(xyz)=k(110)と表せるので、固有ベクトルは1本あり、固有ベクトル p1 は、p1=(110)となる。

(2重解なのに固有ベクトルが1本なので、残り1本は広義固有ベクトル p2 は(3)で求める)

固有値が5のときの固有ベクトルは、(A5E)=(658344668) (030344020) (344020000) (304010000)となる。{3x+4z=0y=0を解くと、任意定数 k を用いて(xyz)=k(403)と表せるので、固有ベクトルは1本あり、固有ベクトル p3 は、p3=(403)となる。

よって固有空間はV(4)={(110)}   V(5)={(403)}となる。

(3) 広義固有ベクトルを求める

固有値4に対する広義固有ベクトルが1本不足しているので、(A4E)p2=p1を満たすようなベクトルを求める。

固有値が4のときの固有ベクトルは、(A4E)=(558133416690) (558133410012) (5501533090012) (110300120000) となる。{xy=3z=2の解の1つに(xyz)=(302)があるので、広義固有ベクトル p2p2=(302)となります。
(広義固有ベクトルは p1, p2, p3 の3つ)

また、広義固有空間はV~(4)={(110),(302)}V~(5)={(403)}となる。

(4) ジョルダン標準形を求める

よって、{(A4E)p1=0Ap1=4p1   {(A4E)p2=p1Ap2=p1+4p2   {(A5E)p3=0Ap3=5p3が成立する。

ここでP=(p1,p2,p3)=(134100023)とすると、AP=A(p1,p2,p3)=(4p1,p1+4p2,5p3)=(p1,p2,p3)(410040005)=P(410040005)となります。

よって、P=(134100023)   J=P1AP=(410040005)とジョルダン標準形にできる。

[検算] AP=PJ=(413204100815)

解答3

(1) 固有値を求める

行列式では列基本変形も行えるというのを忘れないようにしましょう。

固有値を t とする。固有多項式は、|AtE|=|3t3111t1264t|=|2t2t011t1264t|=(2t)|11011t1264t|=(2t)|1100t1044t|=(2t)|t144t|と変形できる。

ここで、|t144t|=t(t4)+4=t24t+4=(t2)2=0となるので、固有方程式は|AtE|=(2t)(t2)2=(t2)2=0を満たす t が固有値となり、固有値は2(3重解)となる。

(2), (3) 固有ベクトル・広義固有ベクトルを求める

固有値2に属する固有ベクトルは、(A2E)= (131131262) (131000000)となる。x+3y+z=0を解くと(要注意)、任意定数 s,t を用いて(xyz)=s(101)+t(310)となります。なので、固有ベクトルを2本求めることができます。
(ただしこの時点では解が出せないので注意)

なので、定義にしたがって、(A2E)2p3=0を満たすベクトルを探します。(A4E)2=(000000000)=Oとなるので、どのような広義固有ベクトル p3 を持ってきてもOKなことがわかりますね。

ここで、p3=(100)とする*5(A2E)p3=p2を満たすような固有ベクトル p2 を探す。

すると、(A4E)p3=(131131262)(100)=(112)=p2となります。これは、s(101)+t(310)s=2, t=1 に相当します。

よって、(A4E)p3=p2を成立させるためには広義固有ベクトル p2, p3p2=(112),   p3=(100)とする必要があります。

残りの固有ベクトル p1 は先ほどの解の (s,t)=(1,0) もしくは (s,t)=(0,1) どちらを選んでもOK。

なので、p1=(101)とします。

よって固有ベクトルは p1, p2 となり、固有空間はV(2)={(101),(112)}となる。

また、広義固有ベクトルは p1, p2, p3 となり、広義固有空間はV(2)={(101),(112),(100)}となる。

(4) ジョルダン標準形を求める

よって、{(A2E)p1=0Ap1=2p1   {(A2E)p2=0Ap2=2p2   {(A2E)p3=p2Ap3=p2+2p3が成立する。

ここでP=(p1,p2,p3)=(111010120)とすると、AP=A(p1,p2,p3)=(2p1,2p2,p2+2p3)=(p1,p2,p3)(200021002)=P(200021002)となります。

よって、P=(111010120)   J=P1AP=(200021002)とジョルダン標準形にできる。

[検算] AP=PJ=(223021242)

8.さいごに

今回は、ジョルダン標準形についての簡単な説明と、3次までのジョルダン標準形の求め方を例題などを用いてまとめました。

ジョルダン標準形を求める際には、固有ベクトルを求めていき、足りなくなったら1つ前の固有ベクトル(もしくは広義固有ベクトル)を借りて(AtE)pi=pi1のようにドミノ倒しのように求めていくものだというのを覚えておきましょう。

また、ドミノ倒し的にもとめていくため、1つ前の固有ベクトル(もしくは広義固有ベクトル)が間違っていくと連鎖的に間違ってしまうため、慎重に計算していきましょう。

次回はジョルダン標準形の応用として行列の n 乗の求め方について簡単にですがまとめていきたいと思います。

*1:2次正方行列の場合、固有ベクトルが2つ求められれば対角化可能。n 次正方行列に対し、n 個の固有ベクトルが求められれば対角化が可能。

*2:例題1のように任意定数を使った書き方だと、(xyz)=(001)+k(100)となる。また、解を1つ選ぶときに他の広義固有ベクトルと1次独立になるように解を選ぶ必要があることに要注意。

*3:(A3E)p3=(A3E)2p2=(A3E)3p1=0を求めるのと同じ計算(上の定義の k=3 と同じ)。

*4:とにかく p1, p2, p3 を満たし、(xyz)=s(110)+t(101)を満たせばどんなベクトルを選んでもOK

*5:s(101)+t(310)s,t をどう頑張ってもつくれないようなものを p3 にすれば p1, p2, p3 を1次独立にできる。

関連広告・スポンサードリンク

おすすめの記事