2
大学数学基礎解説
文献あり

重解のない固有方程式を持つ数列の行列表現とその成分

86
1
$$\newcommand{A}[0]{\boldsymbol{A}} \newcommand{a}[0]{\alpha} \newcommand{aa}[0]{\{a_n\}} \newcommand{ab}[0]{\boldsymbol{a}} \newcommand{c}[0]{\cdots} \newcommand{d}[0]{\ddots} \newcommand{dis}[0]{\displaystyle} \newcommand{v}[0]{\vdots} $$

はじめに

この記事では重解を持たない固有方程式$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)=0$$\aa$の固有方程式$f(x)=0$に等しく、
また$\ab_n=\begin{pmatrix}a_n&a_{n+1}&\cdots&a_{n+k-1}\end{pmatrix}^T$とおくと$\aa$の漸化式は$\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}, D=\begin{pmatrix}\a_1\\&\a_2\\&&\d\\&&&\a_k\end{pmatrix}$

とおくと$\dis 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}$であり$M^n=P^{-1}D^nP$が成り立つ。

そして$P^{-1}D^nP$を実際に計算することで以下の公式を得ます。

$\dis A_n=\sum^k_{m=1}\frac{\a_m^n}{f'(\a_m)}$とおくと$\dis (M^n)_{i,j}=\sum^{i-1}_{l=0}c_{i-1-l}A_{n+j-l-2}$が成り立つ。

この系として apu_yokai 氏が こちらの記事 で提起していた予想が解決されます。

$g(x)$を多項式とし$\dis A'_n=\sum^k_{m=1}\frac{\a_m^n}{f'(\a_m)}g(\a_m)$とおくと$\dis(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$において$\dis (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}$および
$\dis(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$の一般項が求まります。

$\dis 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}$が成り立つ。

ちなみに定理3を$n=0$において適応し$(k,j)$成分を比較すると漸化式から
$\dis (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$の取り方によらず初項が同じというのは扱いやすい事実ですね。

$M$の対角化

$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)$成分は
$\dis\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}\frac{1}{\a^{k-1}}\\\frac{1}{\a^{k-2}}\\\v\\\frac{1}{\a}\\1\end{pmatrix}$ひいては$\begin{pmatrix}1\\\a\\\v\\\a^{k-2}\\\a^{k-1}\end{pmatrix}$
と取れるので
$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}, 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$についてのヴァンデルモンド行列なので この記事 で解説したように逆行列が求められ、定理2を得る。

$P^{-1}D^nP$の計算

いま
$\dis D^nP=(\a_i^{n+j-1})_{1\leq i,j\leq 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}$
であるので
\begin{eqnarray} (P^{-1}D^nP)_{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)}\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(M)$の計算

いま$\dis 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}
と定理3系を得る。

参考文献

投稿日:2021131

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

子葉
子葉
788
149299
主に複素解析、代数学、数論を学んでおります。 私の経験上、その証明が簡単に探しても見つからない、英語の文献を漁らないと載ってない、なんて定理の解説を主にやっていきます。 同じ経験をしている人の助けになれば。最近は自分用のノートになっている節があります。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中