Vector satisfying some of the criteria of an eigenvector
線型代数学 において 、 行列 の 一般化固有ベクトルとは、(通常の) 固有ベクトル よりも緩やかな基準を満たす ベクトル のことである 。 [1] n × n {\displaystyle n\times n} A {\displaystyle A}
を - 次元 ベクトル空間 とし 、をある順序付き 基底 に関するから へ の線型写像の 行列表現 とします 。 V {\displaystyle V} n {\displaystyle n} A {\displaystyle A} V {\displaystyle V} V {\displaystyle V}
の完全な基底を形成するの 線形独立な 固有ベクトル の完全な集合が常に存在するとは限らない 。つまり、行列 は 対角化 できない可能性がある 。 [2] [3] これは、 少なくとも1つの 固有値の代数 的重複度 が、その幾何 学的重複度 ( 行列 の 零値 、または その 零空間の 次元 )よりも大きい場合に発生する。この場合、は 欠陥のある固有値 と呼ばれ 、は 欠陥のある行列 と呼ばれる 。 [4] n {\displaystyle n} A {\displaystyle A} V {\displaystyle V} A {\displaystyle A} λ i {\displaystyle \lambda _{i}} ( A − λ i I ) {\displaystyle (A-\lambda _{i}I)} λ i {\displaystyle \lambda _{i}} A {\displaystyle A}
に対応する 一般化固有ベクトルは 、行列とともに、 不変部分空間 の 基底を形成する線形独立な一般化固有ベクトルのジョルダン連鎖を生成する 。 [5] [6] [7] x i {\displaystyle x_{i}} λ i {\displaystyle \lambda _{i}} ( A − λ i I ) {\displaystyle (A-\lambda _{i}I)} V {\displaystyle V}
一般化固有ベクトルを用いると、 の線形独立な固有ベクトルの集合は 、必要であれば の完全な基底に拡張することができる 。 [8] この基底は、に 類似した ジョルダン標準形 の 「ほぼ対角行列」を決定するために用いることができ 、これは の特定の 行列関数 を計算する際に有用である 。 [9] この行列は、が対角化可能 である必要がない 線形微分方程式系 を解く際にも有用である 。 [10] [11] A {\displaystyle A} V {\displaystyle V} J {\displaystyle J} A {\displaystyle A} A {\displaystyle A} J {\displaystyle J} x ′ = A x , {\displaystyle \mathbf {x} '=A\mathbf {x} ,} A {\displaystyle A}
与えられた固有値に対応する一般化固有空間の次元は、 の代数的重複度である 。 [12] λ {\displaystyle \lambda } λ {\displaystyle \lambda }
概要と定義 通常の固有ベクトル を定義する同等な方法がいくつかあります 。 [13] [14] [15] [16] [17] [18] [19] [20] ここでは、 × 行列 の 固有値に関連付けられた固有ベクトルは の 非ゼロベクトルであり 、 は × 単位行列 、は 長さの ゼロベクトル です 。 [21] つまり、は 変換 の 核 内にあります 。 が 線形独立な固有ベクトル を持つ場合、 は 対角行列 に相似です 。つまり、 のような 可逆行列 が存在し、 は相似変換 により対角化可能です 。 [22] [23] 行列は に対して スペクトル行列 と呼ばれます 。行列は に対して モーダル行列 と呼ばれます 。 [24] 対角化可能な行列は、その行列関数が簡単に計算できるため、特に重要です。 [25] u {\displaystyle \mathbf {u} } λ {\displaystyle \lambda } n {\displaystyle n} n {\displaystyle n} A {\displaystyle A} ( A − λ I ) u = 0 {\displaystyle (A-\lambda I)\mathbf {u} =\mathbf {0} } I {\displaystyle I} n {\displaystyle n} n {\displaystyle n} 0 {\displaystyle \mathbf {0} } n {\displaystyle n} u {\displaystyle \mathbf {u} } ( A − λ I ) {\displaystyle (A-\lambda I)} A {\displaystyle A} n {\displaystyle n} A {\displaystyle A} D {\displaystyle D} M {\displaystyle M} A {\displaystyle A} D = M − 1 A M {\displaystyle D=M^{-1}AM} D {\displaystyle D} A {\displaystyle A} M {\displaystyle M} A {\displaystyle A}
一方、 線形独立な固有ベクトルが関連付けられて いない場合、 対角化できません。 [26] [27] A {\displaystyle A} n {\displaystyle n} A {\displaystyle A}
定義: ベクトルは 行列の 階数 m の一般化固有ベクトルであり 、次の式で表される固有値に対応する 。 x m {\displaystyle \mathbf {x} _{m}} A {\displaystyle A} λ {\displaystyle \lambda }
( A − λ I ) m x m = 0 {\displaystyle (A-\lambda I)^{m}\mathbf {x} _{m}=\mathbf {0} } しかし
( A − λ I ) m − 1 x m ≠ 0 . {\displaystyle (A-\lambda I)^{m-1}\mathbf {x} _{m}\neq \mathbf {0} .} [28] 明らかに、階数1の一般化固有ベクトルは通常の固有ベクトルである。 [29] すべての × 行列に は、線形独立な一般化固有ベクトルが関連付けられており、 ジョルダン標準形の 「ほぼ対角」行列に類似していることが示される。 [30] つまり、 となる可逆行列が存在する 。 [31] この場合の 行列は、 の 一般化様相行列 と呼ばれる。 [32] が代数的重複度 の固有値である 場合 、 は に対応する線形独立な一般化固有ベクトル を持つ 。 [33] これらの結果は、 の特定の行列関数を計算するための簡単な方法を提供する 。 [34] n {\displaystyle n} n {\displaystyle n} A {\displaystyle A} n {\displaystyle n} J {\displaystyle J} M {\displaystyle M} J = M − 1 A M {\displaystyle J=M^{-1}AM} M {\displaystyle M} A {\displaystyle A} λ {\displaystyle \lambda } μ {\displaystyle \mu } A {\displaystyle A} μ {\displaystyle \mu } λ {\displaystyle \lambda } A {\displaystyle A}
注: 体 上の 行列を ジョルダン標準形で表すには、 のすべての固有値 が に含まれていなければならない 。つまり、 特性多項式は 完全に線形因子に因数分解できなければならず、 は代数的に閉体でなければならない。例えば、 が 実数値 要素を持つ場合、固有値と固有ベクトルの成分が 複素数値 を 持つ必要があるかもしれない 。 [35] [36] [37] n × n {\displaystyle n\times n} A {\displaystyle A} F {\displaystyle F} A {\displaystyle A} F {\displaystyle F} f ( x ) {\displaystyle f(x)} F {\displaystyle F} A {\displaystyle A}
与えられたすべての一般化固有ベクトルによって 張られる 集合は、 の 一般化固有空間 を形成する 。 [38] λ {\displaystyle \lambda } λ {\displaystyle \lambda }
例 一般化固有ベクトルの概念を説明するために、いくつかの例を挙げます。詳細については後ほど説明します。
例1 この例は 単純ですが、要点を明確に示しています。このタイプの行列は教科書で頻繁に使用されます。 [39] [40] [41]
A = ( 1 1 0 1 ) . {\displaystyle A={\begin{pmatrix}1&1\\0&1\end{pmatrix}}.} そうすると、固有値は 1 つだけとなり、 その代数的重複度は となります 。 λ = 1 {\displaystyle \lambda =1} m = 2 {\displaystyle m=2}
この行列はジョルダン正規形ですが、 対角行列 ではないことに注意してください。したがって、この行列は対角化できません。 上対角 要素が1つあるため、階数が1より大きい一般化固有ベクトルが1つ存在します(あるいは、ベクトル空間が2次元であるため、階数が1より大きい一般化固有ベクトルは最大で1つしか存在しないことに注意することもできます)。あるいは、 の 零空間 の次元を と計算すると 、 階数が1より大きい一般化固有ベクトルが
存在すること になります。 V {\displaystyle V} A − λ I {\displaystyle A-\lambda I} p = 1 {\displaystyle p=1} m − p = 1 {\displaystyle m-p=1}
通常の固有ベクトル は通常通り計算されます(例については 固有ベクトルの ページを参照してください)。この固有ベクトルを用いて、一般化固有ベクトルを 次のように
計算します。 v 1 = ( 1 0 ) {\displaystyle \mathbf {v} _{1}={\begin{pmatrix}1\\0\end{pmatrix}}} v 2 {\displaystyle \mathbf {v} _{2}}
( A − λ I ) v 2 = v 1 . {\displaystyle (A-\lambda I)\mathbf {v} _{2}=\mathbf {v} _{1}.} 値を書き出す:
( ( 1 1 0 1 ) − 1 ( 1 0 0 1 ) ) ( v 21 v 22 ) = ( 0 1 0 0 ) ( v 21 v 22 ) = ( 1 0 ) . {\displaystyle \left({\begin{pmatrix}1&1\\0&1\end{pmatrix}}-1{\begin{pmatrix}1&0\\0&1\end{pmatrix}}\right){\begin{pmatrix}v_{21}\\v_{22}\end{pmatrix}}={\begin{pmatrix}0&1\\0&0\end{pmatrix}}{\begin{pmatrix}v_{21}\\v_{22}\end{pmatrix}}={\begin{pmatrix}1\\0\end{pmatrix}}.} これを単純化すると
v 22 = 1. {\displaystyle v_{22}=1.} 元に は制約はありません。階数2の一般化固有ベクトルは となり 、ここで a は任意のスカラー値を取ることができます。通常、 a = 0を選択するのが 最も簡単です。 v 21 {\displaystyle v_{21}} v 2 = ( a 1 ) {\displaystyle \mathbf {v} _{2}={\begin{pmatrix}a\\1\end{pmatrix}}}
ご了承ください
( A − λ I ) v 2 = ( 0 1 0 0 ) ( a 1 ) = ( 1 0 ) = v 1 , {\displaystyle (A-\lambda I)\mathbf {v} _{2}={\begin{pmatrix}0&1\\0&0\end{pmatrix}}{\begin{pmatrix}a\\1\end{pmatrix}}={\begin{pmatrix}1\\0\end{pmatrix}}=\mathbf {v} _{1},} これ は一般化固有ベクトルです。 v 2 {\displaystyle \mathbf {v} _{2}}
( A − λ I ) 2 v 2 = ( A − λ I ) [ ( A − λ I ) v 2 ] = ( A − λ I ) v 1 = ( 0 1 0 0 ) ( 1 0 ) = ( 0 0 ) = 0 , {\displaystyle (A-\lambda I)^{2}\mathbf {v} _{2}=(A-\lambda I)[(A-\lambda I)\mathbf {v} _{2}]=(A-\lambda I)\mathbf {v} _{1}={\begin{pmatrix}0&1\\0&0\end{pmatrix}}{\begin{pmatrix}1\\0\end{pmatrix}}={\begin{pmatrix}0\\0\end{pmatrix}}=\mathbf {0} ,} したがって、 は 通常の固有ベクトルであり、 と は 線形独立であるため、ベクトル空間 の基底を構成します 。 v 1 {\displaystyle \mathbf {v} _{1}} v 1 {\displaystyle \mathbf {v} _{1}} v 2 {\displaystyle \mathbf {v} _{2}} V {\displaystyle V}
例2 この例は例1よりも複雑です。残念ながら、低次の興味深い例を構築するのは少し難しいです。 [42] 行列
A = ( 1 0 0 0 0 3 1 0 0 0 6 3 2 0 0 10 6 3 2 0 15 10 6 3 2 ) {\displaystyle A={\begin{pmatrix}1&0&0&0&0\\3&1&0&0&0\\6&3&2&0&0\\10&6&3&2&0\\15&10&6&3&2\end{pmatrix}}} には、 代数的重複度 とを伴う 固有値と があります が 、 幾何学的重複度は 、およびです 。 λ 1 = 1 {\displaystyle \lambda _{1}=1} λ 2 = 2 {\displaystyle \lambda _{2}=2} μ 1 = 2 {\displaystyle \mu _{1}=2} μ 2 = 3 {\displaystyle \mu _{2}=3} γ 1 = 1 {\displaystyle \gamma _{1}=1} γ 2 = 1 {\displaystyle \gamma _{2}=1}
の 一般 化固有空間 は以下のように計算されます。 は に関連付けられた通常の固有ベクトルです 。 は に関連付けられた一般化固有ベクトルです 。 は に関連付けられた通常の固有ベクトルです 。 および は に関連付けられた一般化固有ベクトルです 。 A {\displaystyle A} x 1 {\displaystyle \mathbf {x} _{1}} λ 1 {\displaystyle \lambda _{1}} x 2 {\displaystyle \mathbf {x} _{2}} λ 1 {\displaystyle \lambda _{1}} y 1 {\displaystyle \mathbf {y} _{1}} λ 2 {\displaystyle \lambda _{2}} y 2 {\displaystyle \mathbf {y} _{2}} y 3 {\displaystyle \mathbf {y} _{3}} λ 2 {\displaystyle \lambda _{2}}
( A − 1 I ) x 1 = ( 0 0 0 0 0 3 0 0 0 0 6 3 1 0 0 10 6 3 1 0 15 10 6 3 1 ) ( 0 3 − 9 9 − 3 ) = ( 0 0 0 0 0 ) = 0 , {\displaystyle (A-1I)\mathbf {x} _{1}={\begin{pmatrix}0&0&0&0&0\\3&0&0&0&0\\6&3&1&0&0\\10&6&3&1&0\\15&10&6&3&1\end{pmatrix}}{\begin{pmatrix}0\\3\\-9\\9\\-3\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\0\\0\end{pmatrix}}=\mathbf {0} ,} ( A − 1 I ) x 2 = ( 0 0 0 0 0 3 0 0 0 0 6 3 1 0 0 10 6 3 1 0 15 10 6 3 1 ) ( 1 − 15 30 − 1 − 45 ) = ( 0 3 − 9 9 − 3 ) = x 1 , {\displaystyle (A-1I)\mathbf {x} _{2}={\begin{pmatrix}0&0&0&0&0\\3&0&0&0&0\\6&3&1&0&0\\10&6&3&1&0\\15&10&6&3&1\end{pmatrix}}{\begin{pmatrix}1\\-15\\30\\-1\\-45\end{pmatrix}}={\begin{pmatrix}0\\3\\-9\\9\\-3\end{pmatrix}}=\mathbf {x} _{1},} ( A − 2 I ) y 1 = ( − 1 0 0 0 0 3 − 1 0 0 0 6 3 0 0 0 10 6 3 0 0 15 10 6 3 0 ) ( 0 0 0 0 9 ) = ( 0 0 0 0 0 ) = 0 , {\displaystyle (A-2I)\mathbf {y} _{1}={\begin{pmatrix}-1&0&0&0&0\\3&-1&0&0&0\\6&3&0&0&0\\10&6&3&0&0\\15&10&6&3&0\end{pmatrix}}{\begin{pmatrix}0\\0\\0\\0\\9\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\0\\0\end{pmatrix}}=\mathbf {0} ,} ( A − 2 I ) y 2 = ( − 1 0 0 0 0 3 − 1 0 0 0 6 3 0 0 0 10 6 3 0 0 15 10 6 3 0 ) ( 0 0 0 3 0 ) = ( 0 0 0 0 9 ) = y 1 , {\displaystyle (A-2I)\mathbf {y} _{2}={\begin{pmatrix}-1&0&0&0&0\\3&-1&0&0&0\\6&3&0&0&0\\10&6&3&0&0\\15&10&6&3&0\end{pmatrix}}{\begin{pmatrix}0\\0\\0\\3\\0\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\0\\9\end{pmatrix}}=\mathbf {y} _{1},} ( A − 2 I ) y 3 = ( − 1 0 0 0 0 3 − 1 0 0 0 6 3 0 0 0 10 6 3 0 0 15 10 6 3 0 ) ( 0 0 1 − 2 0 ) = ( 0 0 0 3 0 ) = y 2 . {\displaystyle (A-2I)\mathbf {y} _{3}={\begin{pmatrix}-1&0&0&0&0\\3&-1&0&0&0\\6&3&0&0&0\\10&6&3&0&0\\15&10&6&3&0\end{pmatrix}}{\begin{pmatrix}0\\0\\1\\-2\\0\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\3\\0\end{pmatrix}}=\mathbf {y} _{2}.} これにより、 の一般化固有空間 の それぞれに対する基底が得られます 。2つの一般化固有ベクトルの 鎖 を合わせると、すべての5次元列ベクトルの空間が広がります。 A {\displaystyle A}
{ x 1 , x 2 } = { ( 0 3 − 9 9 − 3 ) , ( 1 − 15 30 − 1 − 45 ) } , { y 1 , y 2 , y 3 } = { ( 0 0 0 0 9 ) , ( 0 0 0 3 0 ) , ( 0 0 1 − 2 0 ) } . {\displaystyle \left\{\mathbf {x} _{1},\mathbf {x} _{2}\right\}=\left\{{\begin{pmatrix}0\\3\\-9\\9\\-3\end{pmatrix}},{\begin{pmatrix}1\\-15\\30\\-1\\-45\end{pmatrix}}\right\},\left\{\mathbf {y} _{1},\mathbf {y} _{2},\mathbf {y} _{3}\right\}=\left\{{\begin{pmatrix}0\\0\\0\\0\\9\end{pmatrix}},{\begin{pmatrix}0\\0\\0\\3\\0\end{pmatrix}},{\begin{pmatrix}0\\0\\1\\-2\\0\end{pmatrix}}\right\}.} ジョルダン正規形 の 「ほぼ対角」行列は 、 次のようにして得られます。 J {\displaystyle J} A {\displaystyle A}
M = ( x 1 x 2 y 1 y 2 y 3 ) = ( 0 1 0 0 0 3 − 15 0 0 0 − 9 30 0 0 1 9 − 1 0 3 − 2 − 3 − 45 9 0 0 ) , {\displaystyle M={\begin{pmatrix}\mathbf {x} _{1}&\mathbf {x} _{2}&\mathbf {y} _{1}&\mathbf {y} _{2}&\mathbf {y} _{3}\end{pmatrix}}={\begin{pmatrix}0&1&0&0&0\\3&-15&0&0&0\\-9&30&0&0&1\\9&-1&0&3&-2\\-3&-45&9&0&0\end{pmatrix}},} J = ( 1 1 0 0 0 0 1 0 0 0 0 0 2 1 0 0 0 0 2 1 0 0 0 0 2 ) , {\displaystyle J={\begin{pmatrix}1&1&0&0&0\\0&1&0&0&0\\0&0&2&1&0\\0&0&0&2&1\\0&0&0&0&2\end{pmatrix}},} ここで は の 一般化モード行列 、 の列は の 標準 基底 、 で ある 。 [43] M {\displaystyle M} A {\displaystyle A} M {\displaystyle M} A {\displaystyle A} A M = M J {\displaystyle AM=MJ}
ジョーダンチェーン 定義: 行列 と固有値に対応する階数 m の一般化固有ベクトル とする 。 によって 生成される連鎖は、 次式で表される ベクトルの集合である。 x m {\displaystyle \mathbf {x} _{m}} A {\displaystyle A} λ {\displaystyle \lambda } x m {\displaystyle \mathbf {x} _{m}} { x m , x m − 1 , … , x 1 } {\displaystyle \left\{\mathbf {x} _{m},\mathbf {x} _{m-1},\dots ,\mathbf {x} _{1}\right\}}
ここで は常に与えられた固有値 を持つ通常の固有ベクトルである 。したがって、一般に、 x 1 {\displaystyle \mathbf {x} _{1}} λ {\displaystyle \lambda }
x j = ( A − λ I ) m − j x m = ( A − λ I ) x j + 1 ( j = 1 , 2 , … , m − 1 ) . {\displaystyle \mathbf {x} _{j}=(A-\lambda I)^{m-j}\mathbf {x} _{m}=(A-\lambda I)\mathbf {x} _{j+1}\qquad (j=1,2,\dots ,m-1).} 2
( 2 )で与えられるベクトルは、 固有値に対応する階数 j の一般化固有ベクトルである 。連鎖とは、線形独立なベクトルの集合である。 [44] x j {\displaystyle \mathbf {x} _{j}} λ {\displaystyle \lambda }
標準的根拠 定義: n 個の線形独立な一般化固有ベクトル の集合は 、完全にジョルダン連鎖で構成されている場合、 標準基底となります。
このように、階数m の一般化固有ベクトルが 標準基底にあると判断すれば、 ジョルダン連鎖で生成される m −1個のベクトル も標準基底にあることがわかります。 [45] x m − 1 , x m − 2 , … , x 1 {\displaystyle \mathbf {x} _{m-1},\mathbf {x} _{m-2},\ldots ,\mathbf {x} _{1}} x m {\displaystyle \mathbf {x} _{m}}
を代数的重複度 の の固有値とします 。 まず、 行列 の 階数 (行列階数)を求めます。整数 は、 の階数を持つ 最初 の整数 として決定されます ( n は の行数または列数 、つまり n × n です )。 λ i {\displaystyle \lambda _{i}} A {\displaystyle A} μ i {\displaystyle \mu _{i}} ( A − λ i I ) , ( A − λ i I ) 2 , … , ( A − λ i I ) m i {\displaystyle (A-\lambda _{i}I),(A-\lambda _{i}I)^{2},\ldots ,(A-\lambda _{i}I)^{m_{i}}} m i {\displaystyle m_{i}} ( A − λ i I ) m i {\displaystyle (A-\lambda _{i}I)^{m_{i}}} n − μ i {\displaystyle n-\mu _{i}} A {\displaystyle A} A {\displaystyle A}
定義する
ρ k = rank ( A − λ i I ) k − 1 − rank ( A − λ i I ) k ( k = 1 , 2 , … , m i ) . {\displaystyle \rho _{k}=\operatorname {rank} (A-\lambda _{i}I)^{k-1}-\operatorname {rank} (A-\lambda _{i}I)^{k}\qquad (k=1,2,\ldots ,m_{i}).} 変数は 、 の標準基底に現れる 固有値に対応する階数 k の線形独立な一般化固有ベクトルの数を表す 。 ρ k {\displaystyle \rho _{k}} λ i {\displaystyle \lambda _{i}} A {\displaystyle A}
rank ( A − λ i I ) 0 = rank ( I ) = n {\displaystyle \operatorname {rank} (A-\lambda _{i}I)^{0}=\operatorname {rank} (I)=n} 。 [46]
一般化固有ベクトルの計算 前のセクションでは、行列 に関連付けられた ベクトル空間の標準基底の線形独立な一般化固有ベクトルを 求める手法を見てきました 。これらの手法は、以下の手順にまとめることができます。 n {\displaystyle n} V {\displaystyle V} n × n {\displaystyle n\times n} A {\displaystyle A}
固有値 とその代数的重複度 についての 特性方程式 を解きます 。 A {\displaystyle A} λ i {\displaystyle \lambda _{i}} μ i {\displaystyle \mu _{i}} それぞれ λ i : {\displaystyle \lambda _{i}:} 決定する ; n − μ i {\displaystyle n-\mu _{i}} 決定する ; m i {\displaystyle m_{i}} を 決定します 。 ρ k {\displaystyle \rho _{k}} ( k = 1 , … , m i ) {\displaystyle (k=1,\ldots ,m_{i})} それぞれのジョーダン連鎖を決定します 。 λ i {\displaystyle \lambda _{i}}
例3 マトリックス
A = ( 5 1 − 2 4 0 5 2 2 0 0 5 3 0 0 0 4 ) {\displaystyle A={\begin{pmatrix}5&1&-2&4\\0&5&2&2\\0&0&5&3\\0&0&0&4\end{pmatrix}}} は代数的重複度 の 固有値と 代数的重複度 の 固有値を持ちます 。また も成り立ちます 。なぜなら が成り立つからです 。 λ 1 = 5 {\displaystyle \lambda _{1}=5} μ 1 = 3 {\displaystyle \mu _{1}=3} λ 2 = 4 {\displaystyle \lambda _{2}=4} μ 2 = 1 {\displaystyle \mu _{2}=1} n = 4 {\displaystyle n=4} λ 1 {\displaystyle \lambda _{1}} n − μ 1 = 4 − 3 = 1 {\displaystyle n-\mu _{1}=4-3=1}
( A − 5 I ) = ( 0 1 − 2 4 0 0 2 2 0 0 0 3 0 0 0 − 1 ) , rank ( A − 5 I ) = 3. {\displaystyle (A-5I)={\begin{pmatrix}0&1&-2&4\\0&0&2&2\\0&0&0&3\\0&0&0&-1\end{pmatrix}},\qquad \operatorname {rank} (A-5I)=3.} ( A − 5 I ) 2 = ( 0 0 2 − 8 0 0 0 4 0 0 0 − 3 0 0 0 1 ) , rank ( A − 5 I ) 2 = 2. {\displaystyle (A-5I)^{2}={\begin{pmatrix}0&0&2&-8\\0&0&0&4\\0&0&0&-3\\0&0&0&1\end{pmatrix}},\qquad \operatorname {rank} (A-5I)^{2}=2.} ( A − 5 I ) 3 = ( 0 0 0 14 0 0 0 − 4 0 0 0 3 0 0 0 − 1 ) , rank ( A − 5 I ) 3 = 1. {\displaystyle (A-5I)^{3}={\begin{pmatrix}0&0&0&14\\0&0&0&-4\\0&0&0&3\\0&0&0&-1\end{pmatrix}},\qquad \operatorname {rank} (A-5I)^{3}=1.} の 順位を持つ 最初の整数 は です 。 m 1 {\displaystyle m_{1}} ( A − 5 I ) m 1 {\displaystyle (A-5I)^{m_{1}}} n − μ 1 = 1 {\displaystyle n-\mu _{1}=1} m 1 = 3 {\displaystyle m_{1}=3}
ここで定義する
ρ 3 = rank ( A − 5 I ) 2 − rank ( A − 5 I ) 3 = 2 − 1 = 1 , {\displaystyle \rho _{3}=\operatorname {rank} (A-5I)^{2}-\operatorname {rank} (A-5I)^{3}=2-1=1,} ρ 2 = rank ( A − 5 I ) 1 − rank ( A − 5 I ) 2 = 3 − 2 = 1 , {\displaystyle \rho _{2}=\operatorname {rank} (A-5I)^{1}-\operatorname {rank} (A-5I)^{2}=3-2=1,} ρ 1 = rank ( A − 5 I ) 0 − rank ( A − 5 I ) 1 = 4 − 3 = 1. {\displaystyle \rho _{1}=\operatorname {rank} (A-5I)^{0}-\operatorname {rank} (A-5I)^{1}=4-3=1.} その結果、階数3、2、1のそれぞれ1つずつ、3つの線形独立な一般化固有ベクトルが存在することになる。 は 3つの線形独立な一般化固有ベクトルの単一の連鎖に対応するので、 階数3の一般化固有ベクトルが存在し 、 λ 1 {\displaystyle \lambda _{1}} x 3 {\displaystyle \mathbf {x} _{3}} λ 1 {\displaystyle \lambda _{1}}
( A − 5 I ) 3 x 3 = 0 {\displaystyle (A-5I)^{3}\mathbf {x} _{3}=\mathbf {0} } 3
しかし
( A − 5 I ) 2 x 3 ≠ 0 . {\displaystyle (A-5I)^{2}\mathbf {x} _{3}\neq \mathbf {0} .} 4
式( 3 )と式( 4 )は 、 について解ける 線形方程式 を表す。 x 3 {\displaystyle \mathbf {x} _{3}}
x 3 = ( x 31 x 32 x 33 x 34 ) . {\displaystyle \mathbf {x} _{3}={\begin{pmatrix}x_{31}\\x_{32}\\x_{33}\\x_{34}\end{pmatrix}}.} それから
( A − 5 I ) 3 x 3 = ( 0 0 0 14 0 0 0 − 4 0 0 0 3 0 0 0 − 1 ) ( x 31 x 32 x 33 x 34 ) = ( 14 x 34 − 4 x 34 3 x 34 − x 34 ) = ( 0 0 0 0 ) {\displaystyle (A-5I)^{3}\mathbf {x} _{3}={\begin{pmatrix}0&0&0&14\\0&0&0&-4\\0&0&0&3\\0&0&0&-1\end{pmatrix}}{\begin{pmatrix}x_{31}\\x_{32}\\x_{33}\\x_{34}\end{pmatrix}}={\begin{pmatrix}14x_{34}\\-4x_{34}\\3x_{34}\\-x_{34}\end{pmatrix}}={\begin{pmatrix}0\\0\\0\\0\end{pmatrix}}} そして
( A − 5 I ) 2 x 3 = ( 0 0 2 − 8 0 0 0 4 0 0 0 − 3 0 0 0 1 ) ( x 31 x 32 x 33 x 34 ) = ( 2 x 33 − 8 x 34 4 x 34 − 3 x 34 x 34 ) ≠ ( 0 0 0 0 ) . {\displaystyle (A-5I)^{2}\mathbf {x} _{3}={\begin{pmatrix}0&0&2&-8\\0&0&0&4\\0&0&0&-3\\0&0&0&1\end{pmatrix}}{\begin{pmatrix}x_{31}\\x_{32}\\x_{33}\\x_{34}\end{pmatrix}}={\begin{pmatrix}2x_{33}-8x_{34}\\4x_{34}\\-3x_{34}\\x_{34}\end{pmatrix}}\neq {\begin{pmatrix}0\\0\\0\\0\end{pmatrix}}.} したがって、条件( 3 )と( 4 )を満たすためには、と が必要となる。 と には制約はない 。を選択すると 、 x 34 = 0 {\displaystyle x_{34}=0} x 33 ≠ 0 {\displaystyle x_{33}\neq 0} x 31 {\displaystyle x_{31}} x 32 {\displaystyle x_{32}} x 31 = x 32 = x 34 = 0 , x 33 = 1 {\displaystyle x_{31}=x_{32}=x_{34}=0,x_{33}=1}
x 3 = ( 0 0 1 0 ) {\displaystyle \mathbf {x} _{3}={\begin{pmatrix}0\\0\\1\\0\end{pmatrix}}} に対応する階数3の一般化固有ベクトルとして定義される。ただし 、、、 および を 異なる値に選ぶことで、階数3の一般化固有ベクトルを無限に得ることが可能であることに注意されたい 。しかしながら、我々の最初の選択は最も単純なものである。 [47] λ 1 = 5 {\displaystyle \lambda _{1}=5} x 31 {\displaystyle x_{31}} x 32 {\displaystyle x_{32}} x 33 {\displaystyle x_{33}} x 33 ≠ 0 {\displaystyle x_{33}\neq 0}
ここで、式( 1 )を用いて、階数2と階数1の一般化固有ベクトルとして、およびを得る。ここ で 、 x 2 {\displaystyle \mathbf {x} _{2}} x 1 {\displaystyle \mathbf {x} _{1}}
x 2 = ( A − 5 I ) x 3 = ( − 2 2 0 0 ) , {\displaystyle \mathbf {x} _{2}=(A-5I)\mathbf {x} _{3}={\begin{pmatrix}-2\\2\\0\\0\end{pmatrix}},} そして
x 1 = ( A − 5 I ) x 2 = ( 2 0 0 0 ) . {\displaystyle \mathbf {x} _{1}=(A-5I)\mathbf {x} _{2}={\begin{pmatrix}2\\0\\0\\0\end{pmatrix}}.} 単純 固有値は 標準的な手法 で扱うことができ 、通常の固有ベクトルを持つ。 λ 2 = 4 {\displaystyle \lambda _{2}=4}
y 1 = ( − 14 4 − 3 1 ) . {\displaystyle \mathbf {y} _{1}={\begin{pmatrix}-14\\4\\-3\\1\end{pmatrix}}.} の標準的な根拠 は A {\displaystyle A}
{ x 3 , x 2 , x 1 , y 1 } = { ( 0 0 1 0 ) ( − 2 2 0 0 ) ( 2 0 0 0 ) ( − 14 4 − 3 1 ) } . {\displaystyle \left\{\mathbf {x} _{3},\mathbf {x} _{2},\mathbf {x} _{1},\mathbf {y} _{1}\right\}=\left\{{\begin{pmatrix}0\\0\\1\\0\end{pmatrix}}{\begin{pmatrix}-2\\2\\0\\0\end{pmatrix}}{\begin{pmatrix}2\\0\\0\\0\end{pmatrix}}{\begin{pmatrix}-14\\4\\-3\\1\end{pmatrix}}\right\}.} x 1 , x 2 {\displaystyle \mathbf {x} _{1},\mathbf {x} _{2}} および はに関連付けられた一般化固有ベクトルであり 、 は に関連付けられた通常の固有ベクトルです 。 x 3 {\displaystyle \mathbf {x} _{3}} λ 1 {\displaystyle \lambda _{1}} y 1 {\displaystyle \mathbf {y} _{1}} λ 2 {\displaystyle \lambda _{2}}
これはかなり単純な例です。一般に、階数を持つ線形独立な一般化固有ベクトルの 数は 必ずしも等しくありません。つまり、特定の固有値に対応する異なる長さの連鎖が複数存在する可能性があります。 [48] ρ k {\displaystyle \rho _{k}} k {\displaystyle k}
一般化モーダル行列 を n × n 行列とする 。 一般化様相行列 は、 n × n 行列であり 、その列をベクトルとみなすと、その標準基底を形成し、 以下の規則に従って 現れる。 A {\displaystyle A} M {\displaystyle M} A {\displaystyle A} A {\displaystyle A} M {\displaystyle M}
1 つのベクトル (つまり、長さが 1 つのベクトル) で構成されるすべてのジョルダン連鎖は、 の最初の列に表示されます 。 M {\displaystyle M} 1 つのチェーンのすべてのベクトルは、 の隣接する列に一緒に表示されます 。 M {\displaystyle M} 各連鎖は昇順で現れる (つまり、ランク1の一般化固有ベクトルは同じ連鎖のランク2の一般化固有ベクトルの前に現れ、ランク2の一般化固有ベクトルは同じ連鎖のランク3の一般化固有ベクトルの前に現れる、など)。 [49] M {\displaystyle M}
[ ⌜ λ 1 1 λ 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌞ λ 1 1 λ 1 λ 1 ⌟ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 n ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌟ ⌞ λ 2 λ 2 ⌟ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 λ n ⌝ ⌜ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌜ λ n 1 n ⌝ ⌞ λ 1 1 λ 1 1 λ 1 ⌝ ⌜ λ 2 1 λ 2 ⌝ [ λ 3 ] ⋱ ⌞ λ n λ n ⌟ ] {\displaystyle {\begin{bmatrix}{\color {red}\ulcorner }\lambda _{1}1{\hphantom {\lambda _{1}\lambda _{1}}}{\color {red}\urcorner }{\hphantom {\ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1}}\lambda _{1}1{\hphantom {\lambda _{1}\urcorner \ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\color {red}\llcorner }{\hphantom {\lambda _{1}1\lambda _{1}}}\lambda _{1}{\color {red}\lrcorner }{\hphantom {\ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\urcorner }}{\color {red}\ulcorner }\lambda _{2}1{\hphantom {n}}{\color {red}\urcorner }{\hphantom {[\lambda _{3}]\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\lrcorner }}{\color {red}\llcorner }{\hphantom {\lambda _{2}}}\lambda _{2}{\color {red}\lrcorner }{\hphantom {[\lambda _{3}]\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\urcorner \ulcorner \lambda _{2}1\lambda _{2}\urcorner }}{\color {red}[}\lambda _{3}{\color {red}]}{\hphantom {\ddots \ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\urcorner \ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]}}\ddots {\hphantom {\ulcorner \lambda _{n}1\lambda _{n}\urcorner }}\\{\hphantom {\ulcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\urcorner \ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]\ddots }}{\color {red}\ulcorner }\lambda _{n}1{\hphantom {n}}{\color {red}\urcorner }\\{\hphantom {\llcorner \lambda _{1}1\lambda _{1}1\lambda _{1}\urcorner \ulcorner \lambda _{2}1\lambda _{2}\urcorner [\lambda _{3}]\ddots }}{\color {red}\llcorner }{\hphantom {\lambda _{n}}}\lambda _{n}{\color {red}\lrcorner }\end{bmatrix}}}
ジョルダン正規形の行列の例。 赤いブロックはジョルダンブロックと呼ばれます。 n次元ベクトル空間を L ( V ) の 線型写像、つまり から 自身へのすべての線型写像の集合をL(V ) の線型写像とする 。また、 をある順序基底に関する 行列表現とする。 の 特性多項式 が 線型因数に因数分解されるとき、 は次のよう な形になる ことが示される。 V {\displaystyle V} ϕ {\displaystyle \phi } V {\displaystyle V} A {\displaystyle A} ϕ {\displaystyle \phi } f ( λ ) {\displaystyle f(\lambda )} A {\displaystyle A} f ( λ ) {\displaystyle f(\lambda )}
f ( λ ) = ± ( λ − λ 1 ) μ 1 ( λ − λ 2 ) μ 2 ⋯ ( λ − λ r ) μ r , {\displaystyle f(\lambda )=\pm (\lambda -\lambda _{1})^{\mu _{1}}(\lambda -\lambda _{2})^{\mu _{2}}\cdots (\lambda -\lambda _{r})^{\mu _{r}},} ここで、 は の異なる固有値であり 、各 は 対応する固有値の代数的重複度であり 、 は ジョルダン標準形 の 行列に類似しており 、各 は対角線上に連続し て現れ 、各 のすぐ上の要素 (つまり、上 対角線 上)は 0 か 1 です。各ブロックでは、各 の最初の出現の上の要素は 常に 0 です(最初のブロックを除く)。上対角線上のその他の要素はすべて 1 です。その他の要素(対角線と上対角線外)はすべて 0 です。(ただし、固有値間、または特定の固有値のブロック間では順序付けは行われません。)この行列は、 の対角化にできるだけ近くなります 。 が対角化可能な場合、対角線より上のすべての要素は 0 です。 [50] 一部の教科書では、 を上 対角線 上ではなく、下対角線上、つまり主対角線のすぐ下に挙げていることに注意してください。固有値は依然として主対角線上にあります。 [51] [52] λ 1 , λ 2 , … , λ r {\displaystyle \lambda _{1},\lambda _{2},\ldots ,\lambda _{r}} A {\displaystyle A} μ i {\displaystyle \mu _{i}} λ i {\displaystyle \lambda _{i}} A {\displaystyle A} J {\displaystyle J} λ i {\displaystyle \lambda _{i}} μ i {\displaystyle \mu _{i}} λ i {\displaystyle \lambda _{i}} λ i {\displaystyle \lambda _{i}} J {\displaystyle J} A {\displaystyle A} A {\displaystyle A}
すべての n × n 行列は 、相似変換によって得られるジョルダン正規形の 行列に相似であり 、ここで は の一般化様相行列である 。 [53] (上記の注を参照。) A {\displaystyle A} J {\displaystyle J} J = M − 1 A M {\displaystyle J=M^{-1}AM} M {\displaystyle M} A {\displaystyle A}
例4 ジョルダン正規形の行列で、次の行列に類似するものを見つけなさい。
A = ( 0 4 2 − 3 8 3 4 − 8 − 2 ) . {\displaystyle A={\begin{pmatrix}0&4&2\\-3&8&3\\4&-8&-2\end{pmatrix}}.} 解答: の特性方程式 は なので、 は 代数的重複度3の固有値である。前節の手順に従うと、 A {\displaystyle A} ( λ − 2 ) 3 = 0 {\displaystyle (\lambda -2)^{3}=0} λ = 2 {\displaystyle \lambda =2}
rank ( A − 2 I ) = 1 {\displaystyle \operatorname {rank} (A-2I)=1} そして
rank ( A − 2 I ) 2 = 0 = n − μ . {\displaystyle \operatorname {rank} (A-2I)^{2}=0=n-\mu .} したがって、 および は、の標準基底が 階数2の線形独立な一般化固有ベクトル1つと階数1の線形独立な一般化固有ベクトル2つ、あるいはそれと同等の、2つのベクトルの鎖1つと 1つのベクトルの鎖1つを含む ことを意味する。をとすれば 、 ρ 2 = 1 {\displaystyle \rho _{2}=1} ρ 1 = 2 {\displaystyle \rho _{1}=2} A {\displaystyle A} { x 2 , x 1 } {\displaystyle \left\{\mathbf {x} _{2},\mathbf {x} _{1}\right\}} { y 1 } {\displaystyle \left\{\mathbf {y} _{1}\right\}} M = ( y 1 x 1 x 2 ) {\displaystyle M={\begin{pmatrix}\mathbf {y} _{1}&\mathbf {x} _{1}&\mathbf {x} _{2}\end{pmatrix}}}
M = ( 2 2 0 1 3 0 0 − 4 1 ) , {\displaystyle M={\begin{pmatrix}2&2&0\\1&3&0\\0&-4&1\end{pmatrix}},} そして
J = ( 2 0 0 0 2 1 0 0 2 ) , {\displaystyle J={\begin{pmatrix}2&0&0\\0&2&1\\0&0&2\end{pmatrix}},} ここで 、 は の一般化モーダル行列であり 、 の列は の標準基底であり 、 である 。 [54] 一般化固有ベクトル自体は一意ではなく、と の列の一部は 交換可能であるため、 と は 一意ではないことに注意すること。 [55] M {\displaystyle M} A {\displaystyle A} M {\displaystyle M} A {\displaystyle A} A M = M J {\displaystyle AM=MJ} M {\displaystyle M} J {\displaystyle J} M {\displaystyle M} J {\displaystyle J}
例5 例3では、行列 の線形独立な一般化固有ベクトルの標準基底を求めました 。 の一般化様相行列 は A {\displaystyle A} A {\displaystyle A}
M = ( y 1 x 1 x 2 x 3 ) = ( − 14 2 − 2 0 4 0 2 0 − 3 0 0 1 1 0 0 0 ) . {\displaystyle M={\begin{pmatrix}\mathbf {y} _{1}&\mathbf {x} _{1}&\mathbf {x} _{2}&\mathbf {x} _{3}\end{pmatrix}}={\begin{pmatrix}-14&2&-2&0\\4&0&2&0\\-3&0&0&1\\1&0&0&0\end{pmatrix}}.} ジョルダン正規形の行列は 、 A {\displaystyle A}
J = ( 4 0 0 0 0 5 1 0 0 0 5 1 0 0 0 5 ) , {\displaystyle J={\begin{pmatrix}4&0&0&0\\0&5&1&0\\0&0&5&1\\0&0&0&5\end{pmatrix}},} となることによって 。 A M = M J {\displaystyle AM=MJ}
アプリケーション
行列関数 正方行列 に対して実行できる最も基本的な演算のうち3つは、 行列の加算、スカラーによる乗算、行列の乗算である。 [56]これらはまさに n × n 行列 の 多項式 関数 を定義するために必要な演算である 。 [57]基本的な 微積分学で多くの関数が マクローリン級数 として表せること を思い出せば 、行列のより一般的な関数をかなり簡単に定義できる。 [58] が対角化可能である 場合 、 A {\displaystyle A} A {\displaystyle A}
D = M − 1 A M , {\displaystyle D=M^{-1}AM,} と
D = ( λ 1 0 ⋯ 0 0 λ 2 ⋯ 0 ⋮ ⋮ ⋱ ⋮ 0 0 ⋯ λ n ) , {\displaystyle D={\begin{pmatrix}\lambda _{1}&0&\cdots &0\\0&\lambda _{2}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &\lambda _{n}\end{pmatrix}},} それから
D k = ( λ 1 k 0 ⋯ 0 0 λ 2 k ⋯ 0 ⋮ ⋮ ⋱ ⋮ 0 0 ⋯ λ n k ) {\displaystyle D^{k}={\begin{pmatrix}\lambda _{1}^{k}&0&\cdots &0\\0&\lambda _{2}^{k}&\cdots &0\\\vdots &\vdots &\ddots &\vdots \\0&0&\cdots &\lambda _{n}^{k}\end{pmatrix}}} の関数に対するマクローリン級数の計算は 大幅に簡素化される。 [59] 例えば、 の任意のべき乗 k を得るには、 を計算し、 を前置乗算し 、その結果に を後置乗算する だけでよい 。 [60] A {\displaystyle A} A {\displaystyle A} D k {\displaystyle D^{k}} D k {\displaystyle D^{k}} M {\displaystyle M} M − 1 {\displaystyle M^{-1}}
一般化固有ベクトルを用いると、のジョルダン正規形を得ることができ 、これらの結果は非対角化行列の関数を計算する簡単な方法に一般化できる。 [61] ( 行列関数#ジョルダン分解を 参照。) A {\displaystyle A}
微分方程式 線形常微分方程式系を解く問題を考えてみましょう
x ′ = A x , {\displaystyle \mathbf {x} '=A\mathbf {x} ,} 5
どこ
x = ( x 1 ( t ) x 2 ( t ) ⋮ x n ( t ) ) , x ′ = ( x 1 ′ ( t ) x 2 ′ ( t ) ⋮ x n ′ ( t ) ) , {\displaystyle \mathbf {x} ={\begin{pmatrix}x_{1}(t)\\x_{2}(t)\\\vdots \\x_{n}(t)\end{pmatrix}},\quad \mathbf {x} '={\begin{pmatrix}x_{1}'(t)\\x_{2}'(t)\\\vdots \\x_{n}'(t)\end{pmatrix}},} そして A = ( a i j ) . {\displaystyle A=(a_{ij}).} 行列が 対角行列で、 に対してとなる場合 、システム( 5 )は、次の形をとる n 個の方程式のシステムになる。 A {\displaystyle A} a i j = 0 {\displaystyle a_{ij}=0} i ≠ j {\displaystyle i\neq j}
この場合、一般解は次のように与えられる。
x 1 = k 1 e a 11 t {\displaystyle x_{1}=k_{1}e^{a_{11}t}} x 2 = k 2 e a 22 t {\displaystyle x_{2}=k_{2}e^{a_{22}t}} ⋮ {\displaystyle \vdots } x n = k n e a n n t . {\displaystyle x_{n}=k_{n}e^{a_{nn}t}.} 一般の場合、 系( 5 )を対角化して、以下のように系( 6 )のような系に縮約する。 が対角化可能な場合、 となる。 ここで は のモード行列である 。 を代入すると 、 式( 5 )は となる。 A {\displaystyle A} A {\displaystyle A} D = M − 1 A M {\displaystyle D=M^{-1}AM} M {\displaystyle M} A {\displaystyle A} A = M D M − 1 {\displaystyle A=MDM^{-1}} M − 1 x ′ = D ( M − 1 x ) {\displaystyle M^{-1}\mathbf {x} '=D(M^{-1}\mathbf {x} )}
y ′ = D y , {\displaystyle \mathbf {y} '=D\mathbf {y} ,} 7
どこ
x = M y . {\displaystyle \mathbf {x} =M\mathbf {y} .} 8
( 7 ) の解は
y 1 = k 1 e λ 1 t {\displaystyle y_{1}=k_{1}e^{\lambda _{1}t}} y 2 = k 2 e λ 2 t {\displaystyle y_{2}=k_{2}e^{\lambda _{2}t}} ⋮ {\displaystyle \vdots } y n = k n e λ n t . {\displaystyle y_{n}=k_{n}e^{\lambda _{n}t}.} ( 5 )の解は 関係式( 8 )を用いて得られる。 [62] x {\displaystyle \mathbf {x} }
一方、 が対角化できない場合、 の一般化様相行列 をとします 。この行列 は のジョルダン正規形となります 。この系は 次の式を持ちます
。 A {\displaystyle A} M {\displaystyle M} A {\displaystyle A} J = M − 1 A M {\displaystyle J=M^{-1}AM} A {\displaystyle A} y ′ = J y {\displaystyle \mathbf {y} '=J\mathbf {y} }
y 1 ′ = λ 1 y 1 + ϵ 1 y 2 ⋮ y n − 1 ′ = λ n − 1 y n − 1 + ϵ n − 1 y n y n ′ = λ n y n , {\displaystyle {\begin{aligned}y_{1}'&=\lambda _{1}y_{1}+\epsilon _{1}y_{2}\\&\vdots \\y_{n-1}'&=\lambda _{n-1}y_{n-1}+\epsilon _{n-1}y_{n}\\y_{n}'&=\lambda _{n}y_{n},\end{aligned}}}
9
ここで、 は の主対角線上の固有値であり 、 は の 上対角線上の 1 と 0 である 。 系 ( 9 ) は ( 5 )より簡単に解けることが多い。 ( 9 ) の最後の方程式を について解くと 、 が得られる。次に、 についてのこの解を( 9 )の最後から2番目の方程式に 代入し 、 について解きます。 この手順を続けて、( 9 )を最後の方程式から最初の方程式まで解き 、 について系全体を解きます。 すると、関係式( 8 )を用いて解が得られます 。 [63] λ i {\displaystyle \lambda _{i}} J {\displaystyle J} ϵ i {\displaystyle \epsilon _{i}} J {\displaystyle J} y n {\displaystyle y_{n}} y n = k n e λ n t {\displaystyle y_{n}=k_{n}e^{\lambda _{n}t}} y n {\displaystyle y_{n}} y n − 1 {\displaystyle y_{n-1}} y {\displaystyle \mathbf {y} } x {\displaystyle \mathbf {x} }
補題:
次のような長さの一般化固有ベクトルの連鎖が与えられる。 r , {\displaystyle r,}
X 1 = v 1 e λ t {\displaystyle X_{1}=v_{1}e^{\lambda t}} X 2 = ( t v 1 + v 2 ) e λ t {\displaystyle X_{2}=(tv_{1}+v_{2})e^{\lambda t}} X 3 = ( t 2 2 v 1 + t v 2 + v 3 ) e λ t {\displaystyle X_{3}=\left({\frac {t^{2}}{2}}v_{1}+tv_{2}+v_{3}\right)e^{\lambda t}} ⋮ {\displaystyle \vdots } X r = ( t r − 1 ( r − 1 ) ! v 1 + . . . + t 2 2 v r − 2 + t v r − 1 + v r ) e λ t {\displaystyle X_{r}=\left({\frac {t^{r-1}}{(r-1)!}}v_{1}+...+{\frac {t^{2}}{2}}v_{r-2}+tv_{r-1}+v_{r}\right)e^{\lambda t}} 、 これらの関数は連立方程式を解きます。
X ′ = A X . {\displaystyle X'=AX.} 証拠:
定義する
v 0 = 0 {\displaystyle v_{0}=0} X j ( t ) = e λ t ∑ i = 1 j t j − i ( j − i ) ! v i . {\displaystyle X_{j}(t)=e^{\lambda t}\sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}v_{i}.} そして、および として 、 t 0 = 1 {\displaystyle {t^{0}}=1} 1 ′ = 0 {\displaystyle 1'=0}
X j ′ ( t ) = e λ t ∑ i = 1 j − 1 t j − i − 1 ( j − i − 1 ) ! v i + e λ t λ ∑ i = 1 j t j − i ( j − i ) ! v i {\displaystyle X'_{j}(t)=e^{\lambda t}\sum _{i=1}^{j-1}{\frac {t^{j-i-1}}{(j-i-1)!}}v_{i}+e^{\lambda t}\lambda \sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}v_{i}} 。 一方、我々は 、 v 0 = 0 {\displaystyle v_{0}=0}
A X j ( t ) = e λ t ∑ i = 1 j t j − i ( j − i ) ! A v i {\displaystyle AX_{j}(t)=e^{\lambda t}\sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}Av_{i}} = e λ t ∑ i = 1 j t j − i ( j − i ) ! ( v i − 1 + λ v i ) {\displaystyle =e^{\lambda t}\sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}(v_{i-1}+\lambda v_{i})} = e λ t ∑ i = 2 j t j − i ( j − i ) ! v i − 1 + e λ t λ ∑ i = 1 j t j − i ( j − i ) ! v i {\displaystyle =e^{\lambda t}\sum _{i=2}^{j}{\frac {t^{j-i}}{(j-i)!}}v_{i-1}+e^{\lambda t}\lambda \sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}v_{i}} = e λ t ∑ i = 1 j − 1 t j − i − 1 ( j − i − 1 ) ! v i + e λ t λ ∑ i = 1 j t j − i ( j − i ) ! v i {\displaystyle =e^{\lambda t}\sum _{i=1}^{j-1}{\frac {t^{j-i-1}}{(j-i-1)!}}v_{i}+e^{\lambda t}\lambda \sum _{i=1}^{j}{\frac {t^{j-i}}{(j-i)!}}v_{i}} = X j ′ ( t ) {\displaystyle =X'_{j}(t)} 必要に応じて。
注記 ^ ブロンソン(1970年、189ページ) ^ Beauregard & Fraleigh (1973、p. 310) ^ ネリング(1970年、118ページ) ^ ゴラブとヴァン・ローン (1996, p. 316) ^ Beauregard & Fraleigh (1973、p. 319) ^ ブロンソン(1970年、194~195ページ) ^ ゴラブとヴァン・ローン (1996, p. 311) ^ ブロンソン(1970年、196ページ) ^ ブロンソン(1970年、189ページ) ^ Beauregard & Fraleigh (1973、pp. 316–318) ^ ネリング(1970年、118ページ) ^ ブロンソン(1970年、196ページ) ^ アントン(1987年、301~302ページ) ^ Beauregard & Fraleigh (1973、p. 266) ^ Burden & Faires (1993, p. 401) ^ ゴラブとヴァン・ローン (1996、pp. 310–311) ^ ハーパー(1976年、58ページ) ^ ハーシュタイン(1964年、225ページ) ^ クライスツィヒ (1972, pp. 273, 684) ^ ネリング(1970年、104ページ) ^ Burden & Faires (1993, p. 401) ^ Beauregard & Fraleigh (1973、pp. 270–274) ^ ブロンソン(1970年、179~183ページ) ^ ブロンソン(1970年、181ページ) ^ ブロンソン(1970年、179ページ) ^ Beauregard & Fraleigh (1973、pp. 270–274) ^ ブロンソン(1970年、179~183ページ) ^ ブロンソン(1970年、189ページ) ^ ブロンソン(1970年、190、202ページ) ^ ブロンソン(1970年、189、203ページ) ^ ブロンソン(1970年、206~207ページ) ^ ブロンソン(1970年、205ページ) ^ ブロンソン(1970年、196ページ) ^ ブロンソン(1970年、189ページ、209~215ページ) ^ ゴラブとヴァン・ローン (1996, p. 316) ^ ハーシュタイン(1964年、259ページ) ^ ネリング(1970年、118ページ) ^ ネリング(1970年、118ページ) ^ ネリング(1970年、118ページ) ^ ハーシュタイン(1964年、261ページ) ^ Beauregard & Fraleigh (1973、p. 310) ^ ネリング(1970年、122、123ページ) ^ ブロンソン(1970年、189~209ページ) ^ ブロンソン(1970年、194~195ページ) ^ ブロンソン(1970年、196、197ページ) ^ ブロンソン(1970年、197、198ページ) ^ ブロンソン(1970年、190~191ページ) ^ ブロンソン(1970年、197~198ページ) ^ ブロンソン(1970年、205ページ) ^ Beauregard & Fraleigh (1973、p. 311) ^ カレン(1966年、114ページ) ^ フランクリン(1968年、122ページ) ^ ブロンソン(1970年、207ページ) ^ ブロンソン(1970年、208ページ) ^ ブロンソン(1970年、206ページ) ^ Beauregard & Fraleigh (1973、pp. 57–61) ^ ブロンソン(1970年、104ページ) ^ ブロンソン(1970年、105ページ) ^ ブロンソン(1970年、184ページ) ^ ブロンソン(1970年、185ページ) ^ ブロンソン(1970年、209~218ページ) ^ Beauregard & Fraleigh (1973、pp. 274–275) ^ Beauregard & Fraleigh (1973、p. 317)
参考文献 アントン・ハワード(1987年)、 初等線形代数 (第5版)、ニューヨーク: ワイリー 、 ISBN 0-471-84819-0 アクラー、シェルドン(1997年) 『線形代数を正しく理解する』 (第2版)シュプリンガー社、 ISBN 978-0-387-98258-8 。 Beauregard, Raymond A.; Fraleigh, John B. (1973), A First Course In Linear Algebra: with Optional Introduction to Groups, Rings, and Fields , Boston: Houghton Mifflin Co. , ISBN 0-395-14017-X ブロンソン、リチャード(1970)、 マトリックス法入門 、ニューヨーク: アカデミックプレス 、 LCCN 70097490 バーデン、リチャード・L.; フェアーズ、J.ダグラス (1993)、『数値解析』(第5版)、ボストン:プリンドル、ウェーバー、シュミット、 ISBN 0-534-93219-3 Cullen, Charles G. (1966), Matrices and Linear Transformations 、Reading: Addison-Wesley 、 LCCN 66021267 フランクリン、ジョエル N. (1968)、 「マトリックス理論」 、エングルウッドクリフス: プレンティスホール 、 LCCN 68016345 ゴルブ、ジーン・H.; ヴァン・ローン、チャールズ・F. (1996)、 『行列計算』 (第3版)、ボルチモア: ジョンズ・ホプキンス大学出版局 、 ISBN 0-8018-5414-8 ハーパー、チャーリー(1976)、 数理物理学入門 、ニュージャージー: プレンティス・ホール 、 ISBN 0-13-487538-9 Herstein, IN (1964), Topics In Algebra , Waltham: Blaisdell Publishing Company, ISBN 978-1114541016 クレイジグ、エルウィン(1972)、Advanced Engineering Mathematics(第3版)、ニューヨーク: Wiley 、 ISBN 0-471-50728-8 Nering, Evar D. (1970), 線形代数と行列理論 (第2版)、ニューヨーク: Wiley 、 LCCN 76091646