目次 / 行列と一次変換 / 旧数学C(課程外)

第5章 ケーリー・ハミルトンの定理と行列の n 乗

—— どんな行列も自分専用の 2 次方程式を背負っている ——

第 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$ を求めるだけの問題に化けるのです。ただし、この定理には**逆が成り立たない**という落とし穴があり、そこが入試でいちばん狙われます。

問題マップ記録を読み込み中…
未回答 25× 0○ 01か月定着 0

行列が満たす 2 次方程式

数の世界では、xx に何か具体的な数を入れれば x2x^2−5x{}- 5x+6{}+ 6 の値が決まります。行列でも同じことができます。AA を 2 次の正方行列とすると、A2A^2 も 5A5A も 6E6E も 2 次の正方行列なので、A2A^2−5A{}- 5A+6E{}+ 6E という式には意味があります。定数項を 66 ではなく 6E6E と書くのは、行列と数は足せないからです。数の 66 を行列の仲間に入れるには、単位行列を使って 6E6E の姿にしてやる必要がありました(第 2 章)。

さて、AA=(3124){}= \begin{pmatrix} 3 & 1 \\ 2 & 4 \end{pmatrix} で試してみます。aa+d{}+ d=7{}= 7,adad−bc{}- bc=12{}= 12−2{}- 2=10{}= 10 です。この 77 と 1010 を係数にして A2A^2−7A{}- 7A+10E{}+ 10E を計算すると

(1171418)\displaystyle \begin{pmatrix} 11 & 7 \\ 14 & 18 \end{pmatrix}−(2171428)\displaystyle {}- \begin{pmatrix} 21 & 7 \\ 14 & 28 \end{pmatrix}+(100010)\displaystyle {}+ \begin{pmatrix} 10 & 0 \\ 0 & 10 \end{pmatrix}=(0000)\displaystyle {}= \begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix}

きれいに零行列になりました。偶然ではありません。どんな 2 次の正方行列でも、対角成分の和 aa+d{}+ d と行列式 Δ\Delta=ad{}= ad−bc{}- bc を係数に選べば、必ずこうなります。

公式1:ケーリー・ハミルトンの定理

AA=(abcd){}= \begin{pmatrix} a & b \\ c & d \end{pmatrix} とし、Δ\Delta=ad{}= ad−bc{}- bc とする。このとき

A2\displaystyle A^2−(a+d)A\displaystyle {}- (a + d)A+ΔE\displaystyle {}+ \Delta E=O\displaystyle {}= O

が常に成り立つ。すなわち

A2\displaystyle A^2=(a+d)A\displaystyle {}= (a + d)A−ΔE\displaystyle {}- \Delta E

aa+d{}+ d(左上から右下へ並ぶ成分の和)を AA の対角和という。

証明は成分を書き下すだけです。

A2\displaystyle A^2=(abcd)(abcd)\displaystyle {}= \begin{pmatrix} a & b \\ c & d \end{pmatrix}\begin{pmatrix} a & b \\ c & d \end{pmatrix}=(a2+bcb(a+d)c(a+d)bc+d2)\displaystyle {}= \begin{pmatrix} a^2 + bc & b(a + d) \\ c(a + d) & bc + d^2 \end{pmatrix}

積の計算で、右上は abab+bd{}+ bd=b(a+d){}= b(a+d)、左下は caca+dc{}+ dc=c(a+d){}= c(a+d) とくくれるのがみそです。ここから (a+d)A(a+d)A を引きます。非対角成分は差し引きゼロ、対角成分は

a2\displaystyle a^2+bc\displaystyle {}+ bc−(a2+ad)\displaystyle {}- (a^2 + ad)=bc\displaystyle {}= bc−ad,\displaystyle {}- ad,bc\displaystyle bc+d2\displaystyle {}+ d^2−(ad+d2)\displaystyle {}- (ad + d^2)=bc\displaystyle {}= bc−ad\displaystyle {}- ad

なので

A2\displaystyle A^2−(a+d)A\displaystyle {}- (a + d)A=(bc−ad00bc−ad)\displaystyle {}= \begin{pmatrix} bc - ad & 0 \\ 0 & bc - ad \end{pmatrix}

対角成分はどちらも bcbc−ad{}- ad=−Δ{}= -\Delta。つまりこれは −ΔE-\Delta E です。移項すれば公式1 が出ます。

この定理のありがたみは、A2A^2 を「AA の何倍か」と「EE の何倍か」の足し算に書き直せるところにあります。2 次以上の項を 1 次以下に引きずり下ろせる、ということです。これを次数下げといいます。

似た経験は数と式の章でしています。α\alpha=1{}= 1+2{}+ \sqrt{2} は x2x^2−2x{}- 2x−1{}- 1=0{}= 0 の解なので α2\alpha^2=2α{}= 2\alpha+1{}+ 1 でした。だから α3\alpha^3 も α10\alpha^{10} も、計算すれば必ず「整数 ++ 整数 ×2\times \sqrt{2}」の形に戻ってきます。2\sqrt{2} が満たす 2 次方程式が、高い累乗を低いところへ押し戻しているのです。行列にも同じ仕組みが備わっている、というのがこの定理の中身です。

どんな 2 次の正方行列も A2A^2=(a+d)A{}= (a + d)A−ΔE{}- \Delta E という自分専用の関係式を持っていて、これを使えば A2A^2 を AA と EE の組み合わせに書き直せるということです。

例題1:定理の確認と A3A^3 の計算

AA=(3124){}= \begin{pmatrix} 3 & 1 \\ 2 & 4 \end{pmatrix} について、次の問いに答えなさい。

(1) ケーリー・ハミルトンの定理の式を書きなさい。

(2) (1) を使って A3A^3 を求めなさい。


【解答】

(1) aa+d{}+ d=3{}= 3+4{}+ 4=7{}= 7,Δ\Delta=3×4{}= 3 \times 4−1×2{}- 1 \times 2=10{}= 10 です。よって

A2‾\displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}A^2}−7A‾\displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}{}- 7A}+10E=O‾\displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}{}+ 10E = O}

(2) (1) より A2A^2=7A{}= 7A−10E{}- 10E です。両辺に右から AA をかけます。

A3\displaystyle A^3=7A2\displaystyle {}= 7A^2−10A\displaystyle {}- 10A

ここで右辺の A2A^2 に、もう一度 7A7A−10E{}- 10E を代入します。

A3\displaystyle A^3=7(7A−10E)\displaystyle {}= 7(7A - 10E)−10A\displaystyle {}- 10A=49A\displaystyle {}= 49A−70E\displaystyle {}- 70E−10A\displaystyle {}- 10A=39A\displaystyle {}= 39A−70E\displaystyle {}- 70E

AA と EE の組み合わせになりました。あとは成分を入れるだけです。

A3\displaystyle A^3=39(3124)\displaystyle {}= 39\begin{pmatrix} 3 & 1 \\ 2 & 4 \end{pmatrix}−70(1001)\displaystyle {}- 70\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}=(117−703978156−70)\displaystyle {}= \begin{pmatrix} 117 - 70 & 39 \\ 78 & 156 - 70 \end{pmatrix}=(47397886)‾\displaystyle {}= \underline{\begin{pmatrix} 47 & 39 \\ 78 & 86 \end{pmatrix}}

素直に A2A^2 を出してから AA をかけても同じ答えになりますが、こちらは行列のかけ算を一度もしていないことに注目してください。やったのは数の計算だけです。A10A^{10} や A100A^{100} になるほど、この差は大きくなります。

定理から逆行列も出てくる

公式1 をもう少し眺めてみます。A2A^2−(a+d)A{}- (a+d)A=−ΔE{}= -\Delta E の左辺は AA でくくれて

A{A−(a+d)E}\displaystyle A\{A - (a + d)E\}=−ΔE\displaystyle {}= -\Delta E

と書けます。両辺を −Δ-\Delta で割れば(Δ\Delta≠0{}\neq 0 のとき)

A⋅1Δ{(a+d)E−A}\displaystyle A \cdot \frac{1}{\Delta}\{(a + d)E - A\}=E\displaystyle {}= E

です。AA にかけて EE になる相手が見つかった、つまりこれが逆行列です。

公式2:定理から得られる逆行列

Δ\Delta≠0{}\neq 0 のとき

A−1\displaystyle A^{-1}=1Δ{(a+d)E−A}\displaystyle {}= \frac{1}{\Delta}\{(a + d)E - A\}

第 3 章で覚えた A−1A^{-1}=1Δ(d−b−ca){}= \dfrac{1}{\Delta}\begin{pmatrix} d & -b \\ -c & a \end{pmatrix} と同じものです。実際

(a+d)E\displaystyle (a + d)E−A\displaystyle {}- A=(a+d00a+d)\displaystyle {}= \begin{pmatrix} a + d & 0 \\ 0 & a + d \end{pmatrix}−(abcd)\displaystyle {}- \begin{pmatrix} a & b \\ c & d \end{pmatrix}=(d−b−ca)\displaystyle {}= \begin{pmatrix} d & -b \\ -c & a \end{pmatrix}

で、第 3 章 公式2 の逆行列の式に出てきた行列そのものが現れます。「aa と dd を入れかえ、bb と cc の符号を変える」というあの操作の正体は、(a+d)E(a+d)E から AA を引くことだったわけです。

例題1 の AA=(3124){}= \begin{pmatrix} 3 & 1 \\ 2 & 4 \end{pmatrix} で確かめます。Δ\Delta=10{}= 10,aa+d{}+ d=7{}= 7 なので

A−1\displaystyle A^{-1}=110{7(1001)−(3124)}\displaystyle {}= \frac{1}{10}\left\{7\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix} - \begin{pmatrix} 3 & 1 \\ 2 & 4 \end{pmatrix}\right\}=110(4−1−23)\displaystyle {}= \frac{1}{10}\begin{pmatrix} 4 & -1 \\ -2 & 3 \end{pmatrix}

第 3 章の公式で求めても同じです。おかげで A−1A^{-1} もまた「AA の何倍か ++ EE の何倍か」の形に書けることがわかりました。負の指数まで含めて、AA の累乗はすべて pApA+qE{}+ qE の形に収まる、ということです。

定理の逆は成り立たない

ここがこの章でいちばん間違えやすいところです。公式1 は「AA を決めれば aa+d{}+ d と Δ\Delta が決まり、その係数で式が成り立つ」と言っています。では逆に、A2A^2−5A{}- 5A+6E{}+ 6E=O{}= O という式が成り立っていたら、aa+d{}+ d=5{}= 5,Δ\Delta=6{}= 6 だと言えるでしょうか。

言えません。 反例はすぐ作れます。AA=2E{}= 2E=(2002){}= \begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix} とすると A2A^2=4E{}= 4E なので

A2\displaystyle A^2−5A\displaystyle {}- 5A+6E\displaystyle {}+ 6E=4E\displaystyle {}= 4E−10E\displaystyle {}- 10E+6E\displaystyle {}+ 6E=O\displaystyle {}= O

たしかに成り立ちますが、この AA の対角和は 22+2{}+ 2=4{}= 4 で、55 ではありません。Δ\Delta=4{}= 4 で、66 でもありません。

種明かしをすると、AA=kE{}= kE の形(数の化身、第 2 章)だけが特別なのです。AA=kE{}= kE なら A2A^2−pA{}- pA+qE{}+ qE=(k2−pk+q)E{}= (k^2 - pk + q)E なので、kk が 2 次方程式 x2x^2−px{}- px+q{}+ q=0{}= 0 の解でありさえすれば、pp と qq がどんな値でも式が成り立ってしまいます。AA 自身の対角和とも行列式とも無関係に、いくらでも式が作れるわけです。

公式3:定理の逆(成り立つのは 1 つの場合を除いたとき)

実数 pp,qq について A2A^2−pA{}- pA+qE{}+ qE=O{}= O が成り立つとき、次のどちらかが起こる。

  • AA=kE{}= kE(kk は実数)の形である
  • AA が kEkE の形でないならば、pp=a{}= a+d{}+ d かつ qq=Δ{}= \Delta である

理由 公式1 の A2A^2−(a+d)A{}- (a + d)A+ΔE{}+ \Delta E=O{}= O から、与えられた式を辺々引くと

{p−(a+d)}A\displaystyle \{p - (a + d)\}A=(q−Δ)E\displaystyle {}= (q - \Delta)E

pp−(a+d){}- (a + d)≠0{}\neq 0 ならば両辺をこれで割って AA=kE{}= kE の形になる。AA が kEkE の形でなければ pp−(a+d){}- (a + d)=0{}= 0、したがって右辺も OO となり qq=Δ{}= \Delta である。

答案では「AA≠kE{}\neq kE のとき」と「AA=kE{}= kE のとき」に分けて書きます。どちらか一方しか書かないと減点されるのが、この手の問題のお約束です。

A2A^2−pA{}- pA+qE{}+ qE=O{}= O という式を見ても、すぐに pp=a{}= a+d{}+ d とはいえず、AA が単位行列の実数倍でないことを確かめてから初めていえるということです。

例題2:定理の逆を使う場合分け

2 次の正方行列 AA が A2A^2−3A{}- 3A+2E{}+ 2E=O{}= O を満たすとき、AA の対角和 aa+d{}+ d と行列式 Δ\Delta の組をすべて求めなさい。


【解答】

(i) AA が kEkE の形でないとき。 公式3 より a‾\underline{\rule[-0.1944em]{0em}{0.8888em}a}+d=3,‾\underline{\rule[-0.1944em]{0em}{0.8888em}{}+ d = 3,} Δ=2‾\underline{\rule[-0.1944em]{0em}{0.8888em}\ \Delta = 2} です。

(ii) AA=kE{}= kE のとき。 A2A^2=k2E{}= k^2E なので

(k2−3k+2)E\displaystyle (k^2 - 3k + 2)E=O\displaystyle {}= O⟹k2\displaystyle {}\Longrightarrow \quad k^2−3k\displaystyle {}- 3k+2\displaystyle {}+ 2=0\displaystyle {}= 0⟹k\displaystyle {}\Longrightarrow \quad k=1,\displaystyle {}= 1, 2\displaystyle \ 2

kk=1{}= 1 のとき AA=E{}= E で、aa+d{}+ d=2{}= 2,Δ\Delta=1{}= 1。

kk=2{}= 2 のとき AA=2E{}= 2E で、aa+d{}+ d=4{}= 4,Δ\Delta=4{}= 4。

以上より、組は (a+d, Δ)=(3, 2),‾\underline{\rule[-0.25em]{0em}{1.0000em}(a+d,\ \Delta) = (3,\ 2),} (2, 1),‾\underline{\rule[-0.25em]{0em}{1.0000em}\ (2,\ 1),} (4, 4)‾\underline{\rule[-0.25em]{0em}{1.0000em}\ (4,\ 4)} の 3 通りです。

(ii) で出てきた kk=1,{}= 1, 2\ 2 が、x2x^2−3x{}- 3x+2{}+ 2=0{}= 0 の解になっていることに注意してください。AA=kE{}= kE のときは、行列の式がまるごと数の 2 次方程式に化けてしまうのです。行列と数の境目がいちばん薄くなる場所が、kEkE という形だといえます。

次数下げで高い累乗を下ろす

道具がそろったので、本題の AnA^n に入ります。公式1 を繰り返し使えば、A3A^3 も A4A^4 も pApA+qE{}+ qE の形に落ちていきます。

公式4:次数下げ

A2A^2=(a+d)A{}= (a + d)A−ΔE{}- \Delta E を使うと、どんな正の整数 nn に対しても

An\displaystyle A^n=pA\displaystyle {}= pA+qE\displaystyle {}+ qE(p, q は実数)\displaystyle (p,\ q \text{ は実数})

の形に書ける。An+1A^{n+1}=An⋅A{}= A^n \cdot A=pA2{}= pA^2+qA{}+ qA の A2A^2 をふたたび下ろす、という手順を繰り返せばよい。

AA の累乗は、何乗であっても「AA の実数倍 ++ EE の実数倍」という 1 次の形に必ず収まるということです。

例題3:次数下げで A3A^3,A4A^4 を求める

AA=(3−120){}= \begin{pmatrix} 3 & -1 \\ 2 & 0 \end{pmatrix} について、A3A^3 と A4A^4 を求めなさい。


【解答】

aa+d{}+ d=3{}= 3+0{}+ 0=3{}= 3,Δ\Delta=3×0{}= 3 \times 0−(−1)×2{}- (-1) \times 2=2{}= 2 なので

A2\displaystyle A^2=3A\displaystyle {}= 3A−2E\displaystyle {}- 2E

両辺に AA をかけ、出てきた A2A^2 をまた下ろします。

A3\displaystyle A^3=3A2\displaystyle {}= 3A^2−2A\displaystyle {}- 2A=3(3A−2E)\displaystyle {}= 3(3A - 2E)−2A\displaystyle {}- 2A=7A\displaystyle {}= 7A−6E\displaystyle {}- 6EA4\displaystyle A^4=7A2\displaystyle {}= 7A^2−6A\displaystyle {}- 6A=7(3A−2E)\displaystyle {}= 7(3A - 2E)−6A\displaystyle {}- 6A=15A\displaystyle {}= 15A−14E\displaystyle {}- 14E

成分に直します。

A3\displaystyle A^3=7(3−120)\displaystyle {}= 7\begin{pmatrix} 3 & -1 \\ 2 & 0 \end{pmatrix}−6E\displaystyle {}- 6E=(15−714−6)‾\displaystyle {}= \underline{\begin{pmatrix} 15 & -7 \\ 14 & -6 \end{pmatrix}}A4\displaystyle A^4=15(3−120)\displaystyle {}= 15\begin{pmatrix} 3 & -1 \\ 2 & 0 \end{pmatrix}−14E\displaystyle {}- 14E=(31−1530−14)‾\displaystyle {}= \underline{\begin{pmatrix} 31 & -15 \\ 30 & -14 \end{pmatrix}}

係数の並びを見ておきましょう。A2A^2 から順に (3, −2)(3,\ -2),(7, −6)(7,\ -6),(15, −14)(15,\ -14)。AA の係数は 3,3, 7,\ 7, 15\ 15、つまり 2n2^n−1{}- 1 です。EE の係数は −2,-2, −6,{}\ -6, −14{}\ -14 で −(2n−2)-(2^n - 2)。予想するなら

An\displaystyle A^n=(2n−1)A\displaystyle {}= (2^n - 1)A−(2n−2)E\displaystyle {}- (2^n - 2)E

ですが、この予想を証明するには帰納法が要ります。次の節では、予想も帰納法も使わずに pp と qq を直接求める方法を身につけます。

割り算の余りで一気に求める

例題3 の手順は確実ですが、nn=100{}= 100 では 99 回繰り返すことになります。もっと直接に pp と qq を出したい。そこで、多項式の割り算を使います。

xnx^n を 2 次式 x2x^2−(a+d)x{}- (a+d)x+Δ{}+ \Delta で割ったときの商を Q(x)Q(x)、余りを pxpx+q{}+ q(1 次以下)とすると

xn\displaystyle x^n={x2−(a+d)x+Δ}Q(x)\displaystyle {}= \{x^2 - (a + d)x + \Delta\}Q(x)+px\displaystyle {}+ px+q\displaystyle {}+ q

この xx に行列 AA を代入します。ここで効くのが公式1 です。{ }\{\ \} の中は AA を入れると OO になるので、第 1 項がまるごと消えて

An\displaystyle A^n=pA\displaystyle {}= pA+qE\displaystyle {}+ qE

が残ります。つまり余りの係数がそのまま答えです。しかも pp,qq を求めるのは、数の計算だけでできます。x2x^2−(a+d)x{}- (a+d)x+Δ{}+ \Delta=0{}= 0 の解を代入すれば、恒等式の第 1 項が 00 になるからです。

整数の世界でも同じことをしています。71007^{100} を 55 で割った余りを知りたいとき、71007^{100} を全部書き下す人はいません。77≡2(mod5){}\equiv 2 \pmod 5 から余りだけを追いかけます(整数 第 5 章)。ここでも、知りたいのは xnx^n そのものではなく、x2x^2−(a+d)x{}- (a+d)x+Δ{}+ \Delta で割った余りだけ。全体を計算せず、必要な部分だけ取り出す発想です。

公式5:AnA^n を割り算の余りで求める

AA の対角和を aa+d{}+ d、行列式を Δ\Delta とし、2 次方程式

x2\displaystyle x^2−(a+d)x\displaystyle {}- (a + d)x+Δ\displaystyle {}+ \Delta=0\displaystyle {}= 0

の 2 つの解を α\alpha,β\beta とする。xnx^n をこの 2 次式で割った余りを pxpx+q{}+ q とすると

An\displaystyle A^n=pA\displaystyle {}= pA+qE\displaystyle {}+ qE

である。α\alpha≠β{}\neq \beta のとき、pp,qq は連立方程式

{αn=pα+qβn=pβ+q\begin{cases} \alpha^n = p\alpha + q \\ \beta^n = p\beta + q \end{cases}

を解いて得られる。

x2x^2−(a+d)x{}- (a+d)x+Δ{}+ \Delta=0{}= 0 という方程式は、この章のあいだ何度も顔を出します。AA が背負っている 2 次方程式そのものです。

AnA^n を求めるには、xnx^n を AA の背負う 2 次式で割った余り pxpx+q{}+ q を求めればよく、その pp,qq は 2 次方程式の解を代入した連立方程式から出るということです。

例題4:AnA^n を求める(解が 2 つあるとき)

AA=(1232){}= \begin{pmatrix} 1 & 2 \\ 3 & 2 \end{pmatrix} について、AnA^n を求めなさい。ただし nn は正の整数とする。


【解答】

aa+d{}+ d=1{}= 1+2{}+ 2=3{}= 3,Δ\Delta=1×2{}= 1 \times 2−2×3{}- 2 \times 3=−4{}= -4 なので、AA が背負う 2 次方程式は

x2\displaystyle x^2−3x\displaystyle {}- 3x−4\displaystyle {}- 4=0\displaystyle {}= 0⟹(x−4)(x+1)\displaystyle {}\Longrightarrow \quad (x - 4)(x + 1)=0\displaystyle {}= 0⟹x\displaystyle {}\Longrightarrow \quad x=4,\displaystyle {}= 4, −1\displaystyle {}\ -1

xnx^n を x2x^2−3x{}- 3x−4{}- 4 で割った商を Q(x)Q(x)、余りを pxpx+q{}+ q とします。

xn\displaystyle x^n=(x2−3x−4)Q(x)\displaystyle {}= (x^2 - 3x - 4)Q(x)+px\displaystyle {}+ px+q\displaystyle {}+ q

xx=4{}= 4 と xx=−1{}= -1 を代入すると、どちらも第 1 項が 00 になります。

{4n=4p+q(−1)n=−p+q\begin{cases} 4^n = 4p + q \\ (-1)^n = -p + q \end{cases}

辺々引いて 4n4^n−(−1)n{}- (-1)^n=5p{}= 5p より

p\displaystyle p=4n−(−1)n5,\displaystyle {}= \frac{4^n - (-1)^n}{5},q\displaystyle q=(−1)n\displaystyle {}= (-1)^n+p\displaystyle {}+ p=4n+4⋅(−1)n5\displaystyle {}= \frac{4^n + 4 \cdot (-1)^n}{5}

ケーリー・ハミルトンの定理より A2A^2−3A{}- 3A−4E{}- 4E=O{}= O なので、上の恒等式に AA を代入すると第 1 項が消えて AnA^n=pA{}= pA+qE{}+ qE です。

An\displaystyle A^n=p(1232)\displaystyle {}= p\begin{pmatrix} 1 & 2 \\ 3 & 2 \end{pmatrix}+q(1001)\displaystyle {}+ q\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}=(p+q2p3p2p+q)\displaystyle {}= \begin{pmatrix} p + q & 2p \\ 3p & 2p + q \end{pmatrix}

pp+q{}+ q=2⋅4n+3(−1)n5{}= \dfrac{2 \cdot 4^n + 3(-1)^n}{5},2p2p+q{}+ q=3⋅4n+2(−1)n5{}= \dfrac{3 \cdot 4^n + 2(-1)^n}{5} を使います。成分が長くなるので、4n4^n=u{}= u,(−1)n(-1)^n=v{}= v と書くと

An\displaystyle A^n=15(2u+3v2(u−v)3(u−v)3u+2v)‾\displaystyle {}= \underline{\frac{1}{5}\begin{pmatrix} 2u + 3v & 2(u - v) \\ 3(u - v) & 3u + 2v \end{pmatrix}}

ただし uu=4n{}= 4^n,vv=(−1)n{}= (-1)^n です。

検算します。nn=1{}= 1 なら uu=4{}= 4,vv=−1{}= -1 で 15(8−32×53×512−2)\dfrac{1}{5}\begin{pmatrix} 8 - 3 & 2 \times 5 \\ 3 \times 5 & 12 - 2 \end{pmatrix}=(1232){}= \begin{pmatrix} 1 & 2 \\ 3 & 2 \end{pmatrix} で AA に戻ります。nn=2{}= 2 なら 15(35304550)\dfrac{1}{5}\begin{pmatrix} 35 & 30 \\ 45 & 50 \end{pmatrix}=(76910){}= \begin{pmatrix} 7 & 6 \\ 9 & 10 \end{pmatrix} で、直接 A2A^2 を計算した値と一致します。nn=1{}= 1 と nn=2{}= 2 で確かめるのが、この種の答えの検算の定石です。

解が重なるとき

公式5 の連立方程式は、α\alpha≠β{}\neq \beta でなければ作れません。解が重なると式が 1 本しか立たないからです。この場合は別の道を通ります。

2 次方程式が (x−α)2(x - \alpha)^2=0{}= 0 の形なら、ケーリー・ハミルトンの定理は

(A−αE)2\displaystyle (A - \alpha E)^2=O\displaystyle {}= O

と書けます。2 乗すると消えてしまう行列が現れるわけです。そこで NN=A{}= A−αE{}- \alpha E とおけば AA=αE{}= \alpha E+N{}+ N で、N2N^2=O{}= O。

αE\alpha E はどんな行列とも交換可能(第 2 章)なので、(αE+N)n(\alpha E + N)^n には二項定理がそのまま使えます。しかも N2N^2 から先はすべて OO なので、展開しても 2 項しか生き残りません。

公式6:解が重なるときの AnA^n

x2x^2−(a+d)x{}- (a + d)x+Δ{}+ \Delta=0{}= 0 が重解 α\alpha をもつとき、NN=A{}= A−αE{}- \alpha E とおくと N2N^2=O{}= O であり

An\displaystyle A^n=(αE+N)n\displaystyle {}= (\alpha E + N)^n=αnE\displaystyle {}= \alpha^n E+nαn−1N\displaystyle {}+ n\alpha^{n-1}N

が成り立つ(nn≧2{}\geqq 2 のとき N2N^2=N3{}= N^3=⋯{}= \cdots=O{}= O だから、二項定理の第 3 項以降は消える)。

2 次方程式が重解をもつときは AA から解の分だけ EE を引くと 2 乗して消える行列になり、二項定理の最初の 2 項だけで AnA^n が書けるということです。

例題5:AnA^n を求める(解が重なるとき)

AA=(41−12){}= \begin{pmatrix} 4 & 1 \\ -1 & 2 \end{pmatrix} について、AnA^n を求めなさい。ただし nn は正の整数とする。


【解答】

aa+d{}+ d=4{}= 4+2{}+ 2=6{}= 6,Δ\Delta=4×2{}= 4 \times 2−1×(−1){}- 1 \times (-1)=9{}= 9 なので

x2\displaystyle x^2−6x\displaystyle {}- 6x+9\displaystyle {}+ 9=0\displaystyle {}= 0⟹(x−3)2\displaystyle {}\Longrightarrow \quad (x - 3)^2=0\displaystyle {}= 0

重解 xx=3{}= 3 です。そこで NN=A{}= A−3E{}- 3E とおきます。

N\displaystyle N=(4−31−12−3)\displaystyle {}= \begin{pmatrix} 4 - 3 & 1 \\ -1 & 2 - 3 \end{pmatrix}=(11−1−1),\displaystyle {}= \begin{pmatrix} 1 & 1 \\ -1 & -1 \end{pmatrix},N2\displaystyle N^2=(1−11−1−1+1−1+1)\displaystyle {}= \begin{pmatrix} 1 - 1 & 1 - 1 \\ -1 + 1 & -1 + 1 \end{pmatrix}=O\displaystyle {}= O

たしかに N2N^2=O{}= O です。AA=3E{}= 3E+N{}+ N で、3E3E と NN は交換可能なので二項定理が使えます。

An\displaystyle A^n=(3E+N)n\displaystyle {}= (3E + N)^n=3nE\displaystyle {}= 3^nE+n⋅3n−1N\displaystyle {}+ n \cdot 3^{n-1}N+nC2⋅3n−2N2\displaystyle {}+ {}_n\mathrm{C}_2 \cdot 3^{n-2}N^2+⋯\displaystyle {}+ \cdots=3nE\displaystyle {}= 3^nE+n⋅3n−1N\displaystyle {}+ n \cdot 3^{n-1}N

第 3 項から先は N2N^2 を含むのですべて OO です。成分に直します。

An\displaystyle A^n=3n(1001)\displaystyle {}= 3^n\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}+n⋅3n−1(11−1−1)\displaystyle {}+ n \cdot 3^{n-1}\begin{pmatrix} 1 & 1 \\ -1 & -1 \end{pmatrix}=3n−1(n+3n−n3−n)‾\displaystyle {}= \underline{3^{n-1}\begin{pmatrix} n + 3 & n \\ -n & 3 - n \end{pmatrix}}

nn=1{}= 1 のとき 30(41−12)3^0\begin{pmatrix} 4 & 1 \\ -1 & 2 \end{pmatrix}=A{}= A、nn=2{}= 2 のとき 3(52−21)3\begin{pmatrix} 5 & 2 \\ -2 & 1 \end{pmatrix}=(156−63){}= \begin{pmatrix} 15 & 6 \\ -6 & 3 \end{pmatrix} で、A2A^2=6A{}= 6A−9E{}- 9E から求めた値と合っています。

N2N^2=O{}= O なのに NN≠O{}\neq O という点は、行列ならではです。数なら x2x^2=0{}= 0 から xx=0{}= 0 が出ますが、行列では出ません。零因子(第 2 章)の話が、こんなところで役に立っています。

対角行列に化ける形が与えられたとき

もう 1 つ、入試でよく出る形を見ておきます。うまい行列 PP を持ってきて P−1APP^{-1}AP を作ると、対角成分だけが残った行列になることがあります。

P−1AP\displaystyle P^{-1}AP=(α00β)\displaystyle {}= \begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix}

対角成分以外が 00 の行列を対角行列といいます。対角行列の累乗はとても簡単で、成分をそれぞれ nn 乗するだけです。

(α00β) n\displaystyle \begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix}^{\,n}=(αn00βn)\displaystyle {}= \begin{pmatrix} \alpha^n & 0 \\ 0 & \beta^n \end{pmatrix}

この簡単さを AA 側に持ち帰れないでしょうか。持ち帰れます。DD=P−1AP{}= P^{-1}AP とおくと

Dn\displaystyle D^n=(P−1AP)(P−1AP)⋯(P−1AP)\displaystyle {}= (P^{-1}AP)(P^{-1}AP)\cdots(P^{-1}AP)

で、となり合う PP−1PP^{-1} がすべて EE になって消え、まん中の AA だけが nn 個残ります。

Dn\displaystyle D^n=P−1AnP\displaystyle {}= P^{-1}A^nP

左から PP、右から P−1P^{-1} をかければ AnA^n が取り出せます。

公式7:P−1APP^{-1}AP が対角行列のとき

PP が正則で P−1APP^{-1}AP=D{}= D=(α00β){}= \begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix} ならば

An\displaystyle A^n=PDnP−1\displaystyle {}= PD^nP^{-1}=P(αn00βn)P−1\displaystyle {}= P\begin{pmatrix} \alpha^n & 0 \\ 0 & \beta^n \end{pmatrix}P^{-1}

P−1APP^{-1}AP という形は、第 3 章の実践問題で一度出ています。det⁡(P−1AP)\det(P^{-1}AP)=det⁡A{}= \det A、つまり PP で挟んでも行列式は変わらない、という問題でした。行列式だけでなく対角和も変わりません。AA と DD は、同じものを別の見方で書いたものだからです。

例えるなら、方眼紙の目盛りを引き直す作業です。ふつうの xyxy 軸で見ると AA は「伸ばしながらねじる」複雑な操作に見えるのに、PP の 2 つの列が指す 2 方向を新しい軸にとると、それぞれの方向に α\alpha 倍・β\beta 倍するだけの単純な操作になっている。目盛りの取り替えが PP、取り替えた世界での姿が DD です。この見方は、第 6 章以降の一次変換でもう一度出てきます。

P−1APP^{-1}AP が対角行列になるなら、AnA^n は「対角行列を nn 乗して PP で挟み直す」だけで求まるということです。

例題6:対角行列を経由して AnA^n を求める

AA=(4213){}= \begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix},PP=(211−1){}= \begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix} とする。P−1APP^{-1}AP を求め、それを使って AnA^n を求めなさい。


【解答】

まず P−1P^{-1} を求めます。PP の行列式は 2×(−1)2 \times (-1)−1×1{}- 1 \times 1=−3{}= -3≠0{}\neq 0 なので

P−1\displaystyle P^{-1}=1−3(−1−1−12)\displaystyle {}= \frac{1}{-3}\begin{pmatrix} -1 & -1 \\ -1 & 2 \end{pmatrix}=13(111−2)\displaystyle {}= \frac{1}{3}\begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix}

APAP を先に計算します。

AP\displaystyle AP=(4213)(211−1)\displaystyle {}= \begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}\begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix}=(1025−2)\displaystyle {}= \begin{pmatrix} 10 & 2 \\ 5 & -2 \end{pmatrix}P−1(AP)\displaystyle P^{-1}(AP)=13(111−2)(1025−2)\displaystyle {}= \frac{1}{3}\begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix}\begin{pmatrix} 10 & 2 \\ 5 & -2 \end{pmatrix}=13(15006)\displaystyle {}= \frac{1}{3}\begin{pmatrix} 15 & 0 \\ 0 & 6 \end{pmatrix}=(5002)‾\displaystyle {}= \underline{\begin{pmatrix} 5 & 0 \\ 0 & 2 \end{pmatrix}}

対角行列になりました。DD=(5002){}= \begin{pmatrix} 5 & 0 \\ 0 & 2 \end{pmatrix} とおくと DnD^n=(5n002n){}= \begin{pmatrix} 5^n & 0 \\ 0 & 2^n \end{pmatrix} なので、公式7 より

An\displaystyle A^n=PDnP−1\displaystyle {}= PD^nP^{-1}=(211−1)(5n002n)\displaystyle {}= \begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix}\begin{pmatrix} 5^n & 0 \\ 0 & 2^n \end{pmatrix}⋅13(111−2)\displaystyle \qquad \cdot \frac{1}{3}\begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix}

左の 2 つを先に計算します。

(211−1)(5n002n)\displaystyle \begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix}\begin{pmatrix} 5^n & 0 \\ 0 & 2^n \end{pmatrix}=(2⋅5n2n5n−2n)\displaystyle {}= \begin{pmatrix} 2 \cdot 5^n & 2^n \\ 5^n & -2^n \end{pmatrix}An\displaystyle A^n=13(2⋅5n2n5n−2n)(111−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}=13(2⋅5n+2n2⋅5n−2⋅2n5n−2n5n+2⋅2n)‾\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}}

nn=1{}= 1 で 13(12639)\dfrac{1}{3}\begin{pmatrix} 12 & 6 \\ 3 & 9 \end{pmatrix}=(4213){}= \begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}、nn=2{}= 2 で 13(54422133)\dfrac{1}{3}\begin{pmatrix} 54 & 42 \\ 21 & 33 \end{pmatrix}=(1814711){}= \begin{pmatrix} 18 & 14 \\ 7 & 11 \end{pmatrix}。どちらも直接計算した値に一致します。

なお DD の対角成分 55,22 は、AA の背負う 2 次方程式 x2x^2−7x{}- 7x+10{}+ 10=0{}= 0 の 2 つの解です。偶然ではありません。PP で挟んでも対角和と行列式は変わらないので、DD も同じ 2 次方程式を背負っているからです。公式5 で α\alpha,β\beta と呼んでいたものの正体が、ここで見えます。

どの方法を使うか

AnA^n の求め方が出そろいました。整理しておきます。

方法使いどころ
A2A^2,A3A^3 を計算して予想し、帰納法成分に 00 が多く、規則が目で見えるとき(第 2 章 例題7)
周期を見つけて nn を割った余りで判断A3=EA^3 = E,A4=EA^4 = E のように単位行列に戻るとき(第 2 章 j13)
次数下げを繰り返すA3A^3,A4A^4 など、具体的な低い累乗を求めるとき
割り算の余り(公式5・公式6)一般の nn で AnA^n を求めるとき。いちばん汎用的
P−1APP^{-1}AP が対角行列(公式7)問題文に PP が与えられているとき

最後に、AnA^n が実際に何かの役に立つ場面を見ておきます。

A 店 B 店 3/4 1/2 1/4 1/2 翌日 B 店へ 翌日 A 店へ
例題7:何日後にどちらの店にいるか

ある人は毎日、昼食を A 店か B 店のどちらかでとる。A 店で食べた翌日にまた A 店を選ぶ確率は 34\dfrac{3}{4}、B 店で食べた翌日に A 店を選ぶ確率は 12\dfrac{1}{2} である。1 日目は A 店で食べた。nn 日目に A 店で食べる確率 pnp_n を求め、nn を限りなく大きくしたときの値を調べなさい。


【解答】

nn 日目に A 店で食べる確率を pnp_n、B 店で食べる確率を qnq_n とします。nn+1{}+ 1 日目に A 店にいるのは「A 店にいて残る」か「B 店にいて移る」かのどちらかなので

{pn+1=34pn+12qnqn+1=14pn+12qn\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}すなわち\displaystyle \text{すなわち}(pn+1qn+1)\displaystyle \begin{pmatrix} p_{n+1} \\ q_{n+1} \end{pmatrix}=A(pnqn),\displaystyle {}= A\begin{pmatrix} p_n \\ q_n \end{pmatrix},A\displaystyle A=(34121412)\displaystyle {}= \begin{pmatrix} \dfrac{3}{4} & \dfrac{1}{2} \\[4pt] \dfrac{1}{4} & \dfrac{1}{2} \end{pmatrix}

1 日目は A 店なので (p1q1)\begin{pmatrix} p_1 \\ q_1 \end{pmatrix}=(10){}= \begin{pmatrix} 1 \\ 0 \end{pmatrix} です。nn 日目の確率は、これに AA を nn−1{}- 1 回かけたもの、つまり

(pnqn)\displaystyle \begin{pmatrix} p_n \\ q_n \end{pmatrix}=An−1(10)\displaystyle {}= A^{n-1}\begin{pmatrix} 1 \\ 0 \end{pmatrix}

An−1A^{n-1} を求めます。aa+d{}+ d=34{}= \dfrac{3}{4}+12{}+ \dfrac{1}{2}=54{}= \dfrac{5}{4},Δ\Delta=34⋅12{}= \dfrac{3}{4} \cdot \dfrac{1}{2}−12⋅14{}- \dfrac{1}{2} \cdot \dfrac{1}{4}=14{}= \dfrac{1}{4} なので、背負う 2 次方程式は

x2\displaystyle x^2−54x\displaystyle {}- \frac{5}{4}x+14\displaystyle {}+ \frac{1}{4}=0\displaystyle {}= 0⟹4x2\displaystyle {}\Longrightarrow \quad 4x^2−5x\displaystyle {}- 5x+1\displaystyle {}+ 1=0\displaystyle {}= 0⟹(4x−1)(x−1)\displaystyle {}\Longrightarrow \quad (4x - 1)(x - 1)=0\displaystyle {}= 0

解は xx=1,{}= 1, 14\ \dfrac{1}{4} です。mm=n{}= n−1{}- 1 とおき、xmx^m をこの 2 次式で割った余りを pxpx+q{}+ q とすると

{1=p+q(14)m=14p+q\begin{cases} 1 = p + q \\[2pt] \left(\dfrac{1}{4}\right)^{m} = \dfrac{1}{4}p + q \end{cases}

辺々引いて 11−(14)m{}- \left(\dfrac{1}{4}\right)^{m}=34p{}= \dfrac{3}{4}p より pp=43{1−(14)m}{}= \dfrac{4}{3}\left\{1 - \left(\dfrac{1}{4}\right)^{m}\right\}、qq=1{}= 1−p{}- p です。

AmA^m=pA{}= pA+qE{}+ qE なので、(10)\begin{pmatrix} 1 \\ 0 \end{pmatrix} にかけると

Am(10)\displaystyle A^m\begin{pmatrix} 1 \\ 0 \end{pmatrix}=p(3/41/4)\displaystyle {}= p\begin{pmatrix} 3/4 \\ 1/4 \end{pmatrix}+q(10)\displaystyle {}+ q\begin{pmatrix} 1 \\ 0 \end{pmatrix}=(34p+q14p)\displaystyle {}= \begin{pmatrix} \dfrac{3}{4}p + q \\[4pt] \dfrac{1}{4}p \end{pmatrix}

上の成分が pnp_n です。qq=1{}= 1−p{}- p を使って

pn\displaystyle p_n=34p\displaystyle {}= \frac{3}{4}p+(1−p)\displaystyle {}+ (1 - p)=1\displaystyle {}= 1−14p\displaystyle {}- \frac{1}{4}p=1\displaystyle {}= 1−13{1−(14)m}\displaystyle {}- \frac{1}{3}\left\{1 - \left(\frac{1}{4}\right)^{m}\right\}=23‾\displaystyle {}= \underline{\rule[-0.95em]{0em}{2.6040em}\frac{2}{3}}+13(14)n−1‾\displaystyle \underline{\rule[-0.95em]{0em}{2.6040em}{}+ \frac{1}{3}\left(\frac{1}{4}\right)^{n-1}}

検算します。nn=1{}= 1 なら 23\dfrac{2}{3}+13{}+ \dfrac{1}{3}=1{}= 1 で、初日は必ず A 店。nn=2{}= 2 なら 23\dfrac{2}{3}+112{}+ \dfrac{1}{12}=34{}= \dfrac{3}{4} で、問題の設定どおりです。

(14)n−1\left(\dfrac{1}{4}\right)^{n-1} は nn が大きくなると 00 に近づくので

pn\displaystyle p_n⟶23\displaystyle {}\longrightarrow \frac{2}{3}

日がたつにつれ、A 店にいる確率は 23\dfrac{2}{3} に落ち着きます。おもしろいのは、初日にどちらの店にいたかが効かなくなることです。q1q_1=1{}= 1(初日は B 店)から始めても、同じ 23\dfrac{2}{3} に近づきます。AnA^n の姿が nn とともに一定の形へ寄っていくためで、この落ち着き先を確率の分野では定常状態と呼びました(数列 第 9 章)。小話では、この考え方が Web 検索の順位づけに使われている話を取り上げます。

ここまでで、行列の計算としての道具はひととおりそろいました。次の第 6 章からは主役を交代させ、行列を平面上の点を動かす仕掛けとして見直します。A(xy)A\begin{pmatrix} x \\ y \end{pmatrix} を「点 (x, y)(x,\ y) の像」と呼んでいたあの見方を正式に一次変換と名づけ、記号 ff で表します。この章で身につけた AnA^n は、そこでは「同じ移動を nn 回くり返す」という意味を持つことになります。

基礎確認問題(全5問)

まずは公式をそのまま使う、ごく簡単な問題で確認しましょう。

問1

AA=(2513){}= \begin{pmatrix} 2 & 5 \\ 1 & 3 \end{pmatrix} について、ケーリー・ハミルトンの定理の式を求めなさい。

答えを見る
答え

aa+d{}+ d=5{}= 5,Δ\Delta=6{}= 6−5{}- 5=1{}= 1 より A2‾\underline{\rule[-0.0833em]{0em}{0.8974em}A^2}−5A‾\underline{\rule[-0.0833em]{0em}{0.8974em}{}- 5A}+E=O‾\underline{\rule[-0.0833em]{0em}{0.8974em}{}+ E = O}

自己採点:
記録を読み込み中…

問2

k01 の AA について、A2A^2 を定理を使って求めなさい。

つまずいたときは:
答えを見る
答え

A2A^2=5A{}= 5A−E{}- E=(925514)‾{}= \underline{\begin{pmatrix} 9 & 25 \\ 5 & 14 \end{pmatrix}}

自己採点:
記録を読み込み中…

問3

k01 の AA について、A−1A^{-1}=1Δ{(a+d)E−A}{}= \dfrac{1}{\Delta}\{(a+d)E - A\} を使って逆行列を求めなさい。

答えを見る
答え

A−1A^{-1}=5E{}= 5E−A{}- A=(3−5−12)‾{}= \underline{\begin{pmatrix} 3 & -5 \\ -1 & 2 \end{pmatrix}}

自己採点:
記録を読み込み中…

問4

AA=(1122){}= \begin{pmatrix} 1 & 1 \\ 2 & 2 \end{pmatrix} について、A2A^2 を AA の実数倍の形で表しなさい。

答えを見る
答え

aa+d{}+ d=3{}= 3,Δ\Delta=2{}= 2−2{}- 2=0{}= 0 より A2=3A‾\underline{A^2 = 3A}

自己採点:
記録を読み込み中…

問5

AA=(2−110){}= \begin{pmatrix} 2 & -1 \\ 1 & 0 \end{pmatrix} について、次数下げを使って A3A^3 を求めなさい。

つまずいたときは:
答えを見る
答え

A2A^2=2A{}= 2A−E{}- E より A3A^3=2A2{}= 2A^2−A{}- A=3A{}= 3A−2E{}- 2E=(4−33−2)‾{}= \underline{\begin{pmatrix} 4 & -3 \\ 3 & -2 \end{pmatrix}}

自己採点:
記録を読み込み中…

実践問題(全20問)

難易度マークは ★=基礎、★★=標準、★★★=入試レベルです。★から順に取り組みましょう。

問1 ★

次の 2 つの行列 (1) AA=(5234){}= \begin{pmatrix} 5 & 2 \\ 3 & 4 \end{pmatrix}、(2) BB=(−1326){}= \begin{pmatrix} -1 & 3 \\ 2 & 6 \end{pmatrix} について、ケーリー・ハミルトンの定理の式をそれぞれ求めなさい。

つまずいたときは:
答えを見る
答え

(1) A2A^2−9A{}- 9A+14E{}+ 14E=O{}= O (2) B2B^2−5B{}- 5B−12E{}- 12E=O{}= O

解説

対角和 aa+d{}+ d と行列式 Δ\Delta=ad{}= ad−bc{}- bc を求め、A2A^2−(a+d)A{}- (a+d)A+ΔE{}+ \Delta E=O{}= O に当てはめるだけです。

(1) aa+d{}+ d=5{}= 5+4{}+ 4=9{}= 9、Δ\Delta=5×4{}= 5 \times 4−2×3{}- 2 \times 3=20{}= 20−6{}- 6=14{}= 14 なので

A2‾\displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}A^2}−9A‾\displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}{}- 9A}+14E=O‾\displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}{}+ 14E = O}

(2) aa+d{}+ d=−1{}= -1+6{}+ 6=5{}= 5、Δ\Delta=(−1)×6{}= (-1) \times 6−3×2{}- 3 \times 2=−6{}= -6−6{}- 6=−12{}= -12 なので

B2\displaystyle B^2−5B\displaystyle {}- 5B+(−12)E\displaystyle {}+ (-12)E=O\displaystyle {}= O⟹B2‾\displaystyle {}\Longrightarrow \quad \underline{\rule[-0.0833em]{0em}{0.9474em}B^2}−5B‾\displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}{}- 5B}−12E=O‾\displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}{}- 12E = O}

定数項の符号に注意します。Δ\Delta が負のときは +ΔE+\Delta E がそのまま −12E-12E になります。Δ\Delta の符号を変えて +12E+12E と書いてしまう誤りが多いところです。

自己採点:
記録を読み込み中…

問2 ★

次の 2 つの行列 (1) AA=(2341){}= \begin{pmatrix} 2 & 3 \\ 4 & 1 \end{pmatrix}、(2) AA=(−2153){}= \begin{pmatrix} -2 & 1 \\ 5 & 3 \end{pmatrix} について、ケーリー・ハミルトンの定理を使って A2A^2 をそれぞれ求めなさい。

つまずいたときは:
答えを見る
答え

(1) (1691213)\begin{pmatrix} 16 & 9 \\ 12 & 13 \end{pmatrix} (2) (91514)\begin{pmatrix} 9 & 1 \\ 5 & 14 \end{pmatrix}

解説

(1) aa+d{}+ d=3{}= 3、Δ\Delta=2×1{}= 2 \times 1−3×4{}- 3 \times 4=−10{}= -10 なので A2A^2=3A{}= 3A+10E{}+ 10E です。

A2\displaystyle A^2=3(2341)\displaystyle {}= 3\begin{pmatrix} 2 & 3 \\ 4 & 1 \end{pmatrix}+10E\displaystyle {}+ 10E=(6+109123+10)\displaystyle {}= \begin{pmatrix} 6 + 10 & 9 \\ 12 & 3 + 10 \end{pmatrix}=(1691213)‾\displaystyle {}= \underline{\begin{pmatrix} 16 & 9 \\ 12 & 13 \end{pmatrix}}

(2) aa+d{}+ d=−2{}= -2+3{}+ 3=1{}= 1、Δ\Delta=(−2)×3{}= (-2) \times 3−1×5{}- 1 \times 5=−11{}= -11 なので A2A^2=A{}= A+11E{}+ 11E です。

A2\displaystyle A^2=(−2153)\displaystyle {}= \begin{pmatrix} -2 & 1 \\ 5 & 3 \end{pmatrix}+11E\displaystyle {}+ 11E=(91514)‾\displaystyle {}= \underline{\begin{pmatrix} 9 & 1 \\ 5 & 14 \end{pmatrix}}

どちらも行列どうしのかけ算をしていません。A2A^2 ひとつなら直接かけたほうが速いくらいですが、AA の実数倍と EE の実数倍を足すだけで済む、という感覚をここでつかんでおきます。

自己採点:
記録を読み込み中…

問3 ★

AA=(1−234){}= \begin{pmatrix} 1 & -2 \\ 3 & 4 \end{pmatrix} について、次数下げを使って A3A^3 を求めなさい。

つまずいたときは:
答えを見る
答え

A3A^3=15A{}= 15A−50E{}- 50E=(−35−304510){}= \begin{pmatrix} -35 & -30 \\ 45 & 10 \end{pmatrix}

解説

aa+d{}+ d=1{}= 1+4{}+ 4=5{}= 5、Δ\Delta=1×4{}= 1 \times 4−(−2)×3{}- (-2) \times 3=4{}= 4+6{}+ 6=10{}= 10 なので

A2\displaystyle A^2=5A\displaystyle {}= 5A−10E\displaystyle {}- 10E

両辺に右から AA をかけ、出てきた A2A^2 をもう一度この式で置きかえます。

A3\displaystyle A^3=5A2\displaystyle {}= 5A^2−10A\displaystyle {}- 10A=5(5A−10E)\displaystyle {}= 5(5A - 10E)−10A\displaystyle {}- 10A=25A\displaystyle {}= 25A−50E\displaystyle {}- 50E−10A\displaystyle {}- 10A=15A\displaystyle {}= 15A−50E\displaystyle {}- 50E

成分を入れます。

A3\displaystyle A^3=15(1−234)\displaystyle {}= 15\begin{pmatrix} 1 & -2 \\ 3 & 4 \end{pmatrix}−50E\displaystyle {}- 50E=(15−50−304560−50)\displaystyle {}= \begin{pmatrix} 15 - 50 & -30 \\ 45 & 60 - 50 \end{pmatrix}=(−35−304510)‾\displaystyle {}= \underline{\begin{pmatrix} -35 & -30 \\ 45 & 10 \end{pmatrix}}

A2A^2 を成分で求めてから AA をかけても同じ答えになります(A2A^2=(−5−101510){}= \begin{pmatrix} -5 & -10 \\ 15 & 10 \end{pmatrix})。確かめておくと安心です。

自己採点:
記録を読み込み中…

問4 ★

A−1A^{-1}=1Δ{(a+d)E−A}{}= \dfrac{1}{\Delta}\{(a+d)E - A\} を使って、次の 2 つの行列 (1) AA=(4312){}= \begin{pmatrix} 4 & 3 \\ 1 & 2 \end{pmatrix}、(2) AA=(3512){}= \begin{pmatrix} 3 & 5 \\ 1 & 2 \end{pmatrix} の逆行列をそれぞれ求めなさい。

つまずいたときは:
答えを見る
答え

(1) 15(2−3−14)\dfrac{1}{5}\begin{pmatrix} 2 & -3 \\ -1 & 4 \end{pmatrix} (2) (2−5−13)\begin{pmatrix} 2 & -5 \\ -1 & 3 \end{pmatrix}

解説

(1) aa+d{}+ d=6{}= 6、Δ\Delta=4×2{}= 4 \times 2−3×1{}- 3 \times 1=5{}= 5≠0{}\neq 0 なので正則です。

A−1\displaystyle A^{-1}=15(6E−A)\displaystyle {}= \frac{1}{5}(6E - A)=15{(6006)−(4312)}\displaystyle {}= \frac{1}{5}\left\{\begin{pmatrix} 6 & 0 \\ 0 & 6 \end{pmatrix} - \begin{pmatrix} 4 & 3 \\ 1 & 2 \end{pmatrix}\right\}=15(2−3−14)‾\displaystyle {}= \underline{\frac{1}{5}\begin{pmatrix} 2 & -3 \\ -1 & 4 \end{pmatrix}}

(2) aa+d{}+ d=5{}= 5、Δ\Delta=3×2{}= 3 \times 2−5×1{}- 5 \times 1=1{}= 1≠0{}\neq 0 なので

A−1\displaystyle A^{-1}=5E\displaystyle {}= 5E−A\displaystyle {}- A=(5−3−5−15−2)\displaystyle {}= \begin{pmatrix} 5 - 3 & -5 \\ -1 & 5 - 2 \end{pmatrix}=(2−5−13)‾\displaystyle {}= \underline{\begin{pmatrix} 2 & -5 \\ -1 & 3 \end{pmatrix}}

第 3 章の公式 A−1A^{-1}=1Δ(d−b−ca){}= \dfrac{1}{\Delta}\begin{pmatrix} d & -b \\ -c & a \end{pmatrix} で求めても、まったく同じ行列が出ます。(a+d)E(a+d)E−A{}- A を計算すると、aa と dd が入れかわり、bb と cc の符号が変わる。あの手順の中身は引き算だったわけです。

自己採点:
記録を読み込み中…

問5 ★

AA=(2412){}= \begin{pmatrix} 2 & 4 \\ 1 & 2 \end{pmatrix} について、A2A^2 と A3A^3 をそれぞれ AA の実数倍の形で表しなさい。

つまずいたときは:
答えを見る
答え

A2A^2=4A{}= 4A,A3A^3=16A{}= 16A

解説

aa+d{}+ d=2{}= 2+2{}+ 2=4{}= 4、Δ\Delta=2×2{}= 2 \times 2−4×1{}- 4 \times 1=0{}= 0 です。Δ\Delta=0{}= 0 なので定理の式の定数項が消えて

A2\displaystyle A^2−4A\displaystyle {}- 4A+0⋅E\displaystyle {}+ 0 \cdot E=O\displaystyle {}= O⟹A2=4A‾\displaystyle {}\Longrightarrow \quad \underline{A^2 = 4A}

これを使えば A3A^3 もすぐです。

A3\displaystyle A^3=A2A\displaystyle {}= A^2 A=(4A)A\displaystyle {}= (4A)A=4A2\displaystyle {}= 4A^2=4×4A\displaystyle {}= 4 \times 4A=16A‾\displaystyle {}= \underline{16A}

成分で確かめると A2A^2=(81648){}= \begin{pmatrix} 8 & 16 \\ 4 & 8 \end{pmatrix}=4A{}= 4A となっています。

Δ\Delta=0{}= 0 の行列はこのように EE の項が現れず、累乗がすべて AA の実数倍で済みます。一般には AnA^n=(a+d)n−1A{}= (a+d)^{n-1}A となり、これは j15 で示します。

自己採点:
記録を読み込み中…

問6 ★

AA=(1−110){}= \begin{pmatrix} 1 & -1 \\ 1 & 0 \end{pmatrix} について、次の問いに答えなさい。

(1) A2A^2 を AA と EE で表しなさい。

(2) A3A^3 を求めなさい。

(3) A2026A^{2026} を求めなさい。

つまずいたときは:
答えを見る
答え

(1) A2A^2=A{}= A−E{}- E (2) A3A^3=−E{}= -E (3) A2026A^{2026}=(−11−10){}= \begin{pmatrix} -1 & 1 \\ -1 & 0 \end{pmatrix}

解説

(1) aa+d{}+ d=1{}= 1+0{}+ 0=1{}= 1、Δ\Delta=1×0{}= 1 \times 0−(−1)×1{}- (-1) \times 1=1{}= 1 なので

A2\displaystyle A^2−A\displaystyle {}- A+E\displaystyle {}+ E=O\displaystyle {}= O⟹A2=A‾\displaystyle {}\Longrightarrow \quad \underline{\rule[-0.0833em]{0em}{0.9474em}A^2 = A}−E‾\displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}{}- E}

(2) 両辺に AA をかけ、A2A^2 をもう一度下ろします。

A3\displaystyle A^3=A2\displaystyle {}= A^2−A\displaystyle {}- A=(A−E)\displaystyle {}= (A - E)−A\displaystyle {}- A=−E‾\displaystyle {}= \underline{-E}

(3) A3A^3=−E{}= -E なので A6A^6=(A3)2{}= (A^3)^2=(−E)2{}= (-E)^2=E{}= E です。6 乗すると単位行列に戻る、つまり周期は 6 です。

20262026=6×337{}= 6 \times 337+4{}+ 4 なので

A2026\displaystyle A^{2026}=(A6)337A4\displaystyle {}= (A^6)^{337}A^4=E⋅A4\displaystyle {}= E \cdot A^4=A3A\displaystyle {}= A^3 A=−A\displaystyle {}= -A=(−11−10)‾\displaystyle {}= \underline{\begin{pmatrix} -1 & 1 \\ -1 & 0 \end{pmatrix}}

EE に戻るところを見つけたら、指数を周期で割った余りだけを見ればよくなります。第 2 章 j13 で A4A^4=E{}= E の行列を扱ったのと同じ考え方です。AnA^n の求め方のうち、いちばん手早く済むのがこの形です。

自己採点:
記録を読み込み中…

問7 ★

A2A^2−5A{}- 5A+6E{}+ 6E=O{}= O を満たす 2 次の正方行列 AA のうち、対角和 aa+d{}+ d が 55 に等しくないものを 1 つ求めなさい。

つまずいたときは:
答えを見る
答え

AA=2E{}= 2E=(2002){}= \begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix}(対角和は 44)。AA=3E{}= 3E でもよい

解説

AA=kE{}= kE の形をためします。A2A^2=k2E{}= k^2E なので

A2\displaystyle A^2−5A\displaystyle {}- 5A+6E\displaystyle {}+ 6E=(k2−5k+6)E\displaystyle {}= (k^2 - 5k + 6)E

これが OO になるのは k2k^2−5k{}- 5k+6{}+ 6=0{}= 0、つまり (k−2)(k−3)(k-2)(k-3)=0{}= 0 より kk=2,{}= 2, 3\ 3 のときです。

kk=2{}= 2 とすると

A\displaystyle A=(2002)‾\displaystyle {}= \underline{\begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix}}

このとき対角和は 22+2{}+ 2=4{}= 4 で、55 ではありません。行列式も 44 で、66 ではありません。kk=3{}= 3 とした AA=3E{}= 3E(対角和 66、行列式 99)も答えになります。

ケーリー・ハミルトンの定理は「AA を決めれば係数が決まる」と言っているだけで、「係数が分かれば AA の対角和が分かる」とは言っていません。AA=kE{}= kE という形は、x2x^2−5x{}- 5x+6{}+ 6=0{}= 0 の解を 1 つでも含んでいれば式が成り立ってしまうので、係数と対角和のつながりが切れます。

自己採点:
記録を読み込み中…

問8 ★

AA=(2134){}= \begin{pmatrix} 2 & 1 \\ 3 & 4 \end{pmatrix} について、次の 2 つの行列 (1) A2A^2−6A{}- 6A+5E{}+ 5E、(2) A2A^2−4A{}- 4A+3E{}+ 3E をそれぞれ求めなさい。

つまずいたときは:
答えを見る
答え

(1) OO (2) (2266)\begin{pmatrix} 2 & 2 \\ 6 & 6 \end{pmatrix}

解説

aa+d{}+ d=2{}= 2+4{}+ 4=6{}= 6、Δ\Delta=2×4{}= 2 \times 4−1×3{}- 1 \times 3=5{}= 5 なので、ケーリー・ハミルトンの定理は

A2\displaystyle A^2−6A\displaystyle {}- 6A+5E\displaystyle {}+ 5E=O\displaystyle {}= O

(1) は定理の式そのものなので O‾\underline{O} です。計算するまでもありません。

(2) 求めたい式を、(1) の式との差でとらえます。

A2\displaystyle A^2−4A\displaystyle {}- 4A+3E\displaystyle {}+ 3E=(A2−6A+5E)\displaystyle {}= (A^2 - 6A + 5E)+2A\displaystyle {}+ 2A−2E\displaystyle {}- 2E=O\displaystyle {}= O+2A\displaystyle {}+ 2A−2E\displaystyle {}- 2E=2A\displaystyle {}= 2A−2E\displaystyle {}- 2E 2A\displaystyle 2A−2E\displaystyle {}- 2E=(4268)\displaystyle {}= \begin{pmatrix} 4 & 2 \\ 6 & 8 \end{pmatrix}−(2002)\displaystyle {}- \begin{pmatrix} 2 & 0 \\ 0 & 2 \end{pmatrix}=(2266)‾\displaystyle {}= \underline{\begin{pmatrix} 2 & 2 \\ 6 & 6 \end{pmatrix}}

A2A^2 を成分で求めると (761819)\begin{pmatrix} 7 & 6 \\ 18 & 19 \end{pmatrix} で、そこから計算しても同じ答えになります。定理の式を「A2A^2 を消すための道具」と見ると、こういう式の値はすべて AA の 1 次式に落ちます。

自己採点:
記録を読み込み中…

問9 ★★

AA=(3122){}= \begin{pmatrix} 3 & 1 \\ 2 & 2 \end{pmatrix} について、正の整数 nn に対する AnA^n を求めなさい。

つまずいたときは:
答えを見る
答え

AnA^n=13(2⋅4n+14n−12⋅4n−24n+2){}= \dfrac{1}{3}\begin{pmatrix} 2 \cdot 4^n + 1 & 4^n - 1 \\ 2 \cdot 4^n - 2 & 4^n + 2 \end{pmatrix}

解説

aa+d{}+ d=3{}= 3+2{}+ 2=5{}= 5、Δ\Delta=3×2{}= 3 \times 2−1×2{}- 1 \times 2=4{}= 4 なので、AA が背負う 2 次方程式は

x2\displaystyle x^2−5x\displaystyle {}- 5x+4\displaystyle {}+ 4=0\displaystyle {}= 0⟹(x−1)(x−4)\displaystyle {}\Longrightarrow \quad (x - 1)(x - 4)=0\displaystyle {}= 0⟹x\displaystyle {}\Longrightarrow \quad x=1,\displaystyle {}= 1, 4\displaystyle \ 4

xnx^n を x2x^2−5x{}- 5x+4{}+ 4 で割った商を Q(x)Q(x)、余りを pxpx+q{}+ q とします。

xn\displaystyle x^n=(x2−5x+4)Q(x)\displaystyle {}= (x^2 - 5x + 4)Q(x)+px\displaystyle {}+ px+q\displaystyle {}+ q

xx=1{}= 1 と xx=4{}= 4 を代入します。

{1=p+q4n=4p+q\begin{cases} 1 = p + q \\ 4^n = 4p + q \end{cases}

辺々引いて 4n4^n−1{}- 1=3p{}= 3p より

p\displaystyle p=4n−13,\displaystyle {}= \frac{4^n - 1}{3},q\displaystyle q=1\displaystyle {}= 1−p\displaystyle {}- p=4−4n3\displaystyle {}= \frac{4 - 4^n}{3}

ケーリー・ハミルトンの定理より A2A^2−5A{}- 5A+4E{}+ 4E=O{}= O なので、上の恒等式に AA を代入すると第 1 項が消えて AnA^n=pA{}= pA+qE{}+ qE です。

An\displaystyle A^n=p(3122)\displaystyle {}= p\begin{pmatrix} 3 & 1 \\ 2 & 2 \end{pmatrix}+qE\displaystyle {}+ qE=(3p+qp2p2p+q)\displaystyle {}= \begin{pmatrix} 3p + q & p \\ 2p & 2p + q \end{pmatrix}

3p3p+q{}+ q=3⋅4n−3+4−4n3{}= \dfrac{3 \cdot 4^n - 3 + 4 - 4^n}{3}=2⋅4n+13{}= \dfrac{2 \cdot 4^n + 1}{3}、2p2p+q{}+ q=2⋅4n−2+4−4n3{}= \dfrac{2 \cdot 4^n - 2 + 4 - 4^n}{3}=4n+23{}= \dfrac{4^n + 2}{3} なので

An\displaystyle A^n=13(2⋅4n+14n−12⋅4n−24n+2)‾\displaystyle {}= \underline{\frac{1}{3}\begin{pmatrix} 2 \cdot 4^n + 1 & 4^n - 1 \\ 2 \cdot 4^n - 2 & 4^n + 2 \end{pmatrix}}

nn=1{}= 1 で 13(9366)\dfrac{1}{3}\begin{pmatrix} 9 & 3 \\ 6 & 6 \end{pmatrix}=A{}= A、nn=2{}= 2 で 13(33153018)\dfrac{1}{3}\begin{pmatrix} 33 & 15 \\ 30 & 18 \end{pmatrix}=(115106){}= \begin{pmatrix} 11 & 5 \\ 10 & 6 \end{pmatrix}。直接計算した A2A^2 と一致します。

自己採点:
記録を読み込み中…

問10 ★★

AA=(5−113){}= \begin{pmatrix} 5 & -1 \\ 1 & 3 \end{pmatrix} について、正の整数 nn に対する AnA^n を求めなさい。

つまずいたときは:
答えを見る
答え

AnA^n=4n−1(n+4−nn4−n){}= 4^{n-1}\begin{pmatrix} n + 4 & -n \\ n & 4 - n \end{pmatrix}

解説

aa+d{}+ d=5{}= 5+3{}+ 3=8{}= 8、Δ\Delta=5×3{}= 5 \times 3−(−1)×1{}- (-1) \times 1=16{}= 16 なので

x2\displaystyle x^2−8x\displaystyle {}- 8x+16\displaystyle {}+ 16=0\displaystyle {}= 0⟹(x−4)2\displaystyle {}\Longrightarrow \quad (x - 4)^2=0\displaystyle {}= 0

重解 xx=4{}= 4 です。解が 1 つしかないので、連立方程式を作る方法は使えません。NN=A{}= A−4E{}- 4E とおきます。

N\displaystyle N=(1−11−1),\displaystyle {}= \begin{pmatrix} 1 & -1 \\ 1 & -1 \end{pmatrix},N2\displaystyle N^2=(1−1−1+11−1−1+1)\displaystyle {}= \begin{pmatrix} 1 - 1 & -1 + 1 \\ 1 - 1 & -1 + 1 \end{pmatrix}=O\displaystyle {}= O

ケーリー・ハミルトンの定理 (A−4E)2(A - 4E)^2=O{}= O が、そのまま N2N^2=O{}= O を保証しています。

AA=4E{}= 4E+N{}+ N で、4E4E はどんな行列とも交換可能なので二項定理が使えます。N2N^2 以降はすべて OO なので

An\displaystyle A^n=(4E+N)n\displaystyle {}= (4E + N)^n=4nE\displaystyle {}= 4^nE+n⋅4n−1N\displaystyle {}+ n \cdot 4^{n-1}N

成分に直します。

An\displaystyle A^n=4n(1001)\displaystyle {}= 4^n\begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}+n⋅4n−1(1−11−1)\displaystyle {}+ n \cdot 4^{n-1}\begin{pmatrix} 1 & -1 \\ 1 & -1 \end{pmatrix}=4n−1(n+4−nn4−n)‾\displaystyle {}= \underline{4^{n-1}\begin{pmatrix} n + 4 & -n \\ n & 4 - n \end{pmatrix}}

nn=1{}= 1 で (5−113)\begin{pmatrix} 5 & -1 \\ 1 & 3 \end{pmatrix}=A{}= A、nn=2{}= 2 で 4(6−222)4\begin{pmatrix} 6 & -2 \\ 2 & 2 \end{pmatrix}=(24−888){}= \begin{pmatrix} 24 & -8 \\ 8 & 8 \end{pmatrix}。A2A^2=8A{}= 8A−16E{}- 16E から求めた値と一致します。

重解かどうかで方法が変わるので、x2x^2−(a+d)x{}- (a+d)x+Δ{}+ \Delta=0{}= 0 を解いた時点で判別式を確かめる習慣をつけておきます。

自己採点:
記録を読み込み中…

問11 ★★

AA=(3011){}= \begin{pmatrix} 3 & 0 \\ 1 & 1 \end{pmatrix} について、A2A^2 と A3A^3 を計算して AnA^n を推定し、それが正しいことを数学的帰納法で示しなさい。

つまずいたときは:
答えを見る
答え

AnA^n=(3n03n−121){}= \begin{pmatrix} 3^n & 0 \\ \dfrac{3^n - 1}{2} & 1 \end{pmatrix}

解説

まず具体的に計算します。

A2\displaystyle A^2=(3011)(3011)\displaystyle {}= \begin{pmatrix} 3 & 0 \\ 1 & 1 \end{pmatrix}\begin{pmatrix} 3 & 0 \\ 1 & 1 \end{pmatrix}=(9041),\displaystyle {}= \begin{pmatrix} 9 & 0 \\ 4 & 1 \end{pmatrix},A3\displaystyle A^3=A2A\displaystyle {}= A^2A=(9041)(3011)\displaystyle {}= \begin{pmatrix} 9 & 0 \\ 4 & 1 \end{pmatrix}\begin{pmatrix} 3 & 0 \\ 1 & 1 \end{pmatrix}=(270131)\displaystyle {}= \begin{pmatrix} 27 & 0 \\ 13 & 1 \end{pmatrix}

左上は 3,3, 9,\ 9, 27\ 27 で 3n3^n、右上は 00、右下は 11 のままです。左下は 1,1, 4,\ 4, 13\ 13。これは 3n−12\dfrac{3^n - 1}{2}(22,\dfrac{2}{2}, 82,\ \dfrac{8}{2}, 262\ \dfrac{26}{2})と読めます。そこで

An\displaystyle A^n=(3n03n−121)\displaystyle {}= \begin{pmatrix} 3^n & 0 \\ \dfrac{3^n - 1}{2} & 1 \end{pmatrix}

と推定します。

数学的帰納法で示します。

nn=1{}= 1 のとき、右辺は (3011)\begin{pmatrix} 3 & 0 \\ 1 & 1 \end{pmatrix}=A{}= A で成り立ちます。

nn=k{}= k で成り立つとします。両辺に右から AA をかけると

Ak+1\displaystyle A^{k+1}=(3k03k−121)(3011)\displaystyle {}= \begin{pmatrix} 3^k & 0 \\ \dfrac{3^k - 1}{2} & 1 \end{pmatrix}\begin{pmatrix} 3 & 0 \\ 1 & 1 \end{pmatrix}=(3k+103(3k−1)2+11)\displaystyle {}= \begin{pmatrix} 3^{k+1} & 0 \\ \dfrac{3(3^k - 1)}{2} + 1 & 1 \end{pmatrix}

左下を整理します。

3k+1−32\displaystyle \frac{3^{k+1} - 3}{2}+1\displaystyle {}+ 1=3k+1−3+22\displaystyle {}= \frac{3^{k+1} - 3 + 2}{2}=3k+1−12\displaystyle {}= \frac{3^{k+1} - 1}{2}

よって nn=k{}= k+1{}+ 1 でも成り立ちます。ゆえにすべての正の整数 nn で正しいといえます。

この AA は x2x^2−4x{}- 4x+3{}+ 3=0{}= 0、つまり解 1,1, 3\ 3 を背負っているので、公式5 の割り算の余りでも求まります(pp=3n−12{}= \dfrac{3^n - 1}{2},qq=3−3n2{}= \dfrac{3 - 3^n}{2})。成分に 00 が多い行列では、規則を目で見つけて帰納法で固めるほうが速いことがあります。

自己採点:
記録を読み込み中…

問12 ★★

A2A^2−4A{}- 4A+3E{}+ 3E=O{}= O を満たす 2 次の正方行列 AA について、対角和 aa+d{}+ d と行列式 Δ\Delta の組をすべて求めなさい。

つまずいたときは:
答えを見る
答え

(a+d, Δ)(a+d,\ \Delta)=(4, 3),{}= (4,\ 3), (2, 1),\ (2,\ 1), (6, 9)\ (6,\ 9)

解説

ケーリー・ハミルトンの定理より A2A^2−(a+d)A{}- (a+d)A+ΔE{}+ \Delta E=O{}= O です。与えられた式との差をとると

{4−(a+d)}A\displaystyle \{4 - (a + d)\}A=(3−Δ)E\displaystyle {}= (3 - \Delta)E

(i) 44−(a+d){}- (a+d)≠0{}\neq 0 のとき。 両辺を 44−(a+d){}- (a+d) で割ると AA=kE{}= kE の形になります。この場合は (ii) で扱います。

(ii) AA が kEkE の形でないとき。 上の式で AA の係数が 00 でなければ AA=kE{}= kE になってしまうので、44−(a+d){}- (a+d)=0{}= 0 でなければなりません。すると右辺も OO となり 33−Δ{}- \Delta=0{}= 0 です。

a‾\displaystyle \underline{\rule[-0.1944em]{0em}{0.8888em}a}+d=4,‾\displaystyle \underline{\rule[-0.1944em]{0em}{0.8888em}{}+ d = 4,}Δ=3‾\displaystyle \underline{\rule[-0.1944em]{0em}{0.8888em}\quad \Delta = 3}

(iii) AA=kE{}= kE のとき。 A2A^2=k2E{}= k^2E なので

(k2−4k+3)E\displaystyle (k^2 - 4k + 3)E=O\displaystyle {}= O⟹(k−1)(k−3)\displaystyle {}\Longrightarrow \quad (k - 1)(k - 3)=0\displaystyle {}= 0⟹k\displaystyle {}\Longrightarrow \quad k=1,\displaystyle {}= 1, 3\displaystyle \ 3

kk=1{}= 1 のとき AA=E{}= E で aa+d{}+ d=2{}= 2、Δ\Delta=1{}= 1。

kk=3{}= 3 のとき AA=3E{}= 3E で aa+d{}+ d=6{}= 6、Δ\Delta=9{}= 9。

以上より、組は (4, 3),‾\underline{\rule[-0.25em]{0em}{1.0000em}(4,\ 3),} (2, 1),‾\underline{\rule[-0.25em]{0em}{1.0000em}\ (2,\ 1),} (6, 9)‾\underline{\rule[-0.25em]{0em}{1.0000em}\ (6,\ 9)} の 3 通りです。

AA=kE{}= kE の場合を落とすと (4, 3) だけを答えてしまいます。A2A^2−pA{}- pA+qE{}+ qE=O{}= O という形の条件を見たら、まず「AA=kE{}= kE かどうか」で分ける、と手が動くようにしておきます。

自己採点:
記録を読み込み中…

問13 ★★

PP=(111−2){}= \begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix} とし、2 次の正方行列 AA が P−1APP^{-1}AP=(200−1){}= \begin{pmatrix} 2 & 0 \\ 0 & -1 \end{pmatrix} を満たすとする。AA と AnA^n を求めなさい。

つまずいたときは:
答えを見る
答え

AA=(1120){}= \begin{pmatrix} 1 & 1 \\ 2 & 0 \end{pmatrix},AnA^n=13(2u+vu−v2u−2vu+2v){}= \dfrac{1}{3}\begin{pmatrix} 2u + v & u - v \\ 2u - 2v & u + 2v \end{pmatrix}(uu=2n{}= 2^n,vv=(−1)n{}= (-1)^n)

解説

PP の行列式は 1×(−2)1 \times (-2)−1×1{}- 1 \times 1=−3{}= -3≠0{}\neq 0 なので

P−1\displaystyle P^{-1}=1−3(−2−1−11)\displaystyle {}= \frac{1}{-3}\begin{pmatrix} -2 & -1 \\ -1 & 1 \end{pmatrix}=13(211−1)\displaystyle {}= \frac{1}{3}\begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix}

DD=(200−1){}= \begin{pmatrix} 2 & 0 \\ 0 & -1 \end{pmatrix} とおくと P−1APP^{-1}AP=D{}= D の両辺を PP と P−1P^{-1} で挟んで AA=PDP−1{}= PDP^{-1} です。

PD\displaystyle PD=(111−2)(200−1)\displaystyle {}= \begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix}\begin{pmatrix} 2 & 0 \\ 0 & -1 \end{pmatrix}=(2−122)\displaystyle {}= \begin{pmatrix} 2 & -1 \\ 2 & 2 \end{pmatrix} A\displaystyle A=(2−122)⋅13(211−1)\displaystyle {}= \begin{pmatrix} 2 & -1 \\ 2 & 2 \end{pmatrix} \cdot \frac{1}{3}\begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix}=13(4−12+14+22−2)\displaystyle {}= \frac{1}{3}\begin{pmatrix} 4 - 1 & 2 + 1 \\ 4 + 2 & 2 - 2 \end{pmatrix}=(1120)‾\displaystyle {}= \underline{\begin{pmatrix} 1 & 1 \\ 2 & 0 \end{pmatrix}}

次に AnA^n です。DnD^n=(2n00(−1)n){}= \begin{pmatrix} 2^n & 0 \\ 0 & (-1)^n \end{pmatrix} で、AnA^n=PDnP−1{}= PD^nP^{-1} が成り立ちます。

PDn\displaystyle PD^n=(111−2)(2n00(−1)n)\displaystyle {}= \begin{pmatrix} 1 & 1 \\ 1 & -2 \end{pmatrix}\begin{pmatrix} 2^n & 0 \\ 0 & (-1)^n \end{pmatrix}=(2n(−1)n2n−2(−1)n)\displaystyle {}= \begin{pmatrix} 2^n & (-1)^n \\ 2^n & -2(-1)^n \end{pmatrix}

成分が長くなるので、2n2^n=u{}= u,(−1)n(-1)^n=v{}= v と書きます。PDnPD^n=(uvu−2v){}= \begin{pmatrix} u & v \\ u & -2v \end{pmatrix} なので

An\displaystyle A^n=(uvu−2v)⋅13(211−1)\displaystyle {}= \begin{pmatrix} u & v \\ u & -2v \end{pmatrix} \cdot \frac{1}{3}\begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix}=13(2u+vu−v2u−2vu+2v)‾\displaystyle {}= \underline{\frac{1}{3}\begin{pmatrix} 2u + v & u - v \\ 2u - 2v & u + 2v \end{pmatrix}}

ただし uu=2n{}= 2^n,vv=(−1)n{}= (-1)^n です。2u2u=2n+1{}= 2^{n+1} なので、たとえば (1, 1)(1,\ 1) 成分は 2n+1+(−1)n3\dfrac{2^{n+1} + (-1)^n}{3} です。

nn=1{}= 1(uu=2{}= 2,vv=−1{}= -1)で 13(3360)\dfrac{1}{3}\begin{pmatrix} 3 & 3 \\ 6 & 0 \end{pmatrix}=A{}= A、nn=2{}= 2 で 13(9366)\dfrac{1}{3}\begin{pmatrix} 9 & 3 \\ 6 & 6 \end{pmatrix}=(3122){}= \begin{pmatrix} 3 & 1 \\ 2 & 2 \end{pmatrix}。直接計算した A2A^2 と一致します。

DD の対角成分 22,−1-1 は、AA の背負う 2 次方程式 x2x^2−x{}- x−2{}- 2=0{}= 0 の解です(AA の対角和は 11、行列式は −2-2)。PP で挟んでも対角和と行列式は変わらないので、こうなるのは当然のことです。

自己採点:
記録を読み込み中…

問14 ★★

2 次の正方行列 AA=(abcd){}= \begin{pmatrix} a & b \\ c & d \end{pmatrix} について、A2A^2=O{}= O となるための必要十分条件が「aa+d{}+ d=0{}= 0 かつ adad−bc{}- bc=0{}= 0」であることを示しなさい。

つまずいたときは:
答えを見る
答え

十分性はケーリー・ハミルトンの定理から直ちに従い、必要性は (a+d)A(a+d)A=ΔE{}= \Delta E を AA=kE{}= kE の場合と分けて調べればよい

解説

ケーリー・ハミルトンの定理より、つねに

A2\displaystyle A^2−(a+d)A\displaystyle {}- (a + d)A+ΔE\displaystyle {}+ \Delta E=O\displaystyle {}= O(Δ=ad−bc)\displaystyle (\Delta = ad - bc)⋯(∗)\displaystyle \cdots (\ast)

が成り立ちます。

【十分性】 aa+d{}+ d=0{}= 0 かつ Δ\Delta=0{}= 0 とすると、(∗)(\ast) は A2A^2=O{}= O となります。

【必要性】 A2A^2=O{}= O とします。(∗)(\ast) に代入すると

−(a+d)A\displaystyle -(a + d)A+ΔE\displaystyle {}+ \Delta E=O\displaystyle {}= O⟹(a+d)A\displaystyle {}\Longrightarrow \quad (a + d)A=ΔE\displaystyle {}= \Delta E

ここで aa+d{}+ d≠0{}\neq 0 と仮定します。両辺を aa+d{}+ d で割ると

A\displaystyle A=Δa+dE\displaystyle {}= \frac{\Delta}{a + d}E

となり、kk=Δa+d{}= \dfrac{\Delta}{a+d} とおけば AA=kE{}= kE です。このとき A2A^2=k2E{}= k^2E=O{}= O より kk=0{}= 0、つまり AA=O{}= O となります。しかし AA=O{}= O なら aa+d{}+ d=0{}= 0 であり、aa+d{}+ d≠0{}\neq 0 という仮定に反します。

よって aa+d{}+ d=0{}= 0 です。これを (a+d)A(a+d)A=ΔE{}= \Delta E に入れると ΔE\Delta E=O{}= O、すなわち Δ\Delta=0{}= 0 を得ます。

以上より、A2A^2=O{}= O と「aa+d{}+ d=0{}= 0 かつ adad−bc{}- bc=0{}= 0」は同値です。

この条件を満たす行列は、たとえば (1−11−1)\begin{pmatrix} 1 & -1 \\ 1 & -1 \end{pmatrix}(j10 で使った NN)や (0100)\begin{pmatrix} 0 & 1 \\ 0 & 0 \end{pmatrix} です。AA≠O{}\neq O なのに A2A^2=O{}= O になるのは、数の世界では起こらないこと。第 2 章で見た零因子の、いちばん極端な形にあたります。

自己採点:
記録を読み込み中…

問15 ★★

2 次の正方行列 AA=(abcd){}= \begin{pmatrix} a & b \\ c & d \end{pmatrix} が adad−bc{}- bc=0{}= 0 を満たすとき、すべての正の整数 nn に対して AnA^n=(a+d)n−1A{}= (a + d)^{n-1}A が成り立つことを示しなさい。

つまずいたときは:
答えを見る
答え

Δ\Delta=0{}= 0 よりケーリー・ハミルトンの定理は A2A^2=(a+d)A{}= (a+d)A となる。これを使って数学的帰納法で示す

解説

ss=a{}= a+d{}+ d とおきます。Δ\Delta=ad{}= ad−bc{}- bc=0{}= 0 なので、ケーリー・ハミルトンの定理は

A2\displaystyle A^2−sA\displaystyle {}- sA+0⋅E\displaystyle {}+ 0 \cdot E=O\displaystyle {}= O⟹A2\displaystyle {}\Longrightarrow \quad A^2=sA\displaystyle {}= sA⋯(∗)\displaystyle \cdots (\ast)

となります。示したい式は AnA^n=sn−1A{}= s^{n-1}A です。

nn=1{}= 1 のとき。 右辺は s0As^0A=A{}= A で、左辺と一致します。

nn=k{}= k で成り立つとする。 AkA^k=sk−1A{}= s^{k-1}A の両辺に右から AA をかけます。

Ak+1\displaystyle A^{k+1}=sk−1A2\displaystyle {}= s^{k-1}A^2

ここで (∗)(\ast) より A2A^2=sA{}= sA なので

Ak+1\displaystyle A^{k+1}=sk−1⋅sA\displaystyle {}= s^{k-1} \cdot sA=skA\displaystyle {}= s^{k}A=s(k+1)−1A\displaystyle {}= s^{(k+1)-1}A

となり、nn=k{}= k+1{}+ 1 でも成り立ちます。

よって数学的帰納法により、すべての正の整数 nn で AnA^n=(a+d)n−1A{}= (a+d)^{n-1}A が成り立ちます。

j05 の AA=(2412){}= \begin{pmatrix} 2 & 4 \\ 1 & 2 \end{pmatrix} は ss=4{}= 4、Δ\Delta=0{}= 0 だったので、AnA^n=4n−1A{}= 4^{n-1}A です。A2A^2=4A{}= 4A、A3A^3=16A{}= 16A と合っています。

なお ss=0{}= 0 の場合、この式は nn≧2{}\geqq 2 で AnA^n=O{}= O を意味します。j14 で調べた A2A^2=O{}= O の条件(aa+d{}+ d=0{}= 0 かつ Δ\Delta=0{}= 0)と、ちょうど重なります。

自己採点:
記録を読み込み中…

問16 ★★

数列 {an}\{a_n\},{bn}\{b_n\} が a1a_1=1{}= 1,b1b_1=0{}= 0 と

{an+1=4an+bnbn+1=2an+3bn\begin{cases} a_{n+1} = 4a_n + b_n \\ b_{n+1} = 2a_n + 3b_n \end{cases}

を満たすとする。行列の nn 乗を使って、ana_n と bnb_n を nn の式で求めなさい。

つまずいたときは:
答えを見る
答え

ana_n=2⋅5n−1+2n−13{}= \dfrac{2 \cdot 5^{n-1} + 2^{n-1}}{3},bnb_n=2(5n−1−2n−1)3{}= \dfrac{2(5^{n-1} - 2^{n-1})}{3}

解説

2 本の漸化式を行列で 1 本にまとめます(第 4 章 公式1 と同じ書き直しです)。

(an+1bn+1)\displaystyle \begin{pmatrix} a_{n+1} \\ b_{n+1} \end{pmatrix}=A(anbn),\displaystyle {}= A\begin{pmatrix} a_n \\ b_n \end{pmatrix},A\displaystyle A=(4123)\displaystyle {}= \begin{pmatrix} 4 & 1 \\ 2 & 3 \end{pmatrix}

これを繰り返し使うと

(anbn)\displaystyle \begin{pmatrix} a_n \\ b_n \end{pmatrix}=An−1(a1b1)\displaystyle {}= A^{n-1}\begin{pmatrix} a_1 \\ b_1 \end{pmatrix}=An−1(10)\displaystyle {}= A^{n-1}\begin{pmatrix} 1 \\ 0 \end{pmatrix}

An−1A^{n-1} を求めます。aa+d{}+ d=7{}= 7、Δ\Delta=4×3{}= 4 \times 3−1×2{}- 1 \times 2=10{}= 10 なので

x2\displaystyle x^2−7x\displaystyle {}- 7x+10\displaystyle {}+ 10=0\displaystyle {}= 0⟹(x−2)(x−5)\displaystyle {}\Longrightarrow \quad (x - 2)(x - 5)=0\displaystyle {}= 0⟹x\displaystyle {}\Longrightarrow \quad x=2,\displaystyle {}= 2, 5\displaystyle \ 5

mm=n{}= n−1{}- 1 とおき、xmx^m を x2x^2−7x{}- 7x+10{}+ 10 で割った余りを pxpx+q{}+ q とすると

{2m=2p+q5m=5p+q\displaystyle \begin{cases} 2^m = 2p + q \\ 5^m = 5p + q \end{cases}⟹p\displaystyle {}\Longrightarrow \quad p=5m−2m3,\displaystyle {}= \frac{5^m - 2^m}{3},q\displaystyle q=2m\displaystyle {}= 2^m−2p\displaystyle {}- 2p=5⋅2m−2⋅5m3\displaystyle {}= \frac{5 \cdot 2^m - 2 \cdot 5^m}{3}

AmA^m=pA{}= pA+qE{}+ qE を (10)\begin{pmatrix} 1 \\ 0 \end{pmatrix} にかけます。AA の第 1 列は (42)\begin{pmatrix} 4 \\ 2 \end{pmatrix}、EE の第 1 列は (10)\begin{pmatrix} 1 \\ 0 \end{pmatrix} なので

(anbn)\displaystyle \begin{pmatrix} a_n \\ b_n \end{pmatrix}=p(42)\displaystyle {}= p\begin{pmatrix} 4 \\ 2 \end{pmatrix}+q(10)\displaystyle {}+ q\begin{pmatrix} 1 \\ 0 \end{pmatrix}=(4p+q2p)\displaystyle {}= \begin{pmatrix} 4p + q \\ 2p \end{pmatrix}

4p4p と qq の分子(分母はどちらも 33)を足すと

4⋅5m\displaystyle 4 \cdot 5^m−4⋅2m\displaystyle {}- 4 \cdot 2^m+5⋅2m\displaystyle {}+ 5 \cdot 2^m−2⋅5m\displaystyle {}- 2 \cdot 5^m=2⋅5m\displaystyle {}= 2 \cdot 5^m+2m\displaystyle {}+ 2^m

なので

4p\displaystyle 4p+q\displaystyle {}+ q=2⋅5m+2m3\displaystyle {}= \frac{2 \cdot 5^m + 2^m}{3}

mm=n{}= n−1{}- 1 に戻して

an\displaystyle a_n=2⋅5n−1+2n−13‾,\displaystyle {}= \underline{\frac{2 \cdot 5^{n-1} + 2^{n-1}}{3}},bn\displaystyle b_n=2p\displaystyle {}= 2p=2(5n−1−2n−1)3‾\displaystyle {}= \underline{\frac{2(5^{n-1} - 2^{n-1})}{3}}

検算します。nn=1{}= 1 で a1a_1=2+13{}= \dfrac{2 + 1}{3}=1{}= 1,b1b_1=0{}= 0。nn=2{}= 2 で a2a_2=10+23{}= \dfrac{10 + 2}{3}=4{}= 4,b2b_2=2×33{}= \dfrac{2 \times 3}{3}=2{}= 2。漸化式から直接求めた a2a_2=4{}= 4,b2b_2=2{}= 2 と一致します。

数列の第 7 章では、2 本の漸化式を扱うのに「足す・引く」でうまい組み合わせを見つける必要がありました。行列で書けば、やることは AnA^n を求めるという 1 つの作業に固定されます。

自己採点:
記録を読み込み中…

問17 ★★★

A2A^2=A{}= A を満たす 2 次の正方行列 AA をすべて求めなさい。

つまずいたときは:
答えを見る
答え

AA=O{}= O、AA=E{}= E、および「aa+d{}+ d=1{}= 1 かつ adad−bc{}- bc=0{}= 0」を満たす行列のすべて

解説

AA=(abcd){}= \begin{pmatrix} a & b \\ c & d \end{pmatrix}、Δ\Delta=ad{}= ad−bc{}- bc とします。ケーリー・ハミルトンの定理より

A2\displaystyle A^2=(a+d)A\displaystyle {}= (a + d)A−ΔE\displaystyle {}- \Delta E

これが AA に等しいので

A\displaystyle A=(a+d)A\displaystyle {}= (a + d)A−ΔE\displaystyle {}- \Delta E⟹{(a+d)−1}A\displaystyle {}\Longrightarrow \quad \{(a + d) - 1\}A=ΔE\displaystyle {}= \Delta E⋯(∗)\displaystyle \cdots (\ast)

(i) aa+d{}+ d≠1{}\neq 1 のとき。 (∗)(\ast) の両辺を (a+d)(a+d)−1{}- 1 で割って AA=kE{}= kE の形になります。このとき A2A^2=k2E{}= k^2E なので、A2A^2=A{}= A より

k2E\displaystyle k^2E=kE\displaystyle {}= kE⟹k2\displaystyle {}\Longrightarrow \quad k^2=k\displaystyle {}= k⟹k\displaystyle {}\Longrightarrow \quad k=0,\displaystyle {}= 0, 1\displaystyle \ 1

kk=0{}= 0 なら AA=O{}= O(対角和 00)、kk=1{}= 1 なら AA=E{}= E(対角和 22)。どちらも aa+d{}+ d≠1{}\neq 1 を満たしており、条件に合います。

(ii) aa+d{}+ d=1{}= 1 のとき。 (∗)(\ast) の左辺が OO になるので ΔE\Delta E=O{}= O、すなわち Δ\Delta=0{}= 0 です。

逆に aa+d{}+ d=1{}= 1 かつ Δ\Delta=0{}= 0 なら、定理より A2A^2=1⋅A{}= 1 \cdot A−0⋅E{}- 0 \cdot E=A{}= A となって条件を満たします。

以上より、求める行列は

A=O,‾\displaystyle \underline{\rule[-0.1944em]{0em}{0.8889em}A = O,}A=E,‾\displaystyle \underline{\rule[-0.1944em]{0em}{0.8889em}\quad A = E,}および a‾\displaystyle \underline{\rule[-0.1944em]{0em}{0.8889em}\quad \text{および } a}+d=1 かつ ad‾\displaystyle \underline{\rule[-0.1944em]{0em}{0.8889em}{}+ d = 1 \text{ かつ } ad}−bc=0 を満たす A‾\displaystyle \underline{\rule[-0.1944em]{0em}{0.8889em}{}- bc = 0 \text{ を満たす } A}

(ii) の例を挙げると、(1000)\begin{pmatrix} 1 & 0 \\ 0 & 0 \end{pmatrix}、(2−21−1)\begin{pmatrix} 2 & -2 \\ 1 & -1 \end{pmatrix}、(1/21/21/21/2)\begin{pmatrix} 1/2 & 1/2 \\ 1/2 & 1/2 \end{pmatrix} などがあり、無数にあります。実際 (2−21−1)\begin{pmatrix} 2 & -2 \\ 1 & -1 \end{pmatrix} を 2 乗すると (4−2−4+22−1−2+1)\begin{pmatrix} 4 - 2 & -4 + 2 \\ 2 - 1 & -2 + 1 \end{pmatrix}=(2−21−1){}= \begin{pmatrix} 2 & -2 \\ 1 & -1 \end{pmatrix} でもとに戻ります。

数の世界で x2x^2=x{}= x を満たすのは xx=0,{}= 0, 1\ 1 の 2 つだけでした。行列では無数にある。OO と EE の「あいだ」に、2 乗しても変わらない行列がぎっしり詰まっているわけです(第 6 章以降の言葉でいえば、これは平面を 1 本の直線に押しつぶす操作にあたります。一度つぶしたものをもう一度つぶしても、何も変わりません)。

自己採点:
記録を読み込み中…

問18 ★★★

数列 {an}\{a_n\} が a1a_1=1{}= 1,a2a_2=4{}= 4、an+2a_{n+2}=5an+1{}= 5a_{n+1}−6an{}- 6a_n を満たすとする。(an+1an)\begin{pmatrix} a_{n+1} \\ a_n \end{pmatrix}=A(anan−1){}= A\begin{pmatrix} a_n \\ a_{n-1} \end{pmatrix} となる行列 AA を定め、AA の累乗を使って一般項 ana_n を求めなさい。

つまずいたときは:
答えを見る
答え

AA=(5−610){}= \begin{pmatrix} 5 & -6 \\ 1 & 0 \end{pmatrix},ana_n=2⋅3n−1{}= 2 \cdot 3^{n-1}−2n−1{}- 2^{n-1}

解説

行列を定める。 求めたいのは (an+1an)\begin{pmatrix} a_{n+1} \\ a_n \end{pmatrix} を (anan−1)\begin{pmatrix} a_n \\ a_{n-1} \end{pmatrix} から作る行列です。第 1 行は an+1a_{n+1}=5an{}= 5a_n−6an−1{}- 6a_{n-1}、第 2 行は ana_n=1⋅an{}= 1 \cdot a_n+0⋅an−1{}+ 0 \cdot a_{n-1} なので

A\displaystyle A=(5−610)‾\displaystyle {}= \underline{\begin{pmatrix} 5 & -6 \\ 1 & 0 \end{pmatrix}}

第 2 行の (10)\begin{pmatrix} 1 & 0 \end{pmatrix} は「ana_n をそのまま下にずらす」ための行です。これを入れておくことで、3 項間の漸化式が 2 成分のベクトルの 1 歩ずつの前進に変わります。

一般項を求める。 繰り返し使うと

(an+1an)\displaystyle \begin{pmatrix} a_{n+1} \\ a_n \end{pmatrix}=An−1(a2a1)\displaystyle {}= A^{n-1}\begin{pmatrix} a_2 \\ a_1 \end{pmatrix}=An−1(41)\displaystyle {}= A^{n-1}\begin{pmatrix} 4 \\ 1 \end{pmatrix}

AA の対角和は 55+0{}+ 0=5{}= 5、行列式は 5×05 \times 0−(−6)×1{}- (-6) \times 1=6{}= 6 なので

x2\displaystyle x^2−5x\displaystyle {}- 5x+6\displaystyle {}+ 6=0\displaystyle {}= 0⟹(x−2)(x−3)\displaystyle {}\Longrightarrow \quad (x - 2)(x - 3)=0\displaystyle {}= 0⟹x\displaystyle {}\Longrightarrow \quad x=2,\displaystyle {}= 2, 3\displaystyle \ 3

mm=n{}= n−1{}- 1 とおき、xmx^m を x2x^2−5x{}- 5x+6{}+ 6 で割った余りを pxpx+q{}+ q とすると

{2m=2p+q3m=3p+q\displaystyle \begin{cases} 2^m = 2p + q \\ 3^m = 3p + q \end{cases}⟹p\displaystyle {}\Longrightarrow \quad p=3m\displaystyle {}= 3^m−2m,\displaystyle {}- 2^m,q\displaystyle q=2m\displaystyle {}= 2^m−2p\displaystyle {}- 2p=3⋅2m\displaystyle {}= 3 \cdot 2^m−2⋅3m\displaystyle {}- 2 \cdot 3^m

AmA^m=pA{}= pA+qE{}+ qE なので、(41)\begin{pmatrix} 4 \\ 1 \end{pmatrix} にかけた結果の第 2 成分が ana_n です。AA の第 2 行は (10)\begin{pmatrix} 1 & 0 \end{pmatrix}、EE の第 2 行は (01)\begin{pmatrix} 0 & 1 \end{pmatrix} なので

an\displaystyle a_n=p(1×4+0×1)\displaystyle {}= p(1 \times 4 + 0 \times 1)+q(0×4+1×1)\displaystyle {}+ q(0 \times 4 + 1 \times 1)=4p\displaystyle {}= 4p+q\displaystyle {}+ q an\displaystyle a_n=4(3m−2m)\displaystyle {}= 4(3^m - 2^m)+3⋅2m\displaystyle {}+ 3 \cdot 2^m−2⋅3m\displaystyle {}- 2 \cdot 3^m=2⋅3m\displaystyle {}= 2 \cdot 3^m−2m\displaystyle {}- 2^m

mm=n{}= n−1{}- 1 に戻して

an\displaystyle a_n=2⋅3n−1‾\displaystyle {}= \underline{\rule[-0.0833em]{0em}{0.9474em}2 \cdot 3^{n-1}}−2n−1‾\displaystyle \underline{\rule[-0.0833em]{0em}{0.9474em}{}- 2^{n-1}}

検算します。a1a_1=2{}= 2−1{}- 1=1{}= 1、a2a_2=6{}= 6−2{}- 2=4{}= 4、a3a_3=18{}= 18−4{}- 4=14{}= 14。漸化式から a3a_3=5×4{}= 5 \times 4−6×1{}- 6 \times 1=14{}= 14 で一致します。

数列 第 7 章では、3 項間漸化式を x2x^2=5x{}= 5x−6{}- 6(特性方程式)の解を使って解きました。その特性方程式が、ここでは AA が背負う 2 次方程式として現れています。同じ方程式が、数列側では解法の道具、行列側では行列そのものの性質として出てくるわけです。偶然ではなく、AA の作り方が漸化式そのままだからです。

自己採点:
記録を読み込み中…

問19 ★★★

2 次の正方行列 AA が A3A^3=O{}= O を満たすならば、A2A^2=O{}= O であることを示しなさい。

つまずいたときは:
答えを見る
答え

行列式の積の法則から Δ\Delta=0{}= 0 を導き、ケーリー・ハミルトンの定理で A2A^2=(a+d)A{}= (a+d)A とすればよい

解説

AA=(abcd){}= \begin{pmatrix} a & b \\ c & d \end{pmatrix}、ss=a{}= a+d{}+ d、Δ\Delta=ad{}= ad−bc{}- bc とします。

第 1 段:Δ\Delta=0{}= 0 を示す。 行列式の積の法則(第 3 章)より det⁡(A3)\det(A^3)=(det⁡A)3{}= (\det A)^3 です。A3A^3=O{}= O なので det⁡(A3)\det(A^3)=0{}= 0、したがって

Δ3\displaystyle \Delta^3=0\displaystyle {}= 0⟹Δ\displaystyle {}\Longrightarrow \quad \Delta=0\displaystyle {}= 0

第 2 段:A2A^2 を AA で表す。 Δ\Delta=0{}= 0 なので、ケーリー・ハミルトンの定理は

A2\displaystyle A^2=sA\displaystyle {}= sA

となります(j15 の (∗)(\ast) と同じです)。

第 3 段:ss で場合分けする。 両辺に AA をかけると

A3\displaystyle A^3=sA2\displaystyle {}= sA^2=s(sA)\displaystyle {}= s(sA)=s2A\displaystyle {}= s^2A

仮定より A3A^3=O{}= O なので s2As^2A=O{}= O です。

ss≠0{}\neq 0 のとき、s2s^2≠0{}\neq 0 なので両辺を s2s^2 で割って AA=O{}= O。このとき A2A^2=O{}= O です。

ss=0{}= 0 のとき、第 2 段の式が A2A^2=0⋅A{}= 0 \cdot A=O{}= O となります。

どちらの場合も A2A^2=O{}= O が成り立ちます。

数の世界なら x3x^3=0{}= 0 から xx=0{}= 0、したがって x2x^2=0{}= 0 は当たり前です。行列では AA≠O{}\neq O でも A2A^2=O{}= O が起こる(j14)ので当たり前ではありませんが、2 次の行列では 2 乗より先まで粘ることはできない、というのがこの問題の結論です。一般に nn 次の正方行列では、AmA^m=O{}= O となる AA は必ず AnA^n=O{}= O を満たします。

自己採点:
記録を読み込み中…

問20 ★★★

AA,BB を 2 次の正方行列とし、CC=AB{}= AB−BA{}- BA とおく。

(1) CC の対角成分の和が 00 であることを示しなさい。

(2) C2C^2 が単位行列 EE の実数倍であることを示しなさい。

(3) 00 でない実数 kk に対して ABAB−BA{}- BA=kE{}= kE を満たす 2 次の正方行列 AA,BB は存在しないことを示しなさい。

つまずいたときは:
答えを見る
答え

(1) 成分計算で (1,1)(1,1) 成分と (2,2)(2,2) 成分の和が消える (2) (1) とケーリー・ハミルトンの定理より C2C^2=−(det⁡C)E{}= -(\det C)E (3) kEkE の対角成分の和は 2k2k≠0{}\neq 0 で (1) に矛盾する

解説

AA=(abcd){}= \begin{pmatrix} a & b \\ c & d \end{pmatrix}、BB=(pqrs){}= \begin{pmatrix} p & q \\ r & s \end{pmatrix} とおきます。

(1) 積を成分で書きます。

AB\displaystyle AB=(ap+braq+bscp+drcq+ds),\displaystyle {}= \begin{pmatrix} ap + br & aq + bs \\ cp + dr & cq + ds \end{pmatrix},BA\displaystyle BA=(pa+qcpb+qdra+scrb+sd)\displaystyle {}= \begin{pmatrix} pa + qc & pb + qd \\ ra + sc & rb + sd \end{pmatrix}

CC=AB{}= AB−BA{}- BA の対角成分は

C11\displaystyle C_{11}=(ap+br)\displaystyle {}= (ap + br)−(pa+qc)\displaystyle {}- (pa + qc)=br\displaystyle {}= br−qc,\displaystyle {}- qc,C22\displaystyle C_{22}=(cq+ds)\displaystyle {}= (cq + ds)−(rb+sd)\displaystyle {}- (rb + sd)=cq\displaystyle {}= cq−rb\displaystyle {}- rb

よって

C11\displaystyle C_{11}+C22\displaystyle {}+ C_{22}=(br−qc)\displaystyle {}= (br - qc)+(cq−rb)\displaystyle {}+ (cq - rb)=0\displaystyle {}= 0

が示されました。apap と papa、dsds と sdsd は数の積なので順番を変えても等しく、打ち消し合います。残るのは brbr と rbrb、qcqc と cqcq で、これらも符号が逆になって消えます。

(2) (1) より CC の対角和は 00 です。ケーリー・ハミルトンの定理を CC に適用すると

C2\displaystyle C^2−0⋅C\displaystyle {}- 0 \cdot C+(det⁡C)E\displaystyle {}+ (\det C)E=O\displaystyle {}= O⟹C2\displaystyle {}\Longrightarrow \quad C^2=−(det⁡C)E\displaystyle {}= -(\det C)E

det⁡C\det C は実数なので、C2C^2 は EE の実数倍です。

(3) ABAB−BA{}- BA=kE{}= kE とすると、左辺の対角成分の和は (1) より 00 です。一方、右辺 kEkE=(k00k){}= \begin{pmatrix} k & 0 \\ 0 & k \end{pmatrix} の対角成分の和は 2k2k です。等しい行列の対角成分の和は等しいので

0\displaystyle 0=2k\displaystyle {}= 2k⟹k\displaystyle {}\Longrightarrow \quad k=0\displaystyle {}= 0

これは kk≠0{}\neq 0 に反します。よってそのような AA,BB は存在しません。

この (3) には、大きな物語がつながっています。第 2 章の小話で、量子力学の行列力学を取り上げました。位置を表す行列 XX と運動量を表す行列 PP が

PX\displaystyle PX−XP\displaystyle {}- XP=h2πiE\displaystyle {}= \frac{h}{2\pi i}E

という関係を満たし、これが不確定性原理の正体だ、という話です。右辺は EE の 00 でない定数倍。いま示したのは、2 次の行列ではこの関係を実現できないということです。

同じ議論は nn 次の正方行列でもそのまま通用します(対角成分の和が 00 になる計算は、サイズによらないからです)。つまり、どんなに大きな有限サイズの行列を持ってきても、この関係は書けません。だから行列力学では、行や列が無限に続く行列が必要になるのです。ハイゼンベルクとボルンが扱っていたのが無限行列だったのは、こうした事情によります。

自己採点:
記録を読み込み中…

数学小話コーナー

「一般の場合を証明する必要は感じない」

この章の定理には 2 人の名前がついていますが、2 人が一緒に研究したわけではありません。それどころか、扱っていた対象すら別物でした。

先に到達したのはアイルランドのウィリアム・ローワン・ハミルトンです。彼は 1843 年に四元数という、複素数をさらに拡張した数を考え出しました。ii のほかに jj,kk という単位があり、ijij=k{}= k、jiji=−k{}= -k のようにかける順番で答えが変わる。人類が真剣に扱った、おそらく最初の非可換な数です。ハミルトンはこの四元数について、それぞれが 2 次方程式を満たすことを示していました。四元数は 2 次の複素行列で表せるので、これは今でいうケーリー・ハミルトンの定理の特別な場合にあたります。

一方、行列の側から同じ式に行き着いたのがイギリスのアーサー・ケイリーです。1858 年の『行列論回想録』で、2×22 \times 2 の行列について A2A^2−(a+d)A{}- (a + d)A+(ad−bc)E{}+ (ad - bc)E=O{}= O を成分計算で確かめました。第 1 章から何度も名前の出ている、行列論の創始者その人です。

おもしろいのはここからで、ケイリーは 3×33 \times 3 の場合も計算して確かめたうえで、こう書き添えています。「一般の nn 次の場合についても、これ以上証明する必要を感じない」。2×22 \times 2 と 3×33 \times 3 で成り立つのだから当然だろう、というわけです。

現代の目から見ると、ずいぶん大胆な態度に映ります。実際、一般の次数に通用する証明が与えられたのは、それから 20 年ほどあとのことでした。ふつうはフロベニウス(1878 年ごろ)の功績とされます(※引用の訳しかたには幅があり、一般の場合を最初に証明したのが誰かにも異説がある。※諸説あり)。

豆知識

この定理には、一見もっともらしい「証明」がよく知られています。det⁡(xE−A)\det(xE - A) を展開すると x2x^2−(a+d)x{}- (a+d)x+Δ{}+ \Delta になるので、この xx に AA を代入すれば det⁡(AE−A)\det(AE - A)=det⁡O{}= \det O=0{}= 0 だから成り立つ——というものです。もっともらしいのですが、これは誤りです。行き着いた 00 は数であって、零行列ではありません。どこがどう間違っているかは、厳密定義のページで解剖します。

フィボナッチ数を一瞬で取り出す

1,1, 1,\ 1, 2,\ 2, 3,\ 3, 5,\ 5, 8,\ 8, 13,\ 13, …\ \ldots と続くフィボナッチ数列は、前の 2 つを足して次を作ります(数列 第 7 章)。F100F_{100} を知りたければ 100 回足せばよい。単純ですが、100 回は 100 回です。

ここで、次の行列に登場してもらいます。

A\displaystyle A=(1110)\displaystyle {}= \begin{pmatrix} 1 & 1 \\ 1 & 0 \end{pmatrix}

これを (Fn+1Fn)\begin{pmatrix} F_{n+1} \\ F_n \end{pmatrix} にかけると (Fn+1+FnFn+1)\begin{pmatrix} F_{n+1} + F_n \\ F_{n+1} \end{pmatrix}=(Fn+2Fn+1){}= \begin{pmatrix} F_{n+2} \\ F_{n+1} \end{pmatrix} となって、数列が 1 つ進みます。つまり AA をかけることが「1 歩進む」ことです。nn 歩進めたければ AnA^n をかければよい。計算してみると、実に整った形が現れます。

An\displaystyle A^n=(Fn+1FnFnFn−1)\displaystyle {}= \begin{pmatrix} F_{n+1} & F_n \\ F_n & F_{n-1} \end{pmatrix}

さて、ここからが本題です。A100A^{100} を求めるのに、AA を 100 回かける必要はありません。A2A^2 を作り、それを 2 乗して A4A^4、また 2 乗して A8A^8……と進めば、7 回の 2 乗で A128A^{128} まで届きます。必要な A100A^{100}=A64A32A4{}= A^{64}A^{32}A^4 も、数回の積で作れる。100 回が 10 回ほどに減るわけです。この考え方を繰り返し 2 乗法といい、かける回数はおよそ log⁡2n\log_2 n で済みます。F1000000F_{1000000} のような桁でも、数十回の行列の積で取り出せてしまいます。

AA の背負う 2 次方程式を書いてみると x2x^2−x{}- x−1{}- 1=0{}= 0。この解 1±52\dfrac{1 \pm \sqrt{5}}{2} は、まぎれもなく黄金比です。公式5 のやり方で AnA^n を書き下せば、フィボナッチ数を nn の式で表すビネの公式がそのまま出てきます。足し算の規則しか与えていないのに、5\sqrt{5} が顔を出す——数列の章で不思議に思った点が、行列の側から見ると「AA が背負う 2 次方程式の解」として説明されるのです。

検索結果の順番は、行列の n 乗で決まった

1998 年、スタンフォード大学の大学院生だったラリー・ペイジとセルゲイ・ブリンが、Web ページの重要さを測る方法を考えました。ページランクです。

発想は、この章の例題7 とまったく同じです。Web をあてもなくリンクづたいに渡り歩く人を想像する。いまページ ii にいる確率を並べた列ベクトルを x⃗\vec{x} とし、リンクのつながり方から作った行列を AA とすると、1 回リンクをたどったあとの確率は Ax⃗A\vec{x}、nn 回たどったあとは Anx⃗A^n\vec{x} です。nn を大きくしていくと、この確率は出発点によらない一定の値に落ち着く。その落ち着き先の大きい順が、検索結果の順番になります。よく見られるページほど「重要」というわけです。

例題7 では店が 2 つでしたから、AA は 2×22 \times 2 でした。Web ページは、当時ですでに数億、いまでは数兆のけたです。AA は数十億行数十億列の、途方もない大きさの行列になります。

ではどうやって AnA^n を求めるのか。逆行列も行列式も、この大きさでは歯が立ちません。使われるのは、いたって素朴な方法です。適当な x⃗\vec{x} から始めて Ax⃗A\vec{x},A2x⃗A^2\vec{x},A3x⃗A^3\vec{x} …… と、ひたすらかけ続ける。値の変化が十分小さくなったところで打ち切る。これをべき乗法といい、数十回も回せば実用になる精度に達します。AnA^n を式として求めるのではなく、AA をかける操作を繰り返して近づける、という割り切りです。

豆知識

実際のページランクにはダンピング係数という工夫が入っていて、ふつうは 0.850.85 という値が使われます。「85 %85\ \% の確率でリンクをたどり、15 %15\ \% の確率でまったく別のページへ飛ぶ」という設定です。これがないと、リンクが外に出ていない袋小路のページに確率がたまってしまい、落ち着き先が求まりません。数列 第 9 章で確率の推移を扱ったときの「落ち着き先」の話を、行列の側から見たのがこの仕組みです(※現在の検索順位の決め方はずっと複雑で、ページランクはその一要素にすぎない)。

厳密定義(発展)

※ここは発展ページです。高校の範囲では「A2A^2=(a+d)A{}= (a+d)A−ΔE{}- \Delta E を使って次数を下げる」と使えれば十分ですが、そもそも aa+d{}+ d と Δ\Delta という係数はどこから来たのか、なぜ定理の逆だけが kEkE のときに壊れるのか、P−1APP^{-1}AP が対角行列になる PP はどうやって見つけるのかが気になる人のために、少し厳密な形にまとめておきます。

本文では公式1 を成分計算で確かめました。正しい証明ですが、これでは「たまたまそうなった」ようにしか見えません。係数 aa+d{}+ d と Δ\Delta の出どころを示すところから始めます。第 3 章の厳密定義(行列式の積の法則、正則の言いかえ)と第 4 章の厳密定義(1 次独立と解集合)を前提にします。

特性多項式

定義1:対角和と特性多項式

nn 次の正方行列 AA=(aij){}= (a_{ij}) に対し、対角成分の和

tr⁡A\displaystyle \operatorname{tr} A=a11\displaystyle {}= a_{11}+a22\displaystyle {}+ a_{22}+⋯\displaystyle {}+ \cdots+ann\displaystyle {}+ a_{nn}

を AA の対角和(トレース)という。また、xx を変数として

φA(x)\displaystyle \varphi_A(x)=det⁡(xE−A)\displaystyle {}= \det(xE - A)

で定まる xx の nn 次多項式を AA の特性多項式という。

2 次の場合に書き下してみます。

xE\displaystyle xE−A\displaystyle {}- A=(x−a−b−cx−d)\displaystyle {}= \begin{pmatrix} x - a & -b \\ -c & x - d \end{pmatrix} φA(x)\displaystyle \varphi_A(x)=(x−a)(x−d)\displaystyle {}= (x - a)(x - d)−(−b)(−c)\displaystyle {}- (-b)(-c)=x2\displaystyle {}= x^2−(a+d)x\displaystyle {}- (a + d)x+(ad−bc)\displaystyle {}+ (ad - bc)

本文で「AA が背負う 2 次方程式」と呼んでいたものの正体です。係数 aa+d{}+ d と Δ\Delta は、どこからともなく降ってきたのではなく、det⁡(xE−A)\det(xE - A) を展開すれば自動的に現れます。xn−1x^{n-1} の係数が −tr⁡A-\operatorname{tr} A、定数項が (−1)ndet⁡A(-1)^n \det A になるのは、一般の nn 次でも変わりません。

定理1:ケーリー・ハミルトンの定理(2 次)

2 次の正方行列 AA に対し φA(A)\varphi_A(A)=O{}= O である。すなわち

A2\displaystyle A^2−(tr⁡A)A\displaystyle {}- (\operatorname{tr} A)A+(det⁡A)E\displaystyle {}+ (\det A)E=O\displaystyle {}= O

証明 AA=(abcd){}= \begin{pmatrix} a & b \\ c & d \end{pmatrix} とする。第 3 章 定理2 の証明で A~\tilde{A}=(d−b−ca){}= \begin{pmatrix} d & -b \\ -c & a \end{pmatrix} が AA~A\tilde{A}=A~A{}= \tilde{A}A=(det⁡A)E{}= (\det A)E を満たすことを示した。ここで

A~\displaystyle \tilde{A}=(a+d00a+d)\displaystyle {}= \begin{pmatrix} a + d & 0 \\ 0 & a + d \end{pmatrix}−(abcd)\displaystyle {}- \begin{pmatrix} a & b \\ c & d \end{pmatrix}=(tr⁡A)E\displaystyle {}= (\operatorname{tr} A)E−A\displaystyle {}- A

であるから、AA~A\tilde{A}=(det⁡A)E{}= (\det A)E に代入して

A{(tr⁡A)E−A}\displaystyle A\{(\operatorname{tr} A)E - A\}=(det⁡A)E\displaystyle {}= (\det A)E

左辺を展開すると (tr⁡A)A(\operatorname{tr} A)A−A2{}- A^2 である。移項して

A2\displaystyle A^2−(tr⁡A)A\displaystyle {}- (\operatorname{tr} A)A+(det⁡A)E\displaystyle {}+ (\det A)E=O\displaystyle {}= O

を得る。(証明終)

本文の成分計算と同じことですが、こう書くと A~\tilde{A} の正体が (tr⁡A)E(\operatorname{tr} A)E−A{}- A であり、逆行列の公式とこの定理がまったく同じ 1 本の等式の言いかえだったことがはっきりします。第 3 章で「aa と dd を入れかえ、bb と cc の符号を変える」と手順で覚えた操作には、こういう意味があったわけです。

有名な誤った証明

この定理には、広く知られた「証明もどき」があります。

φA(x)\varphi_A(x)=det⁡(xE−A){}= \det(xE - A) の xx に AA を代入すると φA(A)\varphi_A(A)=det⁡(AE−A){}= \det(AE - A)=det⁡(A−A){}= \det(A - A)=det⁡O{}= \det O=0{}= 0。よって φA(A)\varphi_A(A)=O{}= O である。

一見すると鮮やかですが、誤りです。どこが壊れているのかを、はっきりさせておきます。

第一に、結論の型が合っていません。det⁡O\det O は数の 00 であって、零行列 OO ではありません。示したいのは「行列が OO に等しい」ことなのに、出てきたのは「数が 00 に等しい」こと。00 と OO を同じ記号のように扱った時点で、主張がすり替わっています。

第二に、そもそも代入が成り立っていません。φA(x)\varphi_A(x)=det⁡(xE−A){}= \det(xE - A) という等式は、xx が数であるときの等式です。左辺の φA\varphi_A は多項式なので、係数を行列の係数と読みかえて xx に行列を入れることができます(x2x^2 は A2A^2、定数項 qq は qEqE と読む)。ところが右辺の xExE は「数 xx と単位行列の積」であって、ここに行列を入れる意味は定められていません。かりに AEAE=A{}= A と読んだとしても、それは右辺の式を勝手に別のものに取りかえたにすぎず、等式が保たれる保証はどこにもありません。

同じ論法が通ってしまうと困る例を挙げれば、あやしさがはっきりします。ψ(x)\psi(x)=det⁡(xE−A){}= \det(xE - A)+1{}+ 1 という多項式を考えると、同じ手つきで ψ(A)\psi(A)=det⁡(A−A){}= \det(A - A)+1{}+ 1=1{}= 1、したがって「ψ(A)\psi(A)=O{}= O」——明らかに間違った結論が出ます。

一般の nn 次では、A~\tilde{A} にあたる余因子行列を使い、その成分が xx の多項式であることを利用して、係数を次数ごとに比べるという筋道で証明します。フロベニウスが与えたのもこの種の議論です。

対角行列に化ける条件

本文の公式7 では、都合のよい PP が問題文で与えられていました。では PP はどうやって見つけるのでしょうか。

定義2:固有値と固有ベクトル

nn 次の正方行列 AA に対し、u⃗\vec{u}≠0⃗{}\neq \vec{0} と実数 λ\lambda が

Au⃗\displaystyle A\vec{u}=λu⃗\displaystyle {}= \lambda\vec{u}

を満たすとき、λ\lambda を AA の固有値、u⃗\vec{u} を λ\lambda に対する固有ベクトルという。

Au⃗A\vec{u}=λu⃗{}= \lambda\vec{u} は (λE−A)u⃗(\lambda E - A)\vec{u}=0⃗{}= \vec{0} と書けるので、λ\lambda が固有値であることと、同次連立一次方程式 (λE−A)x⃗(\lambda E - A)\vec{x}=0⃗{}= \vec{0} が自明でない解をもつこと、すなわち φA(λ)\varphi_A(\lambda)=det⁡(λE−A){}= \det(\lambda E - A)=0{}= 0 は同値である(第 4 章 公式5)。

固有ベクトルとは、AA をかけても向きが変わらず、長さだけが λ\lambda 倍される方向のことです。本文で「方眼紙の目盛りを引き直す」と言ったときの、新しい軸の方向がこれにあたります。なお高校の教科書ではこの語を使わないので、本文では「x2x^2−(a+d)x{}- (a+d)x+Δ{}+ \Delta=0{}= 0 の解」と呼んできました。同じものです。

定理2:相異なる 2 つの実数解をもてば対角化できる

2 次の正方行列 AA の特性多項式 φA(x)\varphi_A(x) が相異なる 2 つの実数解 α,\alpha, β\ \beta をもつとする。このとき正則な行列 PP が存在して

P−1AP\displaystyle P^{-1}AP=(α00β)\displaystyle {}= \begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix}

となる。

証明 φA(α)\varphi_A(\alpha)=det⁡(αE−A){}= \det(\alpha E - A)=0{}= 0 であるから αE\alpha E−A{}- A は正則でない。よって第 4 章 公式5 により (αE−A)u⃗(\alpha E - A)\vec{u}=0⃗{}= \vec{0} を満たす u⃗\vec{u}≠0⃗{}\neq \vec{0} が存在する。同様に (βE−A)v⃗(\beta E - A)\vec{v}=0⃗{}= \vec{0} を満たす v⃗\vec{v}≠0⃗{}\neq \vec{0} が存在する。すなわち Au⃗A\vec{u}=αu⃗{}= \alpha\vec{u},Av⃗A\vec{v}=βv⃗{}= \beta\vec{v} である。

u⃗\vec{u} と v⃗\vec{v} が 1 次独立であることを示す。v⃗\vec{v}=tu⃗{}= t\vec{u}(tt≠0{}\neq 0)と書けたとすると

βv⃗\displaystyle \beta\vec{v}=Av⃗\displaystyle {}= A\vec{v}=tAu⃗\displaystyle {}= tA\vec{u}=tαu⃗\displaystyle {}= t\alpha\vec{u}=αv⃗\displaystyle {}= \alpha\vec{v}

より (β−α)v⃗(\beta - \alpha)\vec{v}=0⃗{}= \vec{0} であり、v⃗\vec{v}≠0⃗{}\neq \vec{0} から β\beta=α{}= \alpha となって仮定に反する。よって 1 次独立である。

u⃗,\vec{u}, v⃗\ \vec{v} をこの順に列として並べた行列を PP=(u⃗  v⃗){}= (\vec{u}\ \ \vec{v}) とおく。2 つの列が 1 次独立だから PP は正則である(第 3 章 定理2)。積を列ごとに計算すると

AP\displaystyle AP=(Au⃗  Av⃗)\displaystyle {}= (A\vec{u}\ \ A\vec{v})=(αu⃗  βv⃗)\displaystyle {}= (\alpha\vec{u}\ \ \beta\vec{v})=P(α00β)\displaystyle {}= P\begin{pmatrix} \alpha & 0 \\ 0 & \beta \end{pmatrix}

左から P−1P^{-1} をかけて結論を得る。(証明終)

PP の作り方まで証明が教えてくれています。φA(x)\varphi_A(x)=0{}= 0 の解ごとに Au⃗A\vec{u}=λu⃗{}= \lambda\vec{u} を解き、出てきたベクトルを列に並べる。本文の例題6 で与えられた PP=(211−1){}= \begin{pmatrix} 2 & 1 \\ 1 & -1 \end{pmatrix} も、この手順で作られたものです。実際 A(21)A\begin{pmatrix} 2 \\ 1 \end{pmatrix}=(105){}= \begin{pmatrix} 10 \\ 5 \end{pmatrix}=5(21){}= 5\begin{pmatrix} 2 \\ 1 \end{pmatrix}、A(1−1)A\begin{pmatrix} 1 \\ -1 \end{pmatrix}=(2−2){}= \begin{pmatrix} 2 \\ -2 \end{pmatrix}=2(1−1){}= 2\begin{pmatrix} 1 \\ -1 \end{pmatrix} となっています。

逆に、特性多項式が重解をもつ場合は、対角行列に化けるとは限りません。本文の例題5 の AA=(41−12){}= \begin{pmatrix} 4 & 1 \\ -1 & 2 \end{pmatrix} がその例です。もし P−1APP^{-1}AP=3E{}= 3E となる PP があったとすれば AA=P(3E)P−1{}= P(3E)P^{-1}=3E{}= 3E ですが、AA≠3E{}\neq 3E なので、そんな PP はありません。あのとき二項定理という別の道を通ったのは、そういう事情です。

対角和と行列式は PP で挟んでも変わらない

定理3:特性多項式の不変性

PP を正則な行列とすると

φP−1AP(x)\displaystyle \varphi_{P^{-1}AP}(x)=φA(x)\displaystyle {}= \varphi_A(x)

である。とくに tr⁡(P−1AP)\operatorname{tr}(P^{-1}AP)=tr⁡A{}= \operatorname{tr} A,det⁡(P−1AP)\det(P^{-1}AP)=det⁡A{}= \det A が成り立つ。

証明 xExE=P−1(xE)P{}= P^{-1}(xE)P であるから

xE\displaystyle xE−P−1AP\displaystyle {}- P^{-1}AP=P−1(xE−A)P\displaystyle {}= P^{-1}(xE - A)P

行列式の積の法則(第 3 章 定理3)より

det⁡(xE−P−1AP)\displaystyle \det(xE - P^{-1}AP)=det⁡(P−1)det⁡(xE−A)det⁡(P)\displaystyle {}= \det(P^{-1})\det(xE - A)\det(P)

det⁡(P−1)det⁡(P)\det(P^{-1})\det(P)=det⁡(P−1P){}= \det(P^{-1}P)=det⁡E{}= \det E=1{}= 1 であるから、右辺は det⁡(xE−A)\det(xE - A) に等しい。

多項式として等しいのだから、各次の係数も等しい。xx の係数を比べて tr⁡(P−1AP)\operatorname{tr}(P^{-1}AP)=tr⁡A{}= \operatorname{tr} A、定数項を比べて det⁡(P−1AP)\det(P^{-1}AP)=det⁡A{}= \det A を得る。(証明終)

第 3 章の実践問題で det⁡(A−1BA)\det(A^{-1}BA)=det⁡B{}= \det B を示したのは、この定理の定数項の部分でした。PP で挟む操作は、本文の言い方でいえば「方眼紙の目盛りを引き直す」こと。目盛りを変えても、対角和と行列式という 2 つの量は動かない——だからこそ、AA と DD が同じ特性多項式を背負うのです。本文の例題6 で DD の対角成分 55,22 が x2x^2−7x{}- 7x+10{}+ 10=0{}= 0 の解になっていたのは、この定理の帰結です。

定理の逆が壊れる本当の理由

本文の公式3 で、AA=kE{}= kE のときだけ定理の逆が成り立たないことを見ました。その理由を、多項式の言葉で言い直しておきます。

定義3:最小多項式

正方行列 AA に対し、f(A)f(A)=O{}= O を満たす 00 でない多項式 ff のうち、最高次の係数が 11 で次数が最小のものを AA の最小多項式という。

ケーリー・ハミルトンの定理により φA(A)\varphi_A(A)=O{}= O ですから、2 次の正方行列の最小多項式の次数は 11 か 22 です。

定理4:最小多項式が 1 次であることと AA=kE{}= kE は同値

2 次の正方行列 AA について、AA の最小多項式が 1 次であることと、ある実数 kk によって AA=kE{}= kE と書けることは同値である。

証明 最小多項式が 1 次ならば、それは xx−k{}- k の形であり、AA−kE{}- kE=O{}= O すなわち AA=kE{}= kE である。逆に AA=kE{}= kE ならば xx−k{}- k が AA を消すから、最小多項式は 1 次である。(証明終)

これで本文の公式3 が言い直せます。A2A^2−pA{}- pA+qE{}+ qE=O{}= O とは「x2x^2−px{}- px+q{}+ q が AA を消す」ということです。AA≠kE{}\neq kE ならば最小多項式は 2 次で、それは φA\varphi_A 自身です。最高次の係数が 11 の 2 次式で AA を消すものは最小多項式に限られるので、x2x^2−px{}- px+q{}+ q=φA(x){}= \varphi_A(x)、つまり pp=tr⁡A{}= \operatorname{tr} A,qq=det⁡A{}= \det A が確定します。

一方 AA=kE{}= kE のときは、最小多項式が 1 次の xx−k{}- k です。xx−k{}- k を因数に含む 2 次式なら何でも AA を消してしまうので、pp,qq は 1 通りに決まりません。本文で AA=2E{}= 2E が A2A^2−5A{}- 5A+6E{}+ 6E=O{}= O を満たしたのは、x2x^2−5x{}- 5x+6{}+ 6=(x−2)(x−3){}= (x - 2)(x - 3) が xx−2{}- 2 を含んでいたからです。AA の対角和が 55 だったからではありません。

「A2A^2−pA{}- pA+qE{}+ qE=O{}= O から pp=a{}= a+d{}+ d を読み取ってはいけない」という注意の中身は、x2x^2−px{}- px+q{}+ q が最小多項式の倍数でありさえすればよく、最小多項式そのものである必要はないということです。

この先

一般の nn 次でも φA(A)\varphi_A(A)=O{}= O が成り立ちます(ケーリー・ハミルトンの定理)。nn 次の行列の nn+1{}+ 1 乗以上はすべて nn−1{}- 1 次以下に落とせる、という意味で、次数下げは次数がいくつでも使えます。

対角行列に化けない行列をどこまで簡単な形にできるか、という問いにはジョルダン標準形という答えが用意されています。本文の例題5 で使った αE\alpha E+N{}+ N(N2N^2=O{}= O)という形は、2 次のジョルダン標準形そのものです。NN のように何乗かすると OO になる行列をべき零行列といい、対角行列とべき零行列の和に分ける、というのが一般の筋書きになります。

もう 1 つ、AnA^n の先には行列の指数関数があります。

eA\displaystyle e^{A}=E\displaystyle {}= E+A\displaystyle {}+ A+12!A2\displaystyle {}+ \frac{1}{2!}A^2+13!A3\displaystyle {}+ \frac{1}{3!}A^3+⋯\displaystyle {}+ \cdots

数の exe^x の級数(微分積分で学ぶ展開)に、そのまま行列を入れたものです。この収束を示すには級数の議論が要りますが、ひとたび定義してしまえば、微分方程式 dx⃗dt\dfrac{d\vec{x}}{dt}=Ax⃗{}= A\vec{x} の解が x⃗(t)\vec{x}(t)=eAtx⃗(0){}= e^{At}\vec{x}(0) と書けます。ばねの振動、電気回路、人口の変化——連立した微分方程式は、行列の累乗の言葉で一気に片づくのです。AnA^n を求める技術は、その入口にあたります。

次の第 6 章からは、行列を計算の対象としてではなく、平面上の点を動かす仕掛けとして見直します。第 2 章の厳密定義で証明した「積は写像の合成」という事実が、そこで主役になります。

この章の学習が終わったら

学習完了テストを受ける