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

An 正方行列せいほうぎょうれつsquare matrixとするとき、v0 かつ

Av=λv

たすスカラー λ固有値こゆうちEigenvaluev対応たいおうする固有こゆうベクトルEigenvectorという。

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

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

A特性多項式とくせいたこうしきCharacteristic polynomialという。n 多項式たこうしきであり、固有値こゆうちeigenvaluep(λ)=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=λvdet(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定義ていぎは、非零ひれいベクトル v0ふくめて

Av=λv(v0)

である。これは

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

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

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

det(A-λI)=0

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

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

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

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

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=00固有値こゆうちeigenvalueA特異行列とくいぎょうれつ)。

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

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

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

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

次元じげんdimensionである。An×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 polynomialdet(A-λI)=(3-λ)(2-λ)=0 より λ=2,3

λ=3たいして (A-3I)v=0v(10)

λ=2たいして (A-2I)v=0v(1-1)

5.55. Cayley-Hamilton の定理ていり

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

特性多項式とくせいたこうしきcharacteristic polynomialq(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表現ひょうげんする(detA0 のとき)。

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

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

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

固有値こゆうちeigenvalue実部じつぶ安定性あんていせいめる(安定あんていせい不安定ふあんてい)。

5Precise explanation

5.11. From the definition to an equation

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

Av=λv(v0).

This is equivalent to

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

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(v0)(A-λI)v=0(v0)

and

(A-λI)x=0hasanontrivialsolutiondet(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=00isaneigenvalueofA.

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 detA0, 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λ1teλ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がすべてことなる → 十分条件じゅうぶんじょうけん
  • trAdetA素早すばやく → 固有値こゆうち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(v0)(A-λI)v=0(v0)
[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(v0)(A-λI)v=0(v0)
[PARSE ERROR: Undefined("Command(\"boxed\")")](A-λI)x=0hasanontrivialsolutiondet(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 となり、vk0 より ak=0 である。

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

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

Let λ1,,λk be distinct eigenvalues of A, and let vi0 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 vk0 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
タブを全て閉じる