リアプノフ指数

出典: フリー百科事典『ウィキペディア(Wikipedia)』
移動: 案内検索

リアプノフ指数(リアプノフしすう、: Lyapunov exponent)とは、力学系においてごく接近した軌道が離れていく度合いを表す量である。リャプノフ指数とも記される[1]

系の相空間上の2つの軌道について考える。2つの軌道上の時刻 t における点の距離をベクトル \boldsymbol{\delta}(t)として、初期状態t=0には、これらの軌道は距離 \boldsymbol{\delta}(0)だけ離れているとする。\boldsymbol{\delta}(t)を近似的に次のように表す[2][3][4]:

 \| \boldsymbol{\delta}(t) \| \approx \|  \boldsymbol{\delta}(0) \| e^{\lambda t}

ここで\| \cdot \|ノルムを意味する。上式で\lambda > 0の場合は軌道は離れていき、\lambda < 0の場合は軌道は近づいていく。よって、軌道が離れていく度合いは λ の値により決定される。この λ がリアプノフ指数である[3][4]。軌道がカオス的であるとき、上式のように軌道は指数関数的に離れていく[5][6]。すなわち、リアプノフ指数が正であることが軌道がカオス的であることの1つの定義とされる[7]

より詳細には、系の状態変数が N 個(N>1)の場合、すなわち相空間が N 次元である場合は各次元ごとに固有のリアプノフ指数を持つ[8]。これらのリアプノフ指数の組をリアプノフスペクトラム: Lyapunov spectrum)と呼び[9]、そのうちの最大のリアプノフ指数を最大リアプノフ指数: maximal Lyapunov exponent, maximum Lyapunov exponent)と呼ぶ[10]。各々のリアプノフ指数を見れば正であったり負であったりするが、最大リアプノフ指数が正であれば、その系はカオスの特徴の1つである初期値鋭敏性を持つといえる[11][10]

1次元離散力学系のリアプノフ指数[編集]

まず、単純な1次元離散力学系の場合のリアプノフ指数について説明する。x \in \Rを系の状態変数、n \in \Nを離散時間としたとき(ここでは\Nは0を含む)、 写像 x_{n+1}=f(x_n)のリアプノフ指数 λ は次のように定義される[12][5][13][14]

 \lambda = \lim_{n \to \infty}  \frac{1}{n} \sum_{i=0}^{n-1}  \ln | f'(x_i)|

ここで、\ln自然対数を意味する。上式は次のように導入される。

初期位置をx_0とする。さらに、x_0からの微小量 \delta_0ずれた点x_0+\delta_0を考える。リアプノフ指数ではx_0から出発する軌道とx_0+\delta_0から出発する軌道がどれだけ離れていくかを定義したい。ずれは時間発展とともに変化していくと考えられるので、時刻 n におけるずれを\delta_nで表す。n=1でのずれは\delta_1=f(x_0+\delta_0) - f(x_0)となり、n=nでのずれも同様に、\delta_n=f^n(x_0+\delta_0) - f^n(x_0)と得られる。ここで、f^n(x)f(x)n反復写像を表す。

本記事の冒頭で定義したように、\delta_nn に指数関数的に比例するとして、

 | \delta_n | = | \delta_0 | e^{\lambda n}

と表す[6]。両辺の自然対数をとると、

 \lambda =\frac{1}{n} \ln \left | \frac{\delta_n}{\delta_0} \right \vert

が得られる。ただし、初期のずれ量\delta_0は微小量としたが、実際にはリアプノフ指数は初期のずれ量を無限小とした\delta_0 \to 0の極限値で定義される[15]。よって、上式は

 \lambda =\frac{1}{n} \ln \left | \lim_{\delta_0 \to 0} \frac{\delta_n}{\delta_0} \right \vert

となる。上式の絶対値の中身に注目すると、

\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)f^n(x)微分を意味する。\prod総乗を意味し、最右辺は合成関数の微分の連鎖律により得ることができる[14]。よって、

 \lambda =\frac{1}{n} \ln \left | \prod_{i=0}^{n-1} f'(x_i) \right \vert =\frac{1}{n}  \sum_{i=0}^{n-1} \left | \ln f'(x_i) \right \vert

となる。さらに上式においてn \to \inftyとした極限値が存在するとき、その極限値を初期値 x0 から出発する軌道のリアプノフ指数と呼ぶ[16][12]

 \lambda = \lim_{n \to \infty}  \frac{1}{n} \sum_{i=0}^{n-1}  \ln | f'(x_i)|

1968年に発表されたロシアのValery Oseledecの多重エルゴード定理により、n \to \inftyの極限値が存在すること、ほとんどすべての初期値 x0λ は同じ値に収束することが証明されている[17]。特に、対象とする力学系のアトラクターの吸引域内の初期値であれば、全ての初期値で同じ λ の値に収束する[18][12]

多次元力学系のリアプノフ指数[編集]

N 次元の相空間で発展方程式 f^t で表される力学系について、各方向のリアプノフ指数 λ が得られる。これにより、次のような N 個のリアプノフ指数の組であるリアプノフスペクトラム[注釈 1]が得られる[20][21]

 \{ \lambda_1, \lambda_2, \cdots , \lambda_N \}

これは一般に出発点x_0に依存する。また、リアプノフスペクトラムでは一般に\lambda_1から値が大きい順に並べ、\lambda_1を最大リアプノフ指数と呼ぶ[9][19][22]。リアプノフ指数は対象相空間内のベクトルの振る舞いを表し、次のヤコビ行列から定義される:

 J^t(x_0) = \left. \frac{ d f^t(x) }{dx} \right|_{x_0}

J^t 行列は x_0 における小さな差がどのように最終的に f^t(x_0) まで伝達されるかを示している。次の極限

  \lim_{t \rightarrow \infty} (J^t \cdot \mathrm{Transpose}(J^t) )^{1/2t}

は、行列 L(x_0) を定義する。極限が存在することは、1次元離散力学系の場合と同様に多重エルゴード定理により証明される[17]L(x_0)固有値 \Lambda_i(x_0) とすると、リアプノフ指数 \lambda_i は次のように定義される。

 \lambda_i(x_0) = \log \Lambda_i(x_0).\,

リアプノフ指数の集合は、その力学系のエルゴード的成分の全初期値とほぼ等しい。

単位[編集]

対数の底に 2 を使用して計算した場合には、ビット/時間を単位として使用することがある。 これは、λ > 0 の場合、単位時間あたりλビットの情報が失われ、λ < 0 の場合λビットの情報が生成することに相当する[23]

基本特性[編集]

保存系の場合、相空間の全エネルギーは保存される。従って全リアプノフ指数の総和はゼロになる。散逸系ではリアプノフ指数の総和は負になる。

力学系が何らかの流れである場合、1つのリアプノフ指数は常にゼロとなる。つまり、流れの方向の固有ベクトルに対応する固有値から得られるリアプノフ指数がゼロになる。

リアプノフスペクトラムの重要性[編集]

リアプノフスペクトラムは、力学系のエントロピープロダクションやフラクタル次元の概算値を求めるのに使われる。特にリアプノフスペクトラムが分かれば、以下のように定義されるリアプノフ次元 D_{L} を計算できる[24]

 D_{L}= k + \sum_{i=1}^k \lambda_i/|\lambda_{k+1}|

ここで k は大きい方から k 個の指数の総和が負にならない最大個数である。リアプノフ次元は D_{L} は系の情報量次元の上限を表している[25]。さらに言えば、Pesin's theorem によれば、正のリアプノフ指数の総和はコルモゴロフ・シナイ・エントロピー(Kolmogorov-Sinai entropy)の近似値を与える。

最大リアプノフ指数の逆数を「リアプノフ時間; Lyapunov time」と呼ぶことがあり、e-folding time の特性を定義する。カオス的軌道ではリアプノフ時間は有限であり、正規の軌道では無限大となる。

脚注[編集]

注釈[編集]

  1. ^ リアプノフスペクトルと表記する場合もある[19]

出典[編集]

  1. ^ アリグッドほか 2012b, p. 1.
  2. ^ Strogatz 2015, pp. 349-350.
  3. ^ a b 合原 2011, p. 158.
  4. ^ a b 下條 1992, p. 86.
  5. ^ a b 高安 2001, p. 82.
  6. ^ a b Strogatz 2015, p. 400.
  7. ^ アリグッドほか 2012b, pp. 3-4.
  8. ^ 下條 1992, p. 91.
  9. ^ a b 高安 2001, p. 83.
  10. ^ a b 船越 2008, p. 146.
  11. ^ 金子、津田 1997, p. 47.
  12. ^ a b c Strogatz 2015, p. 401.
  13. ^ 合原 2011, p. 161.
  14. ^ a b 下條 1992, p. 87.
  15. ^ 船越 2008, p. 157.
  16. ^ アリグッドほか 2012a, pp. 117-118.
  17. ^ a b 合原 2011, p. 165.
  18. ^ 船越 2008, p. 171.
  19. ^ a b 金子、津田 1997, p. 115.
  20. ^ 船越 2008, p. 181.
  21. ^ 合原 2011, p. 163.
  22. ^ 船越 2008, p. 168.
  23. ^ Chaotic oscillators: theory and applications, Tomasz Kapitaniak,pp287
  24. ^ 合原 2011, p. 169.
  25. ^ 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. 
  • Govindan Rangarajan, Salman Habib, and Robert D. Ryne (1998年). "Lyapunov Exponents without Rescaling and Reorthogonalization". Physical Review Letters 80: 3747–3750. 
  • 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. 
  • 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

関連項目[編集]

外部リンク[編集]