logo

固有値と固有ベクトル 📂行列代数

固有値と固有ベクトル

定義1

n×nn\times n 行列 AAが与えられたとしよう。0\mathbf{0}でないn×1n\times 1列ベクトルx\mathbf{x}、そして定数λ\lambdaに対して、次の式を固有値方程式または固有値問題という。

Ax=λx \begin{equation} A \mathbf{x} = \lambda \mathbf{x} \end{equation}

与えられたAAに対して、上のように固有値方程式を満たすλ\lambdaAA固有値と言い、x\mathbf{x}λ\lambdaに対応するAA固有ベクトルという。

説明

上の定義はλR\lambda \in \mathbb{R}xRn\mathbf{x} \in \mathbb{R}^{n}の時だけでなく、λC\lambda \in \mathbb{C}xCn\mathbf{x} \in \mathbb{C}^{n}の時にもそのまま適用される。「0\mathbf{0}でない」という条件がついているのは、下の式から分かるように、x=0\mathbf{x} = \mathbf{0}ならば常に成り立つからだ。

A0=0=λ0 A \mathbf{0} = \mathbf{0} = \lambda \mathbf{0}

幾何学的な動機

ベクトルx\mathbf{x}を行列AAで変換したAxA \mathbf{x}x\mathbf{x}の方向が同じだとすると、何か実数λ\lambdaに対して

Ax=λx A \mathbf{x} = \lambda \mathbf{x}

が成り立つことになる。行列AAは本来、どんな方向の概念も持たないが、AAの固有ベクトルが存在するならば、AAが何か特有の方向を指していると言えるだろう。だから、このようなベクトルx\mathbf{x}固有ベクトルと呼ぶのだ。例えば、以下のような2×22\times 2行列を考えてみよう。

A=[6223] A = \begin{bmatrix} 6 & 2 \\ 2 & 3 \end{bmatrix}

すると、ベクトル[21]\begin{bmatrix} 2 \\ 1 \end{bmatrix}[6223]\begin{bmatrix} 6 & 2 \\ 2 & 3 \end{bmatrix}に変換された時、[147]\begin{bmatrix} 14 \\ 7 \end{bmatrix}となって方向が同じである。ここでベクトル[21]\begin{bmatrix} 2 \\ 1 \end{bmatrix}λ=7\lambda = 7を掛けると、ベクトルの長さも同じになり、固有値方程式

Ax=λx[6223][21]=7[21] \begin{align*} A \mathbf{x} &= \lambda \mathbf{x} \\ \begin{bmatrix} 6 & 2 \\ 2 & 3 \end{bmatrix} \begin{bmatrix} 2 \\ 1 \end{bmatrix} &= 7 \begin{bmatrix} 2 \\ 1 \end{bmatrix} \end{align*}

の形の等式を満たす。このような理由でλ=7\lambda=7固有値と呼ぶのだ。よく見ると、固有ベクトルは[21]\begin{bmatrix} 2 \\ 1 \end{bmatrix}を伸ばしたり縮めたりして無数に見つけることができるが、固有値は変わらないことが分かる。だから、[21]\begin{bmatrix} 2 \\ 1 \end{bmatrix}を固有値77に対応するAAの固有ベクトルと表現するのだ。

このように幾何学的に説明した議論を一般的に拡張すると、固有値は代数的に方程式Ax=λxA \mathbf{x} = \lambda \mathbf{x}を満たすλ\lambdaであり、固有ベクトルは与えられたλ\lambdaに対する方程式の非自明な解である。

固有値方程式の解法

固有値を求めることは、固有値方程式から始まる。(1)(1)の式を整理すると、次のようになる。

Ax=λx    Axλx=0    AxλIx=0    (AλI)x=0 \begin{align*} && A \mathbf{x} &= \lambda \mathbf{x} \\ \implies && A \mathbf{x} - \lambda \mathbf{x} &= \mathbf{0} \\ \implies && A \mathbf{x} - \lambda I \mathbf{x} &= \mathbf{0} \\ \implies && \left( A - \lambda I \right) \mathbf{x} &= \mathbf{0} \end{align*}

この時、固有ベクトルは条件x0\mathbf{x} \ne \mathbf{0}を満たさなければならない。上記線形システム0\mathbf{0}でない解を持つ同値条件は、(AλI)\left( A - \lambda I \right)の逆行列が存在しないことであり、これは次の式が成り立つことと同値である。

det(AλI)=0 \det (A -\lambda I) = 0

したがって、上の式を満たすλ\lambdaAAの固有値になる。上記の式をAA特性方程式と言い、det(AλI)\det (A -\lambda I)AAn×nn\times n行列の時、nn次の多項式になり、これを特性多項式と呼ぶ。

ちなみに、A+BA+Bの固有値はAABBの固有値の和と異なる場合があり、ABABの固有値もAABBの固有値の積と異なる場合がある。また、方程式の解として固有値を求めることから分かるように、必ずしも実数であるという保証は全くない。

固有値を求める

解の例として、再びA=[6223]A = \begin{bmatrix} 6 & 2 \\ 2 & 3 \end{bmatrix}を考えてみよう。AλI=[6λ223λ]A-\lambda I = \begin{bmatrix} 6 - \lambda & 2 \\ 2 & 3 - \lambda \end{bmatrix}なので、AAの特性方程式を解くと次のようになる。

det(AλI)=0    (6λ)(3λ)4=0    λ29λ+184=0    (λ2)(λ7)=0 \begin{align*} && \det (A - \lambda I) &= 0 \\ \implies && (6 - \lambda)(3 - \lambda) - 4 &= 0 \\ \implies && \lambda^2 - 9 \lambda + 18 - 4 &= 0 \\ \implies && (\lambda - 2)(\lambda - 7) &= 0 \end{align*}

したがって、AAの固有値はλ=2\lambda = 2λ=7\lambda = 7である。2277λ\lambdaに代入してみると、それぞれの固有値に対応する固有ベクトルを求めることができる。ここでは、λ=7\lambda = 7の場合のみ紹介する。

λ=7\lambda = 7に対応する固有ベクトルを求める

λ=7\lambda = 7(1)(1)に代入して整理すると、次のようになる。

[6223][x1x2]=7[x1x2]    [6x1+2x22x1+3x2]=[7x17x2]    [x1+2x22x14x2]=[00] \begin{align*} && \begin{bmatrix} 6 & 2 \\ 2 & 3 \end{bmatrix} \begin{bmatrix} x_{1} \\ x_{2} \end{bmatrix} &= 7\begin{bmatrix} x_{1} \\ x_{2} \end{bmatrix} \\ \implies && \begin{bmatrix} 6x_{1} + 2x_{2} \\ 2x_{1} + 3x_{2} \end{bmatrix} &= \begin{bmatrix} 7x_{1} \\ 7x_{2} \end{bmatrix} \\ \implies && \begin{bmatrix} -x_{1} + 2x_{2} \\ 2x_{1} - 4x_{2} \end{bmatrix} &= \begin{bmatrix} 0 \\ 0 \end{bmatrix} \end{align*}

これを解くと、次のようになる。

{x1+2x2=02x14x2=0 \left\{ \begin{align*} -x_{1} + 2x_{2} &= 0 \\ 2x_{1} - 4x_{2} &= 0 \end{align*} \right.

    x1=2x2 \implies x_{1} = 2x_{2}

したがって、00でないすべてのx2x_{2}に対して、ベクトル[2x2x2]\begin{bmatrix} 2x_{2} \\ x_{2} \end{bmatrix}λ=7\lambda = 7に対応する固有ベクトルになる。通常、最も単純な形または大きさが11になる単位ベクトルを選ぶ。x2=1x_{2} = 1を代入すると、以下の固有ベクトルを得る。

A=[21] A = \begin{bmatrix} 2 \\ 1 \end{bmatrix}

性質

  • 正の整数kkに対して、λ\lambdaが行列AAの固有値であり、x\mathbf{x}λ\lambdaに対応する固有ベクトルであれば、λk\lambda ^{k}AkA^{k}の固有値であり、x\mathbf{x}λk\lambda ^{k}に対応する固有ベクトルである。

  1. Howard Anton, Elementary Linear Algebra: Aplications Version (12th Edition, 2019), p291-292 ↩︎