こんにちは、おぐえもん(@oguemon_com)です。
前回の記事では、固有値と固有ベクトルとは何なのかを基礎から解説しました。今回は、固有値と固有ベクトルを手っ取り早く求める方法を扱います!
固有値問題とは
ある正方行列Aについて、Ax=λxを満たすようなλとxの組み合わせを求める問題、言い換えると、Aの固有値とそれに対する固有ベクトルを求める問題のことを固有値問題と呼びます。
固有値と固有ベクトルは行列や線形変換における重要な指標です。しかし、これをノーヒントで探すのは至難の業(というか無理ゲー)。そこで、賢い先人たちは知恵を絞って固有値と固有ベクトルを手取り早く探す(=固有値問題を解く)方法を編み出しました。
固有値と固有ベクトルの求め方
固有値問題を解く方法の 1 つが、固有方程式(特性方程式とも呼びます)というものを解く方法です。解き方は次の通り。
固有値と固有ベクトルの求め方
Step1. 固有方程式を解いて固有値を導く
固有方程式とは、λについての方程式
∣A−λE∣=0
のことです。左辺は、行列(A−λE)の行列式です。これの解λが複数個見つかった場合、その全てがAの固有値です。
Step2. 固有値に対する固有ベクトルを導く
固有方程式の解λの1つ1つに対して、それぞれ連立方程式
(A−λE)x=o
の非自明解(零ベクトル以外の解)を求めます。この非自明解xが固有ベクトルです。
固有ベクトルは、1 つの固有値に対して無数に現れます(ある解のスカラー倍もまた解となるから)。固有値ベクトルを求める問題に遭遇した時は、特別な指示がない限り、任意変数をスカラー倍した網羅的な解答を記述しましょう!どういうことなのか、その具体例は後述します。
固有方程式で固有値がわかるワケ
突然「固有方程式」なんてものを出しました。なぜこれで固有値が出せるのかを説明します。
まず、ある固有値λに対する固有ベクトルを求める問題は、
Ax=λx
を満たす「x (=o)」を求める問題に帰着します。
この式は、① 両辺の左から単位行列Eを掛けて、② 右辺を左辺へ移項することで、次の形に変形できます。(中辺は、左辺からxを括りました)
Ax−λEx=(A−λE)x=o
つまり、ある固有値λに対する固有ベクトルを求める問題は、連立方程式(A−λE)x=oの「非自明解(零ベクトル以外の解)」を求めること、とも言い換えられます。
そもそも固有ベクトルが存在するためには、この連立方程式が非自明解を持たなければなりません(自明解=零ベクトルは固有ベクトルとして扱われないから)。そのための条件は、行列(A−λE)の行列式が 0 であること、つまり逆行列を持たないことです。もし逆行列を持っていたら、連立方程式の左から逆行列を掛けることでx=oが導かれちゃいます。これでは零ベクトル以外の解がありません。
逆に言えば、「∣A−λE∣=0」を満たすλさえ選べば、Ax=λxは非自明解を持ちます。そして、その非自明解を導いたとき、λは固有値、非自明解xはλに対する固有ベクトルそのものになるのです。
実際に固有値問題を解いてみよう!
固有値問題を解く要領を掴むため、簡単な行列の固有値と固有ベクトルを実際に求めてみましょう。
ここでは、前回の記事でも登場した 2 次元の正方行列Aを使用します。
A=(5439)
Step1. 固有方程式を解く
まずは、固有方程式の左辺(固有多項式と呼びます)を整理しましょう。
∣A−λE∣=(5439)−λ(1001)=5−λ439−λ=(5−λ)(9−λ)−3∗4=(λ−3)(λ−11)
よって、固有方程式は次のような式となります。
(λ−3)(λ−11)=0
この解はλ=3,11です。よって、Aの固有値は「3」と「11」です。
Step2. 各固有値に対する固有ベクトルを導く
次は、求めた固有値に対する固有ベクトルを求めましょう。
固有値「3」に対する固有ベクトル
連立方程式(A−3E)x=oの非自明解を求めましょう。
これはつまり、次の式に相当します。
(2436)x=(00)
これの非自明解、つまり固有ベクトルは、以下の通りです。
x=k(−32)
※ ただし、kは 0 以外の任意の実数です。
ちなみに、このベクトルは前回の記事で固有ベクトルとして与えたものです
固有値「11」に対する固有ベクトル
連立方程式(A−11E)x=oの非自明解を求めましょう。
これはつまり、次の式に相当します。
(−643−2)x=(00)
これの非自明解、つまり固有ベクトルは、以下の通りです。
x=k(12)
※ただし、kは 0 以外の任意の実数です。
先述した通り、固有ベクトルの解は、基本的に0 以外の任意の実数倍の形となり、同じ向きの固有ベクトルが無数に生まれることになります。ちなみに、行列によっては、1 つの固有値に対する固有ベクトルが、ある複数のベクトルの 1 次結合の形になることもあります。
おわりに
今回は、固有値問題を解くための方法について説明しました!
次回の記事では、固有方程式の左辺である「固有多項式」を用いて、行列の対角成分の総和がもつ性質を明らかにしていきます。