こんにちは、おぐえもん(@oguemon_com)です。
前回の記事では、ケーリー・ハミルトンの定理について、勘違いしやすいポイントを交えて解説しました。
今回は、ケーリー・ハミルトンの定理と同じく、三角化の概念を用いて成立の理由を追及することができるフロベニウスの定理を解説します!
フロベニウスの定理とは?
ある多項式に行列を代入してできた行列の固有値は、行列の固有値を同じ多項式に代入した時に得られる値だよって旨の定理です。
フロベニウスの定理
n次正方行列Aは、λ1,λ2,⋯,λnを固有値に持つとする。
行列Xの多項式
f(X)=a0Xn+a1Xn−1+⋯+an−1X+anE
に、Aを代入して得られる行列f(A)の固有値は、
f(λ1), f(λ2), ⋯, f(λn)
である。
定理が成り立つ理由
この定理は、三角化を利用することで成り立つ理由を確かめることができます。
まず、任意の正方行列は三角化できます。つまり次の式が成り立つ適当な行列Pを選ぶことができます。
P−1AP=λ10⋮0a12λ2⋮0……⋱…a1na2n⋮λn
ここで、λ1 〜 λnはAの固有値です。
さて、次は、P−1f(A)Pを変形します。
P−1f(A)P=P−1(a0An+a1An−1+⋯+anE)P=a0P−1AnP+a1P−1An−1P+⋯+anE=a0(P−1AP)n+a1(P−1AP)n−1+⋯+anE=f(P−1AP)
2行目から3行目にかけて、P−1AnP=(P−1AP)nを用いた変形がなされています。この等式は次のようなイメージで成立を理解できると思います。
P−1A3P=P−1AAAP=P−1APP−1APP−1AP=(P−1AP)(P−1AP)(P−1AP)=(P−1AP)3
ところで、三角行列をn乗したとき、その対角成分は、元の成分のn乗となります(もちろん右上の成分は全然違う値になります)。実際に計算すると成り立つことがわかります。
λ10⋮0a12λ2⋮0……⋱…a1na2n⋮λnn=λ1n0⋮0b12λ2n⋮0……⋱…b1nb2n⋮λnn
多項式は、行列の累乗にスカラー倍したものを色々足し合わせたものにすぎません。ですので、ある三角行列Wの対角成分の 1 つをwiとした時、Wを多項式に代入してできた行列f(W)における同位置の対角成分はf(wi)となります。
P−1APは、対角成分として各固有値λ1 〜 λnをもつ三角行列です。ですので、f(P−1AP)の対角成分は、多項式f(X)に各固有値を代入したf(λ1) 〜 f(λn)で構成されることがいえます。
f(P−1AP)=f(λ1)0⋮0c12f(λ2)⋮0……⋱…c1nc2n⋮f(λn)
ところで、上で導いた等式P−1f(A)P=f(P−1AP)から、P−1f(A)Pが三角行列であることがいえます。つまり、Pを用いて行列f(A)を三角化したことに他なりません。三角化に関する記事の中で扱いましたが、三角化により生み出した三角行列の対角成分は元の行列の固有値です。
ゆえに、f(A)の固有値は、f(λ1) 〜 f(λn)であることが言えました。
一緒に例題を計算しよう!
簡単な行列を計算してフロベニウスの定理が実際に成立することを確かめてみましょう。
例として、2 次正方行列Aと、多項式f(X)を次のように定めました。
A=[1−124]
f(X)=X2+3X−E
ちなみに、Aの固有値は、「2」と「3」です!(一応理由 ↓)
∣A−tE∣=(1−t)(4−t)−2∗(−1)=t2−5t+6=(t−2)(t−3)
1. 多項式に行列を入れてみる
多項式に行列Aを入れてできた行列f(A)を計算しましょう。
下準備としてA2を計算しておきました。
A2=[−1−51014]
それではf(A)を計算しましょう。
f(A)=A2+3A−E=[−1−51014]+3[1−124]−[1001]=[−1+3∗1−1−5+3∗(−1)+010+3∗2−014+3∗4−1]=[1−81625]
これで、f(A)が求まりました。
2. 固有値を比較する
行列f(A)の固有値を求めましょう。
∣f(A)−tE∣=(1−t)(25−t)−16∗(−8)=t2−26t+153=(t−9)(t−17)
よって、行列f(A)の固有値は「9」と「17」です!
ところで、Aの固有値「2」と「3」をそれぞれ多項式f(X)に入れたらどうなるのでしょうか。スカラーは 1 次の行列とみることができます。f(2)とf(3)はそれぞれ次のようになります。
f(2)f(3)=22+3∗2−1=9=32+3∗3−1=17
まあ不思議、この組み合わせは行列f(A)の固有値と一致するではありませんか。
以上から、f(A)の固有値が、行列Aの固有値を多項式f(X)に通した時の値であることが分かりました。
おわりに
今回はフロベニウスの定理を解説しました!
これでひとまず【固有値編】はおしまいです。お疲れ様でした!