重解とJordan細胞

同義語:double roots and Jordan blocks

概要

重解とJordan細胞(double roots and Jordan blocks)とは、三項間漸化式 $a_{n+2}=pa_{n+1}+qa_n$($q\ne0$)の特性方程式が重解 $\alpha$ をもつとき、一般項に $n$ が掛かって $(A+Bn)\alpha^n$ となる理由を、行列で説明する見方である。このときコンパニオン行列 $C=\begin{pmatrix}p&q\\1&0\end{pmatrix}$ は固有ベクトルが 1 方向しかなく対角化できないが、$P^{-1}CP=\begin{pmatrix}\alpha&1\\0&\alpha\end{pmatrix}$(Jordan 細胞)と変形できる。Jordan 細胞の $n$ 乗の右上の成分 $n\alpha^{n-1}$ が、一般項の $n$ の正体である。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 三項間漸化式の特性方程式, 漸化式の行列表示, 二項定理

この記事で考えること

三項間漸化式 $a_{n+2}=pa_{n+1}+qa_n$($q\ne0$)の特性方程式 $t^2=pt+q$ が重解 $\alpha$ をもつとき、一般項は $(A+Bn)\alpha^n$ の形になる(三項間漸化式の特性方程式)。

重解の漸化式を高校の方法で解く

$a_0=1$、$a_1=4$、$a_{n+2}=4a_{n+1}-4a_n$ とする。特性方程式 $t^2-4t+4=(t-2)^2=0$ は重解 $2$ をもつ。$a_{n+1}-2a_n$ は公比 $2$ の等比数列なので
$$ a_{n+1}-2a_n=2^n(a_1-2a_0)=2^n\cdot2=2^{n+1} $$
である。両辺を $2^{n+1}$ で割ると $\dfrac{a_{n+1}}{2^{n+1}}-\dfrac{a_n}{2^n}=1$ となり、$\dfrac{a_n}{2^n}$ は初項 $1$、公差 $1$ の等差数列である。よって $a_n=(n+1)2^n$ で、$1,4,12,32,80,\dots$ と続く。

答えに $n$ が掛かるのは、特性方程式の解が異なるとき($A\alpha^n+B\beta^n$)には見られなかった現象である。この $n$ はどこから来たのか。漸化式の行列表示 で見たように、漸化式は $v_n:=\begin{pmatrix}a_{n+1}\\a_n\end{pmatrix}$ とコンパニオン行列 $C:=\begin{pmatrix}p&q\\1&0\end{pmatrix}$ により $v_n=C^nv_0$ と書ける。この記事では、重解のとき $C$ が対角化できないことを示し、代わりに使う Jordan 細胞の $n$ 乗から $n$ が現れることを見る。

高校の計算大学の概念ボックス
等比数列の解が $\alpha^n$ の 1 種類しかない固有ベクトルが 1 方向しかなく、対角化できないthm-jb-not-diag
$a_n/\alpha^n$ が等差数列になるJordan 細胞の $n$ 乗 $\begin{pmatrix}\alpha^n&n\alpha^{n-1}\\0&\alpha^n\end{pmatrix}$lem-jb-power
一般項 $(A+Bn)\alpha^n$$P^{-1}CP=J$ と $C^n=PJ^nP^{-1}$thm-jb-general

行列の積・逆行列・固有値・固有ベクトル・固有多項式 $\chi_A(t)=\det(tE-A)$ の定義と計算は 漸化式の行列表示 のとおりとする。以下、数と言えば複素数とする。

対角化できる行列、できない行列

対角化可能と対角化

右上と左下の成分がどちらも $0$ の 2 次正方行列 $\begin{pmatrix}\lambda_1&0\\0&\lambda_2\end{pmatrix}$ を対角行列という。2 次正方行列 $A$ について、逆行列をもつ 2 次正方行列 $P$ と数 $\lambda_1,\lambda_2$ があって
$$ P^{-1}AP=\begin{pmatrix}\lambda_1&0\\0&\lambda_2\end{pmatrix} $$
となるとき、$A$ は対角化可能であるという。このような $P$ と対角行列 $P^{-1}AP$ を求めることを、$A$ の対角化という。

この定義は 漸化式の行列表示 と同じである。
対角化できれば $A^n=P\begin{pmatrix}\lambda_1^n&0\\0&\lambda_2^n\end{pmatrix}P^{-1}$ と $n$ 乗が計算できる。

対角化できる例

$S=\begin{pmatrix}2&1\\1&2\end{pmatrix}$ と、その固有ベクトル $\begin{pmatrix}1\\1\end{pmatrix}$(固有値 $3$)、$\begin{pmatrix}1\\-1\end{pmatrix}$(固有値 $1$)を並べた $P=\begin{pmatrix}1&1\\1&-1\end{pmatrix}$ を考える。$\det P=-1-1=-2$ なので $P^{-1}=\frac1{-2}\begin{pmatrix}-1&-1\\-1&1\end{pmatrix}=\frac12\begin{pmatrix}1&1\\1&-1\end{pmatrix}$ である。
$$ SP=\begin{pmatrix}3&1\\3&-1\end{pmatrix},\qquad P^{-1}SP=\frac12\begin{pmatrix}1&1\\1&-1\end{pmatrix}\begin{pmatrix}3&1\\3&-1\end{pmatrix}=\frac12\begin{pmatrix}6&0\\0&2\end{pmatrix}=\begin{pmatrix}3&0\\0&1\end{pmatrix}. $$

固有ベクトルが 1 方向しかない例

$U=\begin{pmatrix}1&1\\0&1\end{pmatrix}$ の固有多項式は $\chi_U(t)=(t-1)^2$ で、固有値は $1$ だけである。固有ベクトル $\begin{pmatrix}x\\y\end{pmatrix}$ は
$$ (U-E)\begin{pmatrix}x\\y\end{pmatrix}=\begin{pmatrix}0&1\\0&0\end{pmatrix}\begin{pmatrix}x\\y\end{pmatrix}=\begin{pmatrix}y\\0\end{pmatrix}=\mathbf 0 $$
を満たすもの、すなわち $y=0$ のものなので、$\begin{pmatrix}1\\0\end{pmatrix}$ の定数倍しかない。$U$ を順に掛けると $U^2=\begin{pmatrix}1&2\\0&1\end{pmatrix}$、$U^3=\begin{pmatrix}1&3\\0&1\end{pmatrix}$ となり、右上に $n$ が現れる。

対角化できるかどうかは、固有ベクトルが何本とれるかで決まる。まず次を示す。

対角化する行列の列は固有ベクトル

逆行列をもつ $P$ で $P^{-1}AP=\begin{pmatrix}\lambda_1&0\\0&\lambda_2\end{pmatrix}$ となるならば、$P$ の第 1 列は固有値 $\lambda_1$ に、第 2 列は固有値 $\lambda_2$ に属する $A$ の固有ベクトルである。

列ごとに比べる

$P=\begin{pmatrix}x_1&x_2\\y_1&y_2\end{pmatrix}$ とし、列を $u:=\begin{pmatrix}x_1\\y_1\end{pmatrix}$、$w:=\begin{pmatrix}x_2\\y_2\end{pmatrix}$ とおく。仮定の式に左から $P$ を掛けると
$$ AP=P\begin{pmatrix}\lambda_1&0\\0&\lambda_2\end{pmatrix}=\begin{pmatrix}\lambda_1x_1&\lambda_2x_2\\\lambda_1y_1&\lambda_2y_2\end{pmatrix} $$
である。行列の積の定義から、$AP$ の第 1 列は $Au$、第 2 列は $Aw$ である。列を比べて $Au=\lambda_1u$、$Aw=\lambda_2w$ を得る。$P$ は逆行列をもつので $\det P\ne0$ であり、$u=\mathbf 0$ なら $\det P=0\cdot y_2-x_2\cdot0=0$ となってしまうから $u\ne\mathbf 0$ である。同様に $w\ne\mathbf 0$ である。$\square$

主定理 1:重解のコンパニオン行列は対角化できない

重解のコンパニオン行列は対角化できない

特性方程式 $t^2=pt+q$($q\ne0$)が重解 $\alpha$ をもつとき、コンパニオン行列 $C=\begin{pmatrix}p&q\\1&0\end{pmatrix}$ は対角化可能でない。

対角化できたとして矛盾を導く

方針:対角化できたとすると $C$ が $\alpha E$ に等しくなってしまうことを示す。
段 1:逆行列をもつ $P$ で $P^{-1}CP=\begin{pmatrix}\lambda_1&0\\0&\lambda_2\end{pmatrix}$ となったとする。lem-jb-columns により、$\lambda_1,\lambda_2$ は $C$ の固有値である。
段 2:漸化式の行列表示 の主定理 1(特性方程式とコンパニオン行列の固有多項式)により、$C$ の固有値は特性方程式の解であり、ここでは $\alpha$ だけである。よって $\lambda_1=\lambda_2=\alpha$ で、$P^{-1}CP=\alpha E$ である。
段 3:左から $P$、右から $P^{-1}$ を掛けると $C=P(\alpha E)P^{-1}=\alpha PP^{-1}=\alpha E$ となる。しかし $C$ の左下の成分は $1$、$\alpha E=\begin{pmatrix}\alpha&0\\0&\alpha\end{pmatrix}$ の左下の成分は $0$ なので、$C\ne\alpha E$ である。これは矛盾である。$\square$

証明で使ったのは、固有値が $\alpha$ だけ(重解)であること(段 1・2)と、$C$ の左下の成分が $0$ でないこと(段 3)の 2 つである。前者を外すと、異なる 2 解の場合には対角化できる(漸化式の行列表示)。後者を外すと結論は崩れ、重解をもつ行列でも対角化できることがある(ex-jb-scalar)。

Jordan 細胞

対角行列の代わりに、次の形の行列を目標にする。

Jordan 細胞

数 $\alpha$ について、2 次正方行列
$$ J(\alpha):=\begin{pmatrix}\alpha&1\\0&\alpha\end{pmatrix} $$
を固有値 $\alpha$ の(2 次の)Jordan 細胞という。

ex-jb-shear の $U$ は $J(1)$ である。$J(\alpha)$ の固有多項式は $(t-\alpha)^2$ で、thm-jb-not-diag と同じ理由(左下ではなく右上の成分 $1$ を使う)で対角化できない。

Jordan 細胞の冪を計算する

$J(2)=\begin{pmatrix}2&1\\0&2\end{pmatrix}$ を順に掛ける。
$$ J(2)^2=\begin{pmatrix}2\cdot2+1\cdot0&2\cdot1+1\cdot2\\0&2\cdot2\end{pmatrix}=\begin{pmatrix}4&4\\0&4\end{pmatrix},\qquad J(2)^3=\begin{pmatrix}4&4\\0&4\end{pmatrix}\begin{pmatrix}2&1\\0&2\end{pmatrix}=\begin{pmatrix}8&12\\0&8\end{pmatrix},\qquad J(2)^4=\begin{pmatrix}16&32\\0&16\end{pmatrix}. $$
右上の成分 $1,4,12,32$ は $n\cdot2^{n-1}$($n=1,2,3,4$)である。対角成分は $2^n$ のままである。

Jordan 細胞の n 乗

$\alpha\ne0$ のとき、すべての $n\ge0$ で
$$ J(\alpha)^n=\begin{pmatrix}\alpha^n&n\alpha^{n-1}\\0&\alpha^n\end{pmatrix} $$
である。

数学的帰納法

$n=0$ のとき:右辺は $\begin{pmatrix}1&0\cdot\alpha^{-1}\\0&1\end{pmatrix}=E$ で、$J(\alpha)^0=E$ と等しい($\alpha\ne0$ なので $\alpha^{-1}$ は意味をもつ)。
$n$ で正しいとする。冪の定義 $J(\alpha)^{n+1}=J(\alpha)^nJ(\alpha)$ と積の定義により
$$ J(\alpha)^{n+1}=\begin{pmatrix}\alpha^n&n\alpha^{n-1}\\0&\alpha^n\end{pmatrix}\begin{pmatrix}\alpha&1\\0&\alpha\end{pmatrix} =\begin{pmatrix}\alpha^n\cdot\alpha+n\alpha^{n-1}\cdot0&\alpha^n\cdot1+n\alpha^{n-1}\cdot\alpha\\0\cdot\alpha+\alpha^n\cdot0&0\cdot1+\alpha^n\cdot\alpha\end{pmatrix} =\begin{pmatrix}\alpha^{n+1}&(n+1)\alpha^n\\0&\alpha^{n+1}\end{pmatrix} $$
である(右上は $\alpha^n+n\alpha^n=(n+1)\alpha^n$)。これは $n+1$ のときの式である。$\square$

帰納法の 1 段ごとに、右上の成分には「$\alpha$ 倍したもの」に「対角成分 $\alpha^n$」が 1 つずつ足されていく。$n$ 回で $\alpha^{n-1}$ が $n$ 個たまる、というのが $n$ の正体である。

右上の成分は高校の等差数列と同じ計算

$J(\alpha)^n$ の右上の成分を $x_n$ とする。lem-jb-power の証明の積の計算から、右上の成分は $x_{n+1}=\alpha^n\cdot1+x_n\cdot\alpha$、すなわち
$$ x_{n+1}=\alpha x_n+\alpha^n $$
を満たす。両辺を $\alpha^{n+1}$($\ne0$)で割ると
$$ \frac{x_{n+1}}{\alpha^{n+1}}-\frac{x_n}{\alpha^n}=\frac1\alpha $$
となり、$\frac{x_n}{\alpha^n}$ は初項 $\frac{x_0}{\alpha^0}=0$、公差 $\frac1\alpha$ の等差数列である。よって $\frac{x_n}{\alpha^n}=\frac n\alpha$、$x_n=n\alpha^{n-1}$ である。これは、ex-jb-review で $\frac{a_n}{2^n}$ が等差数列になったのと同じ計算である。

同じことを二項定理で見ることもできる。

二項定理で見る

$N:=\begin{pmatrix}0&1\\0&0\end{pmatrix}$ とおくと $J(\alpha)=\alpha E+N$ であり、
$$ N^2=\begin{pmatrix}0\cdot0+1\cdot0&0\cdot1+1\cdot0\\0&0\end{pmatrix}=O $$
である($O$ はすべての成分が $0$ の行列)。$\alpha E$ はどの行列とも交換する($(\alpha E)X=\alpha X=X(\alpha E)$)。二項定理の証明で使うのは分配法則と「掛ける順序を入れ替えてよいこと」だけなので、交換する 2 つの行列 $\alpha E$ と $N$ についても二項定理が成り立ち、
$$ J(\alpha)^n=(\alpha E+N)^n=\sum_{k=0}^n\binom nk\alpha^{n-k}N^k=\alpha^nE+n\alpha^{n-1}N $$
となる($k\ge2$ の項は $N^k=O$ で消える)。重解で現れる $n\alpha^{n-1}$ は、二項展開の 1 次の項である。

負の固有値の Jordan 細胞

$J(-1)=\begin{pmatrix}-1&1\\0&-1\end{pmatrix}$ では、lem-jb-power により $J(-1)^n=\begin{pmatrix}(-1)^n&n(-1)^{n-1}\\0&(-1)^n\end{pmatrix}$ である。$n=2$ で $\begin{pmatrix}1&-2\\0&1\end{pmatrix}$、$n=3$ で $\begin{pmatrix}-1&3\\0&-1\end{pmatrix}$ となり、直接掛けた結果 $J(-1)^2=\begin{pmatrix}1+0&-1-1\\0&1\end{pmatrix}$、$J(-1)^3=\begin{pmatrix}1&-2\\0&1\end{pmatrix}\begin{pmatrix}-1&1\\0&-1\end{pmatrix}=\begin{pmatrix}-1&1+2\\0&-1\end{pmatrix}$ と一致する。

主定理 2:Jordan 細胞への変形と一般項

重解の場合の標準形と一般項

特性方程式 $t^2=pt+q$($q\ne0$)が重解 $\alpha$ をもつとする($p^2+4q=0$、$\alpha=\frac p2$)。$P:=\begin{pmatrix}\alpha&1\\1&0\end{pmatrix}$ とおくと、$P$ は逆行列をもち
$$ P^{-1}CP=J(\alpha),\qquad C^n=PJ(\alpha)^nP^{-1} $$
である。漸化式 $a_{n+2}=pa_{n+1}+qa_n$ の解 $(a_n)$ は、すべての $n\ge0$ で
$$ a_n=a_0\alpha^n+n\,(a_1-\alpha a_0)\,\alpha^{n-1} $$
を満たす。特に $a_n=(A+Bn)\alpha^n$($A=a_0$、$B=\frac{a_1-\alpha a_0}\alpha$)と書ける。

固有ベクトルともう 1 本のベクトルを並べる

方針:$P$ の第 1 列は固有ベクトル、第 2 列は「$C-\alpha E$ を掛けると固有ベクトルになるベクトル」である。$CP=PJ(\alpha)$ を成分で確かめ、lem-jb-power を使う。
段 1(係数を $\alpha$ で表す):重解なので $t^2-pt-q=(t-\alpha)^2=t^2-2\alpha t+\alpha^2$ であり、係数を比べて $p=2\alpha$、$q=-\alpha^2$ である。$q\ne0$ から $\alpha\ne0$ である。したがって $C=\begin{pmatrix}2\alpha&-\alpha^2\\1&0\end{pmatrix}$ である。
段 2($CP=PJ(\alpha)$):
$$ CP=\begin{pmatrix}2\alpha&-\alpha^2\\1&0\end{pmatrix}\begin{pmatrix}\alpha&1\\1&0\end{pmatrix}=\begin{pmatrix}2\alpha^2-\alpha^2&2\alpha\\\alpha&1\end{pmatrix}=\begin{pmatrix}\alpha^2&2\alpha\\\alpha&1\end{pmatrix}, $$
$$ PJ(\alpha)=\begin{pmatrix}\alpha&1\\1&0\end{pmatrix}\begin{pmatrix}\alpha&1\\0&\alpha\end{pmatrix}=\begin{pmatrix}\alpha^2&\alpha+\alpha\\\alpha&1\end{pmatrix}=\begin{pmatrix}\alpha^2&2\alpha\\\alpha&1\end{pmatrix}. $$
両者は等しい。
段 3($P$ の逆行列):$\det P=\alpha\cdot0-1\cdot1=-1\ne0$ なので、2 次正方行列の逆行列の公式により
$$ P^{-1}=\frac1{-1}\begin{pmatrix}0&-1\\-1&\alpha\end{pmatrix}=\begin{pmatrix}0&1\\1&-\alpha\end{pmatrix} $$
である。段 2 の式に左から $P^{-1}$ を掛けて $P^{-1}CP=J(\alpha)$、右から $P^{-1}$ を掛けて $C=PJ(\alpha)P^{-1}$ を得る。
段 4($C^n=PJ(\alpha)^nP^{-1}$):$n=0$ では $PEP^{-1}=E$ である。$C^n=PJ(\alpha)^nP^{-1}$ ならば
$$ C^{n+1}=C^nC=PJ(\alpha)^nP^{-1}PJ(\alpha)P^{-1}=PJ(\alpha)^{n+1}P^{-1} $$
である(真ん中の $P^{-1}P=E$ が消える)。
段 5(一般項):$v_n=C^nv_0=PJ(\alpha)^nP^{-1}v_0$ を右から計算する。
$$ P^{-1}v_0=\begin{pmatrix}0&1\\1&-\alpha\end{pmatrix}\begin{pmatrix}a_1\\a_0\end{pmatrix}=\begin{pmatrix}a_0\\a_1-\alpha a_0\end{pmatrix},\qquad J(\alpha)^nP^{-1}v_0=\begin{pmatrix}\alpha^na_0+n\alpha^{n-1}(a_1-\alpha a_0)\\\alpha^n(a_1-\alpha a_0)\end{pmatrix} $$
(2 つ目は lem-jb-power による)。$v_n$ の下の成分 $a_n$ は、$P$ の下の行 $(1\ \ 0)$ とこのベクトルの積、つまりこのベクトルの上の成分であり、主張の式を得る。最後に $n\alpha^{n-1}(a_1-\alpha a_0)=n\cdot\frac{a_1-\alpha a_0}\alpha\cdot\alpha^n$ と書き直せば $(A+Bn)\alpha^n$ の形になる。$\square$

Jordan 細胞の $n$ 乗を展開するときに使う二項定理は 二項定理と組合せの恒等式 で扱う。

P の 2 本の列の意味

$P$ の第 1 列 $\begin{pmatrix}\alpha\\1\end{pmatrix}$ は固有値 $\alpha$ の固有ベクトルである。第 2 列 $e:=\begin{pmatrix}1\\0\end{pmatrix}$ は固有ベクトルではないが、
$$ (C-\alpha E)\,e=\begin{pmatrix}2\alpha-\alpha&-\alpha^2\\1&-\alpha\end{pmatrix}\begin{pmatrix}1\\0\end{pmatrix}=\begin{pmatrix}\alpha\\1\end{pmatrix} $$
となり、$C-\alpha E$ を 1 回掛けると固有ベクトルになる。言い換えると $Ce=\alpha e+\begin{pmatrix}\alpha\\1\end{pmatrix}$ で、「$\alpha$ 倍」に「固有ベクトルの分だけずれる」が加わる。このずれが $n$ 回たまって $n\alpha^{n-1}$ になる。このような $e$ を一般化固有ベクトルという。

重解の例(1):C の 5 乗

$a_{n+2}=4a_{n+1}-4a_n$ では $C=\begin{pmatrix}4&-4\\1&0\end{pmatrix}$、$\alpha=2$、$P=\begin{pmatrix}2&1\\1&0\end{pmatrix}$、$P^{-1}=\begin{pmatrix}0&1\\1&-2\end{pmatrix}$ である。ex-jb-powers の続きで $J(2)^5=\begin{pmatrix}32&80\\0&32\end{pmatrix}$ なので
$$ C^5=PJ(2)^5P^{-1}=\begin{pmatrix}2&1\\1&0\end{pmatrix}\begin{pmatrix}32&80\\0&32\end{pmatrix}\begin{pmatrix}0&1\\1&-2\end{pmatrix}=\begin{pmatrix}64&192\\32&80\end{pmatrix}\begin{pmatrix}0&1\\1&-2\end{pmatrix}=\begin{pmatrix}192&-320\\80&-128\end{pmatrix} $$
である。$v_0=\begin{pmatrix}4\\1\end{pmatrix}$($a_0=1$、$a_1=4$)に掛けると $C^5v_0=\begin{pmatrix}768-320\\320-128\end{pmatrix}=\begin{pmatrix}448\\192\end{pmatrix}$ で、$a_6=448=7\cdot2^6$、$a_5=192=6\cdot2^5$ となり、ex-jb-review の $(n+1)2^n$ と一致する。

図1:重解の例で v を 2 の冪で割った点は、固有ベクトルと平行な直線の上を毎回同じ幅だけ進む 図1:重解の例で v を 2 の冪で割った点は、固有ベクトルと平行な直線の上を毎回同じ幅だけ進む
図 1 は ex-jb-example-1 の $v_n=\begin{pmatrix}(n+2)2^{n+1}\\(n+1)2^n\end{pmatrix}$ を $2^n$ で割った点 $(2n+4,\ n+1)$ である。異なる 2 解の絶対値が異なる場合なら、絶対値の大きい方の解の冪で割った値は一定の値に近づく(図 2 の青)。2 解の絶対値が等しいときは、周期的にくり返して近づかないこともある(虚数解の周期 6 の例、三項間漸化式の特性方程式)。重解では固有ベクトルの方向 $(2,1)$ に毎回同じ幅だけ進み続け、止まらない。この「等間隔の進み」が、高校の計算の「$a_n/2^n$ は等差数列」である。

重解の例(2):負の重解

$a_{n+2}=-2a_{n+1}-a_n$、$a_0=a_1=1$ とする。特性方程式 $(t+1)^2=0$ の重解は $\alpha=-1$ で、$C=\begin{pmatrix}-2&-1\\1&0\end{pmatrix}$、$P=\begin{pmatrix}-1&1\\1&0\end{pmatrix}$、$P^{-1}=\begin{pmatrix}0&1\\1&1\end{pmatrix}$ である。確かめると
$$ CP=\begin{pmatrix}2-1&-2\\-1&1\end{pmatrix}=\begin{pmatrix}1&-2\\-1&1\end{pmatrix},\qquad PJ(-1)=\begin{pmatrix}-1&1\\1&0\end{pmatrix}\begin{pmatrix}-1&1\\0&-1\end{pmatrix}=\begin{pmatrix}1&-2\\-1&1\end{pmatrix} $$
で等しい。thm-jb-general の式は $a_n=(-1)^n+n\,(1+1)(-1)^{n-1}=(1-2n)(-1)^n$ で、$1,1,-3,5,-7,\dots$ となる。

図2:重解の例では割った値が直線的に増え、異なる 2 解の例では一定の値に近づく 図2:重解の例では割った値が直線的に増え、異なる 2 解の例では一定の値に近づく
図 2 の赤は重解の例の $a_n/2^n=n+1$、青は異なる 2 解の例 $a_{n+2}=2a_{n+1}+3a_n$、$a_0=0$、$a_1=4$ の $a_n/3^n=1-(-\frac13)^n$ である。

例と反例

外した仮定崩れる結論ボックス
特性方程式が異なる 2 解をもつ解が特性方程式の解の冪の和 $A\alpha^n+B\beta^n$ で書けるex-jb-only-geometric
行列がコンパニオン行列(左下が $1$)固有値が重なれば対角化できないex-jb-scalar
反例:等比数列だけでは初期値に合わない

$a_{n+2}=4a_{n+1}-4a_n$、$a_0=1$、$a_1=4$ は、「特性方程式が異なる 2 解をもつ」という仮定だけを満たさない。使える等比数列は $2^n$ だけなので $a_n=A\cdot2^n$ としてみると、$a_0=1$ から $A=1$、すると $a_1=2\ne4$ となる。「解は特性方程式の解の冪の和で書ける」という結論が破れる。行列で言えば、固有ベクトルが $\begin{pmatrix}2\\1\end{pmatrix}$ の方向しかなく、$v_0=\begin{pmatrix}4\\1\end{pmatrix}$ がその方向にないため、固有ベクトルだけでは $v_0$ を表せない。

反例:固有値が重なっても対角化できる行列

$2E=\begin{pmatrix}2&0\\0&2\end{pmatrix}$ の固有多項式は $(t-2)^2$ で、固有値 $2$ が重なっている。しかし $2E$ ははじめから対角行列なので($P=E$)、対角化可能である。どんなベクトル $v$ も $2Ev=2v$ を満たし、固有ベクトルは平面のすべての方向にある。$2E$ は左下の成分が $0$ でコンパニオン行列の形でないので、thm-jb-not-diag の仮定を満たさない。重解であることだけでは、対角化できないとは言えない。

さらに先へ

  • Jordan 標準形:どんな正方行列も、Jordan 細胞を対角線に並べた行列に $P^{-1}AP$ の形で変形できる(Jordan標準形。本記事では証明しない。Axl24 8.46、p. 324)。$m$ 次の Jordan 細胞の $n$ 乗からは $\binom nj\alpha^{n-j}$($j< m$)が現れ、$k$ 項の漸化式で $m$ 重の解から $n\alpha^n,\dots,n^{m-1}\alpha^n$ が現れる理由になる。
  • 微分方程式:$y''=py'+qy$ で特性方程式が重解 $\alpha$ をもつとき、解は $(A+Bx)e^{\alpha x}$ になる(微分方程式としての指数関数・三角関数)。$x$ が掛かる理由も同じ Jordan 細胞である。
  • 余りによる計算:$C^n$ は $t^n$ を $(t-\alpha)^2$ で割った余りからも求まり、重解では余りの係数に微分が現れる(行列のn乗と割り算の余り)。

関連項目

参考文献

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する