広義固有ベクトル

線型代数学において，//n// × //n//}} 行列 A}} の'//広義'//固有ベクトルは，固有ベクトルの定義を緩めたある条件を満たすベクトルであるBronson．

V}} を n}} 次元ベクトル空間とする．φ}} を V}} から V}} への線型写像とする．A}} をある基底についての φ}} の行列表示とする．

V}} の完全な基底をなす A}} の n}} 個の線型独立な固有ベクトルがいつも存在するとは限らない．つまり，行列 A}} は対角化可能とは限らないBeauregardNering．これは少なくとも1つの固有値 λsub の代数的重複度がその幾何学的重複度よりも大きいときに起こる．この場合，λsub は不足固有値と呼ばれ，A}} は不足行列と呼ばれるGolub．

λsub に対応する広義固有ベクトル xsub は，行列 //A// − //λsub とあわせて，V}} の不変部分空間の基底をなす線型独立な広義固有ベクトルのジョルダン鎖を生成するBeauregardBronsonGolub．//

広義固有ベクトルを用いて，A}} の線型独立な固有ベクトルの集合を必要ならば V}} の完全な基底に拡張できるBronson．この基底は A}} に相似なジョルダン標準形にある「ほとんど対角な行列」J}} を決定するのに用いることができ，これは A}} のある行列関数を計算するのに有用であるBronson．行列 J}} は A}} が対角化可能とは限らないときに線形微分方程式系 1=// = mvar を解く際にも有用であるBeauregardNering．

1. 概要と定義
--------

固有ベクトルを定義するいくつかの同値な方法があるAntonBeauregardBurdenGolubHarperHersteinKreyszigNering．//n// × //n//}} 行列 A}} の固有値 λ}} の固有ベクトル x}} とは，1=// = 0}} なる零でないベクトルである，ただし I}} は //n// × //n//}} の単位行列であり，0}} は n}} 次元の零ベクトルであるBurden．つまり，//}} は変換 //A// − //λI//}} の核に属する．A}} が n}} 個の線型独立な固有ベクトルを持てば，A}} は対角行列 D}} に相似である．つまり，ある可逆行列 M}} が存在して，A}} は相似変換 1=//D// = //M//sup により対角化可能であるBeauregardBronson．行列 D}} は A}} のスペクトル行列と呼ばれる．行列 M}} は A}} の モード行列 と呼ばれるBronson．対角化可能な行列は，その行列関数が容易に計算できるなどの特長があるBronson．

一方，//n// × //n//}} 行列 A}} が n}} 個の線型独立な固有ベクトルを持たないとき，A}} は対角化可能ではないBeauregardBronson．

定義: ベクトル //sub が行列 A}} の固有値 λ}} に対応する'//階数 m}} の広義'//固有ベクトル  であるとは， を満たすが， であることをいうBronson．明らかに，階数 1 の広義固有ベクトルは通常の固有ベクトルであるBronson．すべての //n// × //n//}} 行列 A}} は n}} 個の線型独立な広義固有ベクトルを持ち，ジョルダン標準形の「ほとんど対角」な行列 J}} に相似であることを示すことができるBronson．つまり，可逆行列 M}} が存在して，1=//J// = //M//sup となるBronson．このときの行列 M}} は A}} の 広義モード行列 と呼ばれるBronson．λ}} が代数的重複度 μ}} の固有値ならば，A}} は λ}} に対応する μ}} 個の線型独立な広義固有ベクトルを持つBronson．これらの結果は，A}} のある種の行列関数を簡易に計算する際に有用となるBronson．

与えられた λ}} に対するすべての広義固有ベクトルによって張られる集合は λ}} の'//広義'//固有空間  をなすNering．

2. 例
----

広義固有ベクトルの概念を説明するいくつかの例を挙げる．詳細のいくつかは後で記述される．

2.1. 例 1
--------

まず，固有値が重複しても異なる固有ベクトルが得られる例について示す．

とする．A}} の固有値は，1=det = 0}} を満たす λ}} であり，これを解くと， となり，ただ1つの固有値 1=//λ// = 1}} が得られる．その代数的重複度は である．この固有値 1=//λ// = 1}} に対する固有ベクトルを求めてみよう． とおき， を満たすゼロでないベクトルを求めると，

となり，//x//sub, //x//sub とも任意でよい．たとえば， も，もいずれも固有値 1=//λ// = 1}} に対する固有ベクトルとなる．この2つの固有ベクトルは互いに線形独立である．

なお，

であり，幾何学的重複度は 1=//γ// = 2}} である．実際に，固有値 1=//λ// = 1}} に対して，2つの互いに線形独立な固有ベクトルが得られることがわかる．

2.2. 例 2
--------

固有値が重複する場合に異なる固有ベクトルが得られない例について示すNeringHersteinBeauregard．

とする．A}} の固有値は，1=det = 0}} を満たす λ}} であり，これを解くと，1=sup となり，ただ1つの固有値 1=//λ// = 1}} が得られる．その代数的重複度は 1=//μ// = 2}} である．例 1 と異なり，

で，幾何的重複度は 1=//γ// = 1}} である．すなわち，例 1 とは異なり，固有値 1=//λ// = 1}} に対する固有ベクトルは1つしかない．

では，この固有値 1=//λ// = 1}} に対する固有ベクトルを求めよう．今， とおき，1=// = 0}} を満たすゼロでないベクトルを求めると，

となり，//x//sub は任意であるが，1=//x//sub でなくてはならない．したがって，たとえば， は固有値 1=//λ// = 1}} に対する固有ベクトルとなる．なお，//x//sub は任意であるから， もまた固有ベクトルであるが， と は互いに独立ではないことにも注意されたい．

つぎに，この固有ベクトル を //sub として，一般固有ベクトル //sub を求めよう．これは， とおき， を解くことによって求めることができる．具体的には，

を解くと，//x//sub は任意であり，1=//x//sub となる．すると階数 2 の一般固有ベクトルは である，ただし a}} は任意のスカラー値である．1=//a// = 0}} とするのが最も単純である．

//sub と //sub は線型独立であり，ベクトル空間 V}} の基底をなす．

行列 A}} は，対角化可能ではないことに注意されたい．この行列は1つの 優対角成分があるから，階数が 1 よりも大きい一般化固有ベクトルが1つある．あるいは，}} の零空間の次元が 1=//p// = 1}} であることを計算でき，したがって階数が 1 よりも大きい広義固有ベクトルは 1=//m// − //p// = 1}} 個ある．

さて，求めた固有ベクトルと固有ベクトルを並べた行列

に対して，

となり，A}} は対角化はできていないが，J}} には A}} の固有値が対角成分に現れ，右上に“1” が配置されたジョルダン標準形となっていることがわかる．

2.3. 例 3
--------

この例は例 2 よりも複雑である．低い次数のよい例題を構成することはやや少し難しいNering．

行列

の固有値は

の解なので，固有値 1=//λ//sub と 1=//λ//sub を持ち，その代数的重複度はそれぞれ 1=//μ//sub と 1=//μ//sub である．1=//λ//sub に対して，

となるので，幾何学的重複度は 1=//γ//sub である．1=//λ//sub に対して，

となるので，幾何的重複度は 1=//γ//sub である．

はじめに，1=//λ//sub に対する固有ベクトル //sub を求める．幾何学的重複度は 1=//γ//sub なので，残りの一般固有ベクトルは，//sub から逐次的に求められる．具体的には，次のように求めた．

つぎに，1=//λ//sub に対する固有ベクトル //sub を求め，幾何学的重複度は 1=//γ//sub なので，残りの一般固有ベクトル //sub と //sub は，//sub から逐次的に求められる．具体的には，次のように求めた．

。これは A}} の各広義固有空間の基底となる．広義固有ベクトルの2つの鎖と合わせて5次元列ベクトル全体の空間を張る．

A}} に相似な「ほぼ対角」なジョルダン標準形の行列 J}} は以下のようにして得られる：

ただし M}} は A}} の 広義モード行列 であり，M}} の列は A}} の標準基底であり，1=//AM// = //MJ//}} である。

3. ジョルダン鎖
---------

定義: //sub を行列 A}} の固有値 λ}} に対応する階数 m}} の広義固有ベクトルとする．//sub によって生成される鎖とは次で与えられるベクトルの集合 である：

したがって，一般に，

 によって与えられるベクトル xsub は固有値 λ}} に対応する階数 j}} の広義固有ベクトルである．鎖はベクトルの線型独立な集合であるBronson．

4. 標準基底
-------

定義：n}} 個の線型独立な広義固有ベクトルの集合が標準基底  であるとは，ジョルダン鎖の全体からなることをいう．

したがって，階数 m}} の広義固有ベクトルが標準基底に入っていることを一度決定すれば，xsub によって生成されるジョルダン鎖に入っている //m// − 1}} 個のベクトル も標準基底に入っていることが従うBronson．

λsub を A}} の代数的重複度 μsub の固有値とする．まず，行列 の階数を求める．整数 msub は の階数が //n// − //μsub となる「最初の整数」として決定される．//

さて

と定義する．変数 ρsub は A}} の標準基底に現れる固有値 λsub に対応する階数 k}} の線型独立な広義固有ベクトルの個数を表す．

に注意Bronson．

5. 広義固有ベクトルの計算
--------------

これまでの節で //n// × //n//}} 行列 A}} に付随するベクトル空間 V}} に対する標準基底の n}} 個の線型独立な広義固有ベクトルを得る方法を見た．これらの方法を結合して手順を得る：

固有値 λsub と代数的重複度 μsub に対する A}} の特性方程式を解く；各 λsub に対して：//n// − //μsub を決定する；//msub を決定する；1=//k// = 1, ..., //msub に対して ρsub を決定する；//λsub に対して各ジョルダン鎖を決定する；

5.1. 例 4
--------

行列

の固有値は

の解であり，これを解くと 1=//λ//subおよび 1=//λ//subが得られる．また，1=//n// = 4}} である．1=//λ//sub に対して 1=//n// − //μ//sub である．

の階数が 1=//n// − //μ//sub になる最初の整数 //m//sub は 1=//m//sub である．

次のように定義する：

したがって，3つの線型独立な広義固有ベクトルが存在する；階数 3, 2, 1 に1つずつである．//λ//sub は3つの線型独立な広義固有ベクトルのただ1つの鎖に対応するので，//λ//sub に対応する階数 3 の広義固有ベクトルであって

なるものが存在することを知っている．方程式  と  は //sub について解くことができる線型方程式系を表す．

とする．すると

および

である．したがって，条件  と  を満たすためには，1=//x//sub かつ 1=//x//sub でなければならない．//x//sub と //x//sub には何の制約もない．1=//x//sub, 1=//x//sub と選ぶことで，1=//λ//sub に対応する階数 3 の広義固有ベクトルとして

を得る．//x//sub, //x//sub, //x//sub で //x//sub なる異なる値を選ぶことによって階数 3 の他の広義固有ベクトルを無限個得ることができることに注意．しかしながら，我々の最初の選択が最も単純であるBronson．

さて方程式  を用いて，xsub と xsub をそれぞれ階数 2 と 1 の広義固有ベクトルとして得る，ただし

および

である．代数的重複度が 1 の固有値 1=//λ//sub は標準的な手法で扱うことができ，通常の固有ベクトル

を持つ．A}} の標準基底は

である．//sub は //λ//sub に伴う広義固有ベクトルである．//sub は //λ//sub に伴う通常の固有ベクトルである．

これはかなり単純な例であることに注意すべきである．一般に，階数 k}} の線型独立な広義固有ベクトルの個数 ρsub は必ずしも等しくない．つまり，特定の固有値に対応する異なる長さのいくつかの鎖があるかもしれないBronson．

6. 広義モード行列
----------

A}} を //n// × //n//}} 行列とする．A}} の広義モード行列  M}} とは，//n// × //n//}} 行列であって，その列が，ベクトルと考えたときに，A}} の標準基底をなし，M}} において以下の規則に従って現れるものをいう：

1つのベクトルからなるすべてのジョルダン鎖は M}} のはじめの列に現れる．
1つの鎖のすべてのベクトルは M}} の隣接する列に一緒に現れる．
各鎖は M}} において階数が増える順番で現れるBronson．

7. ジョルダン標準形
-----------

V}} を n}} 次元ベクトル空間とする；φ}} を V}} から自身への線型写像全体の集合 End}} の元とする；A}} をある基底に関する φ}} の行列表示とする．次のことを示すことができる．A}} の特性多項式 //f//}} が一次式に分解して

の形，ただし は A}} の相異なる固有値，になれば，各 μsub は対応する固有値 λsub の代数的重複度であり，A}} はジョルダン標準形の行列 J}} に相似である，ただし各 λsub は対角線上連続した μsub 回現れ，各 λsub の上の各成分は 0 または 1 である；各 λsub の最初の出現の上の成分はつねに 0 である．すべての他の成分は 0 である．行列 J}} は A}} の対角化にできるだけ近い．A}} が対角化可能ならば，対角線の上のすべての成分は 0 であるBeauregard．教科書によっては優対角成分ではなく 劣対角成分, すなわち主対角線の直下に 1 たちがあることに注意．固有値はなお主対角線にあるCullenFranklin．

すべての //n// × //n//}} 行列 A}} は相似変換 1=//J// = //M//sup- によって得られるジョルダン標準形の行列 J}} に相似である，ただし M}} は A}} の広義モード行列であるBronson。

7.1. 例 5
--------

に相似なジョルダン標準形の行列を見つけよ．

解：A}} の特性方程式は 1=sup であるので，固有値は 1=//λ// = 2}}である．前の節の手順に従って，

と

が分かる．したがって，1=//ρ//sub と 1=//ρ//sub であり，A}} の標準基底は階数 2 の1つの線型独立な広義固有ベクトルと階数 1 の2つの線型独立な広義固有ベクトルを含むことが分かる，あるいは同じことだが，2つのベクトルの1つの鎖 mset と1つのベクトルの1つの鎖 mset を含む．1=//M// = (//sub と書いて，次が分かる：

および

ただし M}} は A}} の広義モード行列で，M}} の列は A}} の標準基底で，1=//AM// = //MJ//}} であるBronson．広義固有ベクトル自身は一意ではないから，また M}} と J}} の両方の列のいくつかは交換できるから，M}} と J}} はいずれも一意ではないことが従うことに注意Bronson．

7.2. 例 6
--------

例 4 において，行列 A}} に対する線型独立な広義固有ベクトルの標準基底を求めた．A}} の広義モード行列は

である．A}} に相似なジョルダン標準形の行列は

であり，1=//AM// = //MJ//}} である。

8. 応用
-----

8.1. 行列関数
---------

正方行列に実行できる最も基本的な演算の3つは，和とスカラー倍と積であるBeauregard．これらは //n// × //n//}} 行列 A}} の多項式関数を定義するのにちょうど必要な演算であるBronson．多くの関数がマクローリン級数として書けることを基本的な解析学から思い出すと，行列のより一般の関数をきわめて容易に定義できるBronson．A}} が対角化可能ならば，つまり

で

ならば，

であり，A}} の関数のマクローリン級数の計算は大きく単純化されるBronson．例えば，A}} の任意の冪 k}} を得るには，Dsup を計算し，M}} を左から掛け，さらに //M//sup を右から掛けるだけでよいBronson．

広義固有ベクトルを用いて，A}} のジョルダン標準形を得ることができ，これらの結果は対角化可能でない行列の関数を計算する直截的手法に一般化できるBronson．

8.2. 微分方程式
----------

次の線型常微分方程式系を解く問題を考える：

ただし

4}} および 4}} 行列 A}} が対角行列で //i// ≠ //j//}} に対して 1=//asub のとき，系  は次の形の n}} 個の方程式の系に簡約される：//

この場合，一般解は次で与えられる：

一般の場合には，A}} を対角化し系  を  のような系に以下のように簡約しようとする．A}} が対角化可能ならば，M}} を A}} のモード行列として，1=//D// = //M//sup である．1=//A// = //MDM//sup を代入して，方程式  は次の形となる： あるいは

ただし

 の解は

 の解 //x//}} はすると関係式  を用いて得られるBeauregard

一方，A}} が対角化可能でなければ，M}} を A}} の広義モード行列に選び，1=//J// = //M//sup を A}} のジョルダン標準形とする．系 1=//′ = mvar は次の形を持つ：

ただし λsub は J}} の主対角成分にある固有値であり，εsub は J}} の優対角成分にある 1 と 0 である．系  はしばしば  よりも容易に解かれる． の最後の方程式を ysub に対して解いて， を得る．次に ysub のこの解を  の最後から二番目の方程式に代入して，//y//sub に対して解く．この手順を続けて， を最後の方程式から最初までやり，//y//}} に対する全体の系を解く．すると解 //x//}} は関係式  を用いて得られるBeauregard．

9. 脚注
-----

10. 参考文献
--------

11. 外部リンク
---------