重解のない固有方程式を持つ数列の行列表現とその成分
はじめに
この記事では重解を持たない固有方程式$f(x)=0$を持つ数列$a_n$の行列による表現をその各成分まで具体的に求めます。
定義と主張
まず色々と言葉と記号を定めておきます。
定数係数の線形漸化式
$$a_{n+k}=c_{k-1}a_{n+k-1}+c_{k-2}a_{n+k-2}+\cdots+c_0a_n$$
を満たす数列$a_n$に対し固有多項式$f(x)$を
$$f(x)=x^k-c_{k-1}x^{k-1}-c_{k-2}x^{k-2}-\cdots-c_0$$
と定め、固有多項式$f(x)$についての方程式$f(x)=0$を$a_n$の固有方程式と言う。
上のような数列$a_n$に対し$k$次正方行列$M$を
$$M^T=
\begin{pmatrix}
0&1&0&\c&0\\
0&0&1&\c&0\\
\v&\v&\v&\d&\v\\
0&0&0&\c&1\\
c_0&c_1&c_2&\c&c_{k-1}
\end{pmatrix}$$
($A^T$は行列の転置)と定めると$M$の固有多項式$\det(xI-M)$は$a_n$の固有多項式$f(x)$に等しく、また
$$\ab_n=\begin{pmatrix}a_n&a_{n+1}&\cdots&a_{n+k-1}\end{pmatrix}^T$$
とおくと$a_n$の漸化式は$\ab_{n+1}=M^T\ab_n$と表現できる。
いま
$$\ab_n=M^T\ab_{n-1}=(M^T)^n\ab_0$$
より$a_n$は$\ab_0$と$M^n$の値によって定まるので、以下$M^n$の計算に焦点を当てていきます。
ここで$\aa$の固有方程式が重解を持たない、つまりその解$x=\a_1,\a_2,\ldots,\a_k$がそれぞれ異なるものとすると以下が成り立ちます。
$$P=
\begin{pmatrix}
1&\a_1&\c&\a_1^{k-1}\\
1&\a_2&\c&\a_2^{k-1}\\
\v&\v&\d&\v\\
1&\a_k&\c&\a_k^{k-1}\\
\end{pmatrix},\quad
D=\begin{pmatrix}\a_1\\&\a_2\\&&\d\\&&&\a_k\end{pmatrix}$$
とおくと
$$P^{-1}=\left(\frac{\sum^{i-1}_{l=0}c_l\a_j^{l-i}}{f'(\a_j)}\right)_{1\leq i,j\leq k},\quad M^n=P^{-1}D^nP$$
が成り立つ。
そして$P^{-1}D^nP$を実際に計算することで以下の公式が得られます。
$$A_n=\sum^k_{m=1}\frac{\a_m^n}{f'(\a_m)}$$
とおくと
$$(M^n)_{i,j}=\sum^{i-1}_{l=0}c_{i-1-l}A_{n+j-l-2}$$
が成り立つ。
この系として apu_yokai 氏が こちらの記事 で提起していた予想が解決されます。
$g(x)$を多項式とし
$$A'_n=\sum^k_{m=1}\frac{\a_m^n}{f'(\a_m)}g(\a_m)$$
とおくと
$$(g(M))_{i,j}=\sum^{i-1}_{l=0}c_{i-1-l}A'_{j-l-2}$$
が成り立つ。
特に$c_0=c_1=\cdots=c_{k-1}=1$において
$$(M^n)_{i,j}=\sum^{i-1}_{l=0}F^{[k]}_{n+j-l-2},\;(f'(M))_{i,j}=\sum^{i-1}_{l=0}L^{[k]}_{j-l-2}$$
が成り立つ。
$g(x)=x^n$とおくと$A'_{j-l-2}=A_{n+j-l-2}$なので定理3はこの系に埋め込まれることになります。
また$a_n=((M^n)^T\ab_0)_1=\sum^k_{i=1}(M^n)_{i,1}a_{i-1}$および
$$(M^n)_{i,1}=\sum^{i-1}_{l=0}c_{i-1-l}A_{n-l-1}
=A_{n+k-i}-\sum^{k-i-1}_{l=0}c_{i+l}A_{n+l}$$
から以下のように$\aa$の一般項が求まります。
\begin{align}
a_n
&=\sum^{k-1}_{l=0}(\sum^{k-1}_{i=l}c_{i-l}a_i)A_{n-l-1}\\
&=\sum^{k-1}_{l=0}(a_{k-1-l}-\sum^{k-1}_{i=l+1}c_ia_{i-l-1})A_{n+l}
\end{align}
が成り立つ。
ちなみに定理3を$n=0$において適応し$(k,j)$成分を比較すると漸化式から
$$(M^0)_{k,j}=(I)_{k,j}=\sum^{k-1}_{l=0}c_lA_{j-1-k+l}=A_{j-1}$$
となるので$A_n$の一般項を漸化式で表すと以下のようになります。
$\{A_n\}$は漸化式
$$\left\{\begin{array}{ll}
A_n=0&(0\leq n\leq k-2)\\
A_n=1&(n=k-1)\\
\dis A_{n+k}=\sum^{k-1}_{j=0}c_jA_{n+j}&(n\geq0)
\end{array}\right.$$
によって定まる。
$A_n$は定理3で示した定義では扱いづらいですがこうやって漸化式で表すと色々と便利です。特に$c_j$の取り方によらず初項が同じというのは扱いやすい事実ですね。
定理2の証明
$M^T$の固有値$\a=\a_1,\a_2,\ldots,\a_k$に対する固有ベクトルを考える。
$$\a I-M^T=
\begin{pmatrix}
\a&-1&\c&0\\
0&\a&\c&0\\
\v&\v&\d&\v\\
-c_0&-c_1&\c&\a-c_{k-1}
\end{pmatrix}$$
の第$k$行に第$1$行の$\frac{c_0}\a$倍を足し、また第$2$行の$\frac{c_1\a+c_0}{\a^2}$倍を足し、第$3$行の$\frac{c_2\a^2+c_1\a+c_0}{\a^3}$倍を足し、...と掃き出していくことで最終的に$(k,k)$成分は
$$\a-\frac{c_{k-1}\a^{k-1}+c_{k-2}\a^{k-2}+\cdots+c_0}{\a^{k-1}}=\a-\frac{\a^k}{\a^{k-1}}=0$$
となる。そしてその状態からさらに掃き出していくことで$\a I-M^T$の簡約化
$$\begin{pmatrix}
1&0&\c&0&-\frac{1}{\a^{k-1}}\\
0&1&\c&0&-\frac{1}{\a^{k-2}}\\
\v&\v&\d&\v&\v\\
0&0&\c&1&-\frac{1}{\a}\\
0&0&\c&0&0
\end{pmatrix}$$
を得る。すなわち$\a$に対する固有ベクトルは
$$\begin{pmatrix}1&\a&\cdots&\a^{k-2}&\a^{k-1}\end{pmatrix}^T$$
と取れるので
$$P=
\begin{pmatrix}
1&\a_1&\c&\a_1^{k-1}\\
1&\a_2&\c&\a_2^{k-1}\\
\v&\v&\d&\v\\
1&\a_k&\c&\a_k^{k-1}\\
\end{pmatrix},\quad
D=\begin{pmatrix}\a_1\\&\a_2\\&&\d\\&&&\a_k\end{pmatrix}$$
とおくと対角化の理論から
$$(P^T)^{-1}(M^T)^nP^T=D^n$$
ひいては
$$PM^nP^{-1}=D^n$$
が成り立つことになる(ちなみに$\aa$の固有方程式が重解を持つときは$P$が正則にならないので広義固有空間$\mathrm{Ker}(\a I-M)^m$の基底を考える必要が出てくる)。
そして$P$は$\a_1,\a_2,\ldots,\a_k$についてのヴァンデルモンド行列なので
この記事
で解説したようにその逆行列は
$$P^{-1}=\left(\frac{\sum^{i-1}_{l=0}c_l\a_j^{l-i}}{f'(\a_j)}\right)_{1\leq i,j\leq k}$$
と求まることがわかる。
定理3の証明
いま
$$D^nP=(\a_i^{n+j-1})_{1\leq i,j\leq k}$$
に注意すると
\begin{eqnarray}
(M^n)_{i,j}
&=&\sum^k_{m=1}(P^{-1})_{i,m}(D^nP)_{m,j}
\\&=&\sum^k_{m=1}\frac{\sum^{i-1}_{l=0}c_l\a_m^{l-i}}{f'(\a_m)}\cdot\a_m^{n+j-1}
\\&=&\sum^{i-1}_{l=0}c_lA_{n+l-i+j-1}
\\&=&\sum^{i-1}_{l=0}c_{i-1-l}A_{n+j-l-2}
\end{eqnarray}
と計算できる。
定理3系の証明
いま$g(x)=\sum^{d}_{t=0}C_tx^t$とおくと
\begin{eqnarray}
(g(M))_{i,j}&=&\sum^{d}_{t=0}C_t(M^t)_{i,j}
\\&=&\sum^{d}_{t=0}C_t\sum^k_{m=1}\frac{\sum^{i-1}_{l=0}c_{l}\a_m^{l-i}}{f'(\a_m)}\a_m^{t+j-1}
\\&=&\sum^{i-1}_{l=0}c_l\sum^k_{m=1}\frac{\a_m^{j+l-i-1}}{f'(\a_m)}\sum^{d}_{t=0}C_t\a^t_m
\\&=&\sum^{i-1}_{l=0}c_{i-1-l}\sum^k_{m=1}\frac{\a_m^{j-l-2}}{f'(\a_m)}g(\a_m)
\\&=&\sum^{i-1}_{l=0}c_{i-1-l}A'_{j-l-2}
\end{eqnarray}
と計算できる。
参考文献
この記事を高評価した人
この記事に送られたバッジ
投稿者
