Lemma数学ロードマップ

02 線形代数 · 第 5 章

固有値と対角化

目安 12〜16 時間定理など 18演習 10 問実験 1

この章の目標

  • 固有多項式を用いて固有値・固有ベクトル・固有空間を計算できる
  • 対角化可能であるための条件を重複度で判定し、実際に対角化できる
  • 三角化定理とケーリー–ハミルトンの定理を証明できる
  • 最小多項式を定義し、「最小多項式が重根をもたない   ⟺  \iff 対角化可能」を証明して使える
  • 同時対角化を理解し、行列のべき・漸化式・マルコフ連鎖に応用できる

前提:第1章〜第4章。多項式の割り算(余りのある除法、代数学 第5章)と、多項式の素因数分解の一意性(代数学 第6章)を使う。

5.1 動機:よい基底を探す

対角行列は扱いやすい。べき乗は対角成分のべき乗であり、連立方程式は成分ごとに解ける。例 3.28 では、行列 AA(第 1 行 (2,1)(2, 1)、第 2 行 (1,2)(1, 2))が、基底 u1=t(1,1)u_1 = {}^t(1, 1)、u2=t(1,−1)u_2 = {}^t(1, -1) に関しては対角行列 diag⁡(3,1)\operatorname{diag}(3, 1) で表されることを見た。これは Au1=3u1Au_1 = 3u_1、Au2=u2Au_2 = u_2、つまり LAL_A が u1u_1 方向と u2u_2 方向をそれぞれ定数倍するだけの変換だということである。

一般に、線形変換 f ⁣:V→Vf\colon V \to V の表現行列が基底 (v1,…,vn)(v_1, \dots, v_n) に関して対角行列 diag⁡(λ1,…,λn)\operatorname{diag}(\lambda_1, \dots, \lambda_n) になることは、f(vi)=λivif(v_i) = \lambda_i v_i がすべての ii で成り立つことと同じである。そこで、ff によって方向が変わらないベクトルを探すことになる。これが固有ベクトルである。

「行列は線形写像の座標表示」という視点では、対角化とは線形変換の本来の姿が最も見やすくなる座標系を探すことであり、行列 AA を P−1APP^{-1}AP に変える計算はその座標変換にほかならない。

以下、KK は体、VV は KK 上の nn 次元ベクトル空間(n≥1n \geq 1)、f ⁣:V→Vf\colon V \to V は線形変換とする。行列 A∈M⁡n(K)A \in \operatorname{M}_n(K) については LA ⁣:Kn→KnL_A\colon K^n \to K^n を考える。

5.2 固有値と固有ベクトル

定義 5.1(固有値・固有ベクトル・固有空間)λ∈K\lambda \in K に対し、f(v)=λvf(v) = \lambda v を満たす v≠0v \neq 0 が存在するとき、λ\lambda を ff の固有値 (eigenvalue)、vv を λ\lambda に属する固有ベクトル (eigenvector) という。部分空間

V(λ)=Ker⁡(f−λ idV)={v∈V∣f(v)=λv}V(\lambda) = \operatorname{Ker}(f - \lambda\ \mathrm{id}_V) = \lbrace v \in V \mid f(v) = \lambda v \rbrace

を λ\lambda の固有空間 (eigenspace) という。行列 AA の固有値・固有ベクトルは LAL_A のそれとする。

固有ベクトルは 0 でないと約束する(f(0)=λ0f(0) = \lambda 0 はすべての λ\lambda で成り立つので、0 を許すと意味がない)。固有空間は固有ベクトル全体に 0 を加えたものである。

例 5.2

  1. 対角行列 diag⁡(λ1,…,λn)\operatorname{diag}(\lambda_1, \dots, \lambda_n) の固有値は λ1,…,λn\lambda_1, \dots, \lambda_n で、eie_i は λi\lambda_i に属する固有ベクトルである。
  2. 回転行列 R(θ)R(\theta)(例 1.9)は、θ\theta が π\pi の整数倍でなければ実数の固有値をもたない。幾何学的には、原点を通るどの直線も回転で自分自身に写らないからである。ところが C\mathbb{C} 上で考えると、固有値 e±iθe^{\pm i\theta} をもつ(例 5.16)。
  3. 平面の直線 ℓ\ell への射影(例 3.23)の固有値は 11(固有空間 ℓ\ell)と 00(固有空間は射影の方向)である。直線に関する鏡映の固有値は 11 と −1-1 である。
  4. (無限次元の例)微分 D ⁣:C∞(R)→C∞(R)D\colon C^\infty(\mathbb{R}) \to C^\infty(\mathbb{R}) について、任意の λ∈R\lambda \in \mathbb{R} は固有値で、eλxe^{\lambda x} が固有ベクトル(固有関数)である。一方、標数 0 の体上で D ⁣:K[x]≤n→K[x]≤nD\colon K[x]_{\leq n} \to K[x]_{\leq n} の固有値は 00 だけである(次数を下げるので Dp=λpDp = \lambda p, λ≠0\lambda \neq 0 は p=0p = 0 を意味する)。
実験2 次元の線形写像行列で平面をゆがめ、行列式(面積の倍率)と固有ベクトル(向きが変わらない方向)を目で見ます。

この実験は JavaScript を有効にすると動きます。

命題 5.3 λ∈K\lambda \in K について次は同値である:(1) λ\lambda は ff の固有値、(2) f−λ idVf - \lambda\ \mathrm{id}_V は単射でない、(3) det⁡(f−λ idV)=0\det(f - \lambda\ \mathrm{id}_V) = 0。

証明. (1)   ⟺  \iff (2) は定義と命題 3.6 から。(2)   ⟺  \iff (3) は、有限次元では単射と同型が同値であること(系 3.12)と命題 4.37 から。□\square

5.3 固有多項式

命題 5.3 の条件 (3) を λ\lambda の方程式として扱うため、λ\lambda を変数 tt に置き換える。

定義 5.4(固有多項式, characteristic polynomial)A∈M⁡n(K)A \in \operatorname{M}_n(K) に対し、tt を変数として

ΦA(t)=det⁡(tIn−A)\Phi_A(t) = \det(tI_n - A)

を AA の固有多項式(特性多項式)という。方程式 ΦA(t)=0\Phi_A(t) = 0 を固有方程式という。

tIn−AtI_n - A は多項式環 K[t]K[t] に成分をもつ行列であり、その行列式は注意 4.22 により意味をもち、tt の多項式になる。相似な行列の固有多項式は等しい:

det⁡(tI−P−1AP)=det⁡(P−1(tI−A)P)=det⁡P−1det⁡(tI−A)det⁡P=det⁡(tI−A)\det(tI - P^{-1}AP) = \det(P^{-1}(tI - A)P) = \det P^{-1} \det(tI - A) \det P = \det(tI - A)

したがって線形変換 ff の固有多項式 Φf(t):=Φ[f]B(t)\Phi_f(t) := \Phi_{[f]_{\mathcal{B}}}(t) は基底によらずに定まる。

命題 5.5 ΦA(t)\Phi_A(t) は nn 次のモニック多項式(最高次の係数が 1)で

ΦA(t)=tn−(tr⁡A) tn−1+⋯+(−1)ndet⁡A\Phi_A(t) = t^n - (\operatorname{tr} A)\ t^{n-1} + \cdots + (-1)^n \det A

である。λ∈K\lambda \in K が AA の固有値であるための必要十分条件は ΦA(λ)=0\Phi_A(\lambda) = 0 である。

証明. 定義 4.10 で det⁡(tI−A)\det(tI - A) を展開する。σ=id\sigma = \mathrm{id} の項は ∏j(t−ajj)\prod_{j}(t - a_{jj}) である。σ≠id\sigma \neq \mathrm{id} なら σ\sigma は少なくとも 2 点を動かすので、その項が含む対角成分 t−ajjt - a_{jj} は n−2n - 2 個以下であり、項の次数は n−2n - 2 以下である。よって tnt^n と tn−1t^{n-1} の係数は ∏j(t−ajj)\prod_j (t - a_{jj}) から来て、それぞれ 11 と −∑jajj-\sum_j a_{jj} である。定数項は ΦA(0)=det⁡(−A)=(−1)ndet⁡A\Phi_A(0) = \det(-A) = (-1)^n \det A。後半:行列式は成分の多項式なので、tt に λ\lambda を代入してから行列式をとっても同じであり、ΦA(λ)=det⁡(λI−A)\Phi_A(\lambda) = \det(\lambda I - A)。これが 0 であることは命題 5.3 より λ\lambda が固有値であることと同値である。□\square

したがって固有値の計算は固有方程式を解くことに帰着され、固有ベクトルの計算は同次連立方程式 (A−λI)x=0(A - \lambda I)x = 0 を解くことに帰着される。

例 5.6 次の行列の固有値と固有ベクトルを求める。

A=(4123)A = \begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix}

ΦA(t)=(t−4)(t−3)−2=t2−7t+10=(t−2)(t−5)\Phi_A(t) = (t - 4)(t - 3) - 2 = t^2 - 7t + 10 = (t - 2)(t - 5) より固有値は 2,52, 5。λ=5\lambda = 5:(A−5I)x=0(A - 5I)x = 0 は −x1+x2=0-x_1 + x_2 = 0 なので V(5)=span⁡(t(1,1))V(5) = \operatorname{span}({}^t(1, 1))。λ=2\lambda = 2:2x1+x2=02x_1 + x_2 = 0 なので V(2)=span⁡(t(1,−2))V(2) = \operatorname{span}({}^t(1, -2))。検算:A t(1,1)=t(5,5)A\ {}^t(1, 1) = {}^t(5, 5)、A t(1,−2)=t(2,−4)A\ {}^t(1, -2) = {}^t(2, -4)。トレース 7=2+57 = 2 + 5、行列式 10=2⋅510 = 2 \cdot 5 も命題 5.5 と整合する。

例 5.7 次の行列を考える。

A=(122212221)A = \begin{pmatrix} 1 & 2 & 2 \\ 2 & 1 & 2 \\ 2 & 2 & 1 \end{pmatrix}

JJ をすべての成分が 1 の 3 次正方行列とすると A=2J−IA = 2J - I である。ΦA\Phi_A を直接計算してもよいが、命題 5.5 の係数を確かめる形で求めると:tr⁡A=3\operatorname{tr} A = 3、det⁡A=1⋅(1−4)−2⋅(2−4)+2⋅(4−2)=5\det A = 1 \cdot (1 - 4) - 2 \cdot (2 - 4) + 2 \cdot (4 - 2) = 5、tt の係数は 2 次の主小行列式の和 3×(1−4)=−93 \times (1 - 4) = -9(問題 5.2 参照)であり、

ΦA(t)=t3−3t2−9t−5=(t−5)(t+1)2\Phi_A(t) = t^3 - 3t^2 - 9t - 5 = (t - 5)(t + 1)^2

となる(展開して確かめよ)。固有値は 55 と −1-1 である。A+I=2JA + I = 2J なので V(−1)={x∣x1+x2+x3=0}V(-1) = \lbrace x \mid x_1 + x_2 + x_3 = 0 \rbrace は 2 次元で、基底 t(1,−1,0){}^t(1, -1, 0), t(1,0,−1){}^t(1, 0, -1) をもつ。V(5)V(5) は (A−5I)x=0(A - 5I)x = 0 を解いて span⁡(t(1,1,1))\operatorname{span}({}^t(1, 1, 1)) である。

定義 5.8(重複度)λ\lambda を ff の固有値とする。λ\lambda が Φf(t)\Phi_f(t) の根として重複する回数((t−λ)m(t - \lambda)^m が Φf\Phi_f を割り切る最大の mm)を λ\lambda の代数的重複度 m(λ)m(\lambda)、dim⁡V(λ)\dim V(\lambda) を幾何的重複度 d(λ)d(\lambda) という。

命題 5.9 各固有値 λ\lambda について 1≤d(λ)≤m(λ)1 \leq d(\lambda) \leq m(\lambda)。

証明. d=d(λ)≥1d = d(\lambda) \geq 1 は固有ベクトルの存在から。V(λ)V(\lambda) の基底 v1,…,vdv_1, \dots, v_d を VV の基底に延長すると、f(vi)=λvif(v_i) = \lambda v_i (i≤d)(i \leq d) より表現行列はブロック上三角行列

(λIdBOC)\begin{pmatrix} \lambda I_d & B \\ O & C \end{pmatrix}

になる。定理 4.34(K[t]K[t] 成分の行列にも適用できる、注意 4.22)より Φf(t)=(t−λ)dΦC(t)\Phi_f(t) = (t - \lambda)^d \Phi_C(t) なので m(λ)≥dm(\lambda) \geq d。□\square

例 5.10 次の行列を考える。

B=(11−13)B = \begin{pmatrix} 1 & 1 \\ -1 & 3 \end{pmatrix}

ΦB(t)=(t−1)(t−3)+1=(t−2)2\Phi_B(t) = (t - 1)(t - 3) + 1 = (t - 2)^2 なので固有値は 22 だけで m(2)=2m(2) = 2。B−2IB - 2I の第 1 行 (−1,1)(-1, 1)、第 2 行 (−1,1)(-1, 1) より rank⁡(B−2I)=1\operatorname{rank}(B - 2I) = 1、d(2)=2−1=1<m(2)d(2) = 2 - 1 = 1 < m(2) である。例 5.7 では d(−1)=m(−1)=2d(-1) = m(-1) = 2、d(5)=m(5)=1d(5) = m(5) = 1 であった。

代数学の基本定理(複素係数の多項式は C\mathbb{C} の中に根をもつ。証明は複素関数論 第3章)により、K=CK = \mathbb{C} なら ΦA\Phi_A は一次式の積に分解し、すべての複素正方行列は固有値をもつ。重複度を込めて数えると固有値はちょうど nn 個ある。K=RK = \mathbb{R} ではこれは成り立たない(例 5.2 の 2)。

5.4 対角化

定理 5.11 相異なる固有値 λ1,…,λk\lambda_1, \dots, \lambda_k に属する固有ベクトル v1,…,vkv_1, \dots, v_k は一次独立である。したがって、固有空間の和 V(λ1)+⋯+V(λk)V(\lambda_1) + \cdots + V(\lambda_k) は直和である。

証明. kk に関する帰納法。k=1k = 1 は v1≠0v_1 \neq 0 から。∑i=1kcivi=0\sum_{i=1}^k c_i v_i = 0 とし、両辺に f−λk idf - \lambda_k\ \mathrm{id} を施すと ∑i=1k−1ci(λi−λk)vi=0\sum_{i=1}^{k-1} c_i (\lambda_i - \lambda_k) v_i = 0。帰納法の仮定より ci(λi−λk)=0c_i(\lambda_i - \lambda_k) = 0、λi≠λk\lambda_i \neq \lambda_k から ci=0c_i = 0 (i<k)(i < k)。すると ckvk=0c_k v_k = 0 から ck=0c_k = 0。

直和であること:w1+⋯+wk=0w_1 + \cdots + w_k = 0(wi∈V(λi)w_i \in V(\lambda_i))とする。0 でない wiw_i たちは相異なる固有値の固有ベクトルなので、前半より、係数がすべて 1 のこの一次関係は 0 でない wiw_i が一つもない場合しかありえない。命題 2.41 より直和である。□\square

定義 5.12(対角化可能)VV の基底で、それに関する ff の表現行列が対角行列になるものが存在するとき、ff は対角化可能 (diagonalizable) であるという。行列 AA については、P−1APP^{-1}AP が対角行列になる P∈GL⁡n(K)P \in \operatorname{GL}_n(K) が存在することをいう。

表現行列が対角行列 diag⁡(λ1,…,λn)\operatorname{diag}(\lambda_1, \dots, \lambda_n) であることは、基底の各元が固有ベクトルであることと同値である。したがって、ff が対角化可能   ⟺  \iff VV が ff の固有ベクトルからなる基底をもつ。行列の言葉では、PP の第 jj 列を pjp_j とすると AP=Pdiag⁡(λ1,…,λn)  ⟺  Apj=λjpjAP = P\operatorname{diag}(\lambda_1, \dots, \lambda_n) \iff Ap_j = \lambda_j p_j であり、PP の列に固有ベクトルを並べれば、P−1APP^{-1}AP は対応する固有値を並べた対角行列になる。

定理 5.13(対角化可能性の判定)ff の相異なる固有値を λ1,…,λk\lambda_1, \dots, \lambda_k とする。次は同値である。

  1. ff は対角化可能である。
  2. V=V(λ1)⊕⋯⊕V(λk)V = V(\lambda_1) \oplus \cdots \oplus V(\lambda_k)。
  3. d(λ1)+⋯+d(λk)=nd(\lambda_1) + \cdots + d(\lambda_k) = n。
  4. Φf(t)\Phi_f(t) は K[t]K[t] で一次式の積に分解し、かつすべての ii で d(λi)=m(λi)d(\lambda_i) = m(\lambda_i)。

証明. (1 ⇒ 2):固有ベクトルからなる基底を固有値ごとに分けると、V=∑iV(λi)V = \sum_i V(\lambda_i) がわかり、定理 5.11 より直和である。(2 ⇒ 3):直和の次元(命題 2.41)。(3 ⇒ 1):各 V(λi)V(\lambda_i) の基底を合わせたものは、和が直和であること(定理 5.11)から一次独立で、元の個数が nn なので基底である(定理 2.31)。これは固有ベクトルからなる基底である。(3 ⇔ 4):命題 5.9 と、Φf\Phi_f の KK における根の重複度の和が次数 nn 以下であることから

∑id(λi)≤∑im(λi)≤n\sum_i d(\lambda_i) \leq \sum_i m(\lambda_i) \leq n

であり、右の等号は Φf\Phi_f が一次式の積に分解することと同値である。よって左端が nn に等しいことは、両方の等号が成り立つこと、すなわち 4 と同値である。□\square

系 5.14 Φf\Phi_f が KK に nn 個の相異なる根をもつ(ff が nn 個の相異なる固有値をもつ)ならば、ff は対角化可能である。

証明. このとき各 λi\lambda_i で 1≤d(λi)≤m(λi)=11 \leq d(\lambda_i) \leq m(\lambda_i) = 1 である。□\square

例 5.15(対角化の手順)例 5.7 の AA では d(5)+d(−1)=1+2=3d(5) + d(-1) = 1 + 2 = 3 なので対角化可能である。固有ベクトルを並べて

P=(1111−1010−1),P−1AP=(5000−1000−1)P = \begin{pmatrix} 1 & 1 & 1 \\ 1 & -1 & 0 \\ 1 & 0 & -1 \end{pmatrix}, \qquad P^{-1}AP = \begin{pmatrix} 5 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & -1 \end{pmatrix}

となる。検算は AP=PDAP = PD を列ごとに確かめればよい:A t(1,1,1)=t(5,5,5)A\ {}^t(1, 1, 1) = {}^t(5, 5, 5)、A t(1,−1,0)=t(−1,1,0)A\ {}^t(1, -1, 0) = {}^t(-1, 1, 0)、A t(1,0,−1)=t(−1,0,1)A\ {}^t(1, 0, -1) = {}^t(-1, 0, 1)。P−1P^{-1} を計算する必要はない。

例 5.16(対角化可能性は体による)例 5.10 の BB は d(2)<m(2)d(2) < m(2) なので、どんな体の上でも対角化可能でない。一方、回転行列 R(θ)R(\theta)(θ∉πZ\theta \notin \pi\mathbb{Z})は Φ(t)=t2−2cos⁡θ t+1\Phi(t) = t^2 - 2\cos\theta\ t + 1 が実根をもたないので R\mathbb{R} 上では対角化できないが、C\mathbb{C} 上では相異なる 2 根 e±iθe^{\pm i\theta} をもつので対角化可能である。実際、R(θ) t(1,−i)=t(cos⁡θ+isin⁡θ, sin⁡θ−icos⁡θ)=eiθ t(1,−i)R(\theta)\ {}^t(1, -i) = {}^t(\cos\theta + i\sin\theta,\ \sin\theta - i\cos\theta) = e^{i\theta}\ {}^t(1, -i) であり、同様に t(1,i){}^t(1, i) は e−iθe^{-i\theta} に属する固有ベクトルである。

注意

対角化 P−1AP=DP^{-1}AP = D で、PP の列の順序と DD の対角成分の順序は対応していなければならない。また、固有ベクトルは定数倍の自由度があるので PP は一意ではない。答えが教科書と違っても、AP=PDAP = PD と PP の正則性を確かめれば正誤は判定できる。

5.5 不変部分空間と三角化

対角化できない場合でも、表現行列をある程度簡単にすることはできる。鍵になるのは不変部分空間である。

定義 5.17(不変部分空間)部分空間 W⊂VW \subset V が f(W)⊂Wf(W) \subset W を満たすとき、WW は ff-不変 (invariant) であるという。このとき ff の WW への制限 f∣W ⁣:W→Wf\vert_W\colon W \to W が定まる。

Ker⁡f\operatorname{Ker} f、Im⁡f\operatorname{Im} f、固有空間は ff-不変である。より一般に、多項式 p(t)=∑cktkp(t) = \sum c_k t^k に対し p(f)=∑ckfkp(f) = \sum c_k f^k(f0=idf^0 = \mathrm{id})と定めると、p(f)p(f) は ff と可換なので、Ker⁡p(f)\operatorname{Ker} p(f) と Im⁡p(f)\operatorname{Im} p(f) は ff-不変である(v∈Ker⁡p(f)v \in \operatorname{Ker} p(f) なら p(f)(f(v))=f(p(f)(v))=0p(f)(f(v)) = f(p(f)(v)) = 0)。

命題 5.18 WW を ff-不変部分空間とし、WW の基底 w1,…,wkw_1, \dots, w_k を VV の基底 w1,…,wk,uk+1,…,unw_1, \dots, w_k, u_{k+1}, \dots, u_n に延長する。この基底に関する ff の表現行列は

(A1BOA2)\begin{pmatrix} A_1 & B \\ O & A_2 \end{pmatrix}

の形であり、A1=[f∣W]A_1 = [f\vert_W] である。A2A_2 は、ff が商空間に誘導する線形変換 fˉ ⁣:V/W→V/W\bar{f}\colon V/W \to V/W, v+W↦f(v)+Wv + W \mapsto f(v) + W の、基底 uk+1+W,…,un+Wu_{k+1} + W, \dots, u_n + W に関する表現行列であり、Φf=Φf∣WΦfˉ\Phi_f = \Phi_{f\vert_W} \Phi_{\bar{f}} が成り立つ。さらに V=W⊕W′V = W \oplus W' で W′W' も ff-不変なら、W′W' の基底をとることで B=OB = O(ブロック対角)にできる。

証明. f(wj)∈Wf(w_j) \in W は w1,…,wkw_1, \dots, w_k の一次結合なので左下は OO。fˉ\bar{f} は f(W)⊂Wf(W) \subset W より well-defined であり、f(uj)=∑i≤k(⋯ )wi+∑i>k(A2)ijuif(u_j) = \sum_{i \leq k}(\cdots)w_i + \sum_{i > k} (A_2)_{ij} u_i を WW を法として見れば A2A_2 が fˉ\bar{f} の表現行列であることがわかる(uk+1+W,…,un+Wu_{k+1} + W, \dots, u_n + W が V/WV/W の基底であることは、π\pi を span⁡(uk+1,…,un)\operatorname{span}(u_{k+1}, \dots, u_n) に制限したものが同型であることから従う)。固有多項式の式は定理 4.34 から。最後の主張は明らか。□\square

定理 5.19(三角化)ΦA(t)\Phi_A(t) が K[t]K[t] で一次式の積に分解するならば、P−1APP^{-1}AP が上三角行列になる P∈GL⁡n(K)P \in \operatorname{GL}_n(K) が存在する。このとき P−1APP^{-1}AP の対角成分には、AA の固有値が代数的重複度の回数だけ並ぶ。特に、任意の複素正方行列は三角化可能である。

証明. nn に関する帰納法。n=1n = 1 は自明。n≥2n \geq 2 とする。ΦA\Phi_A は一次式の積なので AA は固有値 λ1∈K\lambda_1 \in K をもつ。その固有ベクトル v1v_1 を含む KnK^n の基底をとり、それを列に並べた行列を P1P_1 とすると、命題 5.18(W=Kv1W = Kv_1)より

P1−1AP1=(λ1b0A1)(A1∈M⁡n−1(K), b は行ベクトル)P_1^{-1} A P_1 = \begin{pmatrix} \lambda_1 & b \\ 0 & A_1 \end{pmatrix} \qquad (A_1 \in \operatorname{M}_{n-1}(K),\ b \text{ は行ベクトル})

となる。ΦA(t)=(t−λ1)ΦA1(t)\Phi_A(t) = (t - \lambda_1)\Phi_{A_1}(t) なので ΦA1\Phi_{A_1} も一次式の積に分解する(K[t]K[t] の素因数分解の一意性)。帰納法の仮定により Q−1A1Q=T1Q^{-1}A_1Q = T_1 が上三角行列となる Q∈GL⁡n−1(K)Q \in \operatorname{GL}_{n-1}(K) がある。P=P1diag⁡(1,Q)P = P_1 \operatorname{diag}(1, Q)(左上が 11、右下が QQ のブロック対角行列との積)とおくと、ブロック積により

P−1AP=(100Q−1)(λ1b0A1)(100Q)=(λ1bQ0T1)P^{-1}AP = \begin{pmatrix} 1 & 0 \\ 0 & Q^{-1} \end{pmatrix} \begin{pmatrix} \lambda_1 & b \\ 0 & A_1 \end{pmatrix} \begin{pmatrix} 1 & 0 \\ 0 & Q \end{pmatrix} = \begin{pmatrix} \lambda_1 & bQ \\ 0 & T_1 \end{pmatrix}

は上三角行列である。上三角行列 TT の固有多項式は ∏i(t−tii)\prod_i (t - t_{ii}) なので、対角成分は ΦA\Phi_A の根を重複度込めて並べたものである。□\square

系 5.20 ΦA(t)=(t−λ1)⋯(t−λn)\Phi_A(t) = (t - \lambda_1) \cdots (t - \lambda_n)(重複を許す)と分解するとき、

tr⁡A=λ1+⋯+λn,det⁡A=λ1⋯λn\operatorname{tr} A = \lambda_1 + \cdots + \lambda_n, \qquad \det A = \lambda_1 \cdots \lambda_n

であり、任意の多項式 pp について Φp(A)(t)=(t−p(λ1))⋯(t−p(λn))\Phi_{p(A)}(t) = (t - p(\lambda_1)) \cdots (t - p(\lambda_n))。すなわち p(A)p(A) の固有値は p(λ1),…,p(λn)p(\lambda_1), \dots, p(\lambda_n) である。

証明. 前半は命題 5.5 の係数と根の関係から。後半:T=P−1APT = P^{-1}AP を上三角とすると p(T)=P−1p(A)Pp(T) = P^{-1}p(A)P も上三角で、対角成分は p(tii)=p(λi)p(t_{ii}) = p(\lambda_i)(並べ替えを除いて)。□\square

例えば実行列 AA の固有値を C\mathbb{C} の中で重複度込めて λ1,…,λn\lambda_1, \dots, \lambda_n とすれば、tr⁡Ak=∑iλik\operatorname{tr} A^k = \sum_i \lambda_i^k である。

5.6 ケーリー–ハミルトンの定理

定理 5.21(ケーリー–ハミルトンの定理, Cayley–Hamilton theorem)任意の A∈M⁡n(K)A \in \operatorname{M}_n(K) について ΦA(A)=O\Phi_A(A) = O。

証明. B(t)=adj⁡(tI−A)B(t) = \operatorname{adj}(tI - A) とおく。定理 4.28 は K[t]K[t] 成分の行列でも成り立つ(注意 4.22)ので

(tI−A)B(t)=ΦA(t)I(1)(tI - A) B(t) = \Phi_A(t) I \tag{1}

B(t)B(t) の各成分は tI−AtI - A の (n−1)(n-1) 次小行列式(の符号つき)なので tt について n−1n - 1 次以下の多項式である。よって B(t)=B0+B1t+⋯+Bn−1tn−1B(t) = B_0 + B_1 t + \cdots + B_{n-1} t^{n-1}(Bk∈M⁡n(K)B_k \in \operatorname{M}_n(K))と書ける。ΦA(t)=tn+cn−1tn−1+⋯+c0\Phi_A(t) = t^n + c_{n-1}t^{n-1} + \cdots + c_0 とし、(1) の両辺の tkt^k の係数(行列)を比較すると

tn:Bn−1=Itk (1≤k≤n−1):Bk−1−ABk=ckIt0:−AB0=c0I\begin{aligned} t^n &: \quad B_{n-1} = I \\ t^k \ (1 \leq k \leq n-1) &: \quad B_{k-1} - A B_k = c_k I \\ t^0 &: \quad -A B_0 = c_0 I \end{aligned}

tkt^k の式に左から AkA^k を掛けて k=0,…,nk = 0, \dots, n について足し合わせると、右辺は An+cn−1An−1+⋯+c0I=ΦA(A)A^n + c_{n-1}A^{n-1} + \cdots + c_0 I = \Phi_A(A) となる。左辺は

AnBn−1+∑k=1n−1(AkBk−1−Ak+1Bk)−AB0A^n B_{n-1} + \sum_{k=1}^{n-1} \left( A^k B_{k-1} - A^{k+1} B_k \right) - A B_0

であり、和の部分は ∑j=0n−2Aj+1Bj−∑j=1n−1Aj+1Bj=AB0−AnBn−1\sum_{j=0}^{n-2} A^{j+1}B_j - \sum_{j=1}^{n-1} A^{j+1}B_j = AB_0 - A^nB_{n-1} と打ち消し合うので、左辺全体は OO である。よって ΦA(A)=O\Phi_A(A) = O。□\square

注意

「ΦA(A)=det⁡(AI−A)=det⁡O=0\Phi_A(A) = \det(AI - A) = \det O = 0」は誤った証明である。ΦA(t)=det⁡(tI−A)\Phi_A(t) = \det(tI - A) の tt に行列 AA を代入するというのは、展開して得た多項式 tn+cn−1tn−1+⋯t^n + c_{n-1}t^{n-1} + \cdots に AA を代入することであり、行列式の記号の中の tt にスカラー行列として代入することではない。左辺は行列、右辺はスカラーで、そもそも型が合わない。

線形変換 ff についても、表現行列に適用すれば Φf(f)=0\Phi_f(f) = 0 が成り立つ。K=CK = \mathbb{C} の場合には三角化を用いた別証明がある(問題 5.8)。

例 5.22 例 5.6 の AA では ΦA(t)=t2−7t+10\Phi_A(t) = t^2 - 7t + 10 で、実際

A2=(1871411),A2−7A+10I=(18−28+107−714−1411−21+10)=OA^2 = \begin{pmatrix} 18 & 7 \\ 14 & 11 \end{pmatrix}, \qquad A^2 - 7A + 10I = \begin{pmatrix} 18 - 28 + 10 & 7 - 7 \\ 14 - 14 & 11 - 21 + 10 \end{pmatrix} = O

である。応用として、A(7I−A)=10IA(7I - A) = 10I から A−1=110(7I−A)A^{-1} = \frac{1}{10}(7I - A) が得られる。また t5t^5 を ΦA(t)\Phi_A(t) で割った余りを at+bat + b とすると A5=aA+bIA^5 = aA + bI である。t=2,5t = 2, 5 を代入して 32=2a+b32 = 2a + b、3125=5a+b3125 = 5a + b から a=1031a = 1031、b=−2030b = -2030、よって

A5=1031A−2030I=(2094103120621063)A^5 = 1031A - 2030I = \begin{pmatrix} 2094 & 1031 \\ 2062 & 1063 \end{pmatrix}

(例 5.34 の一般式でも確かめられる)。

5.7 最小多項式

ケーリー–ハミルトンの定理により、AA を代入して OO になる多項式(零化多項式)が存在する。そのうち次数最小のものが、対角化可能性の鍵を握る。

定義 5.23(最小多項式, minimal polynomial)p(A)=Op(A) = O を満たす 00 でない多項式のうち、次数が最小でモニックなものを AA の最小多項式といい、μA(t)\mu_A(t) と書く。

命題 5.24 p∈K[t]p \in K[t] について、p(A)=O  ⟺  μAp(A) = O \iff \mu_A が pp を割り切る。特に μA\mu_A はただ一つに定まり、ΦA\Phi_A を割り切る(したがって deg⁡μA≤n\deg \mu_A \leq n)。相似な行列の最小多項式は等しい。

証明. pp を μA\mu_A で割って p=qμA+rp = q\mu_A + r(r=0r = 0 または deg⁡r<deg⁡μA\deg r < \deg \mu_A)とすると、p(A)=q(A)μA(A)+r(A)=r(A)p(A) = q(A)\mu_A(A) + r(A) = r(A)。p(A)=Op(A) = O なら r(A)=Or(A) = O で、μA\mu_A の次数の最小性から r=0r = 0。逆は明らか。一意性:二つの最小多項式は互いに割り切り合い、ともにモニックなので等しい。ΦA(A)=O\Phi_A(A) = O(定理 5.21)より μA∣ΦA\mu_A \mid \Phi_A。最後の主張は p(P−1AP)=P−1p(A)Pp(P^{-1}AP) = P^{-1}p(A)P から。□\square

命題 5.25 λ∈K\lambda \in K について、μA(λ)=0  ⟺  λ\mu_A(\lambda) = 0 \iff \lambda は AA の固有値。

証明. Av=λvAv = \lambda v (v≠0)(v \neq 0) なら、任意の多項式 pp で p(A)v=p(λ)vp(A)v = p(\lambda)v である。p=μAp = \mu_A とすると 0=μA(λ)v0 = \mu_A(\lambda)v で、v≠0v \neq 0 より μA(λ)=0\mu_A(\lambda) = 0。逆に μA(λ)=0\mu_A(\lambda) = 0 なら、μA∣ΦA\mu_A \mid \Phi_A より ΦA(λ)=0\Phi_A(\lambda) = 0 なので固有値である。□\square

つまり μA\mu_A と ΦA\Phi_A は同じ根をもつが、重複度が違いうる。次の補題は、第6章でも使う基本的な道具である。

補題 5.26(ベズーの等式)p1,p2∈K[t]p_1, p_2 \in K[t] が互いに素(1 以外にモニックな公約元をもたない)ならば、ap1+bp2=1a p_1 + b p_2 = 1 となる a,b∈K[t]a, b \in K[t] が存在する。

証明. p1,p2p_1, p_2 の少なくとも一方は 0 でない。I={ap1+bp2∣a,b∈K[t]}I = \lbrace a p_1 + b p_2 \mid a, b \in K[t] \rbrace の 0 でない元のうち次数最小のものを dd とし、モニックにしておく。p1=qd+rp_1 = qd + r と割ると r=p1−qd∈Ir = p_1 - qd \in I で deg⁡r<deg⁡d\deg r < \deg d なので r=0r = 0、すなわち d∣p1d \mid p_1。同様に d∣p2d \mid p_2 なので d=1d = 1 であり、1∈I1 \in I。□\square

補題 5.27(核の分解)p=p1p2⋯pkp = p_1 p_2 \cdots p_k で、p1,…,pk∈K[t]p_1, \dots, p_k \in K[t] はどの 2 つも互いに素であるとする。このとき

Ker⁡p(f)=Ker⁡p1(f)⊕⋯⊕Ker⁡pk(f)\operatorname{Ker} p(f) = \operatorname{Ker} p_1(f) \oplus \cdots \oplus \operatorname{Ker} p_k(f)

であり、各 Ker⁡pi(f)\operatorname{Ker} p_i(f) は ff-不変である。特に p(f)=0p(f) = 0 なら V=⨁iKer⁡pi(f)V = \bigoplus_i \operatorname{Ker} p_i(f)。

証明. 不変性は前節で述べた。k=2k = 2 の場合を示す。Wi=Ker⁡pi(f)W_i = \operatorname{Ker} p_i(f)、U=Ker⁡p(f)U = \operatorname{Ker} p(f) とする。p(f)=p2(f)p1(f)=p1(f)p2(f)p(f) = p_2(f)p_1(f) = p_1(f)p_2(f) なので W1,W2⊂UW_1, W_2 \subset U。補題 5.26 の ap1+bp2=1a p_1 + b p_2 = 1 に ff を代入して a(f)p1(f)+b(f)p2(f)=idVa(f)p_1(f) + b(f)p_2(f) = \mathrm{id}_V。u∈Uu \in U に対し

u=a(f)p1(f)u⏟=:u2+b(f)p2(f)u⏟=:u1u = \underbrace{a(f)p_1(f)u}_{=: u_2} + \underbrace{b(f)p_2(f)u}_{=: u_1}

と書くと、p2(f)u2=a(f)p(f)u=0p_2(f)u_2 = a(f)p(f)u = 0 より u2∈W2u_2 \in W_2、同様に u1∈W1u_1 \in W_1。よって U=W1+W2U = W_1 + W_2。v∈W1∩W2v \in W_1 \cap W_2 なら v=a(f)p1(f)v+b(f)p2(f)v=0v = a(f)p_1(f)v + b(f)p_2(f)v = 0 なので直和である。

k≥3k \geq 3 は帰納法による。ajp1+bjpj=1a_j p_1 + b_j p_j = 1 (j=2,…,k)(j = 2, \dots, k) を辺々掛けて展開すると、p1p_1 の倍元と (b2⋯bk)(p2⋯pk)(b_2 \cdots b_k)(p_2 \cdots p_k) の和が 1 に等しくなるので、p1p_1 と p2⋯pkp_2 \cdots p_k は互いに素である。k=2k = 2 の場合を適用して U=W1⊕Ker⁡(p2⋯pk)(f)U = W_1 \oplus \operatorname{Ker}(p_2 \cdots p_k)(f) とし、後者に帰納法の仮定を適用すればよい(直和の結合性は命題 2.41 から従う)。□\square

証明からわかるように、UU から各成分への射影は ff の多項式で書ける。これは後で(第6章、第7章)有用になる。

定理 5.28(最小多項式による対角化可能性の判定)次は同値である。

  1. ff は対角化可能である。
  2. μf(t)=(t−λ1)⋯(t−λk)\mu_f(t) = (t - \lambda_1) \cdots (t - \lambda_k)(λi∈K\lambda_i \in K は相異なる)、すなわち μf\mu_f は KK 上で一次式の積に分解し、重根をもたない。
  3. KK 上で一次式の積に分解し重根をもたない多項式 pp で、p(f)=0p(f) = 0 となるものが存在する。

証明. (1 ⇒ 3):固有ベクトルからなる基底 v1,…,vnv_1, \dots, v_n をとり、相異なる固有値を λ1,…,λk\lambda_1, \dots, \lambda_k として p=∏i(t−λi)p = \prod_i (t - \lambda_i) とおく。vjv_j の固有値を λi\lambda_i とすると、p(f)vjp(f)v_j は因子 (f−λi)(f - \lambda_i) を含み(因子の順序は交換できる)0 である。よって p(f)=0p(f) = 0。

(3 ⇒ 2):μf∣p\mu_f \mid p(命題 5.24)であり、重根のない一次式の積の約元はまた重根のない一次式の積である。

(2 ⇒ 1):t−λ1,…,t−λkt - \lambda_1, \dots, t - \lambda_k はどの 2 つも互いに素(λi≠λj\lambda_i \neq \lambda_j なので (t−λi)−(t−λj)(t - \lambda_i) - (t - \lambda_j) は 0 でない定数)なので、補題 5.27 より V=⨁iKer⁡(f−λi)=⨁iV(λi)V = \bigoplus_i \operatorname{Ker}(f - \lambda_i) = \bigoplus_i V(\lambda_i)。定理 5.13 より対角化可能である。□\square

例 5.29 例 5.7 の AA では ΦA=(t−5)(t+1)2\Phi_A = (t - 5)(t + 1)^2 で、(A−5I)(A+I)=(2J−6I)(2J)=4J2−12J=O(A - 5I)(A + I) = (2J - 6I)(2J) = 4J^2 - 12J = O(J2=3JJ^2 = 3J を使った)なので、μA=(t−5)(t+1)\mu_A = (t - 5)(t + 1) であり(A−5I≠OA - 5I \neq O、A+I≠OA + I \neq O だから次数 1 の零化多項式はない)、対角化可能であることが再確認できる。例 5.10 の BB では B−2I≠OB - 2I \neq O、(B−2I)2=O(B - 2I)^2 = O なので μB=(t−2)2\mu_B = (t - 2)^2 で、重根をもつから対角化可能でない。

系 5.30

  1. A2=AA^2 = A を満たす行列(射影)は対角化可能で、固有値は 0,10, 1 のみである。
  2. K=CK = \mathbb{C} で Am=IA^m = I (m∈Nm \in \mathbb{N}) ならば AA は対角化可能で、固有値は 1 の mm 乗根である。
  3. 冪零行列(Am=OA^m = O)が対角化可能なら A=OA = O である。

証明. 1:t2−t=t(t−1)t^2 - t = t(t - 1)。2:tm−1t^m - 1 は C\mathbb{C} で相異なる mm 個の根 e2πik/me^{2\pi i k/m} をもつ。3:固有値はすべて 0(Av=λvAv = \lambda v なら 0=Amv=λmv0 = A^m v = \lambda^m v)なので、対角化すると P−1AP=OP^{-1}AP = O。□\square

命題 5.31 ff が対角化可能で WW が ff-不変部分空間ならば、f∣Wf\vert_W も対角化可能である。

証明. μf(f∣W)=μf(f)∣W=0\mu_f(f\vert_W) = \mu_f(f)\vert_W = 0 なので、μf∣W\mu_{f\vert_W} は μf\mu_f を割り切り、重根のない一次式の積である。定理 5.28 を f∣Wf\vert_W に適用する。□\square

5.8 同時対角化

複数の線形変換を同じ基底で同時に対角化できるのはいつか。

定理 5.32(同時対角化)f,g ⁣:V→Vf, g\colon V \to V がともに対角化可能で f∘g=g∘ff \circ g = g \circ f ならば、ff と gg の両方の固有ベクトルからなる VV の基底が存在する。逆に、そのような基底があれば f∘g=g∘ff \circ g = g \circ f である。

証明. ff の固有値を λ1,…,λk\lambda_1, \dots, \lambda_k とし、V=⨁iV(λi)V = \bigoplus_i V(\lambda_i) と分解する。v∈V(λi)v \in V(\lambda_i) なら f(g(v))=g(f(v))=λig(v)f(g(v)) = g(f(v)) = \lambda_i g(v) なので、各 V(λi)V(\lambda_i) は gg-不変である。命題 5.31 より g∣V(λi)g\vert_{V(\lambda_i)} は対角化可能であり、V(λi)V(\lambda_i) は gg の固有ベクトルからなる基底をもつ。これらは V(λi)V(\lambda_i) の元なので ff の固有ベクトルでもある。ii について合わせれば求める基底が得られる。逆:その基底に関して f,gf, g の表現行列はともに対角行列で、対角行列どうしは可換なので f∘g=g∘ff \circ g = g \circ f(定理 3.21)。□\square

帰納法により、互いに可換な有限個の対角化可能な線形変換も同時対角化できる。

例 5.33 次の二つの行列は可換である(AB=BAAB = BA は第 1 行 (1,2)(1, 2)、第 2 行 (2,1)(2, 1) の行列)。

A=(2112),B=(0110)A = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix}, \qquad B = \begin{pmatrix} 0 & 1 \\ 1 & 0 \end{pmatrix}

t(1,1){}^t(1, 1) は AA の固有値 3、BB の固有値 1 の固有ベクトル、t(1,−1){}^t(1, -1) は AA の固有値 1、BB の固有値 −1-1 の固有ベクトルであり、これらを並べた例 3.28 の行列 PP によって A,BA, B は同時に対角化される。実は A=2I+BA = 2I + B なので、これは当然である。

同時対角化は、量子力学における「同時に観測可能な物理量」や、可換群の表現の分解(リー群と表現論 第1章)の基礎となる考え方である。

5.9 応用

行列のべき

A=PDP−1A = PDP^{-1}(DD は対角行列)なら Ak=PDkP−1A^k = PD^kP^{-1} であり、DkD^k は対角成分の kk 乗である。

例 5.34 例 5.6 の AA で、固有ベクトルを並べた PP とその逆行列は

P=(111−2),P−1=13(211−1),P−1AP=(5002)P = \begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix}, \qquad P^{-1} = \frac{1}{3} \begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix}, \qquad P^{-1}AP = \begin{pmatrix} 5 & 0 \\ 0 & 2 \end{pmatrix}

であるから

Ak=P(5k002k)P−1=13(2⋅5k+2k5k−2k2⋅5k−2k+15k+2k+1)A^k = P \begin{pmatrix} 5^k & 0 \\ 0 & 2^k \end{pmatrix} P^{-1} = \frac{1}{3} \begin{pmatrix} 2 \cdot 5^k + 2^k & 5^k - 2^k \\ 2 \cdot 5^k - 2^{k+1} & 5^k + 2^{k+1} \end{pmatrix}

k=0k = 0 で II、k=1k = 1 で AA、k=5k = 5 で例 5.22 の結果になることを確かめよ。

線形漸化式

例 5.35(フィボナッチ数列)F0=0F_0 = 0、F1=1F_1 = 1、Fn+2=Fn+1+FnF_{n+2} = F_{n+1} + F_n とする。xn=t(Fn+1,Fn)x_n = {}^t(F_{n+1}, F_n) とおくと xn+1=Cxnx_{n+1} = Cx_n、ここで

C=(1110)C = \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}

であり、xn=Cnx0x_n = C^n x_0、x0=t(1,0)x_0 = {}^t(1, 0) となる。ΦC(t)=t2−t−1\Phi_C(t) = t^2 - t - 1 の根は φ=(1+5)/2\varphi = (1 + \sqrt{5})/2、ψ=(1−5)/2\psi = (1 - \sqrt{5})/2 で、相異なるので CC は対角化可能である。C t(λ,1)=t(λ+1,λ)=λ t(λ,1)C\ {}^t(\lambda, 1) = {}^t(\lambda + 1, \lambda) = \lambda\ {}^t(\lambda, 1)(λ2=λ+1\lambda^2 = \lambda + 1 による)なので t(λ,1){}^t(\lambda, 1) が固有ベクトルである。x0=t(1,0)=15(t(φ,1)−t(ψ,1))x_0 = {}^t(1, 0) = \frac{1}{\sqrt{5}}\left( {}^t(\varphi, 1) - {}^t(\psi, 1) \right)(φ−ψ=5\varphi - \psi = \sqrt{5})より

xn=15(φn(φ1)−ψn(ψ1)),Fn=φn−ψn5x_n = \frac{1}{\sqrt{5}} \left( \varphi^n \begin{pmatrix} \varphi \\ 1 \end{pmatrix} - \psi^n \begin{pmatrix} \psi \\ 1 \end{pmatrix} \right), \qquad F_n = \frac{\varphi^n - \psi^n}{\sqrt{5}}

(ビネの公式)。∣ψ∣<1\lvert \psi \rvert < 1 なので、nn が大きいとき FnF_n は φn/5\varphi^n/\sqrt{5} に最も近い整数であり、Fn+1/Fn→φF_{n+1}/F_n \to \varphi(黄金比)となる。

一般に漸化式 an+2=pan+1+qana_{n+2} = pa_{n+1} + qa_n の解全体は 2 次元のベクトル空間であり(a0,a1a_0, a_1 で決まる)、特性方程式 t2−pt−q=0t^2 - pt - q = 0 が相異なる根 α,β\alpha, \beta をもつなら (αn)(\alpha^n)、(βn)(\beta^n) がその基底となる:an=c1αn+c2βna_n = c_1\alpha^n + c_2\beta^n。重根の場合は対角化できず、(nαn)(n\alpha^n) のような解が現れる。これは第6章のジョルダン標準形で説明される。

マルコフ連鎖

例 5.36 ある町の天気が、晴れの翌日は確率 0.90.9 で晴れ・0.10.1 で雨、雨の翌日は確率 0.50.5 で晴れ・0.50.5 で雨になるとする。kk 日目に晴れ・雨である確率を並べた xk=t(pk,qk)x_k = {}^t(p_k, q_k) は xk+1=Pxkx_{k+1} = Px_k を満たす。ここで

P=(0.90.50.10.5)P = \begin{pmatrix} 0.9 & 0.5 \\ 0.1 & 0.5 \end{pmatrix}

は各列の和が 1 の確率行列(推移行列)である。ΦP(t)=t2−1.4t+0.4=(t−1)(t−0.4)\Phi_P(t) = t^2 - 1.4t + 0.4 = (t - 1)(t - 0.4) で、固有値 1 の固有ベクトルは t(5,1){}^t(5, 1)、固有値 0.40.4 の固有ベクトルは t(1,−1){}^t(1, -1) である。初日が晴れ(x0=t(1,0)x_0 = {}^t(1, 0))とすると x0=16 t(5,1)+16 t(1,−1)x_0 = \frac{1}{6}\ {}^t(5, 1) + \frac{1}{6}\ {}^t(1, -1) なので

xk=16(51)+(0.4)k6(1−1)⟶(5/61/6)(k→∞)x_k = \frac{1}{6}\begin{pmatrix} 5 \\ 1 \end{pmatrix} + \frac{(0.4)^k}{6} \begin{pmatrix} 1 \\ -1 \end{pmatrix} \longrightarrow \begin{pmatrix} 5/6 \\ 1/6 \end{pmatrix} \quad (k \to \infty)

となる。長期的には晴れの確率が 5/65/6 の定常分布に近づき、その速さは第 2 固有値の絶対値 0.40.4 で決まる。一般の確率行列が固有値 1 をもつこと(tP{}^t P が固有ベクトル t(1,…,1){}^t(1, \dots, 1) をもつことから)や、定常分布への収束の条件(ペロン–フロベニウスの定理)は確率論 第6章で扱う。

補足

確率論の本では、確率を行ベクトルに並べ、各行の和が 1 の行列を右から掛ける流儀が多い。ここでの PP はその転置にあたる。どちらの流儀でも、固有値の計算は同じである。

まとめ

  • 固有値 λ\lambda と固有ベクトル v≠0v \neq 0 は f(v)=λvf(v) = \lambda v で定義され、固有値は固有多項式 ΦA(t)=det⁡(tI−A)\Phi_A(t) = \det(tI - A) の根である。ΦA\Phi_A は相似不変で、係数に tr⁡A\operatorname{tr} A と det⁡A\det A が現れる。
  • 幾何的重複度 ≤\le 代数的重複度。相異なる固有値の固有ベクトルは一次独立で、固有空間の和は直和である。
  • 対角化可能   ⟺  \iff 固有ベクトルからなる基底がある   ⟺  \iff Φf\Phi_f が一次式の積に分解し、各固有値で d(λ)=m(λ)d(\lambda) = m(\lambda)。相異なる固有値が nn 個あれば対角化可能。
  • 固有多項式が一次式の積に分解すれば三角化できる(特に複素行列は常に三角化可能)。
  • ケーリー–ハミルトンの定理 ΦA(A)=O\Phi_A(A) = O は余因子行列の等式から証明できる。
  • 最小多項式 μA\mu_A は零化多項式をすべて割り切り、固有多項式と同じ根をもつ。ff が対角化可能   ⟺  \iff μf\mu_f が重根をもたない一次式の積。その証明の鍵は、互いに素な因子への核の分解(ベズーの等式)である。
  • 可換な対角化可能変換は同時対角化できる。
  • 対角化により、行列のべき、線形漸化式、マルコフ連鎖の長期挙動が計算できる。

演習問題

問題 5.1 ★ 次の行列 AA を対角化せよ(P−1APP^{-1}AP が対角行列となる PP と、その対角行列を求めよ)。

A=(200130111)A = \begin{pmatrix} 2 & 0 & 0 \\ 1 & 3 & 0 \\ 1 & 1 & 1 \end{pmatrix}
解答

下三角行列なので固有値は 2,3,12, 3, 1(相異なるので対角化可能)。

  • λ=2\lambda = 2:(A−2I)x=0(A - 2I)x = 0 は x1+x2=0x_1 + x_2 = 0、x1+x2−x3=0x_1 + x_2 - x_3 = 0 で、t(1,−1,0){}^t(1, -1, 0)。
  • λ=3\lambda = 3:−x1=0-x_1 = 0、x1=0x_1 = 0、x1+x2−2x3=0x_1 + x_2 - 2x_3 = 0 で、t(0,2,1){}^t(0, 2, 1)。
  • λ=1\lambda = 1:x1=0x_1 = 0、x1+2x2=0x_1 + 2x_2 = 0 で、t(0,0,1){}^t(0, 0, 1)。
P=(100−120011),P−1AP=(200030001)P = \begin{pmatrix} 1 & 0 & 0 \\ -1 & 2 & 0 \\ 0 & 1 & 1 \end{pmatrix}, \qquad P^{-1}AP = \begin{pmatrix} 2 & 0 & 0 \\ 0 & 3 & 0 \\ 0 & 0 & 1 \end{pmatrix}

検算:A t(1,−1,0)=t(2,−2,0)A\ {}^t(1, -1, 0) = {}^t(2, -2, 0)、A t(0,2,1)=t(0,6,3)A\ {}^t(0, 2, 1) = {}^t(0, 6, 3)、A t(0,0,1)=t(0,0,1)A\ {}^t(0, 0, 1) = {}^t(0, 0, 1)。

問題 5.2 ★ A∈M⁡3(K)A \in \operatorname{M}_3(K) の固有多項式は t3−(tr⁡A)t2+c1t−det⁡At^3 - (\operatorname{tr} A)t^2 + c_1 t - \det A の形で、c1c_1 は 3 つの 2 次主小行列式({1,2}\lbrace 1,2 \rbrace, {1,3}\lbrace 1,3 \rbrace, {2,3}\lbrace 2,3 \rbrace 行・列を選んだ小行列式)の和であることを示せ。また AA と tA{}^t A は同じ固有多項式をもつが、固有ベクトルは一般に異なることを例で示せ。

解答

det⁡(tI−A)\det(tI - A) の各列は tej−ajte_j - a_j である。多重線形性で展開すると、tt の 1 次の項は、2 つの列で −aj-a_j を、残り 1 つの列 kk で tekte_k を選んだ項の和である。この項を第 kk 列で展開すると、tt と「AA から第 kk 行・第 kk 列を除いた 2 次行列の行列式」との積になる(他の 2 列が −1-1 倍されているので符号は (−1)2=1(-1)^2 = 1)。よって c1c_1 は 2 次主小行列式の和である。

後半:det⁡(tI−tA)=det⁡t(tI−A)=det⁡(tI−A)\det(tI - {}^t A) = \det {}^t(tI - A) = \det(tI - A)。例:第 1 行 (1,1)(1, 1)、第 2 行 (0,1)(0, 1) の行列の固有ベクトルは e1e_1 の定数倍のみだが、その転置(第 1 行 (1,0)(1, 0)、第 2 行 (1,1)(1, 1))の固有ベクトルは e2e_2 の定数倍のみである。

問題 5.3 ★ 次の行列の nn 乗を求めよ。

A=(1221)A = \begin{pmatrix} 1 & 2 \\ 2 & 1 \end{pmatrix}
解答

ΦA=(t−1)2−4=(t−3)(t+1)\Phi_A = (t - 1)^2 - 4 = (t - 3)(t + 1)。固有ベクトルは t(1,1){}^t(1, 1)(固有値 3)、t(1,−1){}^t(1, -1)(固有値 −1-1)。PP を例 3.28 のものとすると P−1=12PP^{-1} = \frac{1}{2}P で

An=P(3n00(−1)n)12P=12(3n+(−1)n3n−(−1)n3n−(−1)n3n+(−1)n)A^n = P \begin{pmatrix} 3^n & 0 \\ 0 & (-1)^n \end{pmatrix} \frac{1}{2} P = \frac{1}{2}\begin{pmatrix} 3^n + (-1)^n & 3^n - (-1)^n \\ 3^n - (-1)^n & 3^n + (-1)^n \end{pmatrix}

n=1n = 1 で AA に戻ることを確かめよ。

問題 5.4 ★★ a∈Ra \in \mathbb{R} とする。次の行列が対角化可能であるための aa の条件を求めよ。

A=(2a1020001)A = \begin{pmatrix} 2 & a & 1 \\ 0 & 2 & 0 \\ 0 & 0 & 1 \end{pmatrix}
解答

上三角なので ΦA=(t−2)2(t−1)\Phi_A = (t - 2)^2(t - 1)、m(2)=2m(2) = 2、m(1)=1m(1) = 1。対角化可能   ⟺  d(2)=2  ⟺  rank⁡(A−2I)=1\iff d(2) = 2 \iff \operatorname{rank}(A - 2I) = 1。A−2IA - 2I の零でない行は (0,a,1)(0, a, 1) と (0,0,−1)(0, 0, -1) で、a≠0a \neq 0 なら一次独立(階数 2)、a=0a = 0 なら (0,0,1)(0, 0, 1) と (0,0,−1)(0, 0, -1) で階数 1。よって対角化可能   ⟺  a=0\iff a = 0。

問題 5.5 ★★ A,B∈M⁡n(K)A, B \in \operatorname{M}_n(K) に対し、ABAB と BABA の固有多項式は等しいことを示せ。(ヒント:次の二つの 2n2n 次行列が相似であることを示す。)

(ABOBO),(OOBBA)\begin{pmatrix} AB & O \\ B & O \end{pmatrix}, \qquad \begin{pmatrix} O & O \\ B & BA \end{pmatrix}
解答

SS を左上・右下が InI_n、右上が AA、左下が OO のブロック行列とする(正則、例 1.38)。ブロック積で

(ABOBO)S=(ABABABBA)=S(OOBBA)\begin{pmatrix} AB & O \\ B & O \end{pmatrix} S = \begin{pmatrix} AB & ABA \\ B & BA \end{pmatrix} = S \begin{pmatrix} O & O \\ B & BA \end{pmatrix}

が確かめられるので、二つの行列は相似であり、固有多項式が等しい。定理 4.34(ブロック三角)より前者の固有多項式は ΦAB(t)⋅tn\Phi_{AB}(t) \cdot t^n、後者は tn⋅ΦBA(t)t^n \cdot \Phi_{BA}(t)。K[t]K[t] では tnt^n で約せる(整域なので)ので ΦAB=ΦBA\Phi_{AB} = \Phi_{BA}。

(注:ABAB と BABA は相似とは限らない。A=E11A = E_{11}、B=E12B = E_{12} なら AB=E12≠OAB = E_{12} \neq O、BA=OBA = O。)

問題 5.6 ★★ 次の行列の最小多項式を求め、対角化可能かどうか判定せよ。

A=(210020002),B=(200020003),C=(210021002)A = \begin{pmatrix} 2 & 1 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 2 \end{pmatrix}, \qquad B = \begin{pmatrix} 2 & 0 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 3 \end{pmatrix}, \qquad C = \begin{pmatrix} 2 & 1 & 0 \\ 0 & 2 & 1 \\ 0 & 0 & 2 \end{pmatrix}
解答
  • AA:ΦA=(t−2)3\Phi_A = (t - 2)^3。A−2I=E12≠OA - 2I = E_{12} \neq O、(A−2I)2=E122=O(A - 2I)^2 = E_{12}^2 = O より μA=(t−2)2\mu_A = (t - 2)^2。重根があるので対角化不可能。
  • BB:対角行列なので対角化可能で、μB=(t−2)(t−3)\mu_B = (t - 2)(t - 3)(ΦB=(t−2)2(t−3)\Phi_B = (t-2)^2(t-3) と異なることに注意)。
  • CC:C−2I=E12+E23C - 2I = E_{12} + E_{23}、(C−2I)2=E13≠O(C - 2I)^2 = E_{13} \neq O、(C−2I)3=O(C - 2I)^3 = O より μC=(t−2)3=ΦC\mu_C = (t - 2)^3 = \Phi_C。対角化不可能。

AA と CC は固有多項式が同じだが最小多項式が異なるので相似でない。

問題 5.7 ★★

  1. A∈M⁡n(C)A \in \operatorname{M}_n(\mathbb{C}) が A3=AA^3 = A を満たすならば、AA は対角化可能で固有値は 0,1,−10, 1, -1 のいずれかであることを示せ。
  2. A∈M⁡n(R)A \in \operatorname{M}_n(\mathbb{R}) が A2=IA^2 = I を満たすならば、AA は R\mathbb{R} 上で対角化可能であり、tr⁡A\operatorname{tr} A は nn と偶奇が一致する −n-n 以上 nn 以下の整数であることを示せ。
解答
  1. t3−t=t(t−1)(t+1)t^3 - t = t(t - 1)(t + 1) は重根のない一次式の積なので、定理 5.28 より対角化可能。固有値は μA\mu_A の根で、μA∣t3−t\mu_A \mid t^3 - t より {0,±1}\lbrace 0, \pm 1 \rbrace に含まれる。
  2. t2−1=(t−1)(t+1)t^2 - 1 = (t - 1)(t + 1) より R\mathbb{R} 上対角化可能で、固有値 11 が pp 個、−1-1 が qq 個(p+q=np + q = n)並ぶ。tr⁡A=p−q=n−2q\operatorname{tr} A = p - q = n - 2q。

問題 5.8 ★★ K=CK = \mathbb{C} のとき、三角化(定理 5.19)を用いてケーリー–ハミルトンの定理を証明せよ。

解答

T=P−1APT = P^{-1}AP を上三角、対角成分を λ1,…,λn\lambda_1, \dots, \lambda_n とすると ΦA(t)=∏i(t−λi)\Phi_A(t) = \prod_i (t - \lambda_i) で、ΦA(A)=PΦA(T)P−1\Phi_A(A) = P\Phi_A(T)P^{-1} なので ΦA(T)=(T−λ1I)⋯(T−λnI)=O\Phi_A(T) = (T - \lambda_1 I) \cdots (T - \lambda_n I) = O を示せばよい。Wk=span⁡(e1,…,ek)W_k = \operatorname{span}(e_1, \dots, e_k)(W0={0}W_0 = \lbrace 0 \rbrace)とおく。TT は上三角なので T−λkIT - \lambda_k I は WkW_k を WkW_k に写し、しかも (T−λkI)ek(T - \lambda_k I)e_k の第 kk 成分は tkk−λk=0t_{kk} - \lambda_k = 0 なので (T−λkI)Wk⊂Wk−1(T - \lambda_k I)W_k \subset W_{k-1} である。因子は可換なので ΦA(T)=(T−λ1I)(T−λ2I)⋯(T−λnI)\Phi_A(T) = (T - \lambda_1 I)(T - \lambda_2 I) \cdots (T - \lambda_n I) を右から順に Wn=CnW_n = \mathbb{C}^n に施すと、Wn→Wn−1→⋯→W0={0}W_n \to W_{n-1} \to \cdots \to W_0 = \lbrace 0 \rbrace となり、ΦA(T)=O\Phi_A(T) = O。

問題 5.9 ★★★ ff が対角化可能で WW が ff-不変部分空間ならば、W=⨁λ(W∩V(λ))W = \bigoplus_{\lambda} (W \cap V(\lambda)) であることを示せ。

解答

命題 5.31 より f∣Wf\vert_W は対角化可能なので、WW は f∣Wf\vert_W の固有空間の直和である。f∣Wf\vert_W の固有値 λ\lambda の固有空間は {w∈W∣f(w)=λw}=W∩V(λ)\lbrace w \in W \mid f(w) = \lambda w \rbrace = W \cap V(\lambda) である。ff の固有値で f∣Wf\vert_W の固有値でないものについては W∩V(λ)={0}W \cap V(\lambda) = \lbrace 0 \rbrace なので、和に加えても変わらない。

(別証:補題 5.27 の証明から、V(λi)V(\lambda_i) への射影 πi\pi_i は ff の多項式である。WW は ff-不変なので πi(W)⊂W\pi_i(W) \subset W、よって w=∑iπi(w)w = \sum_i \pi_i(w) の各成分は W∩V(λi)W \cap V(\lambda_i) に属する。)

問題 5.10 ★★★ A,B∈M⁡n(C)A, B \in \operatorname{M}_n(\mathbb{C}) が AB=BAAB = BA を満たすとする。

  1. AA と BB は共通の固有ベクトルをもつことを示せ。
  2. P−1APP^{-1}AP と P−1BPP^{-1}BP がともに上三角行列になる P∈GL⁡n(C)P \in \operatorname{GL}_n(\mathbb{C}) が存在することを示せ。
解答
  1. AA の固有値 λ\lambda をとる(代数学の基本定理)。v∈V(λ)v \in V(\lambda) なら A(Bv)=BAv=λBvA(Bv) = BAv = \lambda Bv なので V(λ)V(\lambda) は BB-不変である。B∣V(λ)B\vert_{V(\lambda)} は(V(λ)≠{0}V(\lambda) \neq \lbrace 0 \rbrace 上の複素線形変換なので)固有ベクトル vv をもち、これは AA の固有ベクトルでもある。
  2. nn に関する帰納法。1 の共通固有ベクトルを含む基底をとると、両者は同時に
(λa0A1),(μb0B1)\begin{pmatrix} \lambda & a \\ 0 & A_1 \end{pmatrix}, \qquad \begin{pmatrix} \mu & b \\ 0 & B_1 \end{pmatrix}

の形になる。積の右下ブロックを比べると A1B1=B1A1A_1B_1 = B_1A_1 なので、帰納法の仮定から A1,B1A_1, B_1 を同時に上三角化する QQ があり、定理 5.19 の証明と同様に PP を作ればよい。

この章を読み終えたら

「読了」にすると学習記録とロードマップに反映されます。演習の自己採点もお忘れなく。

この章の誤りを報告GitHub で見る