大学・線形代数・ベクトル・行列

固有値・固有ベクトルと対角化の練習問題

固有方程式から固有値・固有ベクトルを求め、行列を対角化できる。

考え方

行列の変換で向きが変わらない特別な方向(固有ベクトル)と、その方向の伸び縮みの倍率(固有値)を求め、行列を対角化できるようになります。主成分分析、レコメンド、振動の解析、長期の推移の予測などの土台です。

固有ベクトル
Av=λvA\boldsymbol{v} = \lambda\boldsymbol{v}(v≠0\boldsymbol{v} \neq \boldsymbol{0})となるベクトル。変換しても同じ直線の上にある方向。
固有値
そのときの倍率 λ\lambda。
固有方程式
det⁡(A−λE)=0\det(A - \lambda E) = 0。固有値はこの方程式の解。
対角化
固有ベクトルを並べた行列 PP で、P−1APP^{-1}AP を対角行列(対角成分以外が 0)にすること。

行列 AA でベクトルを変換すると、ふつうは向きも長さも変わります。ところが特別な方向のベクトルは、向き(直線)が変わらず、長さだけが λ\lambda 倍になります。この方向が固有ベクトル、倍率が固有値です。

Av=λvA\boldsymbol{v} = \lambda\boldsymbol{v} は (A−λE)v=0(A - \lambda E)\boldsymbol{v} = \boldsymbol{0} と書けます。0\boldsymbol{0} でない解 v\boldsymbol{v} があるのは、A−λEA - \lambda E が平面をつぶすとき、つまり det⁡(A−λE)=0\det(A - \lambda E) = 0 のときです。2×2 なら λ\lambda の2次方程式になります。

det⁡(A−λE)\det(A - \lambda E) は、対角成分の両方から λ\lambda を引いた行列の行列式 (a−λ)(d−λ)−bc(a - \lambda)(d - \lambda) - bc です。展開すると λ2−(a+d)λ+(ad−bc)\lambda^2 - (a + d)\lambda + (ad - bc) になり、固有値の和は a+da + d、積は det⁡A\det A になります(検算に便利です)。

固有値が求まったら、(A−λE)v=0(A - \lambda E)\boldsymbol{v} = \boldsymbol{0} を解いて固有ベクトルを求めます。固有ベクトルを列に並べた PP を使うと P−1APP^{-1}AP が固有値を対角に並べた行列になり、AnA^n の計算などが簡単になります。

v = (1, 1) は A で (3, 3) に移り、同じ直線(点線)の上で3倍に伸びる。u = (1, 0) は (2, 1) に移り、向きが変わる123123OxyvAv = 3vuAu
A=(2112)A = \begin{pmatrix} 2 & 1 \\ 1 & 2 \end{pmatrix} で変換すると、v=(1,1)\boldsymbol{v} = (1, 1) は同じ直線(点線)上の (3,3)(3, 3) に移り、長さだけが 3 倍になります(固有値 3)。u=(1,0)\boldsymbol{u} = (1, 0) は (2,1)(2, 1) に移り、向きが変わります。

ポイント

Av=λvA\boldsymbol{v} = \lambda\boldsymbol{v}(v≠0\boldsymbol{v} \neq \boldsymbol{0})

固有方程式: det⁡(A−λE)=(a−λ)(d−λ)−bc=λ2−(a+d)λ+(ad−bc)=0\det(A - \lambda E) = (a - \lambda)(d - \lambda) - bc = \lambda^2 - (a + d)\lambda + (ad - bc) = 0

対角化: P−1AP=(λ100λ2)P^{-1}AP = \begin{pmatrix} \lambda_1 & 0 \\ 0 & \lambda_2 \end{pmatrix}(PP は固有ベクトルを列に並べた行列)

解き方の手順

  1. 対角成分の両方から λ\lambda を引いた行列の行列式 (a−λ)(d−λ)−bc(a - \lambda)(d - \lambda) - bc を作る。
  2. 展開して λ2−(a+d)λ+(ad−bc)=0\lambda^2 - (a + d)\lambda + (ad - bc) = 0 の形にし、因数分解して固有値を求める(λ=1,6\lambda = 1, 6 のように答える)。
  3. 固有値の和が a+da + d、積が ad−bcad - bc になっているか確かめる。
  4. 固有ベクトルは、(A−λE)v=0(A - \lambda E)\boldsymbol{v} = \boldsymbol{0} の1行目から v=(1,t)\boldsymbol{v} = (1, t) の形で求める。

例

A=(5222)A = \begin{pmatrix} 5 & 2 \\ 2 & 2 \end{pmatrix} の固有値と、大きいほうの固有値の固有ベクトルを求めてください。

  1. det⁡(A−λE)=(5−λ)(2−λ)−2⋅2=λ2−7λ+6\det(A - \lambda E) = (5 - \lambda)(2 - \lambda) - 2 \cdot 2 = \lambda^2 - 7\lambda + 6
  2. (λ−1)(λ−6)=0(\lambda - 1)(\lambda - 6) = 0 より λ=1,6\lambda = 1, 6(和 7=5+27 = 5 + 2、積 6=det⁡A6 = \det A で確かめ)
  3. λ=6\lambda = 6 の固有ベクトル: (5−6)v1+2v2=0(5 - 6)v_1 + 2v_2 = 0 より v2=12v1v_2 = \dfrac{1}{2}v_1。v=(1,12)\boldsymbol{v} = \left(1, \dfrac{1}{2}\right) の定数倍
  4. 確かめ: A(21)=(126)=6(21)A\begin{pmatrix} 2 \\ 1 \end{pmatrix} = \begin{pmatrix} 12 \\ 6 \end{pmatrix} = 6\begin{pmatrix} 2 \\ 1 \end{pmatrix}

よくある間違い

練習問題でこの間違い方をしたときは、その場でお伝えします。

  • 固有方程式の立て方の誤りdet(A − λI) の対角成分から λ を引き忘れる、符号を間違える。固有方程式は det(A − λI) = 0 です。対角成分すべてから λ を引きましょう。

見る

計算例

主成分分析やレコメンドでは、行列の固有値・固有ベクトルを使って、データのいちばん大きな傾向(よく広がっている方向)を取り出します。

問いA=(12−36)A = \begin{pmatrix}1 & 2 \\ -3 & 6\end{pmatrix} の固有値をすべて求めてください。

途中の式

  1. 固有値は固有方程式 det⁡(A−λE)=0\det(A - \lambda E) = 0 の解です(対角成分の両方から λ\lambda を引く)。
  2. det⁡(A−λE)=∣1−λ2−36−λ∣=(1−λ)(6−λ)−2⋅(−3)\det(A - \lambda E) = \begin{vmatrix} 1 - \lambda & 2 \\ -3 & 6 - \lambda \end{vmatrix} = (1 - \lambda)(6 - \lambda) - 2\cdot (-3)
  3. =λ2−7λ+12=(λ−4)(λ−3)= \lambda ^{2} - 7\lambda + 12 = (\lambda-4)(\lambda-3)
  4. よって λ=4,3\lambda = 4, 3
  5. λ=4\lambda = 4 の固有ベクトル: (A−4E)v=0(A - 4E)\boldsymbol{v} = \boldsymbol{0} の1行目 −3v1+2v2=0-3v_1 + 2v_2 = 0 より v=(1,32)\boldsymbol{v} = \left(1, \frac{3}{2}\right) の定数倍

答え固有値(例: λ = 2, 5)は λ=4,3\lambda = 4, 3

一緒にやる

小さく分けて、一緒に解く

1つずつ答えを確かめながら進みます。間違えても大丈夫です。

主成分分析やレコメンドでは、行列の固有値・固有ベクトルを使って、データのいちばん大きな傾向(よく広がっている方向)を取り出します。

A=(524−2)A = \begin{pmatrix}5 & 2 \\ 4 & -2\end{pmatrix} の固有値をすべて求めてください。

小さな問いに分けて、1つずつ一緒に進めます。

ステップ 1 / 3

固有方程式 det⁡(A−λE)=0\det(A - \lambda E) = 0 の左辺を展開した、λ\lambda の2次式は?

λ などの文字を使って式で答えます。累乗は λ^2、かけ算の記号は省略できます(2λ)。

読み取り:(入力すると、ここに数式で表示されます)

自分でやる

練習問題

数字は毎回変わります。2〜3問続けて解けたら、次の単元や元のレッスンに進みましょう。

問題を準備しています…

学び直しマップ

前後の単元

分からないところがあったら前提まで戻り、分かったら使う先へ進みましょう。

学び直しマップで、この単元のまわりを見る