リアプノフ指数

リアプノフ指数（リアプノフしすう、英: Lyapunov exponent）とは、力学系においてごく接近した軌道が離れていく度合いを表す量である。リャプノフ指数とも表記される。ロシア人科学者 Алекса́ндр Ляпуно́в（アレクサンドル・リプノーフ、Aleksandr Lyapunov）にその名をちなむ。 系の相空間上の2つの軌道について考える。2つの軌道上の時刻 t における点の距離をベクトル δ(t) として、初期状態 t = 0 には、これらの軌道は距離 δ(0) だけ離れているとする。δ(t) を近似的に次のように表す。 ‖ δ ( t ) ‖ ≈ ‖ δ ( 0 ) ‖ e λ t {\displaystyle \|{\boldsymbol {\delta }}(t)\|\approx \|{\boldsymbol {\delta }}(0)\|e^{\lambda t}} ここで ‖ ⋅ ‖ {\displaystyle \|\cdot \|} はユークリッドノルムを意味する。上式で λ > 0 の場合は軌道は離れていき、 λ < 0 の場合は軌道は近づいていく。よって、軌道が離れていく度合いは λ の値により決定される。この λ がリアプノフ指数である。軌道がカオス的であるとき、上式のように軌道は指数関数的に離れていく。すなわち、リアプノフ指数が正であることが軌道がカオス的であることの1つの定義とされる。 より詳細には、系の状態変数が k 個（k > 1）の場合、すなわち相空間が k 次元である場合は各次元ごとに固有のリアプノフ指数を持つ。これらのリアプノフ指数の組をリアプノフスペクトラムと呼び、そのうちの最大のリアプノフ指数を最大リアプノフ指数と呼ぶ。各々のリアプノフ指数を見れば正であったり負であったりするが、最大リアプノフ指数が正であれば、その系はカオスの特徴の1つである初期値鋭敏性を持つといえる。 まず、単純な1次元離散力学系の場合のリアプノフ指数について説明する。 x ∈ R {\displaystyle x\in \mathbb {R} } を系の状態変数、 n ∈ N {\displaystyle n\in \mathbb {N} } を離散時間としたとき（ここでは N {\displaystyle \mathbb {N} } は0を含む）、 写像 xn+1 = f(xn) のリアプノフ指数 λ は次のように定義される。 λ = lim n → ∞ 1 n ∑ i = 0 n − 1 ln ⁡ | f ′ ( x i ) | {\displaystyle \lambda =\lim _{n\to \infty }{\frac {1}{n}}\sum _{i=0}^{n-1}\ln |f'(x_{i})|} ここで、ln は自然対数を意味する。上式は次のように導入される。 初期位置を x0 とする。さらに、x0 からの微小量 λ0 ずれた点 x0 + λ0 を考える。リアプノフ指数では x0 から出発する軌道と x0 + λ0 から出発する軌道がどれだけ離れていくかを定義したい。ずれは時間発展とともに変化していくと考えられるので、時刻 n におけるずれを λn で表す。n = 1 でのずれは δ 1 = f ( x 0 + δ 0 ) − f ( x 0 ) {\displaystyle \delta _{1}=f(x_{0}+\delta _{0})-f(x_{0})} となり、n = n でのずれも同様に、 δ n = f n ( x 0 + δ 0 ) − f n ( x 0 ) {\displaystyle \delta _{n}=f^{n}(x_{0}+\delta _{0})-f^{n}(x_{0})} と得られる。ここで、fn(x)は f(x)の n 回反復写像を表す。 本記事の冒頭で定義したように、λn が n に指数関数的に比例するとして、 | δ n | = | δ 0 | e λ n {\displaystyle |\delta _{n}|=|\delta _{0}|e^{\lambda n}} と表す。両辺の自然対数をとると、 λ = 1 n ln ⁡ | δ n δ 0 | {\displaystyle \lambda ={\frac {1}{n}}\ln \left|{\frac {\delta _{n}}{\delta _{0}}}\right\vert } が得られる。ただし、初期のずれ量 λ0 は微小量としたが、実際にはリアプノフ指数は初期のずれ量を無限小とした λ0 → 0 の極限値で定義される。よって、上式は λ = 1 n ln ⁡ | lim δ 0 → 0 δ n δ 0 | {\displaystyle \lambda ={\frac {1}{n}}\ln \left|\lim _{\delta _{0}\to 0}{\frac {\delta _{n}}{\delta _{0}}}\right\vert } となる。上式の絶対値の中身に注目すると、 lim δ 0 → 0 δ n δ 0 = lim δ 0 → 0 f n ( x 0 + δ 0 ) − f n ( x 0 ) δ 0 = ( f n ) ′ ( x 0 ) = ∏ i = 0 n − 1 f ′ ( x i ) {\displaystyle \lim _{\delta _{0}\to 0}{\frac {\delta _{n}}{\delta _{0}}}=\lim _{\delta _{0}\to 0}{\frac {f^{n}(x_{0}+\delta _{0})-f^{n}(x_{0})}{\delta _{0}}}=(f^{n})'(x_{0})=\prod _{i=0}^{n-1}f'(x_{i})} とできる。ここで (f n)'(x) は、fn(x) の微分を意味する。 ∏ {\displaystyle \prod } は総乗を意味し、最右辺は合成関数の微分の連鎖律により得ることができる。よって、 λ = 1 n ln ⁡ | ∏ i = 0 n − 1 f ′ ( x i ) | = 1 n ∑ i = 0 n − 1 ln ⁡ | f ′ ( x i ) | {\displaystyle \lambda ={\frac {1}{n}}\ln \left|\prod _{i=0}^{n-1}f'(x_{i})\right\vert ={\frac {1}{n}}\sum _{i=0}^{n-1}\ln \left|f'(x_{i})\right\vert } となる。さらに上式において n → ∞ とした極限値が存在するとき、その極限値を初期値 x0 から出発する軌道のリアプノフ指数と呼ぶ。 λ = lim n → ∞ 1 n ∑ i = 0 n − 1 ln ⁡ | f ′ ( x i ) | {\displaystyle \lambda =\lim _{n\to \infty }{\frac {1}{n}}\sum _{i=0}^{n-1}\ln |f'(x_{i})|} 1968年に発表されたValery Oseledecの多重エルゴード定理により、n → ∞ の極限値が存在すること、ほとんどすべての初期値 x0 で λ は同じ値に収束することが証明されている。対象とする力学系のアトラクターの吸引域内の初期値であれば、全ての初期値で同じ λ の値に収束する。 力学系が k 次元の相空間を持つ高次元力学系の場合、各方向に別々のリアプノフ指数が存在する。すなわち高次元力学系であれば、軌道のずれは、ある方向には離れていくが、別の方向では縮まっていく状況がありえる。よって k 個のリアプノフ指数を得ることができる。このような k 個のリアプノフ指数の組を、リアプノフスペクトラム（英: Lyapunov spectrum）と呼ぶ。 λ i = { λ 1 , λ 2 , ⋯ , λ k } {\displaystyle \lambda _{i}=\{\lambda _{1},\lambda _{2},\cdots ,\lambda _{k}\}} リアプノフスペクトラムでは一般に λ1 から値が大きい順に並べる。最大値である λ1 を、特に最大リアプノフ指数（英: maximal Lyapunov exponent, maximum Lyapunov exponent）と呼ぶ。記事冒頭で述べたように、相空間上の2つの軌道上の時刻 t における点の間の距離、すなわちずれを δ(t) とする。リアプノフスペクトラム λi は以下のように定義される。 λ i = lim t → ∞ 1 t ln ⁡ α i ( i = 1 , 2 , ⋯ , k ) {\displaystyle \lambda _{i}=\lim _{t\to \infty }{\frac {1}{t}}\ln \alpha _{i}\quad (i=1,2,\cdots ,k)} 一般に λi は初期値 x(0) に依存する。しかし1次元離散力学系の場合と同様に、ほとんどすべての初期位置 x0 から同一の λi を得ることができる。λi の定義式にある αi は、次式で定義される k × k 正定値行列 Λ の固有値である。 Λ = ( M T M ) 1 2 t {\displaystyle {\boldsymbol {\Lambda }}=({\boldsymbol {M}}^{T}{\boldsymbol {M}})^{\frac {1}{2t}}} さらに M は、δ(t) の解を次の形式で表したときの δ(0) に対する乗数として得られる。 δ ( t ) = M δ ( 0 ) {\displaystyle {\boldsymbol {\delta }}(t)={\boldsymbol {M}}{\boldsymbol {\delta }}(0)} 系が連続力学系の場合、 k 個の状態変数 {x1, x2, ..., xk}、常微分方程式 {f1, f2, ..., fk} から成る常微分方程式系 d x d t = f ( t , x ) {\displaystyle {\frac {d{\boldsymbol {x}}}{dt}}={\boldsymbol {f}}(t,{\boldsymbol {x}})} が与えられる。f が線形近似可能な場合、f のヤコビ行列を用いて、 d δ ( t ) d t = J δ ( t ) {\displaystyle {\frac {d{\boldsymbol {\delta }}(t)}{dt}}={\boldsymbol {J}}{\boldsymbol {\delta }}(t)} と表すことができる。また、系が離散力学系の場合、 k 個の状態変数、常差分方程式から成る差分方程式系 x(t + 1) = f(x(t))が与えられる。同じく、差分方程式系 f が線形近似可能な場合、f のヤコビ行列を用いて、 δ ( t + 1 ) = J δ ( t ) {\displaystyle {\boldsymbol {\delta }}(t+1)={\boldsymbol {J}}{\boldsymbol {\delta }}(t)} と表すことができる。ここに、J は以下に示すようなヤコビ行列による線形写像で、軌道 x(t) に依存し、すなわち初期値 x(0)、時間 t に依存して変化する。 J = ( ∂ f 1 ∂ x 1 ⋯ ∂ f 1 ∂ x k ⋮ ⋱ ⋮ ∂ f k ∂ x 1 ⋯ ∂ f k ∂ x k ) {\displaystyle {\boldsymbol {J}}={\begin{pmatrix}{\cfrac {\partial f_{1}}{\partial x_{1}}}&\cdots &{\cfrac {\partial f_{1}}{\partial x_{k}}}\\\vdots &\ddots &\vdots \\{\cfrac {\partial f_{k}}{\partial x_{1}}}&\cdots &{\cfrac {\partial f_{k}}{\partial x_{k}}}\end{pmatrix}}} 常微分方程式系の場合は、 d δ d t = J ( t ) δ {\displaystyle {\frac {d{\boldsymbol {\delta }}}{dt}}={\boldsymbol {J}}(t){\boldsymbol {\delta }}} を解いて δ ( t ) = M δ ( 0 ) {\displaystyle {\boldsymbol {\delta }}(t)={\boldsymbol {M}}{\boldsymbol {\delta }}(0)} を得ることで、上記の定義で出てきた正方行列 M を得ることができる。差分方程式系の場合は、J を t 回繰り返し適用することで次のような δ(t) と δ(0) の関係式で書き表すことできるので、M は Jn の 0 から t − 1 までの総乗として得ることができる。 δ ( t ) = J ( t − 1 ) J ( t − 2 ) ⋯ J ( 0 ) δ ( 0 ) = M δ ( 0 ) {\displaystyle {\boldsymbol {\delta }}(t)={\boldsymbol {J}}(t-1){\boldsymbol {J}}(t-2)\cdots {\boldsymbol {J}}(0){\boldsymbol {\delta }}(0)={\boldsymbol {M}}{\boldsymbol {\delta }}(0)} 対数の底に 2 を使用して計算した場合には、ビット/時間を単位として使用することがある。 これは、λ > 0 の場合、単位時間あたりλビットの情報が失われ、λ < 0 の場合λビットの情報が生成することに相当する。 保存系の場合、相空間の全エネルギーは保存される。従って全リアプノフ指数の総和はゼロになる。散逸系ではリアプノフ指数の総和は負になる。 力学系が何らかの流れである場合、1つのリアプノフ指数は常にゼロとなる。つまり、流れの方向の固有ベクトルに対応する固有値から得られるリアプノフ指数がゼロになる。 Pesin's theorem によれば、正のリアプノフ指数の総和はコルモゴロフ・シナイ・エントロピー（Kolmogorov-Sinai entropy）の近似値を与える。 最大リアプノフ指数の逆数を「リアプノフ時間; Lyapunov time」と呼ぶことがあり、e-folding time の特性を定義する。カオス的軌道ではリアプノフ時間は有限であり、正規の軌道では無限大となる。 軌道がカオス的振る舞いをみせるストレンジアトラクターは、フラクタル構造を持つことが多い。このようなアトラクターのフラクタル次元とリアプノフスペクトラムの間には関係が存在する。アトラクターのリアプノフスペクトラムが得られたとして、その各リアプノフ指数 λi がそれらの値の大きさで降順に並べられているとしたとき、次のようなフラクタル次元の一種 DL が定義される。 D L = j + ξ j | λ j + 1 | ( j < k ) {\displaystyle D_{L}=j+{\frac {\xi _{j}}{|\lambda _{j+1}|}}\quad (j<k)} ここで、 ξj は、 ξ j = ∑ i = 1 j λ i {\displaystyle \xi _{j}=\sum _{i=1}^{j}\lambda _{i}} であり、j は ξj が負にならない最大値、すなわち、リアプノフ指数を λ1 + λ2 + λ3 ...と順に足していったときに、その総和が負となる直前における、足し合わせたリアプノフ指数の個数である。ξj ≥ 0 を満たす j が存在しないときは、DL = 0、系の次元数 k と j が一致する場合は、DL = k とする。 このように定義されたフラクタル次元 DL は、リアプノフ次元 （英: Lyapunov dimension）と呼ばれる。リアプノフ次元は、James L. KaplanとJames A. Yorkeにより1979年に提案された。そのためカプラン・ヨーク次元（英: Kaplan–Yorke dimension）とも呼ばれ、記号 DKY とも記される。 上式で定義されるリアプノフ次元は、フラクタル次元の1つである容量次元の考え方をもとにして、次のように導入される。リアプノフ指数の総和は相空間内の k 次元の体積要素の体積変化率を与えるので、同様に部分和である ξj から、j 次元までの体積要素は拡大するが、j + 1 次元以上の体積要素は縮小することになる。そのため、アトラクタを収めるためには少なくとも j 次元は必要で、考えられるフラクタル次元の下限を与えている。一方で、j + 1 次元は考えられるフラクタル次元の上限といえる。k 次元相空間上の一辺が d の立方体は、時間発展により各辺 exp(λ1t), exp(λ2t), exp(λ3t),... の直方体に写像される。ここで、１辺の長さ ε が、 ϵ = d exp ⁡ ( λ j + 1 t ) {\displaystyle \epsilon =d\exp(\lambda _{j+1}t)} である j + 1 次元立方体の箱を考え、容量次元と同じように、ストレンジアトラクタの不変集合をこの箱を何個も当てながら全体を覆うことを考える。このとき、覆うのに必要な箱の数 N は、 N = exp ⁡ [ ( λ 1 − λ j + 1 ) t ] × exp ⁡ [ ( λ 2 − λ j + 1 ) t ] × ⋯ × exp ⁡ [ ( λ j − λ j + 1 ) t ] {\displaystyle N=\exp \left[(\lambda _{1}-\lambda _{j+1})t\right]\times \exp \left[(\lambda _{2}-\lambda _{j+1})t\right]\times \cdots \times \exp \left[(\lambda _{j}-\lambda _{j+1})t\right]} と推論できる。よって、容量次元と同様の定義から、 lim ϵ → 0 ln ⁡ N ln ⁡ ( 1 / ϵ ) ≈ lim d → 0 ( − j λ j + 1 + λ 1 + λ 2 + λ 3 + ⋯ + λ j ) t − λ j + 1 t − ln ⁡ d = j + λ 1 + λ 2 + λ 3 + ⋯ + λ j − λ j + 1 {\displaystyle \lim _{\epsilon \to 0}{\frac {\ln N}{\ln(1/\epsilon )}}\approx \lim _{d\to 0}{\frac {(-j\lambda _{j+1}+\lambda _{1}+\lambda _{2}+\lambda _{3}+\cdots +\lambda _{j})t}{-\lambda _{j+1}t-\ln d}}=j+{\frac {\lambda _{1}+\lambda _{2}+\lambda _{3}+\cdots +\lambda _{j}}{-\lambda _{j+1}}}} となり、j の定義より −λj + 1 = |λj + 1| なので、リアプノフ次元の定義 D L = j + λ 1 + λ 2 + λ 3 + ⋯ + λ j | λ j + 1 | {\displaystyle D_{L}=j+{\frac {\lambda _{1}+\lambda _{2}+\lambda _{3}+\cdots +\lambda _{j}}{|\lambda _{j+1}|}}} を得ることができる。 他のフラクタル次元と比較したリアプノフ次元の利点は、リアプノフスペクトラムさえ得ることができれば簡単に計算可能なことである。また、リアプノフ次元は系の情報量次元の上限を表している。 ^ リアプノフスペクトルと表記する場合もある。 ^ a b アリグッドほか 2012b, p. 1. ^ イアン・スチュアート、須田不二夫・三村和男（訳）、1998、『カオス的世界像 ―非定形の理論から複雑系の科学へ』第一版、白揚社 ISBN 4-8269-0085-6. p. 377 ^ Strogatz 2015, pp. 349–350. ^ a b 合原 2011, p. 158. ^ a b 下條 1992, p. 86. ^ a b 高安 2001, p. 82. ^ a b Strogatz 2015, p. 400. ^ アリグッドほか 2012b, pp. 3–4. ^ a b 下條 1992, p. 91. ^ a b c 高安 2001, p. 83. ^ a b 船越 2008, p. 146. ^ 金子、津田 1997, p. 47. ^ a b c Strogatz 2015, p. 401. ^ 合原 2011, p. 161. ^ a b 下條 1992, p. 87. ^ 船越 2008, p. 157. ^ アリグッドほか 2012a, pp. 117–118. ^ a b c d e 合原 2011, p. 165. ^ 船越 2008, p. 171. ^ 船越 2008, p. 181. ^ a b 合原 2011, p. 163. ^ a b 金子、津田 1997, p. 115. ^ 船越 2008, p. 168. ^ a b c d e 高安 2001, p. 84. ^ a b ベルジュほか 1992, p. 265. ^ 合原 2011, pp. 164–165. ^ 小室 2005, p. 17. ^ 合原 2011, p. 167. ^ 小室 2005, p. 22. ^ a b 合原 2011, p. 164. ^ 高安 2001, pp. 83–84. ^ ベルジュほか 1992, pp. 261–262. ^ アリグッドほか 2012b, p. 8. ^ Chaotic oscillators: theory and applications, Tomasz Kapitaniak,pp287 ^ 高安 2001, p. 93. ^ アリグッドほか 2012b, p. 11. ^ 合原 2011, p. 169. ^ アリグッドほか 2012b, p. 12. ^ a b c d Chlouverakis, Konstantinos E.; Sprott, J.C. (January 2005). “A comparison of correlation and Lyapunov dimensions”. Physica D: Nonlinear Phenomena (Elsevier) 200 (1–2): 156–164. doi:10.1016/j.physd.2004.10.006.  ^ 合原一幸ほか 著、合原一幸 編『カオス―カオス理論の基礎と応用』（初版）サイエンス社、1990年、97頁。ISBN 4-7819-0592-7。  ^ アリグッドほか 2012b, p. 15. ^ a b c アリグッドほか 2012b, p. 16. ^ アリグッドほか 2012b, pp. 11–12. ^ J. Kaplan and J. Yorke Chaotic behavior of multidimensional difference equations In Peitgen, H. O. & Walther, H. O., editors, Functional Differential Equations and Approximation of Fixed Points Springer, New York (1987) Cvitanovi? P., Artuso R., Mainieri R. , Tanner G. and Vattay G.Chaos: Classical and Quantum Niels Bohr Institute, Copenhagen 2005 Freddy Christiansen and Hans Henrik Rugh (1997). “Computing Lyapunov spectra with continuous Gram-Schmidt orthonormalization”. Nonlinearity 10: 1063–1072. オリジナルの2006年4月25日時点におけるアーカイブ。. https://web.archive.org/web/20060425194442/http://www.mpipks-dresden.mpg.de/eprint/freddy/9702017/9702017.ps.  Govindan Rangarajan, Salman Habib, and Robert D. Ryne (1998). “Lyapunov Exponents without Rescaling and Reorthogonalization”. Physical Review Letters 80: 3747–3750. https://arxiv.org/pdf/chao-dyn/9803017.  X. Zeng, R. Eykholt, and R. A. Pielke (1991). “Estimating the Lyapunov-exponent spectrum from short time series of low precision”. Physical Review Letters 66: 3229. http://link.aps.org/abstract/PRL/v66/p3229.  K.T.アリグッド・T.D.サウアー・J.A.ヨーク、シュプリンガー・ジャパン（編）、津田一郎（監訳）、星野高志・阿部巨仁・黒田拓・松本和宏（訳）、2012、『カオス　第1巻　力学系入門』、丸善出版 ISBN 978-4-621-06223-4 K.T.アリグッド・T.D.サウアー・J.A.ヨーク、シュプリンガー・ジャパン（編）、津田一郎（監訳）、星野高志・阿部巨仁・黒田拓・松本和宏（訳）、2012、『カオス　第2巻　力学系入門』、丸善出版 ISBN 978-4-621-06279-1 池口徹・山田泰司・小室元政、合原一幸（編）、2011、『カオス時系列解析の基礎と応用』初版第4刷、産業図書 ISBN 978-4-7828-1010-1 Steven H. Strogatz、田中久陽・中尾裕也・千葉逸人（訳）、2015、『ストロガッツ　非線形ダイナミクスとカオス―数学的基礎から物理・生物・化学・工学への応用まで』、丸善出版 ISBN 978-4-621-08580-6 下條隆嗣、1998、『カオス力学入門―古典力学からカオス力学へ』初版第4刷、近代科学社〈シミュレーション物理学6〉 ISBN 4-7649-2005-0 高安秀樹・本田勝也・佐野雅己・田崎睛明・村山和郎・伊藤敬祐、2001、『フラクタル科学』初版第11刷、朝倉書店 ISBN 4-254-10063-9 船越満明、2008、『カオス』初版、朝倉書店〈シリーズ 非線形科学入門3〉 ISBN 978-4-254-11613-7 金子邦彦・津田一郎、1997、『複雑系のカオス的シナリオ』初版第4刷、朝倉書店〈複雑系双書1〉 ISBN 4-254-10514-2 ピエール・ベルジュ、イヴェ・ポモウ、クリスチャン・ビダル、相澤洋二（訳）、1992、『カオスの中の秩序―乱流の理解に向けて』初版、産業図書 ISBN 4-7828-0068-1 小室元政、2005、『基礎からの力学系―分岐解析からカオス的遍歴へ』新版、サイエンス社 ISBN 4-7819-1118-8 エルゴード理論 リウヴィルの定理 (物理学) フロケ理論 Lyapunov exponent （英語） - スカラーペディア百科事典「リアプノフ指数」の項目。 Weisstein, Eric W. "Lyapunov Characteristic Exponent". mathworld.wolfram.com (英語).