リアプノフスペクトルの定義
定義
空間 $X = \mathbb{R}^{n}$ と 函数 $f : X \to X$ に対して次のような ベクトル場 が 微分方程式 として与えられているとする。
$$
\dot{x} = f(x)
$$
簡単な定義
フロー $F_{T} (v)$ の時間-$1$ マップ に対する 多次元マップのリャプノフ数とリャプノフ指数 をそれぞれ $F_{T} (v)$ の リャプノフ数Lyapunov number、リャプノフ指数Lyapunov exponent と定義する1。
難しい定義
変分方程式: $f$ の ヤコビアン行列 $J$ に関して次を 変分方程式variational equation とする。
$$ \dot{Y} = J Y $$ ここで 行列函数 $Y = Y(t) \in \mathbb{R}^{n \times n}$ の初期条件は 単位行列 $Y(0) = I$ とする。
…
幾何学的に、$Y$ は元のシステムの$x(0)$ から少し動かした $x(t)$ に変わる間に、その 接線ベクトル 自体がどのように作用するかを示していると言える。
$$
\lambda_{k} := \lim_{t \to \infty} \log \left[ \left( {\frac{ \left\| Y(t, v) v \right\| }{ \left\| v \right\| }} \right)^{1/t} \right]
$$
このように定義された $\left\{ \lambda_{1} , \cdots , \lambda_{n} \right\}$ を リャプノフスペクトルLyapunov spectrum とするか、
$$
\Lambda_{v} := \lim_{t \to \infty} \left[ Y(t)^{\ast} Y(t) \right]^{1/2t}
$$
このように定義された 行列 $\Lambda_{v}$ の 固有値 $\mu_{1} , \cdots , \mu_{n}$ に ログ を取った $\lambda_{k} := \log \mu_{k}$ をリャプノフスペクトルとする2。
- $Y^{\ast}$ は行列 $Y$ の 共役転置行列 である。
説明
実際、両方の定義ともそれほど簡単ではない。実際に連続的なシステムでリャプノフスペクトルというものを理解し扱うことはかなり難しい。
一次元マップのリャプノフ数 と同様に、リャプノフスペクトルのそのモチーフは $x_{0}$ と $x_{0} + \delta_{0}$ の間の小さな差である $\delta_{0}$ と $t$ 分だけ時間が経った後の差である $\delta_{t}$ を次のようにある $\lambda$ に関連した式で表現しようということから始まる。
$$
\left| \delta_{t} \right| \approx \left| \delta_{0} \right| e^{t \lambda}
$$
$t = N$ 時点で $T_{N}$ という オペレーター が $T_{N} : v_{N} \mapsto v_{N+1}$ のようにマッピングの役割を果たすとすると、$T_{N}$ が空間を拡張または縮小させる割合の 幾何平均 は次のようになる。
$$
\left( {\frac{ \left\| T_{1} v \right\| }{ \left\| v \right\| }} \cdot {\frac{ \left\| T_{2} v \right\| }{ \left\| T_{1} v \right\| }} \cdot \cdots \cdot \cdot {\frac{ \left\| T_{N} v \right\| }{ \left\| T_{N-1} v \right\| }} \right)^{1/N} = \left( {\frac{ \left\| T_{N} v \right\| }{ \left\| v \right\| }} \right)^{1/N}
$$
連続的なシステムでは変分章程式の $Y$ が $T_{n}$ の役割を果たすが、これに従い 正規直交集合 $\left\{ v_{1} , \cdots , v_{n} \right\}$ の $v_{k}$ に対する $k$ 番目のリャプノフ指数 $\lambda_{k}$ は次のように定義される3。
$$ \begin{align*} \lambda_{k} :=& \lim_{t \to \infty} \log \left[ \left( {\frac{ \left\| Y(t, v) v \right\| }{ \left\| v \right\| }} \right)^{1/t} \right] \\ =& \lim_{t \to \infty} {\frac{ 1 }{ t }} \log \left( \left\| Y(t, v) v \right\| - \left\| v \right\| \right) \\ =& \lim_{t \to \infty} {\frac{ 1 }{ t }} \log \left\| Y(t, v) v \right\| \end{align*} $$
一方で、$Y(t,u)$ の 特異値分解 $Y = U \Sigma V^{\ast}$ に従い $k$ 番目の特異値 $\sigma_{k} (t)$ は $U, V$ の $k$ 番目のカラムベクトル $u_{k}, v_{k}$ に対して次のように表すことができる4。
$$
Y V = U \Sigma \implies Y v_{k} = \sigma_{k} (t) u_{k}
$$
これは $Y$ の特異値、すなわち 固有値 に相当する概念がリャプノフスペクトルと関連している手がかりを提供する。実際に $Y$ の左側に $Y^{\ast} = V \Sigma^{\ast} U^{\ast}$ を掛けた
$$
Y^{\ast} Y = V \Sigma^{2} V^{\ast}
$$
の固有値は $\sigma_{k}^{2} (t)$ であることは難しく推測できないし、
$$
\Lambda_{v} := \lim_{t \to \infty} \left[ Y(t)^{\ast} Y(t) \right]^{1/2t}
$$
と定義された 行列 $\Lambda_{v}$ の 固有値 $\mu_{k}$ に ログ を取った $\log \mu_{k}$ はリャプノフスペクトルとなる。前述の $Y v_{k} = \sigma_{k} (t) u_{k}$ と関連付けて考えてみると、
$$
\begin{align*}
\log \mu_{k} =& \log \lim_{t \to \infty} \left[ \sigma_{k}^{2} (t) \right]^{1/2t}
\\ =& \lim_{t \to \infty} \log \left[ \sigma_{k} (t) \right]^{1/t}
\\ =& \lim_{t \to \infty} {\frac{ 1 }{ t }} \log \sigma_{k} (t)
\\ =& \lim_{t \to \infty} {\frac{ 1 }{ t }} \log \left\| \sigma_{k} (t) u \right\|
\\ =& \lim_{t \to \infty} {\frac{ 1 }{ t }} \log \left\| Y(t, v) v \right\|
\\ =& \lambda_{k}
\end{align*}
$$
のように直感として両定義が同じであることを受け入れることができる。
Yorke. (1996). CHAOS: An Introduction to Dynamical Systems: p381. ↩︎
http://crossgroup.caltech.edu/Chaos_Course/Lesson7/Lyapunov.pdf ↩︎
Karlheinz Geist, Ulrich Parlitz, Werner Lauterborn, Comparison of Different Methods for Computing Lyapunov Exponents, Progress of Theoretical Physics, Volume 83, Issue 5, May 1990, Pages 875–893, https://doi.org/10.1143/PTP.83.875 ↩︎