markdown
固有値と固有ベクトルmd 759c783
lecture/math/linear-algebra/eigenvalues-and-eigenvectors.lecture.n.md
Download PDF

固有値こゆうちeigenvalueと固有ベクトルeigenvector

Eigenvalues固有値こゆうち and eigenvectors固有こゆうベクトル

1導入どうにゅう

この講義こうぎでは、固有値こゆうちeigenvalue・固有ベクトルeigenvectorとは「線型変換せんけいへんかんlinear transformationの中なかで向むきが変かわらない方向ほうこう(固有ベクトルeigenvector)とその伸縮率しんしゅくりつ(固有値こゆうちeigenvalue)」であり、変換へんかんの本質的ほんしつてきな構造こうぞうstructureをもっとも単純たんじゅんな基底きていbasisで表あらわすという観点かんてんである。

「特性方程式とくせいほうていしきを解とく技術ぎじゅつ」として扱あつかうと、「なぜ det(A-λI)=0 なのか」「代数的重複度だいすうてきちょうふくどと幾何的重複度きかてきちょうふくどの違ちがいとは何なにか」が不明瞭ふめいりょうになる。変換へんかんが「向むきを変かえない方向ほうこう」を持もつ条件じょうけんから出発しゅっぱつすると、すべてが自然しぜんに導みちびかれる。

1Introduction

The central idea of this lecture is that an eigenvalue固有値こゆうち and an eigenvector固有こゆうベクトル describe a direction that a linear transformation線型変換せんけいへんかん does not turn, together with the scale factor on that direction. This is the starting point for expressing the essential structure構造こうぞう of a transformation in the simplest possible basis基底きてい.

If eigenvalues are treated only as roots of a characteristic equation, the reasons behind det(A-λI)=0 and the distinction between algebraic and geometric multiplicity become unclear. Starting from the condition that a transformation preserves a direction makes the determinant condition, multiplicities, diagonalization, and applications follow naturally.

2用語ようごと定義ていぎ

2.1固有値こゆうちEigenvalueと固有こゆうベクトルEigenvector

A を n 次じ正方行列せいほうぎょうれつsquare matrixとするとき、v≠0 かつ

Av=λv

を満みたすスカラー λ を固有値こゆうちEigenvalue、v を対応たいおうする固有こゆうベクトルEigenvectorという。

2.2特性多項式とくせいたこうしきCharacteristic polynomial

p(λ)=det(A-λI)

を A の特性多項式とくせいたこうしきCharacteristic polynomialという。n 次じ多項式たこうしきであり、固有値こゆうちeigenvalueは p(λ)=0 の解かいである。

2Terms and definitions

2.1Eigenvalue固有値こゆうち and eigenvector固有こゆうベクトル

Let A be an n by n square matrix正方行列せいほうぎょうれつ. A scalar λ is an eigenvalue固有値こゆうち of A if there exists a nonzero vector v such that

Av=λv.

The nonzero vector v is an eigenvector固有こゆうベクトル corresponding to λ.

2.2Characteristic polynomial特性多項式とくせいたこうしき

The polynomial

p(λ)=det(A-λI)

is called the characteristic polynomial特性多項式とくせいたこうしき of A. It has degree n, and the eigenvalues固有値こゆうち are precisely the roots of p(λ)=0 over the field being considered.

3方針ほうしん

  1. Av=λv を det(A-λI)=0 へ翻訳ほんやくする論理ろんりを確認かくにんする
  2. 代数的重複度だいすうてきちょうふくどと幾何的重複度きかてきちょうふくどの差さが対角化可能性たいかくかかのうせいを決きめる
  3. Cayley-Hamilton の定理ていりと微分方程式びぶんほうていしきへの接続せつぞくを確認かくにんする
data/lecture/math/linear-algebra/linear-maps-and-matrices.lecture.n.md

固有値問題こゆうちもんだいeigenvalue problemとしての具体的ぐたいてきな解法手順かいほうてじゅんは、次つぎの講義こうぎで独立どくりつに扱あつかう。

data/lecture/math/linear-algebra/eigenvalue-problem-basics.lecture.n.md

3Plan

  1. Translate Av=λv into det(A-λI)=0.
  2. Explain how algebraic multiplicity and geometric multiplicity control diagonalizability.
  3. Connect the Cayley-Hamilton theorem with powers of matrices and systems of differential equations.
data/lecture/math/linear-algebra/linear-maps-and-matrices.lecture.n.md

The detailed solving workflow for the eigenvalue problem固有値問題こゆうちもんだい is treated in the next lecture.

data/lecture/math/linear-algebra/eigenvalue-problem-basics.lecture.n.md

4直感的ちょっかんてきな説明せつめい

一般いっぱんの線型変換せんけいへんかんlinear transformationは、ベクトルの向むきを変かえる。たとえば平面へいめんを斜ななめに押おし潰つぶしたり、回転かいてんさせたりすると、ほとんどの矢印やじるしは別べつの方向ほうこうを向むく。

その中なかで、変換後へんかんごも同おなじ直線上ちょくせんじょうに残のこる方向ほうこうがあるなら、その方向ほうこうが固有ベクトルeigenvectorの方向ほうこうである。変換へんかんはその方向ほうこうでは、回転かいてんやせん断だんではなく、ただの伸縮しんしゅくとして見みえる。その伸縮率しんしゅくりつが固有値こゆうちeigenvalueである。

つまり固有値こゆうちeigenvalue・固有ベクトルeigenvectorは、複雑ふくざつな線型変換せんけいへんかんlinear transformationを「方向ほうこうごとの倍率ばいりつ」として読よめる場所ばしょを探さがす道具どうぐである。この見方みかたが対角化たいかくかdiagonalizationへ接続せつぞくする。

data/lecture/math/linear-algebra/diagonalization-basics.lecture.n.md

4Intuitive explanation

A general linear transformation線型変換せんけいへんかん changes the directions of most vectors. A shear, rotation, or non-uniform stretch usually sends an arrow to a different direction.

An eigenvector固有こゆうベクトル is a direction that remains on the same line after the transformation. On that line the transformation looks like pure stretching or shrinking, not rotation or shearing. The scale factor is the eigenvalue固有値こゆうち.

Thus eigenvalues and eigenvectors are tools for finding the places where a complicated transformation can be read as "one scale factor per direction." This viewpoint leads directly to diagonalization対角化たいかくか.

data/lecture/math/linear-algebra/diagonalization-basics.lecture.n.md

5厳密げんみつな説明せつめい

5.11. 方程式ほうていしきへの翻訳ほんやく

固有ベクトルeigenvectorの定義ていぎは、非零ひれいベクトル v≠0 を含ふくめて

Av=λv(v≠0)

である。これは

(A-λI)v=0(v≠0)

と同値どうちequivalentである。ここまでは固有値こゆうちeigenvalueの定義ていぎを移項いこうしただけであり、まだ行列式ぎょうれつしきdeterminantは登場とうじょうしていない。

次つぎに、同次方程式どうじほうていしき (A-λI)x=0 が非自明解ひじめいかいを持もつ条件じょうけんを用もちいる。A-λI が可逆かぎゃくinvertibleなら解かいは x=0 のみであるため、非自明解ひじめいかいが存在そんざいすることは

det(A-λI)=0

と同値どうちequivalentである。つまり、定義ていぎと特性方程式とくせいほうていしきは次つぎの 2 段階だんかいで接続せつぞくされる。

Av=λv(v≠0)⟺(A-λI)v=0(v≠0)
(A-λI)x=0が非自明解を持つ⟺det(A-λI)=0

5.22. 特性多項式とくせいたこうしきcharacteristic polynomialの性質せいしつ

n×n 行列ぎょうれつmatrix A の特性多項式とくせいたこうしきcharacteristic polynomialは n 次じであり:

det(A-λI)=(-1)nλn+(-1)n-1(trA)λn-1+⋯+detA

重要じゅうような関係かんけい:

trA=∑i=1nλi,detA=∏i=1nλi

(λ1,…,λn は重複ちょうふくを込こみにした固有値こゆうちeigenvalue全体ぜんたい)

ここでいう固有値こゆうちeigenvalue全体ぜんたいとは、特性多項式とくせいたこうしきcharacteristic polynomialが考かんがえている体たいの上うえで一次式いちじしきに分解ぶんかいする場合ばあいの根こんを、代数的重複度だいすうてきちょうふくどを込こめて数かぞえたものである。実行列じつぎょうれつでも実数じっすうreal numberの範囲はんいで固有値こゆうちeigenvalueが足たりないことがあるため、必要ひつように応おうじて複素数体ふくそすうたいまで拡張かくちょうして読よむ。

応用おうよう:detA=0⟺0 が固有値こゆうちeigenvalue(A が特異行列とくいぎょうれつ)。

5.33. 重複度ちょうふくどの二種類にしゅるい

代数的重複度だいすうてきちょうふくど(AM):固有値こゆうちeigenvalue λ0 が特性多項式とくせいたこうしきcharacteristic polynomialの根こんとして現あらわれる回数かいすう。

幾何的重複度きかてきちょうふくど(GM):固有空間こゆうくうかんeigenspace

Eλ0=ker(A-λ0I)

の次元じげんdimensionである。A が n×n 行列ぎょうれつmatrixなら、階数かいすうrank・退化次数たいかじすうnullity定理ていりにより

GM(λ0)=dimker(A-λ0I)=n-rank(A-λ0I)

である。

固有値こゆうちeigenvalue λ0 について、常つねに

1[PARSE ERROR: Undefined("Command(\"le\")")]GM(λ0)[PARSE ERROR: Undefined("Command(\"le\")")]AM(λ0)

である。また、特性多項式とくせいたこうしきcharacteristic polynomialが対象たいしょうの体たいの上うえで一次式いちじしきに分解ぶんかいしている場合ばあい、すべての固有値こゆうちeigenvalueについて幾何的重複度きかてきちょうふくどと代数的重複度だいすうてきちょうふくどが一致いっちすることは、対角化可能性たいかくかかのうせいと同値どうちequivalentである。

ここで Spec(A) は A の固有値こゆうちeigenvalue全体ぜんたいの集合しゅうごうである。A が対角化可能たいかくかかのうdiagonalizableであるとは、固有ベクトルeigenvectorからなる基底きていbasisが存在そんざいし、ある可逆行列かぎゃくぎょうれつinvertible matrix P と対角行列たいかくぎょうれつdiagonal matrix D により A=PDP-1 と書かけることである。この同値性どうちせいの証明しょうめいは、次つぎの対角化たいかくかdiagonalizationの講義こうぎで扱あつかう。

data/lecture/math/linear-algebra/diagonalization-basics.lecture.n.md
∀λ∈Spec(A),GM(λ)=AM(λ)⟺Aisdiagonalizable

例れい(対角化不可能たいかくかふかのうな場合ばあい):A=(1101) は λ=1(AM = 2)だが GM = 1。

5.44. 具体例ぐたいれい

A=(3102)

特性多項式とくせいたこうしきcharacteristic polynomial:det(A-λI)=(3-λ)(2-λ)=0 より λ=2,3。

λ=3 に対たいして (A-3I)v=0 → v∝(10)

λ=2 に対たいして (A-2I)v=0 → v∝(1-1)

5.55. Cayley-Hamilton の定理ていり

行列ぎょうれつmatrix A は自分自身じぶんじしんの特性多項式とくせいたこうしきcharacteristic polynomialを満みたす:

特性多項式とくせいたこうしきcharacteristic polynomialを q(t)=det(tI-A) と置おくと、

q(A)=0

である。ここで 0 は零行列れいぎょうれつzero matrixを表あらわす。記号きごうとして p(λ)=det(A-λI) を用もちいる流儀りゅうぎでは最高次係数さいこうじけいすうの符号ふごうが変かわるが、Cayley-Hamilton の内容ないようは「行列ぎょうれつmatrixが自分自身じぶんじしんの特性多項式とくせいたこうしきcharacteristic polynomialを満みたす」という命題めいだいである。

例れい:A=(abcd) の特性多項式とくせいたこうしきcharacteristic polynomialは λ2-(a+d)λ+(ad-bc)。Cayley-Hamilton より A2-(a+d)A+(ad-bc)I=0。

応用おうよう:A-1 を特性多項式とくせいたこうしきcharacteristic polynomialで表現ひょうげんする(detA≠0 のとき)。

5.66. 連立れんりつ微分方程式びぶんほうていしきへの応用おうよう

x′=Ax の解かいは x(t)=eAtx(0)。A が対角化可能たいかくかかのうで A=PDP-1 なら:

eAt=PeDtP-1,eDt=(eλ1t⋱eλnt)

固有値こゆうちeigenvalueの実部じつぶが安定性あんていせいを決きめる(負ふ → 安定あんてい、正せい → 不安定ふあんてい)。

5Precise explanation

5.11. From the definition to an equation

The definition of an eigenvector固有こゆうベクトル includes the nonzero condition:

Av=λv(v≠0).

This is equivalent to

(A-λI)v=0(v≠0).

At this point no determinant has been used; we have only moved all terms to one side. The determinant appears when we ask when the homogeneous system

(A-λI)x=0

has a nontrivial solution. If A-λI is invertible, the only solution is x=0. Therefore a nontrivial solution exists exactly when

det(A-λI)=0.

The two logical steps are

Av=λv(v≠0)⟺(A-λI)v=0(v≠0)

and

(A-λI)x=0hasanontrivialsolution⟺det(A-λI)=0.

5.22. Properties of the characteristic polynomial特性多項式とくせいたこうしき

For an n×n matrix行列ぎょうれつ A,

det(A-λI)=(-1)nλn+(-1)n-1(trA)λn-1+⋯+detA.

When the characteristic polynomial splits into linear factors and the roots are counted with algebraic multiplicity,

trA=∑i=1nλi,detA=∏i=1nλi.

For real matrices, real eigenvalues may not be enough. A rotation matrix can have complex eigenvalues, so one often extends scalars to the complex numbers when discussing all roots of the characteristic polynomial.

An important consequence is

detA=0⟺0isaneigenvalueofA.

Equivalently, A is singular exactly when 0 belongs to its spectrum.

5.33. Algebraic and geometric multiplicity

The algebraic multiplicity代数的重複度だいすうてきちょうふくど of an eigenvalue固有値こゆうち λ0 is its multiplicity as a root of the characteristic polynomial特性多項式とくせいたこうしき.

The geometric multiplicity幾何的重複度きかてきちょうふくど is the dimension of the eigenspace固有空間こゆうくうかん

Eλ0=ker(A-λ0I).

For an n×n matrix,

GM(λ0)=dimker(A-λ0I)=n-rank(A-λ0I).

For every eigenvalue,

1[PARSE ERROR: Undefined("Command(\"le\")")]GM(λ0)[PARSE ERROR: Undefined("Command(\"le\")")]AM(λ0).

If the characteristic polynomial splits over the field under consideration, then A is diagonalizable exactly when every eigenvalue has equal algebraic and geometric multiplicity:

Here Spec(A) is the set of all eigenvalues固有値こゆうち of A. A matrix A is diagonalizable対角化可能たいかくかかのう when there is a basis of eigenvectors, equivalently when A=PDP-1 for some invertible matrix P and diagonal matrix D. The proof of this criterion is treated in the next diagonalization lecture.

data/lecture/math/linear-algebra/diagonalization-basics.lecture.n.md
∀λ∈Spec(A),GM(λ)=AM(λ)⟺Aisdiagonalizable.

Example of failure: for

A=(1101),

λ=1 has AM =2 but GM =1, so A is not diagonalizable.

5.44. Concrete example

Let

A=(3102).

The characteristic polynomial is

det(A-λI)=(3-λ)(2-λ),

so the eigenvalues are λ=2,3. For λ=3,

(A-3I)v=0

gives v∝(10). For λ=2,

(A-2I)v=0

gives v∝(1-1).

5.55. The Cayley-Hamilton theorem

A matrix satisfies its own characteristic polynomial. If

q(t)=det(tI-A),

then

q(A)=0,

where 0 denotes the zero matrix. If one uses p(λ)=det(A-λI) instead, the leading coefficient differs by a sign, but the theorem still says that the matrix annihilates itself through its characteristic polynomial.

For

A=(abcd),

the characteristic polynomial is λ2-(a+d)λ+(ad-bc). Hence

A2-(a+d)A+(ad-bc)I=0.

When detA≠0, this relation can be rearranged to express A-1 using powers of A.

5.66. Application to systems of differential equations

For

x′=Ax,

the solution is

x(t)=eAtx(0).

If A=PDP-1 is diagonalizable, then

eAt=PeDtP-1,eDt=(eλ1t⋱eλnt).

The real parts of the eigenvalues固有値こゆうち determine stability: negative real parts indicate decay, while positive real parts indicate growth.

6判定基準はんていきじゅん

  • Av=λv → 特性多項式とくせいたこうしきcharacteristic polynomial det(A-λI)=0 を解とく
  • 対角化可能たいかくかかのうか → 固有値こゆうちeigenvalueがすべて異ことなる → 可か(十分条件じゅうぶんじょうけん)
  • trA と detA を素早すばやく → 固有値こゆうちeigenvalueの和わ・積せきで検算けんざん
  • 微分方程式びぶんほうていしき x′=Ax の安定性あんていせい → 固有値こゆうちeigenvalueの実部じつぶの符号ふごう

6Criteria

  • To find eigenvalues, solve det(A-λI)=0 and then solve (A-λI)v=0 for nonzero v.
  • If all eigenvalues are distinct, A is diagonalizable; this is sufficient but not necessary.
  • Use trA and detA as checks: they match the sum and product of eigenvalues counted with multiplicity.
  • For x′=Ax, stability is controlled by the signs of the real parts of eigenvalues.

7どこまで成立せいりつするか

実数じっすうreal number行列ぎょうれつmatrixにおいても固有値こゆうちeigenvalueは複素数ふくそすうcomplex numberになりうる(回転かいてん行列ぎょうれつmatrixの λ=e±iθ)。対角化たいかくかdiagonalizationの問題もんだいはスペクトル定理ていり(対称行列たいしょうぎょうれつsymmetric matrixは直交対角化可能ちょっこうたいかくかかのう)やジョルダン標準形ひょうじゅんけい(一般いっぱんの場合ばあい)へ拡張かくちょうされる。スペクトル定理ていりは対称行列たいしょうぎょうれつsymmetric matrixに限かぎらず正規行列せいきぎょうれつまで拡張かくちょうでき、作用素さようそ理論りろんでは微分作用素びぶんさようその固有関数こゆうかんすう(eλx)が同おなじ発想はっそうの無限むげん次元じげんdimensionへの拡張かくちょうである。

7Scope and limitations

Even real matrices can have complex eigenvalues固有値こゆうち. The diagonalization problem extends to the spectral theorem for symmetric or normal matrices and to Jordan normal form in the general case. In operator theory, eigenfunctions of differential operators, such as eλx, are infinite-dimensional versions of the same idea.

8最終形さいしゅうけい

[PARSE ERROR: Undefined("Command(\"boxed\")")]Av=λv(v≠0)⟺(A-λI)v=0(v≠0)
[PARSE ERROR: Undefined("Command(\"boxed\")")](A-λI)x=0が非自明解を持つ⟺det(A-λI)=0
[PARSE ERROR: Undefined("Command(\"boxed\")")]trA=∑λi,detA=∏λi
[PARSE ERROR: Undefined("Command(\"boxed\")")]p(A)=0(CayleyHamilton)[PARSE ERROR: Undefined("RBrace")]
[PARSE ERROR: Undefined("Command(\"boxed\")")]GM(λ)=dimker(A-λI),1[PARSE ERROR: Undefined("Command(\"le\")")]GM(λ)[PARSE ERROR: Undefined("Command(\"le\")")]AM(λ)

8Final forms

[PARSE ERROR: Undefined("Command(\"boxed\")")]Av=λv(v≠0)⟺(A-λI)v=0(v≠0)
[PARSE ERROR: Undefined("Command(\"boxed\")")](A-λI)x=0hasanontrivialsolution⟺det(A-λI)=0
[PARSE ERROR: Undefined("Command(\"boxed\")")]trA=∑λi,detA=∏λi
[PARSE ERROR: Undefined("Command(\"boxed\")")]p(A)=0Cayley-Hamilton)[PARSE ERROR: Undefined("RBrace")]
[PARSE ERROR: Undefined("Command(\"boxed\")")]GM(λ)=dimker(A-λI),1[PARSE ERROR: Undefined("Command(\"le\")")]GM(λ)[PARSE ERROR: Undefined("Command(\"le\")")]AM(λ)

9一言ひとことでいうと

固有値こゆうちeigenvalueは変換へんかんの「向むきを変かえない方向ほうこうでの倍率ばいりつ」であり、特性多項式とくせいたこうしきcharacteristic polynomialはその情報じょうほうを n 次じ多項式たこうしきに圧縮あっしゅくする。Cayley-Hamilton の定理ていりは行列ぎょうれつmatrixが自分自身じぶんじしんの特性多項式とくせいたこうしきcharacteristic polynomialを満みたすという対称性たいしょうせいを示しめし、対角化たいかくかdiagonalizationは固有ベクトルeigenvectorの「向むきが変かわらない軸じく」を基底きていbasisにした変換へんかんのもっとも単純たんじゅんな記述きじゅつである。

9In one sentence

An eigenvalue固有値こゆうち is the scale factor along a direction that the transformation does not turn, the characteristic polynomial特性多項式とくせいたこうしき compresses that information into a degree-n polynomial, and diagonalization対角化たいかくか uses eigenvectors as the coordinate axes of the transformation.

10定理ていりと証明しょうめい:異ことなる固有値こゆうちの固有ベクトルeigenvectorは線型独立せんけいどくりつである

A の固有値こゆうち λ1,…,λk が互たがいに異ことなり、vi が λi に対応たいおうする 0 でない固有ベクトルeigenvectorだとする。このとき v1,…,vk は線型独立せんけいどくりつである。

証明しょうめいは k についての帰納法きのうほうで行おこなう。k=1 は 0 でないベクトルvector一ひとつなので明あきらかである。

a1v1+⋯+akvk=0

とする。両辺りょうへんに A を作用さようさせると

a1λ1v1+⋯+akλkvk=0

である。最初さいしょの式しきに λk を掛かけて引ひくと

a1(λ1-λk)v1+⋯+ak-1(λk-1-λk)vk-1=0

である。帰納法きのうほうの仮定かていより v1,…,vk-1 は線型独立せんけいどくりつなので、ai(λi-λk)=0 である。λi≠λk だから ai=0 である。すると最初さいしょの式しきから akvk=0 となり、vk≠0 より ak=0 である。

したがって、すべての係数けいすうが 0 であり、線型独立せんけいどくりつである。

10Theorem and proof: eigenvectors for distinct eigenvalues are linearly independent

Let λ1,…,λk be distinct eigenvalues of A, and let vi≠0 be an eigenvector for λi. Then v1,…,vk are linearly independent線型独立せんけいどくりつ.

We prove this by induction on k. The case k=1 is immediate because one nonzero vector is linearly independent.

Assume

a1v1+⋯+akvk=0.

Applying A gives

a1λ1v1+⋯+akλkvk=0.

Subtract λk times the first equation:

a1(λ1-λk)v1+⋯+ak-1(λk-1-λk)vk-1=0.

By the induction hypothesis, v1,…,vk-1 are linearly independent, so ai(λi-λk)=0 for i<k. Since λi≠λk, we get ai=0 for i<k. The original equation then gives akvk=0, and vk≠0 implies ak=0.

All coefficients are zero, so the vectors are linearly independent.

raw .n.md をコピー
loc をコピー (filepath:line ~ line)
copy share link
copy encoded share link
path をコピー
copy share link
copy encoded share link
copy share link
copy encoded share link
タブを全て閉じる