代数 > 線形代数 > 行列 > 行列 固有値 >
代数 > 線形代数 > 行列 > 行列 分解 >

EXPLORE THIS TOPIC IN MathWorld Classroom

Eigenvectors are special set of vectors associated with a linear system of equations (i.)…e., 782>

系の固有ベクトルおよび固有値を決定することは、物理学および工学において非常に重要であり、それは行列対角化と等価であり、安定性解析、回転体の物理、振動系の小振動など、いくつかの例を挙げるとよくある応用で発生するものである。 各固有ベクトルは、対応する固有値と呼ばれるものと対になっている。 数学的には、固有ベクトルを左固有ベクトルと右固有ベクトルの2種類に区別する必要がある。 しかし、物理学や工学の多くの問題では、右固有ベクトルのみを考慮すれば十分です。 782>

正方行列Aを固有値と固有ベクトルに分解することをこの作業では固有分解と呼び、Aの固有ベクトルからなる行列が正方である限り常にこの分解が可能であることを固有分解定理と呼んでいる。

右固有ベクトルを

 AX_R=lambda_RX_R,
(1)

ここにAは行列と定義される。 ということは

 (A-lambda_RI)X_R=0,
(2)

つまり右の固有値はゼロ行列式を持たねばならない、即ち。e.,

 det(A-lambda_RI)=0.
(3)

同様に、左固有ベクトルを

X_LA=lambda_LX_Lを満たす行ベクトル <img src= と定義します。 “>
(4)

各辺の転置を取ると

 (X_LA)^(T)=λ_LX_L^(T) が得られます。
(5)

これは

 A^(T)X_L^(T)=lambda_LX_L^(T) に書き換え可能である。
(6)

もう一度整理すると、

 (A^(T)-λ_LI)X_L^(T)=0 が得られる。
(7)

つまり

 det(A^(T)-lambda_LI)=0 ですね。
(8)

書き換えると

<8254>

0=det(A^(T)-λ_LI) =
(9)
= det(A-)lambda_LI)^(T)
(10)
= det(A-lambda_LI),
(11)

ここで最後のステップは恒等式

 det(A)=det(A^(T)) から導かれます。
(12)

任意のAXに対して共に0になる式(◇)と(11)を等化するには、したがってラムダ_R=ラムダ_L、即ち。

ここで、X_Rを右固有ベクトルの列で形成される行列とし、X_Lを左固有ベクトルの行で形成される行列とすると、左と右固有値は同等である。

 D=とする。
(13)

次に

とすると

AX_R = x_rd
(14)
x_la = dx_l
(15)

and

X_LAX_R = x_lx_rd
(16)
x_lax_r = dx_lx_r.Lx_R,
(17)

 X_LX_RD=DX_LX_Rです。
(18)

しかしこの式は

という形になっているのです。

 CD=DC
(19)

ここで、Dは対角行列である。 ということは、C=X_LX_Rも対角であることが真である。 特にAが対称行列であれば、左右の固有ベクトルは単に互いの転置であり、Aが自己共役行列(すなわちエルミート)であれば、左右の固有ベクトルは共役行列となる。

固有ベクトルはゼロベクトルと等しくないことがある。 固有ベクトルの0でないスカラー倍は元の固有ベクトルと等価である。 したがって、一般性を損なうことなく、固有ベクトルはしばしば単位長に正規化されます。

n×n行列は常に n 固有値を持ち、その一部または全部は縮退しているかもしれませんが、そのような行列は 0 から n 線形独立固有ベクトルを持つかもしれません。 例えば、行列は単一の固有ベクトル(1,0)のみを持ちます。

固有ベクトルはWolfram LanguageでEigenvectorsを使って計算することができます。 このコマンドは常に長さnのリストを返すので,線形独立でない固有ベクトルはすべて0ベクトルとして返されます. 固有ベクトルと固有値はEigensystemというコマンドで一緒に返すことができます.

固有ベクトルx_1, x_2, x_3 と対応する固有値 lambda_1 を持つ 3×3 行列 A が与えられたとき、次のようになる。 lambda_2, lambda_3 とすると、任意のベクトル y

 y=b_1x_1+b_2x_2+b_3x_3 と書くことができる。
(20)

行列Aを適用する。

)

Ay = b_1Ax_1+b_2Ax_2+b_3Ax_3
(21)
= lambda_1(b_1x_1+(lambda_2)/(lambda_1)b_2x_2+(lambda_3)/(lambda_1)b_3x_3),
(22)

so

 A^ny=lambda_1^n となる。
(23)

もしlambda_1lambda_2,lambda_3、そしてb_1!=0であるから、

 lim_(n-infty)A^ny=lambda_1^nb_1x_1 が成り立つ。
(24)

従って、任意のベクトルに行列を繰り返し驚くほど適用すると、最大固有値を持つ固有ベクトルに比例したベクトルになります。