前回の記事では,固有値と固有ベクトルが正方行列の対角化と密接に関わっていることを説明し
- 対角化を用いた正方行列$A$の冪$A^n$の計算
- 対角化可能であるための基本定理
を説明しました.
今回の記事では,この固有値と固有ベクトルの求め方を説明します.
結論から言えば
- 固有方程式から固有値を求める
- 連立方程式を解いて固有ベクトルを求める
という順番で固有値,固有ベクトルを求めることができます.
この記事では,
- 固有多項式と固有方程式とは何か
- 固有値と固有ベクトルの求め方
を説明します.
「線形代数学の基本」の一連の記事はこちら
【線形代数の初学者のための道案内|線形代数のイメージを知る】
・行列と列ベクトル 【線形代数1|線形代数は「多変数バージョンの比例」という話】 【線形代数2|行列の計算の基本!行列の積はなぜこうなる?】 【線形代数3|逆行列を考えると何が嬉しいのか?】 【線形代数4|連立1次方程式の掃き出し法と行列の基本変形】 【線形代数5|正則の条件を簡単に!基本変形と行列の積の話】 【線形代数6|行列のランクと,行列が逆行列をもつための条件】 【線形代数7|連立1次方程式が解をもつ条件と解の自由度】 【線形代数8|線形独立のイメージと線形独立であるための条件】
・行列式 【線形代数9|行列の正則性を判定できる行列式のイメージ】 【線形代数10|行列式を定義するための置換の性質を理解する】 【線形代数11|「行列式」は線形代数の要!定義と性質を解説】 【線形代数12|[余因子展開]と[行列式による正則条件]を解説】
・部分空間と基底 【線形代数13|列ベクトル空間の部分空間の定義と具体例】 【線形代数14|部分空間の基底を具体例から理解する(準備中)】 【線形代数15|部分空間の同型と部分空間の次元(準備中)】 【線形代数16|線形写像の像と核と次元定理(準備中)】 【線形代数17|部分空間の和空間と共通部分の空間(準備中)】
・固有値と固有ベクトル 【線形代数18|「固有値」「固有ベクトル」「対角化」とは?】←今の記事 【線形代数19|固有値と固有ベクトルは2ステップで求める!】 【線形代数20|固有値・固有ベクトルの基本性質のまとめ】 【線形代数21|固有空間はなぜ大切か?対角化の必要十分条件】
目次
オススメの入門書
以下は初学者向けのオススメの教科書です.この記事の著者も数学教室の集団授業で使っています.
手を動かしてまなぶ 線形代数
[藤岡敦 著/裳華房]
線形代数の入門書で,説明も非常に丁寧なので初学者にも読み進めやすい教科書です.
数学で初めて出会った概念で詰まった時には,具体例を考えることで理解できるようになることはよくあります.
特に線形代数は高校数学で扱ってきた数学よりも抽象度がやや増すので,いきなり抽象的に理解するよりも「具体例を理解→抽象化」という学び方が効果的です.
本書は具体例と例題が豊富で,実際に手を動かしながらイメージを掴んで抽象的に理解することを目指しています.
また,続巻も発行されていますが,この第1巻だけでも正方行列の対角化(固有値・固有ベクトル)まで学ぶことができます.
固有値と固有ベクトルの復習
念のため,前回の記事で説明した固有値と固有ベクトルの定義を確認します.
この記事では
- 実数成分の$n$次正方行列の集合を$\Mat_{n}(\R)$
- 複素数成分の$n$次正方行列の集合を$\Mat_{n}(\C)$
- 実数成分の$n$次ベクトルを$\R^n$
- 複素数成分の$n$次ベクトルを$\C^n$
と表します.
[固有値・固有ベクトル] $A\in\Mat_{n}(\C)$に対して,$\lambda\in\C$と$\m{a}\in\C^n\setminus\{\m{0}\}$が存在して$\lambda\m{a}=A\m{a}$が成り立つとき,$\lambda$を$A$の固有値,$\m{a}$を$A$の固有値$\lambda$に属する固有ベクトルという.
定義で零ベクトル$\m{0}$を除いているように「固有ベクトルは$\m{0}$でない」というのは忘れがちですが,実は大事なのでしっかり意識してください.
図形的には零ベクトルでない$\m{a}\in\C^n$を勝手にとってきたとき,$\m{a}$は$A\m{a}\in\C^n$と平行とは限りませんが,うまく$\m{a}$をとれば$A\m{a}$と$\m{a}$が平行になることがあります.
つまり,$\lambda\in\C$をうまくとって$A\m{a}=\lambda\m{a}$となるような$\m{a}$をとれることがあります.
このときの$\m{a}$を固有値ベクトル,伸び率$\lambda$を固有値というわけですね.
$A\m{a}=\m{0}$のときは$\lambda=0$で$A\m{a}=\lambda\m{a}$を満たすので,$0$は$A$の固有値となりますね.
2つのベクトルの片方が零ベクトルのときは平行とは言わないので,$A\m{a}=\m{0}$の場合は$A\m{a}$と$\m{a}$は平行ではありませんが,ここはイメージの説明なのでざっくり「平行」と説明している次第です.
【前回の記事:線形代数13|「固有値」「固有ベクトル」「対角化」とは?】
正方行列$A$に対して,$A$の冪$A^n$を直接計算するのは一般には非常に面倒です.そこで,$A$の冪$A^n$を簡単に計算する方法として対角化があります.前回の記事では,対角化を考えるために固有値・固有ベクトルが重要であることを説明しています.
固有値と固有ベクトルの求め方
冒頭でも説明したように,固有値と固有ベクトルは
- 固有方程式から固有値を求める
- 連立方程式を解いて固有ベクトルを求める
という順番で求めることができます.
固有多項式と固有方程式
結論から言えば,正方行列$A$の固有値は$A$の固有方程式の解に一致します.
ということで,まずは固有多項式と固有方程式の定義を述べます.
[固有多項式,固有方程式] $A\in\Mat_{n}(\C)$に対して
- $x$の多項式$|xI-A|$を$A$の固有多項式 (characteristic polynomial)
- $x$の方程式$|xI-A|=0$を$A$の固有方程式 (characteristic equation)
という.ただし,$I$は$n$次正方行列であり,$|xI-A|$は正方行列$xI-A$の行列式である.
具体例を考えると,固有多項式と固有方程式がどういうものかすぐに分かると思います.
次の正方行列の固有多項式,固有方程式を求めよ.また,固有方程式を解け.
- $A=\bmat{1&2\\2&1}$
- $B=\bmat{1&1&0\\2&0&0\\0&0&-1}$
(1) $A=\bmat{1&2\\2&1}$の固有多項式は
である.また,$A$の固有方程式は$(x-3)(x+1)=0$であり,解くと$x=-1,3$となる.
(2) $B=\bmat{1&1&0\\2&0&0\\0&0&-1}$の固有多項式は
である.また,$B$の固有方程式は$(x+1)^2(x-2)=0$であり,解くと$x=-1,2$となる.
$A$が$n$次正方行列なら,行列式として対角成分の積に$n$次の項が現れるので,固有多項式$|xI-A|$は$n$次式となりますね.
固有値の求め方
先程書いたように固有値と固有方程式の解は一致します.
[固有値] $\lambda\in\C$, $A\in\Mat_{n}(\C)$に対して,次は同値である.
- $\lambda$は$A$の固有値である.
- $\lambda$は固有方程式$|xI-A|=0$の解(固有多項式$|xI-A|$の根)である.
2つ目の条件は「$\lambda$は$|\lambda I-A|=0$を満たす」と書いても同じことですね.
[$(1)\Ra(2)$の証明] $A$の固有値$\lambda$に属する固有ベクトルを$\m{a}$とする.
このとき,$A\m{a}=\lambda\m{a}$が成り立つから,移項して
が成り立つ.もし$\lambda I-A$が正則なら,両辺に左から$(\lambda I-A)^{-1}$をかけて$\m{0}=\m{a}$となるが,$\m{a}$は固有ベクトルだから$\m{a}\neq\m{0}$なので矛盾する.
よって,$\lambda I-A$は正則でないから,$|\lambda I-A|=0$を満たす.
[$(2)\Ra(1)$の証明] $|\lambda I-A|=0$が成り立つなら,$\m{x}$の連立方程式$(\lambda I-A)\m{x}=\m{0}$は非自明解($\m{x}\neq\m{0}$なる解)$\m{a}$をもつ.
よって,$\m{a}$は非自明解だったから$\m{a}\neq\m{0}$であり,
が成り立つから,$\lambda$は$A$の固有値(で,$\m{a}$は$A$の固有値$\lambda$に属する固有ベクトル)である.
この定理より,具体的に先程の問題の行列について
- $A=\bmat{1&2\\2&1}$の固有値は,固有方程式の解と一致して$-1$, $3$
- $A=\bmat{1&2&0\\2&0&0\\0&0&-1}$の固有値は,固有方程式の解と一致して$-1$, $2$
となりますね.
固有ベクトルの求め方
$\lambda$が正方行列$A$の固有値と分かっていれば,零ベクトルでない$\m{v}$が固有値$\lambda$に属する固有ベクトルであることは
を満たすことに他なりません.
これは$\m{x}$の連立方程式$(A-\lambda I)\m{x}=\m{0}$の解が$\m{v}$であることと同じことですから,この連立方程式の解が固有ベクトルとなりますね.
正方行列$A=\bmat{1&2\\2&1}$の固有ベクトルを求めよ.
$\m{x}$の連立方程式$P\m{x}=\m{q}$は,拡大係数行列$[P,\m{q}]$に行基本変形を施すことで解けることを思い出しておきましょう.
中学校以来,連立1次方程式は消去法によって解くことがよくありましたが,消去法の式変形は線形代数では行列の行基本変形で考えることになり,この解法を掃き出し法といいます.この記事では,掃き出し法の考え方を具体例から考えています.
$A=\bmat{1&2\\2&1}$の固有値が$-1$ ,$3$であることは上で説明したので
- 固有値$-1$に属する固有ベクトル
- 固有値$3$に属する固有ベクトル
を求めれば良い.
(i) $\m{v}$が「$A$の固有値$-1$の固有ベクトル」であることと
を満たすことは同値である.すなわち,$\m{v}$の連立方程式$(A+I)\m{v}=\m{0}$の解が「$A$の固有値$1$の固有ベクトル」である.
$A+I$は行基本変形により
となるから,掃き出し法の考え方から$\m{x}=c\bmat{1\\-1}$ ($c\neq0$)が$(A+I)\m{v}=\m{0}$の解である.
よって,「$A$の固有値$-1$の固有ベクトル」は$c\bmat{1\\-1}$ ($c\neq0$)である.
(ii) $\m{v}$が「$A$の固有値$3$の固有ベクトル」であることと
を満たすことは同値である.すなわち,$\m{v}$の連立方程式$(A-3I)\m{v}=\m{0}$の解が「$A$の固有値$3$の固有ベクトル」である.
$A-3I$は行基本変形により
となるから,掃き出し法の考え方から$\m{x}=c\bmat{1\\1}$ ($c\neq0$)が$(A-3I)\m{v}=\m{0}$の解である.
よって,「$A$の固有値$3$の固有ベクトル」は$d\bmat{1\\1}$ ($d\neq0$)である.
この例からも分かることですが,
- 固有ベクトルは固有値ごとに異なる
- 固有ベクトルは1つではない($c$, $d$のとりかたによって変わる)
ということに注意してください.
2つ目の注意について$A$の固有値$-1$の固有ベクトルは$c\bmat{1\\-1}$と求まりましたが,これは$c=1$とした$\bmat{1\\-1}$や,$c=-2$とした$\bmat{-2\\2}$なども固有ベクトルです(実際に$A\m{v}=(-1)\m{v}$を満たすことを確かめてみてください).
また,$c\neq0$としているのは,固有ベクトルは定義で$\m{0}$でないとしているためです.
固有値・固有ベクトルの性質
前回の記事と今のこの記事で
- 固有値・固有ベクトルと対角化の関係
- 固有値・固有ベクトルの求め方
を説明しましたが,固有値・固有ベクトルの性質の説明はしていませんでした.
次の記事では,固有値・固有ベクトルの性質を説明します.
「線形代数学の基本」の一連の記事はこちら
【線形代数の初学者のための道案内|線形代数のイメージを知る】
・行列と列ベクトル 【線形代数1|線形代数は「多変数バージョンの比例」という話】 【線形代数2|行列の計算の基本!行列の積はなぜこうなる?】 【線形代数3|逆行列を考えると何が嬉しいのか?】 【線形代数4|連立1次方程式の掃き出し法と行列の基本変形】 【線形代数5|正則の条件を簡単に!基本変形と行列の積の話】 【線形代数6|行列のランクと,行列が逆行列をもつための条件】 【線形代数7|連立1次方程式が解をもつ条件と解の自由度】 【線形代数8|線形独立のイメージと線形独立であるための条件】
・行列式 【線形代数9|行列の正則性を判定できる行列式のイメージ】 【線形代数10|行列式を定義するための置換の性質を理解する】 【線形代数11|「行列式」は線形代数の要!定義と性質を解説】 【線形代数12|[余因子展開]と[行列式による正則条件]を解説】
・部分空間と基底 【線形代数13|列ベクトル空間の部分空間の定義と具体例】 【線形代数14|部分空間の基底を具体例から理解する(準備中)】 【線形代数15|部分空間の同型と部分空間の次元(準備中)】 【線形代数16|線形写像の像と核と次元定理(準備中)】 【線形代数17|部分空間の和空間と共通部分の空間(準備中)】
・固有値と固有ベクトル 【線形代数18|「固有値」「固有ベクトル」「対角化」とは?】 【線形代数19|固有値と固有ベクトルは2ステップで求める!】←次の記事 【線形代数20|固有値・固有ベクトルの基本性質のまとめ】 【線形代数21|固有空間はなぜ大切か?対角化の必要十分条件】
参考文献
以下は参考文献です.
手を動かしてまなぶ 線形代数
[藤岡敦 著/裳華房]
線形代数の入門書で,説明も非常に丁寧なので初学者にも読み進めやすい教科書です.
数学で初めて出会った概念で詰まった時には,具体例を考えることで理解できるようになることはよくあります.
特に線形代数は高校数学で扱ってきた数学よりも抽象度がやや増すので,いきなり抽象的に理解するよりも「具体例を理解→抽象化」という学び方が効果的です.
本書は具体例と例題が豊富で,実際に手を動かしながらイメージを掴んで抽象的に理解することを目指しています.
また,続巻も発行されていますが,この第1巻だけでも正方行列の対角化(固有値・固有ベクトル)まで学ぶことができます.
線型代数入門
[齋藤正彦 著/東京大学出版会]
線形代数の教科書として半世紀に渡って売れ続けている超ロングセラーの教科書です.
発刊されてから本書の内容の流れが線形代数の教科書のスタンダードとなったほど,日本の線形代数の指導にインパクトを与えた名著です.
その証拠に,著者の齋藤正彦氏は本書で日本数学会出版賞を受賞しています.
「線形代数をとりあえず使えるようにするための教科書」ではなく「線形代数を理解するための教科書」のため,論理的に非常に詳しく書かれているのが特徴です.
また,テキストのレベルとしては少なくとも理論系(特に数学系)の学部生であれば,確実に理解しておきたい程度のものとなっています.
なお,本書については,以下の記事で書評としてまとめています.
本書の目次・必要な知識・良い点と気になる点・オススメの使い方などをレビューしています.