行列差分方程式

力学系の理論における行列差分方程式（ぎょうれつさぶんほうていしき、テンプレート:Lang-en-short）は、各時点においてベクトルとして与えられている変数に関する系で、各時点の値がそれ以前の時点における変数の値と行列を用いた関係性で結ばれているような差分方程式を言う。

そのような方程式の階数 (order) とは、変数ベクトルの値を指定するために必要となる任意の二時点の間隔の最大値を言う。例えば、テンプレート:Mvar を テンプレート:Mvar 成分列ベクトル、テンプレート:Mvar を テンプレート:Mvar-次正方行列として、$x_t=A{x}_{t-1}+B{x}_{t-2}$ は二階方程式の例ということになる。特にこれは右辺に定ベクトル項を持たないから斉次方程式の例でもある。添字をずらしたり、別の文字を用いたりして $x_{t+2}=A{x}_{t+1}+B{x}_t$ や $x_n=A{x}_{n-1}+B{x}_{n-2}$ のように書いても方程式としては同じものである。

最もよく遭遇する行列差分方程式は、一階のものであろう。

一階方程式

特に非斉次の一階方程式を考え、その定常状態について確認しよう。すなわち、定ベクトル テンプレート:Mvar を持つ方程式 $$x_t=A{x}_{t-1}+b$$ に関して、それが定める系の定常状態 テンプレート:Mvar とは、それ以降に新たな値を導くことがない状態に到達したベクトル テンプレート:Mvar の値を言う。定義により テンプレート:Mvar は当該の方程式において $x_t=x_{t-1}=x^{*}$ と置いたものを満足するものでなくてはならず、それを テンプレート:Mvar に関して解いて $x^{*}=(I-A{)}^{-1}b$ を得るが、ここで テンプレート:Math は可逆であるものと仮定している（ただし、テンプレート:Mvar は テンプレート:Mvar-次単位行列）。

定常状態が存在するとき、上記の非斉次方程式は、$$[x_t-x^{*}]=A[x_{t-1}-x^{*}]$$ なる形の斉次方程式に書き直すことができる。この斉次方程式が安定である—すなわち、テンプレート:Mvar が定常状態 テンプレート:Mvar に漸近収束する—ための必要十分条件は、（実または複素）変換行列 テンプレート:Mvar の固有値が テンプレート:Math より小さい絶対値を持つことである。

方程式を斉次形で $y_t=A{y}_{t-1}$ と書くことにすれば、よく知られた仕方により、初期条件すなわち テンプレート:Mvar の初期値 テンプレート:Mvar から反復適用と置換から直ちに $$\begin{array}{cc}{y}_1& =A{y}_0,\\ y_2& =A{y}_1=AA{y}_0=A^2{y}_0,\\ y_3& =A{y}_2=A{A}^2{y}_0=A^3{y}_0\end{array}$$ のように各項の値は順次求まり、そして テンプレート:Mvar に関する数学的帰納法によって一般に $y_t=A^t{y}_0$ が成立することが証明される。さらに テンプレート:Mvar が対角化可能ならば、テンプレート:Mvar はその固有値と固有ベクトルを用いて書けるから、その解は $$y_t=P{D}^t{P}^{-1}{y}_0$$ の形に与えられる。ここに テンプレート:Mvar はその列ベクトルが テンプレート:Mvar の固有ベクトルの全体からなる テンプレート:Mvar-次正方行列（固有値が全て異なると仮定した場合には必ず可逆）であり、テンプレート:Mvar は テンプレート:Mvar の固有値が主対角線上に並ぶ対角行列である。

この解の様子が上記の安定性について得ることを動機づける。すなわち テンプレート:Mvar が十分な時間経過後に零行列となるための必要十分条件は テンプレート:Mvar の任意の固有値が絶対値に関して テンプレート:Math より小さいことである。

一階行列系から高階スカラー系へ

テンプレート:Mvar-次元ベクトルの成す系 $y_t=A{y}_{t-1}$ から一つの状態変数、具体的に テンプレート:Math に関する力学系に展開することができる。上で見た方程式のベクトル解 テンプレート:Mvar を テンプレート:Mvar の テンプレート:Mvar 個の固有値に関して展開して、それを テンプレート:Math と書けば、それは テンプレート:Math の発展を記述するものであって、それ自身もまた同じ固有値を含む解を持つ。この記述は、直観的には テンプレート:Math の発展の方程式 $$y_{1,t}=a_1{y}_{1,t-1}+a_2{y}_{1,t-2}+\dots +a_n{y}_{1,t-n}$$ を動機づける。ここに各パラメータ テンプレート:Mvar は テンプレート:Mvar の固有方程式 $${\lambda }^n-a_1{\lambda }^{n-1}-a_2{\lambda }^{n-2}-\dots -a_n{\lambda }^0=0$$ の係数である。

すなわち、テンプレート:Mvar-次元一階線型系の変数ベクトルの各成分スカラーは、各々個別の一スカラー変数 テンプレート:Mvar-階差分方程式に従って発展する。そのような一変数方程式の安定性はもとの行列差分方程式の安定性と同じである。

高階方程式の解と安定性

高階—すなわち最大の時間間隔が一周期よりも長い—行列方程式も、区分行列に関する一階の形の方程式に変換すれば、解くことができて、それらの安定性を調べることができる。

例えば、テンプレート:Mvar-成分変数ベクトル テンプレート:Mvar と テンプレート:Mvar-次係数行列 テンプレート:Mvar を持つ二階方程式 $x_t=A{x}_{t-1}+B{x}_{t-2}$ は テンプレート:Mvar-次単位行列 テンプレート:Mvar と テンプレート:Mvar-次零行列 テンプレート:Mvar を用いて $$\left(\begin{array}{c}{x}_t\\ x_{t-1}\end{array}\right)=\left(\begin{array}{cc}A& B\\ I& O\end{array}\right)\left(\begin{array}{c}{x}_{t-1}\\ x_{t-2}\end{array}\right)$$ の形にまとめることができる。これを、現在時点とその一つ前の時点における変数をひとまとめにした テンプレート:Math 成分変数ベクトル テンプレート:Mvar と テンプレート:Math-次正方区分行列 テンプレート:Mvar に関する一階方程式と見れば、既にみたように解 $$z_t=L^t{z}_0$$ が得られるというわけである。ゆえに以前と同じように、このまとめられた方程式、したがってもとの二階方程式が安定となるための必要十分条件として、行列 テンプレート:Mvar の任意の固有値が絶対値に関して テンプレート:Math より小さいこと、を挙げることができる。

関連項目

• テンプレート:Ill2
• 差分方程式
• 線型差分方程式
• 力学系
• テンプレート:Ill2