第 2 章で $A^n$ を求めたときは、$A^2$,$A^3$ と地道に計算して規則を予想し、帰納法で確かめました。行き当たりばったりの方法です。この章では、**どんな 2 次の正方行列も $A^2 - (a + d)A + (ad - bc)E = O$ という関係を必ず満たす**というケーリー・ハミルトンの定理を手に入れます。この 1 本の式があれば $A^2$ を $A$ と $E$ の組み合わせに書き直せるので、$A^3$ も $A^{100}$ も $A^{2026}$ も、けっきょく $pA + qE$ という形に収まってしまいます。$p$ と $q$ を求めるだけの問題に化けるのです。ただし、この定理には**逆が成り立たない**という落とし穴があり、そこが入試でいちばん狙われます。
行列が満たす 2 次方程式
数の世界では、x x x に何か具体的な数を入れれば x 2 x^2 x 2 − 5 x {}- 5x − 5 x + 6 {}+ 6 + 6 の値が決まります。行列でも同じことができます。A A A を 2 次の正方行列とすると、A 2 A^2 A 2 も 5 A 5A 5 A も 6 E 6E 6 E も 2 次の正方行列なので、A 2 A^2 A 2 − 5 A {}- 5A − 5 A + 6 E {}+ 6E + 6 E という式には意味があります。定数項を 6 6 6 ではなく 6 E 6E 6 E と書くのは、行列と数は足せないからです。数の 6 6 6 を行列の仲間に入れるには、単位行列を使って 6 E 6E 6 E の姿にしてやる必要がありました(第 2 章)。
さて、A A A = ( 3 1 2 4 ) {}= \begin{pmatrix} 3 & 1 \\ 2 & 4 \end{pmatrix} = ( 3 2 1 4 ) で試してみます。a a a + d {}+ d + d = 7 {}= 7 = 7 ,a d ad a d − b c {}- bc − b c = 12 {}= 12 = 12 − 2 {}- 2 − 2 = 10 {}= 10 = 10 です。この 7 7 7 と 10 10 10 を係数にして A 2 A^2 A 2 − 7 A {}- 7A − 7 A + 10 E {}+ 10E + 10 E を計算すると
( 11 7 14 18 ) \displaystyle \begin{pmatrix} 11 & 7 \\ 14 & 18 \end{pmatrix} ( 11 14 7 18 ) − ( 21 7 14 28 ) \displaystyle {}- \begin{pmatrix} 21 & 7 \\ 14 & 28 \end{pmatrix} − ( 21 14 7 28 ) + ( 10 0 0 10 ) \displaystyle {}+ \begin{pmatrix} 10 & 0 \\ 0 & 10 \end{pmatrix} + ( 10 0 0 10 ) = ( 0 0 0 0 ) \displaystyle {}= \begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix} = ( 0 0 0 0 )
きれいに零行列になりました。偶然ではありません。どんな 2 次の正方行列でも 、対角成分の和 a a a + d {}+ d + d と行列式 Δ \Delta Δ = a d {}= ad = a d − b c {}- bc − b c を係数に選べば、必ずこうなります。
証明は成分を書き下すだけです。
A 2 \displaystyle A^2 A 2 = ( a b c d ) ( a b c d ) \displaystyle {}= \begin{pmatrix} a & b \\ c & d \end{pmatrix}\begin{pmatrix} a & b \\ c & d \end{pmatrix} = ( a c b d ) ( a c b d ) = ( a 2 + b c b ( a + d ) c ( a + d ) b c + d 2 ) \displaystyle {}= \begin{pmatrix} a^2 + bc & b(a + d) \\ c(a + d) & bc + d^2 \end{pmatrix} = ( a 2 + b c c ( a + d ) b ( a + d ) b c + d 2 )
積の計算で、右上は a b ab ab + b d {}+ bd + b d = b ( a + d ) {}= b(a+d) = b ( a + d ) 、 左下は c a ca c a + d c {}+ dc + d c = c ( a + d ) {}= c(a+d) = c ( a + d ) とくくれるのがみそです。ここから ( a + d ) A (a+d)A ( a + d ) A を引きます。非対角成分は差し引きゼロ、対角成分は
a 2 \displaystyle a^2 a 2 + b c \displaystyle {}+ bc + b c − ( a 2 + a d ) \displaystyle {}- (a^2 + ad) − ( a 2 + a d ) = b c \displaystyle {}= bc = b c − a d , \displaystyle {}- ad, − a d , b c \displaystyle bc b c + d 2 \displaystyle {}+ d^2 + d 2 − ( a d + d 2 ) \displaystyle {}- (ad + d^2) − ( a d + d 2 ) = b c \displaystyle {}= bc = b c − a d \displaystyle {}- ad − a d
なので
A 2 \displaystyle A^2 A 2 − ( a + d ) A \displaystyle {}- (a + d)A − ( a + d ) A = ( b c − a d 0 0 b c − a d ) \displaystyle {}= \begin{pmatrix} bc - ad & 0 \\ 0 & bc - ad \end{pmatrix} = ( b c − a d 0 0 b c − a d )
対角成分はどちらも b c bc b c − a d {}- ad − a d = − Δ {}= -\Delta = − Δ 。 つまりこれは − Δ E -\Delta E − Δ E です。移項すれば公式1 が出ます。
この定理のありがたみは、A 2 A^2 A 2 を「A A A の何倍か」と「E E E の何倍か」の足し算に書き直せる ところにあります。2 次以上の項を 1 次以下に引きずり下ろせる、ということです。これを次数下げ といいます。
似た経験は数と式の章でしています。α \alpha α = 1 {}= 1 = 1 + 2 {}+ \sqrt{2} + 2 は x 2 x^2 x 2 − 2 x {}- 2x − 2 x − 1 {}- 1 − 1 = 0 {}= 0 = 0 の解なので α 2 \alpha^2 α 2 = 2 α {}= 2\alpha = 2 α + 1 {}+ 1 + 1 でした。だから α 3 \alpha^3 α 3 も α 10 \alpha^{10} α 10 も、計算すれば必ず「整数 + + + 整数 × 2 \times \sqrt{2} × 2 」 の形に戻ってきます。2 \sqrt{2} 2 が満たす 2 次方程式が、高い累乗を低いところへ押し戻しているのです。行列にも同じ仕組みが備わっている、というのがこの定理の中身です。
どんな 2 次の正方行列も A 2 A^2 A 2 = ( a + d ) A {}= (a + d)A = ( a + d ) A − Δ E {}- \Delta E − Δ E という自分専用の関係式を持っていて、これを使えば A 2 A^2 A 2 を A A A と E E E の組み合わせに書き直せるということです。
A A A = ( 3 1 2 4 ) {}= \begin{pmatrix} 3 & 1 \\ 2 & 4 \end{pmatrix} = ( 3 2 1 4 ) について、次の問いに答えなさい。
(1) ケーリー・ハミルトンの定理の式を書きなさい。
(2) (1) を使って A 3 A^3 A 3 を求めなさい。
【解答】
(1) a a a + d {}+ d + d = 3 {}= 3 = 3 + 4 {}+ 4 + 4 = 7 {}= 7 = 7 ,Δ \Delta Δ = 3 × 4 {}= 3 \times 4 = 3 × 4 − 1 × 2 {}- 1 \times 2 − 1 × 2 = 10 {}= 10 = 10 です。よって
A 2 ‾ \displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}A^2} A 2 − 7 A ‾ \displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}{}- 7A} − 7 A + 10 E = O ‾ \displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}{}+ 10E = O} + 10 E = O (2) (1) より A 2 A^2 A 2 = 7 A {}= 7A = 7 A − 10 E {}- 10E − 10 E です。両辺に右から A A A をかけます。
A 3 \displaystyle A^3 A 3 = 7 A 2 \displaystyle {}= 7A^2 = 7 A 2 − 10 A \displaystyle {}- 10A − 10 A ここで右辺の A 2 A^2 A 2 に、もう一度 7 A 7A 7 A − 10 E {}- 10E − 10 E を代入します。
A 3 \displaystyle A^3 A 3 = 7 ( 7 A − 10 E ) \displaystyle {}= 7(7A - 10E) = 7 ( 7 A − 10 E ) − 10 A \displaystyle {}- 10A − 10 A = 49 A \displaystyle {}= 49A = 49 A − 70 E \displaystyle {}- 70E − 70 E − 10 A \displaystyle {}- 10A − 10 A = 39 A \displaystyle {}= 39A = 39 A − 70 E \displaystyle {}- 70E − 70 E A A A と E E E の組み合わせになりました。あとは成分を入れるだけです。
A 3 \displaystyle A^3 A 3 = 39 ( 3 1 2 4 ) \displaystyle {}= 39\begin{pmatrix} 3 & 1 \\ 2 & 4 \end{pmatrix} = 39 ( 3 2 1 4 ) − 70 ( 1 0 0 1 ) \displaystyle {}- 70\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} − 70 ( 1 0 0 1 ) = ( 117 − 70 39 78 156 − 70 ) \displaystyle {}= \begin{pmatrix} 117 - 70 & 39 \\ 78 & 156 - 70 \end{pmatrix} = ( 117 − 70 78 39 156 − 70 ) = ( 47 39 78 86 ) ‾ \displaystyle {}= \underline{\begin{pmatrix} 47 & 39 \\ 78 & 86 \end{pmatrix}} = ( 47 78 39 86 ) 素直に A 2 A^2 A 2 を出してから A A A をかけても同じ答えになりますが、こちらは行列のかけ算を一度もしていない ことに注目してください。やったのは数の計算だけです。A 10 A^{10} A 10 や A 100 A^{100} A 100 になるほど、この差は大きくなります。
定理から逆行列も出てくる
公式1 をもう少し眺めてみます。A 2 A^2 A 2 − ( a + d ) A {}- (a+d)A − ( a + d ) A = − Δ E {}= -\Delta E = − Δ E の左辺は A A A でくくれて
A { A − ( a + d ) E } \displaystyle A\{A - (a + d)E\} A { A − ( a + d ) E } = − Δ E \displaystyle {}= -\Delta E = − Δ E
と書けます。両辺を − Δ -\Delta − Δ で割れば(Δ \Delta Δ ≠ 0 {}\neq 0 = 0 のとき)
A ⋅ 1 Δ { ( a + d ) E − A } \displaystyle A \cdot \frac{1}{\Delta}\{(a + d)E - A\} A ⋅ Δ 1 {( a + d ) E − A } = E \displaystyle {}= E = E
です。A A A にかけて E E E になる相手が見つかった、つまりこれが逆行列です。
第 3 章で覚えた A − 1 A^{-1} A − 1 = 1 Δ ( d − b − c a ) {}= \dfrac{1}{\Delta}\begin{pmatrix} d & -b \\ -c & a \end{pmatrix} = Δ 1 ( d − c − b a ) と同じものです。実際
( a + d ) E \displaystyle (a + d)E ( a + d ) E − A \displaystyle {}- A − A = ( a + d 0 0 a + d ) \displaystyle {}= \begin{pmatrix} a + d & 0 \\ 0 & a + d \end{pmatrix} = ( a + d 0 0 a + d ) − ( a b c d ) \displaystyle {}- \begin{pmatrix} a & b \\ c & d \end{pmatrix} − ( a c b d ) = ( d − b − c a ) \displaystyle {}= \begin{pmatrix} d & -b \\ -c & a \end{pmatrix} = ( d − c − b a )
で、第 3 章 公式2 の逆行列の式に出てきた行列そのものが現れます。「a a a と d d d を入れかえ、b b b と c c c の符号を変える」というあの操作の正体は、( a + d ) E (a+d)E ( a + d ) E から A A A を引くことだったわけです。
例題1 の A A A = ( 3 1 2 4 ) {}= \begin{pmatrix} 3 & 1 \\ 2 & 4 \end{pmatrix} = ( 3 2 1 4 ) で確かめます。Δ \Delta Δ = 10 {}= 10 = 10 ,a a a + d {}+ d + d = 7 {}= 7 = 7 なので
A − 1 \displaystyle A^{-1} A − 1 = 1 10 { 7 ( 1 0 0 1 ) − ( 3 1 2 4 ) } \displaystyle {}= \frac{1}{10}\left\{7\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} - \begin{pmatrix} 3 & 1 \\ 2 & 4 \end{pmatrix}\right\} = 10 1 { 7 ( 1 0 0 1 ) − ( 3 2 1 4 ) } = 1 10 ( 4 − 1 − 2 3 ) \displaystyle {}= \frac{1}{10}\begin{pmatrix} 4 & -1 \\ -2 & 3 \end{pmatrix} = 10 1 ( 4 − 2 − 1 3 )
第 3 章の公式で求めても同じです。おかげで A − 1 A^{-1} A − 1 もまた「A A A の何倍か + + + E E E の何倍か」の形に書けることがわかりました。負の指数まで含めて、A A A の累乗はすべて p A pA p A + q E {}+ qE + q E の形に収まる、ということです。
定理の逆は成り立たない
ここがこの章でいちばん間違えやすいところです。公式1 は「A A A を決めれば a a a + d {}+ d + d と Δ \Delta Δ が決まり、その係数で式が成り立つ」と言っています。では逆に、A 2 A^2 A 2 − 5 A {}- 5A − 5 A + 6 E {}+ 6E + 6 E = O {}= O = O という式が成り立っていたら、a a a + d {}+ d + d = 5 {}= 5 = 5 ,Δ \Delta Δ = 6 {}= 6 = 6 だと言えるでしょうか。
言えません。 反例はすぐ作れます。A A A = 2 E {}= 2E = 2 E = ( 2 0 0 2 ) {}= \begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix} = ( 2 0 0 2 ) とすると A 2 A^2 A 2 = 4 E {}= 4E = 4 E なので
A 2 \displaystyle A^2 A 2 − 5 A \displaystyle {}- 5A − 5 A + 6 E \displaystyle {}+ 6E + 6 E = 4 E \displaystyle {}= 4E = 4 E − 10 E \displaystyle {}- 10E − 10 E + 6 E \displaystyle {}+ 6E + 6 E = O \displaystyle {}= O = O
たしかに成り立ちますが、この A A A の対角和は 2 2 2 + 2 {}+ 2 + 2 = 4 {}= 4 = 4 で、5 5 5 ではありません。Δ \Delta Δ = 4 {}= 4 = 4 で、6 6 6 でもありません。
種明かしをすると、A A A = k E {}= kE = k E の形(数の化身 、第 2 章)だけが特別なのです。A A A = k E {}= kE = k E なら A 2 A^2 A 2 − p A {}- pA − p A + q E {}+ qE + q E = ( k 2 − p k + q ) E {}= (k^2 - pk + q)E = ( k 2 − p k + q ) E なので、k k k が 2 次方程式 x 2 x^2 x 2 − p x {}- px − p x + q {}+ q + q = 0 {}= 0 = 0 の解でありさえすれば、p p p と q q q がどんな値でも式が成り立ってしまいます。A A A 自身の対角和とも行列式とも無関係に、いくらでも式が作れるわけです。
答案では「A A A ≠ k E {}\neq kE = k E のとき」と「A A A = k E {}= kE = k E のとき」に分けて書きます。どちらか一方しか書かないと減点されるのが、この手の問題のお約束です。
A 2 A^2 A 2 − p A {}- pA − p A + q E {}+ qE + q E = O {}= O = O という式を見ても、すぐに p p p = a {}= a = a + d {}+ d + d とはいえず、A A A が単位行列の実数倍でないことを確かめてから初めていえるということです。
例題2:定理の逆を使う場合分け
2 次の正方行列 A A A が A 2 A^2 A 2 − 3 A {}- 3A − 3 A + 2 E {}+ 2E + 2 E = O {}= O = O を満たすとき、A A A の対角和 a a a + d {}+ d + d と行列式 Δ \Delta Δ の組をすべて求めなさい。
【解答】
(i) A A A が k E kE k E の形でないとき。 公式3 より a ‾ \underline{\rule[-0.1944em]{0em}{0.8888em}a} a + d = 3 , ‾ \underline{\rule[-0.1944em]{0em}{0.8888em}{}+ d = 3,} + d = 3 , Δ = 2 ‾ \underline{\rule[-0.1944em]{0em}{0.8888em}\ \Delta = 2} Δ = 2 です。
(ii) A A A = k E {}= kE = k E のとき。 A 2 A^2 A 2 = k 2 E {}= k^2E = k 2 E なので
( k 2 − 3 k + 2 ) E \displaystyle (k^2 - 3k + 2)E ( k 2 − 3 k + 2 ) E = O \displaystyle {}= O = O ⟹ k 2 \displaystyle {}\Longrightarrow \quad k^2 ⟹ k 2 − 3 k \displaystyle {}- 3k − 3 k + 2 \displaystyle {}+ 2 + 2 = 0 \displaystyle {}= 0 = 0 ⟹ k \displaystyle {}\Longrightarrow \quad k ⟹ k = 1 , \displaystyle {}= 1, = 1 , 2 \displaystyle \ 2 2 k k k = 1 {}= 1 = 1 のとき A A A = E {}= E = E で、a a a + d {}+ d + d = 2 {}= 2 = 2 ,Δ \Delta Δ = 1 {}= 1 = 1 。
k k k = 2 {}= 2 = 2 のとき A A A = 2 E {}= 2E = 2 E で、a a a + d {}+ d + d = 4 {}= 4 = 4 ,Δ \Delta Δ = 4 {}= 4 = 4 。
以上より、組は ( a + d , Δ ) = ( 3 , 2 ) , ‾ \underline{\rule[-0.25em]{0em}{1.0000em}(a+d,\ \Delta) = (3,\ 2),} ( a + d , Δ ) = ( 3 , 2 ) , ( 2 , 1 ) , ‾ \underline{\rule[-0.25em]{0em}{1.0000em}\ (2,\ 1),} ( 2 , 1 ) , ( 4 , 4 ) ‾ \underline{\rule[-0.25em]{0em}{1.0000em}\ (4,\ 4)} ( 4 , 4 ) の 3 通りです。
(ii) で出てきた k k k = 1 , {}= 1, = 1 , 2 \ 2 2 が、x 2 x^2 x 2 − 3 x {}- 3x − 3 x + 2 {}+ 2 + 2 = 0 {}= 0 = 0 の解になっていることに注意してください。A A A = k E {}= kE = k E のときは、行列の式がまるごと数の 2 次方程式に化けてしまうのです。行列と数の境目がいちばん薄くなる場所が、k E kE k E という形だといえます。
次数下げで高い累乗を下ろす
道具がそろったので、本題の A n A^n A n に入ります。公式1 を繰り返し使えば、A 3 A^3 A 3 も A 4 A^4 A 4 も p A pA p A + q E {}+ qE + q E の形に落ちていきます。
A A A の累乗は、何乗であっても「A A A の実数倍 + + + E E E の実数倍」という 1 次の形に必ず収まるということです。
例題3:次数下げで
A 3 A^3 A 3 ,A 4 A^4 A 4 を求める
A A A = ( 3 − 1 2 0 ) {}= \begin{pmatrix} 3 & -1 \\ 2 & 0 \end{pmatrix} = ( 3 2 − 1 0 ) について、A 3 A^3 A 3 と A 4 A^4 A 4 を求めなさい。
【解答】
a a a + d {}+ d + d = 3 {}= 3 = 3 + 0 {}+ 0 + 0 = 3 {}= 3 = 3 ,Δ \Delta Δ = 3 × 0 {}= 3 \times 0 = 3 × 0 − ( − 1 ) × 2 {}- (-1) \times 2 − ( − 1 ) × 2 = 2 {}= 2 = 2 なので
A 2 \displaystyle A^2 A 2 = 3 A \displaystyle {}= 3A = 3 A − 2 E \displaystyle {}- 2E − 2 E 両辺に A A A をかけ、出てきた A 2 A^2 A 2 をまた下ろします。
A 3 \displaystyle A^3 A 3 = 3 A 2 \displaystyle {}= 3A^2 = 3 A 2 − 2 A \displaystyle {}- 2A − 2 A = 3 ( 3 A − 2 E ) \displaystyle {}= 3(3A - 2E) = 3 ( 3 A − 2 E ) − 2 A \displaystyle {}- 2A − 2 A = 7 A \displaystyle {}= 7A = 7 A − 6 E \displaystyle {}- 6E − 6 E A 4 \displaystyle A^4 A 4 = 7 A 2 \displaystyle {}= 7A^2 = 7 A 2 − 6 A \displaystyle {}- 6A − 6 A = 7 ( 3 A − 2 E ) \displaystyle {}= 7(3A - 2E) = 7 ( 3 A − 2 E ) − 6 A \displaystyle {}- 6A − 6 A = 15 A \displaystyle {}= 15A = 15 A − 14 E \displaystyle {}- 14E − 14 E 成分に直します。
A 3 \displaystyle A^3 A 3 = 7 ( 3 − 1 2 0 ) \displaystyle {}= 7\begin{pmatrix} 3 & -1 \\ 2 & 0 \end{pmatrix} = 7 ( 3 2 − 1 0 ) − 6 E \displaystyle {}- 6E − 6 E = ( 15 − 7 14 − 6 ) ‾ \displaystyle {}= \underline{\begin{pmatrix} 15 & -7 \\ 14 & -6 \end{pmatrix}} = ( 15 14 − 7 − 6 ) A 4 \displaystyle A^4 A 4 = 15 ( 3 − 1 2 0 ) \displaystyle {}= 15\begin{pmatrix} 3 & -1 \\ 2 & 0 \end{pmatrix} = 15 ( 3 2 − 1 0 ) − 14 E \displaystyle {}- 14E − 14 E = ( 31 − 15 30 − 14 ) ‾ \displaystyle {}= \underline{\begin{pmatrix} 31 & -15 \\ 30 & -14 \end{pmatrix}} = ( 31 30 − 15 − 14 ) 係数の並びを見ておきましょう。A 2 A^2 A 2 から順に ( 3 , − 2 ) (3,\ -2) ( 3 , − 2 ) ,( 7 , − 6 ) (7,\ -6) ( 7 , − 6 ) ,( 15 , − 14 ) (15,\ -14) ( 15 , − 14 ) 。A A A の係数は 3 , 3, 3 , 7 , \ 7, 7 , 15 \ 15 15 、 つまり 2 n 2^n 2 n − 1 {}- 1 − 1 です。E E E の係数は − 2 , -2, − 2 , − 6 , {}\ -6, − 6 , − 14 {}\ -14 − 14 で − ( 2 n − 2 ) -(2^n - 2) − ( 2 n − 2 ) 。 予想するなら
A n \displaystyle A^n A n = ( 2 n − 1 ) A \displaystyle {}= (2^n - 1)A = ( 2 n − 1 ) A − ( 2 n − 2 ) E \displaystyle {}- (2^n - 2)E − ( 2 n − 2 ) E ですが、この予想を証明するには帰納法が要ります。次の節では、予想も帰納法も使わずに p p p と q q q を直接求める方法を身につけます。
割り算の余りで一気に求める
例題3 の手順は確実ですが、n n n = 100 {}= 100 = 100 では 99 回繰り返すことになります。もっと直接に p p p と q q q を出したい。そこで、多項式の割り算を使います。
x n x^n x n を 2 次式 x 2 x^2 x 2 − ( a + d ) x {}- (a+d)x − ( a + d ) x + Δ {}+ \Delta + Δ で割ったときの商を Q ( x ) Q(x) Q ( x ) 、 余りを p x px p x + q {}+ q + q (1 次以下)とすると
x n \displaystyle x^n x n = { x 2 − ( a + d ) x + Δ } Q ( x ) \displaystyle {}= \{x^2 - (a + d)x + \Delta\}Q(x) = { x 2 − ( a + d ) x + Δ } Q ( x ) + p x \displaystyle {}+ px + p x + q \displaystyle {}+ q + q
この x x x に行列 A A A を代入します。ここで効くのが公式1 です。{ } \{\ \} { } の中は A A A を入れると O O O になるので、第 1 項がまるごと消えて
A n \displaystyle A^n A n = p A \displaystyle {}= pA = p A + q E \displaystyle {}+ qE + q E
が残ります。つまり余りの係数がそのまま答え です。しかも p p p ,q q q を求めるのは、数の計算だけでできます。x 2 x^2 x 2 − ( a + d ) x {}- (a+d)x − ( a + d ) x + Δ {}+ \Delta + Δ = 0 {}= 0 = 0 の解を代入すれば、恒等式の第 1 項が 0 0 0 になるからです。
整数の世界でも同じことをしています。7 100 7^{100} 7 100 を 5 5 5 で割った余りを知りたいとき、7 100 7^{100} 7 100 を全部書き下す人はいません。7 7 7 ≡ 2 ( m o d 5 ) {}\equiv 2 \pmod 5 ≡ 2 ( mod 5 ) から余りだけを追いかけます(整数 第 5 章)。ここでも、知りたいのは x n x^n x n そのものではなく、x 2 x^2 x 2 − ( a + d ) x {}- (a+d)x − ( a + d ) x + Δ {}+ \Delta + Δ で割った余りだけ。全体を計算せず、必要な部分だけ取り出す発想です。
x 2 x^2 x 2 − ( a + d ) x {}- (a+d)x − ( a + d ) x + Δ {}+ \Delta + Δ = 0 {}= 0 = 0 という方程式は、この章のあいだ何度も顔を出します。A A A が背負っている 2 次方程式そのものです。
A n A^n A n を求めるには、x n x^n x n を A A A の背負う 2 次式で割った余り p x px p x + q {}+ q + q を求めればよく、その p p p ,q q q は 2 次方程式の解を代入した連立方程式から出るということです。
例題4:
A n A^n A n を求める(解が 2 つあるとき)
A A A = ( 1 2 3 2 ) {}= \begin{pmatrix} 1 & 2 \\ 3 & 2 \end{pmatrix} = ( 1 3 2 2 ) について、A n A^n A n を求めなさい。ただし n n n は正の整数とする。
【解答】
a a a + d {}+ d + d = 1 {}= 1 = 1 + 2 {}+ 2 + 2 = 3 {}= 3 = 3 ,Δ \Delta Δ = 1 × 2 {}= 1 \times 2 = 1 × 2 − 2 × 3 {}- 2 \times 3 − 2 × 3 = − 4 {}= -4 = − 4 なので、A A A が背負う 2 次方程式は
x 2 \displaystyle x^2 x 2 − 3 x \displaystyle {}- 3x − 3 x − 4 \displaystyle {}- 4 − 4 = 0 \displaystyle {}= 0 = 0 ⟹ ( x − 4 ) ( x + 1 ) \displaystyle {}\Longrightarrow \quad (x - 4)(x + 1) ⟹ ( x − 4 ) ( x + 1 ) = 0 \displaystyle {}= 0 = 0 ⟹ x \displaystyle {}\Longrightarrow \quad x ⟹ x = 4 , \displaystyle {}= 4, = 4 , − 1 \displaystyle {}\ -1 − 1 x n x^n x n を x 2 x^2 x 2 − 3 x {}- 3x − 3 x − 4 {}- 4 − 4 で割った商を Q ( x ) Q(x) Q ( x ) 、 余りを p x px p x + q {}+ q + q とします。
x n \displaystyle x^n x n = ( x 2 − 3 x − 4 ) Q ( x ) \displaystyle {}= (x^2 - 3x - 4)Q(x) = ( x 2 − 3 x − 4 ) Q ( x ) + p x \displaystyle {}+ px + p x + q \displaystyle {}+ q + q x x x = 4 {}= 4 = 4 と x x x = − 1 {}= -1 = − 1 を代入すると、どちらも第 1 項が 0 0 0 になります。
{ 4 n = 4 p + q ( − 1 ) n = − p + q \begin{cases} 4^n = 4p + q \\ (-1)^n = -p + q \end{cases} { 4 n = 4 p + q ( − 1 ) n = − p + q 辺々引いて 4 n 4^n 4 n − ( − 1 ) n {}- (-1)^n − ( − 1 ) n = 5 p {}= 5p = 5 p より
p \displaystyle p p = 4 n − ( − 1 ) n 5 , \displaystyle {}= \frac{4^n - (-1)^n}{5}, = 5 4 n − ( − 1 ) n , q \displaystyle q q = ( − 1 ) n \displaystyle {}= (-1)^n = ( − 1 ) n + p \displaystyle {}+ p + p = 4 n + 4 ⋅ ( − 1 ) n 5 \displaystyle {}= \frac{4^n + 4 \cdot (-1)^n}{5} = 5 4 n + 4 ⋅ ( − 1 ) n ケーリー・ハミルトンの定理より A 2 A^2 A 2 − 3 A {}- 3A − 3 A − 4 E {}- 4E − 4 E = O {}= O = O なので、上の恒等式に A A A を代入すると第 1 項が消えて A n A^n A n = p A {}= pA = p A + q E {}+ qE + q E です。
A n \displaystyle A^n A n = p ( 1 2 3 2 ) \displaystyle {}= p\begin{pmatrix} 1 & 2 \\ 3 & 2 \end{pmatrix} = p ( 1 3 2 2 ) + q ( 1 0 0 1 ) \displaystyle {}+ q\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} + q ( 1 0 0 1 ) = ( p + q 2 p 3 p 2 p + q ) \displaystyle {}= \begin{pmatrix} p + q & 2p \\ 3p & 2p + q \end{pmatrix} = ( p + q 3 p 2 p 2 p + q ) p p p + q {}+ q + q = 2 ⋅ 4 n + 3 ( − 1 ) n 5 {}= \dfrac{2 \cdot 4^n + 3(-1)^n}{5} = 5 2 ⋅ 4 n + 3 ( − 1 ) n ,2 p 2p 2 p + q {}+ q + q = 3 ⋅ 4 n + 2 ( − 1 ) n 5 {}= \dfrac{3 \cdot 4^n + 2(-1)^n}{5} = 5 3 ⋅ 4 n + 2 ( − 1 ) n を使います。成分が長くなるので、4 n 4^n 4 n = u {}= u = u ,( − 1 ) n (-1)^n ( − 1 ) n = v {}= v = v と書くと
A n \displaystyle A^n A n = 1 5 ( 2 u + 3 v 2 ( u − v ) 3 ( u − v ) 3 u + 2 v ) ‾ \displaystyle {}= \underline{\frac{1}{5}\begin{pmatrix} 2u + 3v & 2(u - v) \\ 3(u - v) & 3u + 2v \end{pmatrix}} = 5 1 ( 2 u + 3 v 3 ( u − v ) 2 ( u − v ) 3 u + 2 v ) ただし u u u = 4 n {}= 4^n = 4 n ,v v v = ( − 1 ) n {}= (-1)^n = ( − 1 ) n です。
検算します。n n n = 1 {}= 1 = 1 なら u u u = 4 {}= 4 = 4 ,v v v = − 1 {}= -1 = − 1 で 1 5 ( 8 − 3 2 × 5 3 × 5 12 − 2 ) \dfrac{1}{5}\begin{pmatrix} 8 - 3 & 2 \times 5 \\ 3 \times 5 & 12 - 2 \end{pmatrix} 5 1 ( 8 − 3 3 × 5 2 × 5 12 − 2 ) = ( 1 2 3 2 ) {}= \begin{pmatrix} 1 & 2 \\ 3 & 2 \end{pmatrix} = ( 1 3 2 2 ) で A A A に戻ります。n n n = 2 {}= 2 = 2 なら 1 5 ( 35 30 45 50 ) \dfrac{1}{5}\begin{pmatrix} 35 & 30 \\ 45 & 50 \end{pmatrix} 5 1 ( 35 45 30 50 ) = ( 7 6 9 10 ) {}= \begin{pmatrix} 7 & 6 \\ 9 & 10 \end{pmatrix} = ( 7 9 6 10 ) で、直接 A 2 A^2 A 2 を計算した値と一致します。n n n = 1 {}= 1 = 1 と n n n = 2 {}= 2 = 2 で確かめるのが、この種の答えの検算の定石です。
解が重なるとき
公式5 の連立方程式は、α \alpha α ≠ β {}\neq \beta = β でなければ作れません。解が重なると式が 1 本しか立たないからです。この場合は別の道を通ります。
2 次方程式が ( x − α ) 2 (x - \alpha)^2 ( x − α ) 2 = 0 {}= 0 = 0 の形なら、ケーリー・ハミルトンの定理は
( A − α E ) 2 \displaystyle (A - \alpha E)^2 ( A − α E ) 2 = O \displaystyle {}= O = O
と書けます。2 乗すると消えてしまう行列が現れるわけです。そこで N N N = A {}= A = A − α E {}- \alpha E − α E とおけば A A A = α E {}= \alpha E = α E + N {}+ N + N で、N 2 N^2 N 2 = O {}= O = O 。
α E \alpha E α E はどんな行列とも交換可能(第 2 章)なので、( α E + N ) n (\alpha E + N)^n ( α E + N ) n には二項定理がそのまま使えます 。しかも N 2 N^2 N 2 から先はすべて O O O なので、展開しても 2 項しか生き残りません。
2 次方程式が重解をもつときは A A A から解の分だけ E E E を引くと 2 乗して消える行列になり、二項定理の最初の 2 項だけで A n A^n A n が書けるということです。
例題5:
A n A^n A n を求める(解が重なるとき)
A A A = ( 4 1 − 1 2 ) {}= \begin{pmatrix} 4 & 1 \\ -1 & 2 \end{pmatrix} = ( 4 − 1 1 2 ) について、A n A^n A n を求めなさい。ただし n n n は正の整数とする。
【解答】
a a a + d {}+ d + d = 4 {}= 4 = 4 + 2 {}+ 2 + 2 = 6 {}= 6 = 6 ,Δ \Delta Δ = 4 × 2 {}= 4 \times 2 = 4 × 2 − 1 × ( − 1 ) {}- 1 \times (-1) − 1 × ( − 1 ) = 9 {}= 9 = 9 なので
x 2 \displaystyle x^2 x 2 − 6 x \displaystyle {}- 6x − 6 x + 9 \displaystyle {}+ 9 + 9 = 0 \displaystyle {}= 0 = 0 ⟹ ( x − 3 ) 2 \displaystyle {}\Longrightarrow \quad (x - 3)^2 ⟹ ( x − 3 ) 2 = 0 \displaystyle {}= 0 = 0 重解 x x x = 3 {}= 3 = 3 です。そこで N N N = A {}= A = A − 3 E {}- 3E − 3 E とおきます。
N \displaystyle N N = ( 4 − 3 1 − 1 2 − 3 ) \displaystyle {}= \begin{pmatrix} 4 - 3 & 1 \\ -1 & 2 - 3 \end{pmatrix} = ( 4 − 3 − 1 1 2 − 3 ) = ( 1 1 − 1 − 1 ) , \displaystyle {}= \begin{pmatrix} 1 & 1 \\ -1 & -1 \end{pmatrix}, = ( 1 − 1 1 − 1 ) , N 2 \displaystyle N^2 N 2 = ( 1 − 1 1 − 1 − 1 + 1 − 1 + 1 ) \displaystyle {}= \begin{pmatrix} 1 - 1 & 1 - 1 \\ -1 + 1 & -1 + 1 \end{pmatrix} = ( 1 − 1 − 1 + 1 1 − 1 − 1 + 1 ) = O \displaystyle {}= O = O たしかに N 2 N^2 N 2 = O {}= O = O です。A A A = 3 E {}= 3E = 3 E + N {}+ N + N で、3 E 3E 3 E と N N N は交換可能なので二項定理が使えます。
A n \displaystyle A^n A n = ( 3 E + N ) n \displaystyle {}= (3E + N)^n = ( 3 E + N ) n = 3 n E \displaystyle {}= 3^nE = 3 n E + n ⋅ 3 n − 1 N \displaystyle {}+ n \cdot 3^{n-1}N + n ⋅ 3 n − 1 N + n C 2 ⋅ 3 n − 2 N 2 \displaystyle {}+ {}_n\mathrm{C}_2 \cdot 3^{n-2}N^2 + n C 2 ⋅ 3 n − 2 N 2 + ⋯ \displaystyle {}+ \cdots + ⋯ = 3 n E \displaystyle {}= 3^nE = 3 n E + n ⋅ 3 n − 1 N \displaystyle {}+ n \cdot 3^{n-1}N + n ⋅ 3 n − 1 N 第 3 項から先は N 2 N^2 N 2 を含むのですべて O O O です。成分に直します。
A n \displaystyle A^n A n = 3 n ( 1 0 0 1 ) \displaystyle {}= 3^n\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} = 3 n ( 1 0 0 1 ) + n ⋅ 3 n − 1 ( 1 1 − 1 − 1 ) \displaystyle {}+ n \cdot 3^{n-1}\begin{pmatrix} 1 & 1 \\ -1 & -1 \end{pmatrix} + n ⋅ 3 n − 1 ( 1 − 1 1 − 1 ) = 3 n − 1 ( n + 3 n − n 3 − n ) ‾ \displaystyle {}= \underline{3^{n-1}\begin{pmatrix} n + 3 & n \\ -n & 3 - n \end{pmatrix}} = 3 n − 1 ( n + 3 − n n 3 − n ) n n n = 1 {}= 1 = 1 のとき 3 0 ( 4 1 − 1 2 ) 3^0\begin{pmatrix} 4 & 1 \\ -1 & 2 \end{pmatrix} 3 0 ( 4 − 1 1 2 ) = A {}= A = A 、n n n = 2 {}= 2 = 2 のとき 3 ( 5 2 − 2 1 ) 3\begin{pmatrix} 5 & 2 \\ -2 & 1 \end{pmatrix} 3 ( 5 − 2 2 1 ) = ( 15 6 − 6 3 ) {}= \begin{pmatrix} 15 & 6 \\ -6 & 3 \end{pmatrix} = ( 15 − 6 6 3 ) で、A 2 A^2 A 2 = 6 A {}= 6A = 6 A − 9 E {}- 9E − 9 E から求めた値と合っています。
N 2 N^2 N 2 = O {}= O = O なのに N N N ≠ O {}\neq O = O という点は、行列ならではです。数なら x 2 x^2 x 2 = 0 {}= 0 = 0 から x x x = 0 {}= 0 = 0 が出ますが、行列では出ません。零因子(第 2 章)の話が、こんなところで役に立っています。
対角行列に化ける形が与えられたとき
もう 1 つ、入試でよく出る形を見ておきます。うまい行列 P P P を持ってきて P − 1 A P P^{-1}AP P − 1 A P を作ると、対角成分だけが残った行列になることがあります。
P − 1 A P \displaystyle P^{-1}AP P − 1 A P = ( α 0 0 β ) \displaystyle {}= \begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix} = ( α 0 0 β )
対角成分以外が 0 0 0 の行列を対角行列 といいます。対角行列の累乗はとても簡単で、成分をそれぞれ n n n 乗するだけです。
( α 0 0 β ) n \displaystyle \begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix}^{\,n} ( α 0 0 β ) n = ( α n 0 0 β n ) \displaystyle {}= \begin{pmatrix} \alpha^n & 0 \\ 0 & \beta^n \end{pmatrix} = ( α n 0 0 β n )
この簡単さを A A A 側に持ち帰れないでしょうか。持ち帰れます。D D D = P − 1 A P {}= P^{-1}AP = P − 1 A P とおくと
D n \displaystyle D^n D n = ( P − 1 A P ) ( P − 1 A P ) ⋯ ( P − 1 A P ) \displaystyle {}= (P^{-1}AP)(P^{-1}AP)\cdots(P^{-1}AP) = ( P − 1 A P ) ( P − 1 A P ) ⋯ ( P − 1 A P )
で、となり合う P P − 1 PP^{-1} P P − 1 がすべて E E E になって消え、まん中の A A A だけが n n n 個残ります。
D n \displaystyle D^n D n = P − 1 A n P \displaystyle {}= P^{-1}A^nP = P − 1 A n P
左から P P P 、 右から P − 1 P^{-1} P − 1 をかければ A n A^n A n が取り出せます。
P − 1 A P P^{-1}AP P − 1 A P という形は、第 3 章の実践問題で一度出ています。det ( P − 1 A P ) \det(P^{-1}AP) det ( P − 1 A P ) = det A {}= \det A = det A 、 つまり P P P で挟んでも行列式は変わらない、という問題でした。行列式だけでなく対角和も変わりません。A A A と D D D は、同じものを別の見方で書いたものだからです。
例えるなら、方眼紙の目盛りを引き直す作業です。ふつうの x y xy x y 軸で見ると A A A は「伸ばしながらねじる」複雑な操作に見えるのに、P P P の 2 つの列が指す 2 方向を新しい軸にとると、それぞれの方向に α \alpha α 倍・β \beta β 倍するだけの単純な操作になっている。目盛りの取り替えが P P P 、 取り替えた世界での姿が D D D です。この見方は、第 6 章以降の一次変換でもう一度出てきます。
P − 1 A P P^{-1}AP P − 1 A P が対角行列になるなら、A n A^n A n は「対角行列を n n n 乗して P P P で挟み直す」だけで求まるということです。
例題6:対角行列を経由して
A n A^n A n を求める
A A A = ( 4 2 1 3 ) {}= \begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix} = ( 4 1 2 3 ) ,P P P = ( 2 1 1 − 1 ) {}= \begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix} = ( 2 1 1 − 1 ) とする。P − 1 A P P^{-1}AP P − 1 A P を求め、それを使って A n A^n A n を求めなさい。
【解答】
まず P − 1 P^{-1} P − 1 を求めます。P P P の行列式は 2 × ( − 1 ) 2 \times (-1) 2 × ( − 1 ) − 1 × 1 {}- 1 \times 1 − 1 × 1 = − 3 {}= -3 = − 3 ≠ 0 {}\neq 0 = 0 なので
P − 1 \displaystyle P^{-1} P − 1 = 1 − 3 ( − 1 − 1 − 1 2 ) \displaystyle {}= \frac{1}{-3}\begin{pmatrix} -1 & -1 \\ -1 & 2 \end{pmatrix} = − 3 1 ( − 1 − 1 − 1 2 ) = 1 3 ( 1 1 1 − 2 ) \displaystyle {}= \frac{1}{3}\begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix} = 3 1 ( 1 1 1 − 2 ) A P AP A P を先に計算します。
A P \displaystyle AP A P = ( 4 2 1 3 ) ( 2 1 1 − 1 ) \displaystyle {}= \begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}\begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix} = ( 4 1 2 3 ) ( 2 1 1 − 1 ) = ( 10 2 5 − 2 ) \displaystyle {}= \begin{pmatrix} 10 & 2 \\ 5 & -2 \end{pmatrix} = ( 10 5 2 − 2 ) P − 1 ( A P ) \displaystyle P^{-1}(AP) P − 1 ( A P ) = 1 3 ( 1 1 1 − 2 ) ( 10 2 5 − 2 ) \displaystyle {}= \frac{1}{3}\begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix}\begin{pmatrix} 10 & 2 \\ 5 & -2 \end{pmatrix} = 3 1 ( 1 1 1 − 2 ) ( 10 5 2 − 2 ) = 1 3 ( 15 0 0 6 ) \displaystyle {}= \frac{1}{3}\begin{pmatrix} 15 & 0 \\ 0 & 6 \end{pmatrix} = 3 1 ( 15 0 0 6 ) = ( 5 0 0 2 ) ‾ \displaystyle {}= \underline{\begin{pmatrix} 5 & 0 \\ 0 & 2 \end{pmatrix}} = ( 5 0 0 2 ) 対角行列になりました。D D D = ( 5 0 0 2 ) {}= \begin{pmatrix} 5 & 0 \\ 0 & 2 \end{pmatrix} = ( 5 0 0 2 ) とおくと D n D^n D n = ( 5 n 0 0 2 n ) {}= \begin{pmatrix} 5^n & 0 \\ 0 & 2^n \end{pmatrix} = ( 5 n 0 0 2 n ) なので、公式7 より
A n \displaystyle A^n A n = P D n P − 1 \displaystyle {}= PD^nP^{-1} = P D n P − 1 = ( 2 1 1 − 1 ) ( 5 n 0 0 2 n ) \displaystyle {}= \begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix}\begin{pmatrix} 5^n & 0 \\ 0 & 2^n \end{pmatrix} = ( 2 1 1 − 1 ) ( 5 n 0 0 2 n ) ⋅ 1 3 ( 1 1 1 − 2 ) \displaystyle \qquad \cdot \frac{1}{3}\begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix} ⋅ 3 1 ( 1 1 1 − 2 ) 左の 2 つを先に計算します。
( 2 1 1 − 1 ) ( 5 n 0 0 2 n ) \displaystyle \begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix}\begin{pmatrix} 5^n & 0 \\ 0 & 2^n \end{pmatrix} ( 2 1 1 − 1 ) ( 5 n 0 0 2 n ) = ( 2 ⋅ 5 n 2 n 5 n − 2 n ) \displaystyle {}= \begin{pmatrix} 2 \cdot 5^n & 2^n \\ 5^n & -2^n \end{pmatrix} = ( 2 ⋅ 5 n 5 n 2 n − 2 n ) A n \displaystyle A^n A n = 1 3 ( 2 ⋅ 5 n 2 n 5 n − 2 n ) ( 1 1 1 − 2 ) \displaystyle {}= \frac{1}{3}\begin{pmatrix} 2 \cdot 5^n & 2^n \\ 5^n & -2^n \end{pmatrix}\begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix} = 3 1 ( 2 ⋅ 5 n 5 n 2 n − 2 n ) ( 1 1 1 − 2 ) = 1 3 ( 2 ⋅ 5 n + 2 n 2 ⋅ 5 n − 2 ⋅ 2 n 5 n − 2 n 5 n + 2 ⋅ 2 n ) ‾ \displaystyle {}= \underline{\frac{1}{3}\begin{pmatrix} 2 \cdot 5^n + 2^n & 2 \cdot 5^n - 2 \cdot 2^n \\ 5^n - 2^n & 5^n + 2 \cdot 2^n \end{pmatrix}} = 3 1 ( 2 ⋅ 5 n + 2 n 5 n − 2 n 2 ⋅ 5 n − 2 ⋅ 2 n 5 n + 2 ⋅ 2 n ) n n n = 1 {}= 1 = 1 で 1 3 ( 12 6 3 9 ) \dfrac{1}{3}\begin{pmatrix} 12 & 6 \\ 3 & 9 \end{pmatrix} 3 1 ( 12 3 6 9 ) = ( 4 2 1 3 ) {}= \begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix} = ( 4 1 2 3 ) 、n n n = 2 {}= 2 = 2 で 1 3 ( 54 42 21 33 ) \dfrac{1}{3}\begin{pmatrix} 54 & 42 \\ 21 & 33 \end{pmatrix} 3 1 ( 54 21 42 33 ) = ( 18 14 7 11 ) {}= \begin{pmatrix} 18 & 14 \\ 7 & 11 \end{pmatrix} = ( 18 7 14 11 ) 。 どちらも直接計算した値に一致します。
なお D D D の対角成分 5 5 5 ,2 2 2 は、A A A の背負う 2 次方程式 x 2 x^2 x 2 − 7 x {}- 7x − 7 x + 10 {}+ 10 + 10 = 0 {}= 0 = 0 の 2 つの解です。偶然ではありません。P P P で挟んでも対角和と行列式は変わらないので、D D D も同じ 2 次方程式を背負っているからです。公式5 で α \alpha α ,β \beta β と呼んでいたものの正体が、ここで見えます。
どの方法を使うか
A n A^n A n の求め方が出そろいました。整理しておきます。
最後に、A n A^n A n が実際に何かの役に立つ場面を見ておきます。
A 店
B 店
3/4
1/2
1/4
1/2
翌日 B 店へ
翌日 A 店へ
例題7:何日後にどちらの店にいるか
ある人は毎日、昼食を A 店か B 店のどちらかでとる。A 店で食べた翌日にまた A 店を選ぶ確率は 3 4 \dfrac{3}{4} 4 3 、 B 店で食べた翌日に A 店を選ぶ確率は 1 2 \dfrac{1}{2} 2 1 である。1 日目は A 店で食べた。n n n 日目に A 店で食べる確率 p n p_n p n を求め、n n n を限りなく大きくしたときの値を調べなさい。
【解答】
n n n 日目に A 店で食べる確率を p n p_n p n 、 B 店で食べる確率を q n q_n q n とします。n n n + 1 {}+ 1 + 1 日目に A 店にいるのは「A 店にいて残る」か「B 店にいて移る」かのどちらかなので
{ p n + 1 = 3 4 p n + 1 2 q n q n + 1 = 1 4 p n + 1 2 q n \displaystyle \begin{cases} p_{n+1} = \dfrac{3}{4}p_n + \dfrac{1}{2}q_n \\[4pt] q_{n+1} = \dfrac{1}{4}p_n + \dfrac{1}{2}q_n \end{cases} ⎩ ⎨ ⎧ p n + 1 = 4 3 p n + 2 1 q n q n + 1 = 4 1 p n + 2 1 q n すなわち \displaystyle \text{すなわち} すなわち ( p n + 1 q n + 1 ) \displaystyle \begin{pmatrix} p_{n+1} \\ q_{n+1} \end{pmatrix} ( p n + 1 q n + 1 ) = A ( p n q n ) , \displaystyle {}= A\begin{pmatrix} p_n \\ q_n \end{pmatrix}, = A ( p n q n ) , A \displaystyle A A = ( 3 4 1 2 1 4 1 2 ) \displaystyle {}= \begin{pmatrix} \dfrac{3}{4} & \dfrac{1}{2} \\[4pt] \dfrac{1}{4} & \dfrac{1}{2} \end{pmatrix} = 4 3 4 1 2 1 2 1 1 日目は A 店なので ( p 1 q 1 ) \begin{pmatrix} p_1 \\ q_1 \end{pmatrix} ( p 1 q 1 ) = ( 1 0 ) {}= \begin{pmatrix} 1 \\ 0 \end{pmatrix} = ( 1 0 ) です。n n n 日目の確率は、これに A A A を n n n − 1 {}- 1 − 1 回かけたもの、つまり
( p n q n ) \displaystyle \begin{pmatrix} p_n \\ q_n \end{pmatrix} ( p n q n ) = A n − 1 ( 1 0 ) \displaystyle {}= A^{n-1}\begin{pmatrix} 1 \\ 0 \end{pmatrix} = A n − 1 ( 1 0 ) A n − 1 A^{n-1} A n − 1 を求めます。a a a + d {}+ d + d = 3 4 {}= \dfrac{3}{4} = 4 3 + 1 2 {}+ \dfrac{1}{2} + 2 1 = 5 4 {}= \dfrac{5}{4} = 4 5 ,Δ \Delta Δ = 3 4 ⋅ 1 2 {}= \dfrac{3}{4} \cdot \dfrac{1}{2} = 4 3 ⋅ 2 1 − 1 2 ⋅ 1 4 {}- \dfrac{1}{2} \cdot \dfrac{1}{4} − 2 1 ⋅ 4 1 = 1 4 {}= \dfrac{1}{4} = 4 1 なので、背負う 2 次方程式は
x 2 \displaystyle x^2 x 2 − 5 4 x \displaystyle {}- \frac{5}{4}x − 4 5 x + 1 4 \displaystyle {}+ \frac{1}{4} + 4 1 = 0 \displaystyle {}= 0 = 0 ⟹ 4 x 2 \displaystyle {}\Longrightarrow \quad 4x^2 ⟹ 4 x 2 − 5 x \displaystyle {}- 5x − 5 x + 1 \displaystyle {}+ 1 + 1 = 0 \displaystyle {}= 0 = 0 ⟹ ( 4 x − 1 ) ( x − 1 ) \displaystyle {}\Longrightarrow \quad (4x - 1)(x - 1) ⟹ ( 4 x − 1 ) ( x − 1 ) = 0 \displaystyle {}= 0 = 0 解は x x x = 1 , {}= 1, = 1 , 1 4 \ \dfrac{1}{4} 4 1 です。m m m = n {}= n = n − 1 {}- 1 − 1 とおき、x m x^m x m をこの 2 次式で割った余りを p x px p x + q {}+ q + q とすると
{ 1 = p + q ( 1 4 ) m = 1 4 p + q \begin{cases} 1 = p + q \\[2pt] \left(\dfrac{1}{4}\right)^{m} = \dfrac{1}{4}p + q \end{cases} ⎩ ⎨ ⎧ 1 = p + q ( 4 1 ) m = 4 1 p + q 辺々引いて 1 1 1 − ( 1 4 ) m {}- \left(\dfrac{1}{4}\right)^{m} − ( 4 1 ) m = 3 4 p {}= \dfrac{3}{4}p = 4 3 p より p p p = 4 3 { 1 − ( 1 4 ) m } {}= \dfrac{4}{3}\left\{1 - \left(\dfrac{1}{4}\right)^{m}\right\} = 3 4 { 1 − ( 4 1 ) m } 、q q q = 1 {}= 1 = 1 − p {}- p − p です。
A m A^m A m = p A {}= pA = p A + q E {}+ qE + q E なので、( 1 0 ) \begin{pmatrix} 1 \\ 0 \end{pmatrix} ( 1 0 ) にかけると
A m ( 1 0 ) \displaystyle A^m\begin{pmatrix} 1 \\ 0 \end{pmatrix} A m ( 1 0 ) = p ( 3 / 4 1 / 4 ) \displaystyle {}= p\begin{pmatrix} 3/4 \\ 1/4 \end{pmatrix} = p ( 3/4 1/4 ) + q ( 1 0 ) \displaystyle {}+ q\begin{pmatrix} 1 \\ 0 \end{pmatrix} + q ( 1 0 ) = ( 3 4 p + q 1 4 p ) \displaystyle {}= \begin{pmatrix} \dfrac{3}{4}p + q \\[4pt] \dfrac{1}{4}p \end{pmatrix} = 4 3 p + q 4 1 p 上の成分が p n p_n p n です。q q q = 1 {}= 1 = 1 − p {}- p − p を使って
p n \displaystyle p_n p n = 3 4 p \displaystyle {}= \frac{3}{4}p = 4 3 p + ( 1 − p ) \displaystyle {}+ (1 - p) + ( 1 − p ) = 1 \displaystyle {}= 1 = 1 − 1 4 p \displaystyle {}- \frac{1}{4}p − 4 1 p = 1 \displaystyle {}= 1 = 1 − 1 3 { 1 − ( 1 4 ) m } \displaystyle {}- \frac{1}{3}\left\{1 - \left(\frac{1}{4}\right)^{m}\right\} − 3 1 { 1 − ( 4 1 ) m } = 2 3 ‾ \displaystyle {}= \underline{\rule[-0.95em]{0em}{2.6040em}\frac{2}{3}} = 3 2 + 1 3 ( 1 4 ) n − 1 ‾ \displaystyle \underline{\rule[-0.95em]{0em}{2.6040em}{}+ \frac{1}{3}\left(\frac{1}{4}\right)^{n-1}} + 3 1 ( 4 1 ) n − 1 検算します。n n n = 1 {}= 1 = 1 なら 2 3 \dfrac{2}{3} 3 2 + 1 3 {}+ \dfrac{1}{3} + 3 1 = 1 {}= 1 = 1 で、初日は必ず A 店。n n n = 2 {}= 2 = 2 なら 2 3 \dfrac{2}{3} 3 2 + 1 12 {}+ \dfrac{1}{12} + 12 1 = 3 4 {}= \dfrac{3}{4} = 4 3 で、問題の設定どおりです。
( 1 4 ) n − 1 \left(\dfrac{1}{4}\right)^{n-1} ( 4 1 ) n − 1 は n n n が大きくなると 0 0 0 に近づくので
p n \displaystyle p_n p n ⟶ 2 3 \displaystyle {}\longrightarrow \frac{2}{3} ⟶ 3 2 日がたつにつれ、A 店にいる確率は 2 3 \dfrac{2}{3} 3 2 に落ち着きます。おもしろいのは、初日にどちらの店にいたかが効かなくなる ことです。q 1 q_1 q 1 = 1 {}= 1 = 1 (初日は B 店)から始めても、同じ 2 3 \dfrac{2}{3} 3 2 に近づきます。A n A^n A n の姿が n n n とともに一定の形へ寄っていくためで、この落ち着き先を確率の分野では定常状態と呼びました(数列 第 9 章)。小話では、この考え方が Web 検索の順位づけに使われている話を取り上げます。
ここまでで、行列の計算としての道具はひととおりそろいました。次の第 6 章からは主役を交代させ、行列を平面上の点を動かす仕掛け として見直します。A ( x y ) A\begin{pmatrix} x \\ y \end{pmatrix} A ( x y ) を「点 ( x , y ) (x,\ y) ( x , y ) の像」と呼んでいたあの見方を正式に一次変換 と名づけ、記号 f f f で表します。この章で身につけた A n A^n A n は、そこでは「同じ移動を n n n 回くり返す」という意味を持つことになります。
※ここは発展ページ です。高校の範囲では「A 2 A^2 A 2 = ( a + d ) A {}= (a+d)A = ( a + d ) A − Δ E {}- \Delta E − Δ E を使って次数を下げる」と使えれば十分ですが、そもそも a a a + d {}+ d + d と Δ \Delta Δ という係数はどこから来たのか、なぜ定理の逆だけが k E kE k E のときに壊れるのか、P − 1 A P P^{-1}AP P − 1 A P が対角行列になる P P P はどうやって見つけるのかが気になる人のために、少し厳密な形にまとめておきます。
本文では公式1 を成分計算で確かめました。正しい証明ですが、これでは「たまたまそうなった」ようにしか見えません。係数 a a a + d {}+ d + d と Δ \Delta Δ の出どころを示すところから始めます。第 3 章の厳密定義(行列式の積の法則、正則の言いかえ)と第 4 章の厳密定義(1 次独立と解集合)を前提にします。
特性多項式
定義1:対角和と特性多項式
n n n 次の正方行列 A A A = ( a i j ) {}= (a_{ij}) = ( a ij ) に対し、対角成分の和
tr A \displaystyle \operatorname{tr} A tr A = a 11 \displaystyle {}= a_{11} = a 11 + a 22 \displaystyle {}+ a_{22} + a 22 + ⋯ \displaystyle {}+ \cdots + ⋯ + a n n \displaystyle {}+ a_{nn} + a nn を A A A の対角和 (トレース)という。また、x x x を変数として
φ A ( x ) \displaystyle \varphi_A(x) φ A ( x ) = det ( x E − A ) \displaystyle {}= \det(xE - A) = det ( x E − A ) で定まる x x x の n n n 次多項式を A A A の特性多項式 という。
2 次の場合に書き下してみます。
x E \displaystyle xE x E − A \displaystyle {}- A − A = ( x − a − b − c x − d ) \displaystyle {}= \begin{pmatrix} x - a & -b \\ -c & x - d \end{pmatrix} = ( x − a − c − b x − d )
φ A ( x ) \displaystyle \varphi_A(x) φ A ( x ) = ( x − a ) ( x − d ) \displaystyle {}= (x - a)(x - d) = ( x − a ) ( x − d ) − ( − b ) ( − c ) \displaystyle {}- (-b)(-c) − ( − b ) ( − c ) = x 2 \displaystyle {}= x^2 = x 2 − ( a + d ) x \displaystyle {}- (a + d)x − ( a + d ) x + ( a d − b c ) \displaystyle {}+ (ad - bc) + ( a d − b c )
本文で「A A A が背負う 2 次方程式」と呼んでいたものの正体です。係数 a a a + d {}+ d + d と Δ \Delta Δ は、どこからともなく降ってきたのではなく、det ( x E − A ) \det(xE - A) det ( x E − A ) を展開すれば自動的に現れます。x n − 1 x^{n-1} x n − 1 の係数が − tr A -\operatorname{tr} A − tr A 、 定数項が ( − 1 ) n det A (-1)^n \det A ( − 1 ) n det A になるのは、一般の n n n 次でも変わりません。
定理1:ケーリー・ハミルトンの定理(2 次)
2 次の正方行列 A A A に対し φ A ( A ) \varphi_A(A) φ A ( A ) = O {}= O = O である。すなわち
A 2 \displaystyle A^2 A 2 − ( tr A ) A \displaystyle {}- (\operatorname{tr} A)A − ( tr A ) A + ( det A ) E \displaystyle {}+ (\det A)E + ( det A ) E = O \displaystyle {}= O = O 証明 A A A = ( a b c d ) {}= \begin{pmatrix} a & b \\ c & d \end{pmatrix} = ( a c b d ) とする。第 3 章 定理2 の証明で A ~ \tilde{A} A ~ = ( d − b − c a ) {}= \begin{pmatrix} d & -b \\ -c & a \end{pmatrix} = ( d − c − b a ) が A A ~ A\tilde{A} A A ~ = A ~ A {}= \tilde{A}A = A ~ A = ( det A ) E {}= (\det A)E = ( det A ) E を満たすことを示した。ここで
A ~ \displaystyle \tilde{A} A ~ = ( a + d 0 0 a + d ) \displaystyle {}= \begin{pmatrix} a + d & 0 \\ 0 & a + d \end{pmatrix} = ( a + d 0 0 a + d ) − ( a b c d ) \displaystyle {}- \begin{pmatrix} a & b \\ c & d \end{pmatrix} − ( a c b d ) = ( tr A ) E \displaystyle {}= (\operatorname{tr} A)E = ( tr A ) E − A \displaystyle {}- A − A であるから、A A ~ A\tilde{A} A A ~ = ( det A ) E {}= (\det A)E = ( det A ) E に代入して
A { ( tr A ) E − A } \displaystyle A\{(\operatorname{tr} A)E - A\} A {( tr A ) E − A } = ( det A ) E \displaystyle {}= (\det A)E = ( det A ) E 左辺を展開すると ( tr A ) A (\operatorname{tr} A)A ( tr A ) A − A 2 {}- A^2 − A 2 である。移項して
A 2 \displaystyle A^2 A 2 − ( tr A ) A \displaystyle {}- (\operatorname{tr} A)A − ( tr A ) A + ( det A ) E \displaystyle {}+ (\det A)E + ( det A ) E = O \displaystyle {}= O = O を得る。(証明終)
本文の成分計算と同じことですが、こう書くと A ~ \tilde{A} A ~ の正体が ( tr A ) E (\operatorname{tr} A)E ( tr A ) E − A {}- A − A であり、逆行列の公式とこの定理がまったく同じ 1 本の等式の言いかえだったことがはっきりします。第 3 章で「a a a と d d d を入れかえ、b b b と c c c の符号を変える」と手順で覚えた操作には、こういう意味があったわけです。
有名な誤った証明
この定理には、広く知られた「証明もどき」があります。
φ A ( x ) \varphi_A(x) φ A ( x ) = det ( x E − A ) {}= \det(xE - A) = det ( x E − A ) の x x x に A A A を代入すると φ A ( A ) \varphi_A(A) φ A ( A ) = det ( A E − A ) {}= \det(AE - A) = det ( A E − A ) = det ( A − A ) {}= \det(A - A) = det ( A − A ) = det O {}= \det O = det O = 0 {}= 0 = 0 。 よって φ A ( A ) \varphi_A(A) φ A ( A ) = O {}= O = O である。
一見すると鮮やかですが、誤りです。どこが壊れているのかを、はっきりさせておきます。
第一に、結論の型が合っていません 。det O \det O det O は数の 0 0 0 であって、零行列 O O O ではありません。示したいのは「行列が O O O に等しい」ことなのに、出てきたのは「数が 0 0 0 に等しい」こと。0 0 0 と O O O を同じ記号のように扱った時点で、主張がすり替わっています。
第二に、そもそも代入が成り立っていません 。φ A ( x ) \varphi_A(x) φ A ( x ) = det ( x E − A ) {}= \det(xE - A) = det ( x E − A ) という等式は、x x x が数であるときの等式です。左辺の φ A \varphi_A φ A は多項式なので、係数を行列の係数と読みかえて x x x に行列を入れることができます(x 2 x^2 x 2 は A 2 A^2 A 2 、 定数項 q q q は q E qE q E と読む)。ところが右辺の x E xE x E は「数 x x x と単位行列の積」であって、ここに行列を入れる意味は定められていません。かりに A E AE A E = A {}= A = A と読んだとしても、それは右辺の式を勝手に別のものに取りかえたにすぎず、等式が保たれる保証はどこにもありません。
同じ論法が通ってしまうと困る例を挙げれば、あやしさがはっきりします。ψ ( x ) \psi(x) ψ ( x ) = det ( x E − A ) {}= \det(xE - A) = det ( x E − A ) + 1 {}+ 1 + 1 という多項式を考えると、同じ手つきで ψ ( A ) \psi(A) ψ ( A ) = det ( A − A ) {}= \det(A - A) = det ( A − A ) + 1 {}+ 1 + 1 = 1 {}= 1 = 1 、 したがって「ψ ( A ) \psi(A) ψ ( A ) = O {}= O = O 」 ——明らかに間違った結論が出ます。
一般の n n n 次では、A ~ \tilde{A} A ~ にあたる余因子行列 を使い、その成分が x x x の多項式であることを利用して、係数を次数ごとに比べるという筋道で証明します。フロベニウスが与えたのもこの種の議論です。
対角行列に化ける条件
本文の公式7 では、都合のよい P P P が問題文で与えられていました。では P P P はどうやって見つけるのでしょうか。
定義2:固有値と固有ベクトル
n n n 次の正方行列 A A A に対し、u ⃗ \vec{u} u ≠ 0 ⃗ {}\neq \vec{0} = 0 と実数 λ \lambda λ が
A u ⃗ \displaystyle A\vec{u} A u = λ u ⃗ \displaystyle {}= \lambda\vec{u} = λ u を満たすとき、λ \lambda λ を A A A の固有値 、u ⃗ \vec{u} u を λ \lambda λ に対する固有ベクトル という。
A u ⃗ A\vec{u} A u = λ u ⃗ {}= \lambda\vec{u} = λ u は ( λ E − A ) u ⃗ (\lambda E - A)\vec{u} ( λ E − A ) u = 0 ⃗ {}= \vec{0} = 0 と書けるので、λ \lambda λ が固有値であることと、同次連立一次方程式 ( λ E − A ) x ⃗ (\lambda E - A)\vec{x} ( λ E − A ) x = 0 ⃗ {}= \vec{0} = 0 が自明でない解をもつこと、すなわち φ A ( λ ) \varphi_A(\lambda) φ A ( λ ) = det ( λ E − A ) {}= \det(\lambda E - A) = det ( λ E − A ) = 0 {}= 0 = 0 は同値である(第 4 章 公式5)。
固有ベクトルとは、A A A をかけても向きが変わらず、長さだけが λ \lambda λ 倍される方向のことです。本文で「方眼紙の目盛りを引き直す」と言ったときの、新しい軸の方向がこれにあたります。なお高校の教科書ではこの語を使わないので、本文では「x 2 x^2 x 2 − ( a + d ) x {}- (a+d)x − ( a + d ) x + Δ {}+ \Delta + Δ = 0 {}= 0 = 0 の解」と呼んできました。同じものです。
定理2:相異なる 2 つの実数解をもてば対角化できる
2 次の正方行列 A A A の特性多項式 φ A ( x ) \varphi_A(x) φ A ( x ) が相異なる 2 つの実数解 α , \alpha, α , β \ \beta β をもつとする。このとき正則な行列 P P P が存在して
P − 1 A P \displaystyle P^{-1}AP P − 1 A P = ( α 0 0 β ) \displaystyle {}= \begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix} = ( α 0 0 β ) となる。
証明 φ A ( α ) \varphi_A(\alpha) φ A ( α ) = det ( α E − A ) {}= \det(\alpha E - A) = det ( α E − A ) = 0 {}= 0 = 0 であるから α E \alpha E α E − A {}- A − A は正則でない。よって第 4 章 公式5 により ( α E − A ) u ⃗ (\alpha E - A)\vec{u} ( α E − A ) u = 0 ⃗ {}= \vec{0} = 0 を満たす u ⃗ \vec{u} u ≠ 0 ⃗ {}\neq \vec{0} = 0 が存在する。同様に ( β E − A ) v ⃗ (\beta E - A)\vec{v} ( β E − A ) v = 0 ⃗ {}= \vec{0} = 0 を満たす v ⃗ \vec{v} v ≠ 0 ⃗ {}\neq \vec{0} = 0 が存在する。すなわち A u ⃗ A\vec{u} A u = α u ⃗ {}= \alpha\vec{u} = α u ,A v ⃗ A\vec{v} A v = β v ⃗ {}= \beta\vec{v} = β v である。
u ⃗ \vec{u} u と v ⃗ \vec{v} v が 1 次独立であることを示す。v ⃗ \vec{v} v = t u ⃗ {}= t\vec{u} = t u (t t t ≠ 0 {}\neq 0 = 0 ) と書けたとすると
β v ⃗ \displaystyle \beta\vec{v} β v = A v ⃗ \displaystyle {}= A\vec{v} = A v = t A u ⃗ \displaystyle {}= tA\vec{u} = t A u = t α u ⃗ \displaystyle {}= t\alpha\vec{u} = t α u = α v ⃗ \displaystyle {}= \alpha\vec{v} = α v より ( β − α ) v ⃗ (\beta - \alpha)\vec{v} ( β − α ) v = 0 ⃗ {}= \vec{0} = 0 であり、v ⃗ \vec{v} v ≠ 0 ⃗ {}\neq \vec{0} = 0 から β \beta β = α {}= \alpha = α となって仮定に反する。よって 1 次独立である。
u ⃗ , \vec{u}, u , v ⃗ \ \vec{v} v をこの順に列として並べた行列を P P P = ( u ⃗ v ⃗ ) {}= (\vec{u}\ \ \vec{v}) = ( u v ) とおく。2 つの列が 1 次独立だから P P P は正則である(第 3 章 定理2)。積を列ごとに計算すると
A P \displaystyle AP A P = ( A u ⃗ A v ⃗ ) \displaystyle {}= (A\vec{u}\ \ A\vec{v}) = ( A u A v ) = ( α u ⃗ β v ⃗ ) \displaystyle {}= (\alpha\vec{u}\ \ \beta\vec{v}) = ( α u β v ) = P ( α 0 0 β ) \displaystyle {}= P\begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix} = P ( α 0 0 β ) 左から P − 1 P^{-1} P − 1 をかけて結論を得る。(証明終)
P P P の作り方まで証明が教えてくれています。φ A ( x ) \varphi_A(x) φ A ( x ) = 0 {}= 0 = 0 の解ごとに A u ⃗ A\vec{u} A u = λ u ⃗ {}= \lambda\vec{u} = λ u を解き、出てきたベクトルを列に並べる 。本文の例題6 で与えられた P P P = ( 2 1 1 − 1 ) {}= \begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix} = ( 2 1 1 − 1 ) も、この手順で作られたものです。実際 A ( 2 1 ) A\begin{pmatrix} 2 \\ 1 \end{pmatrix} A ( 2 1 ) = ( 10 5 ) {}= \begin{pmatrix} 10 \\ 5 \end{pmatrix} = ( 10 5 ) = 5 ( 2 1 ) {}= 5\begin{pmatrix} 2 \\ 1 \end{pmatrix} = 5 ( 2 1 ) 、A ( 1 − 1 ) A\begin{pmatrix} 1 \\ -1 \end{pmatrix} A ( 1 − 1 ) = ( 2 − 2 ) {}= \begin{pmatrix} 2 \\ -2 \end{pmatrix} = ( 2 − 2 ) = 2 ( 1 − 1 ) {}= 2\begin{pmatrix} 1 \\ -1 \end{pmatrix} = 2 ( 1 − 1 ) となっています。
逆に、特性多項式が重解をもつ場合は、対角行列に化けるとは限りません。本文の例題5 の A A A = ( 4 1 − 1 2 ) {}= \begin{pmatrix} 4 & 1 \\ -1 & 2 \end{pmatrix} = ( 4 − 1 1 2 ) がその例です。もし P − 1 A P P^{-1}AP P − 1 A P = 3 E {}= 3E = 3 E となる P P P があったとすれば A A A = P ( 3 E ) P − 1 {}= P(3E)P^{-1} = P ( 3 E ) P − 1 = 3 E {}= 3E = 3 E ですが、A A A ≠ 3 E {}\neq 3E = 3 E なので、そんな P P P はありません。あのとき二項定理という別の道を通ったのは、そういう事情です。
対角和と行列式は P P P で挟んでも変わらない
定理3:特性多項式の不変性
P P P を正則な行列とすると
φ P − 1 A P ( x ) \displaystyle \varphi_{P^{-1}AP}(x) φ P − 1 A P ( x ) = φ A ( x ) \displaystyle {}= \varphi_A(x) = φ A ( x ) である。とくに tr ( P − 1 A P ) \operatorname{tr}(P^{-1}AP) tr ( P − 1 A P ) = tr A {}= \operatorname{tr} A = tr A ,det ( P − 1 A P ) \det(P^{-1}AP) det ( P − 1 A P ) = det A {}= \det A = det A が成り立つ。
証明 x E xE x E = P − 1 ( x E ) P {}= P^{-1}(xE)P = P − 1 ( x E ) P であるから
x E \displaystyle xE x E − P − 1 A P \displaystyle {}- P^{-1}AP − P − 1 A P = P − 1 ( x E − A ) P \displaystyle {}= P^{-1}(xE - A)P = P − 1 ( x E − A ) P 行列式の積の法則(第 3 章 定理3)より
det ( x E − P − 1 A P ) \displaystyle \det(xE - P^{-1}AP) det ( x E − P − 1 A P ) = det ( P − 1 ) det ( x E − A ) det ( P ) \displaystyle {}= \det(P^{-1})\det(xE - A)\det(P) = det ( P − 1 ) det ( x E − A ) det ( P ) det ( P − 1 ) det ( P ) \det(P^{-1})\det(P) det ( P − 1 ) det ( P ) = det ( P − 1 P ) {}= \det(P^{-1}P) = det ( P − 1 P ) = det E {}= \det E = det E = 1 {}= 1 = 1 であるから、右辺は det ( x E − A ) \det(xE - A) det ( x E − A ) に等しい。
多項式として等しいのだから、各次の係数も等しい。x x x の係数を比べて tr ( P − 1 A P ) \operatorname{tr}(P^{-1}AP) tr ( P − 1 A P ) = tr A {}= \operatorname{tr} A = tr A 、 定数項を比べて det ( P − 1 A P ) \det(P^{-1}AP) det ( P − 1 A P ) = det A {}= \det A = det A を得る。(証明終)
第 3 章の実践問題で det ( A − 1 B A ) \det(A^{-1}BA) det ( A − 1 B A ) = det B {}= \det B = det B を示したのは、この定理の定数項の部分でした。P P P で挟む操作は、本文の言い方でいえば「方眼紙の目盛りを引き直す」こと。目盛りを変えても、対角和と行列式という 2 つの量は動かない——だからこそ、A A A と D D D が同じ特性多項式を背負うのです。本文の例題6 で D D D の対角成分 5 5 5 ,2 2 2 が x 2 x^2 x 2 − 7 x {}- 7x − 7 x + 10 {}+ 10 + 10 = 0 {}= 0 = 0 の解になっていたのは、この定理の帰結です。
定理の逆が壊れる本当の理由
本文の公式3 で、A A A = k E {}= kE = k E のときだけ定理の逆が成り立たないことを見ました。その理由を、多項式の言葉で言い直しておきます。
定義3:最小多項式
正方行列 A A A に対し、f ( A ) f(A) f ( A ) = O {}= O = O を満たす 0 0 0 でない多項式 f f f のうち、最高次の係数が 1 1 1 で次数が最小のものを A A A の最小多項式 という。
ケーリー・ハミルトンの定理により φ A ( A ) \varphi_A(A) φ A ( A ) = O {}= O = O ですから、2 次の正方行列の最小多項式の次数は 1 1 1 か 2 2 2 です。
定理4:最小多項式が 1 次であることと
A A A = k E {}= kE = k E は同値
2 次の正方行列 A A A について、A A A の最小多項式が 1 次であることと、ある実数 k k k によって A A A = k E {}= kE = k E と書けることは同値である。
証明 最小多項式が 1 次ならば、それは x x x − k {}- k − k の形であり、A A A − k E {}- kE − k E = O {}= O = O すなわち A A A = k E {}= kE = k E である。逆に A A A = k E {}= kE = k E ならば x x x − k {}- k − k が A A A を消すから、最小多項式は 1 次である。(証明終)
これで本文の公式3 が言い直せます。A 2 A^2 A 2 − p A {}- pA − p A + q E {}+ qE + q E = O {}= O = O とは「x 2 x^2 x 2 − p x {}- px − p x + q {}+ q + q が A A A を消す」ということです。A A A ≠ k E {}\neq kE = k E ならば最小多項式は 2 次で、それは φ A \varphi_A φ A 自身です。最高次の係数が 1 1 1 の 2 次式で A A A を消すものは最小多項式に限られるので、x 2 x^2 x 2 − p x {}- px − p x + q {}+ q + q = φ A ( x ) {}= \varphi_A(x) = φ A ( x ) 、 つまり p p p = tr A {}= \operatorname{tr} A = tr A ,q q q = det A {}= \det A = det A が確定します。
一方 A A A = k E {}= kE = k E のときは、最小多項式が 1 次の x x x − k {}- k − k です。x x x − k {}- k − k を因数に含む 2 次式なら何でも A A A を消してしまうので、p p p ,q q q は 1 通りに決まりません。本文で A A A = 2 E {}= 2E = 2 E が A 2 A^2 A 2 − 5 A {}- 5A − 5 A + 6 E {}+ 6E + 6 E = O {}= O = O を満たしたのは、x 2 x^2 x 2 − 5 x {}- 5x − 5 x + 6 {}+ 6 + 6 = ( x − 2 ) ( x − 3 ) {}= (x - 2)(x - 3) = ( x − 2 ) ( x − 3 ) が x x x − 2 {}- 2 − 2 を含んでいたからです。A A A の対角和が 5 5 5 だったからではありません。
「A 2 A^2 A 2 − p A {}- pA − p A + q E {}+ qE + q E = O {}= O = O から p p p = a {}= a = a + d {}+ d + d を読み取ってはいけない」という注意の中身は、x 2 x^2 x 2 − p x {}- px − p x + q {}+ q + q が最小多項式の倍数でありさえすればよく、最小多項式そのものである必要はない ということです。
この先
一般の n n n 次でも φ A ( A ) \varphi_A(A) φ A ( A ) = O {}= O = O が成り立ちます(ケーリー・ハミルトンの定理)。n n n 次の行列の n n n + 1 {}+ 1 + 1 乗以上はすべて n n n − 1 {}- 1 − 1 次以下に落とせる、という意味で、次数下げは次数がいくつでも使えます。
対角行列に化けない行列をどこまで簡単な形にできるか、という問いにはジョルダン標準形 という答えが用意されています。本文の例題5 で使った α E \alpha E α E + N {}+ N + N (N 2 N^2 N 2 = O {}= O = O ) という形は、2 次のジョルダン標準形そのものです。N N N のように何乗かすると O O O になる行列をべき零行列 といい、対角行列とべき零行列の和に分ける、というのが一般の筋書きになります。
もう 1 つ、A n A^n A n の先には行列の指数関数 があります。
e A \displaystyle e^{A} e A = E \displaystyle {}= E = E + A \displaystyle {}+ A + A + 1 2 ! A 2 \displaystyle {}+ \frac{1}{2!}A^2 + 2 ! 1 A 2 + 1 3 ! A 3 \displaystyle {}+ \frac{1}{3!}A^3 + 3 ! 1 A 3 + ⋯ \displaystyle {}+ \cdots + ⋯
数の e x e^x e x の級数(微分積分で学ぶ展開)に、そのまま行列を入れたものです。この収束を示すには級数の議論が要りますが、ひとたび定義してしまえば、微分方程式 d x ⃗ d t \dfrac{d\vec{x}}{dt} d t d x = A x ⃗ {}= A\vec{x} = A x の解が x ⃗ ( t ) \vec{x}(t) x ( t ) = e A t x ⃗ ( 0 ) {}= e^{At}\vec{x}(0) = e A t x ( 0 ) と書けます。ばねの振動、電気回路、人口の変化——連立した微分方程式は、行列の累乗の言葉で一気に片づくのです。A n A^n A n を求める技術は、その入口にあたります。
次の第 6 章からは、行列を計算の対象としてではなく、平面上の点を動かす仕掛けとして見直します。第 2 章の厳密定義で証明した「積は写像の合成」という事実が、そこで主役になります。