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

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

1導入どうにゅう

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

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

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 の解かいである。

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

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

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

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

つまり固有値こゆうちeigenvalue・固有ベクトルeigenvectorは、複雑ふくざつな線型変換せんけいへんかんlinear transformationを「方向ほうこうごとの倍率ばいりつ」として読よめる場所ばしょを探さがす道具どうぐである。この見方みかたが対角化たいかくか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の実部じつぶが安定性あんていせいを決きめる(負ふ → 安定あんてい、正せい → 不安定ふあんてい)。

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

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

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

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

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(λ)

9一言ひとことでいうと

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

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 であり、線型独立せんけいどくりつである。

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
タブを全て閉じる