リアプノフ指数

リアプノフ指数とは、力学系においてごく接近した軌道が離れていく度合いを表す量である。リャプノフ指数とも表記されるアリグッドほか。ロシア人科学者 Алекса́ндр Ляпуно́в}}にその名をちなむ。

系の相空間上の2つの軌道について考える。2つの軌道上の時刻 //t// における点の距離をベクトル // として、初期状態 //t// = 0 には、これらの軌道は距離 // だけ離れているとする。// を近似的に次のように表すStrogatz合原下條。

ここで はユークリッドノルムを意味する。上式で //λ// > 0 の場合は軌道は離れていき、 //λ// < 0 の場合は軌道は近づいていく。よって、軌道が離れていく度合いは //λ// の値により決定される。この //λ// がリアプノフ指数である合原下條。軌道がカオス的であるとき、上式のように軌道は指数関数的に離れていく高安Strogatz。すなわち、リアプノフ指数が正であることが軌道がカオス的であることの1つの定義とされるアリグッドほか。

より詳細には、系の状態変数が //k// 個の場合、すなわち相空間が //k// 次元である場合は各次元ごとに固有のリアプノフ指数を持つ下條。これらのリアプノフ指数の組をリアプノフスペクトラムと呼び高安、そのうちの最大のリアプノフ指数を最大リアプノフ指数と呼ぶ船越。各々のリアプノフ指数を見れば正であったり負であったりするが、最大リアプノフ指数が正であれば、その系はカオスの特徴の1つである初期値鋭敏性を持つといえる金子、津田船越。

1. 1次元離散時間力学系のリアプノフ指数
---------------------

まず、単純な1次元離散力学系の場合のリアプノフ指数について説明する。を系の状態変数、を離散時間としたとき、 写像 //x////n//+1 = //f// のリアプノフ指数 //λ// は次のように定義されるStrogatz高安合原下條。

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

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

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

と表すStrogatz。両辺の自然対数をとると、

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

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

とできる。ここで ' は、//fn// の微分を意味する。は総乗を意味し、最右辺は合成関数の微分の連鎖律により得ることができる下條。よって、

となる。さらに上式において //n// → ∞ とした極限値が存在するとき、その極限値を初期値 //x//0 から出発する軌道のリアプノフ指数と呼ぶアリグッドほかStrogatz。

1968年に発表されたValery Oseledecの多重エルゴード定理により、//n// → ∞ の極限値が存在すること、ほとんどすべての初期値 //x//0 で //λ// は同じ値に収束することが証明されている合原。対象とする力学系のアトラクターの吸引域内の初期値であれば、全ての初期値で同じ //λ// の値に収束する船越Strogatz。

2. 高次元力学系のリアプノフ指数
-----------------

力学系が //k// 次元の相空間を持つ高次元力学系の場合、各方向に別々のリアプノフ指数が存在する。すなわち高次元力学系であれば、軌道のずれは、ある方向には離れていくが、別の方向では縮まっていく状況がありえるアリグッドほか。よって //k// 個のリアプノフ指数を得ることができる。このような //k// 個のリアプノフ指数の組を、リアプノフスペクトラムと呼ぶ船越合原リアプノフスペクトルと表記する場合もあるSfn。

リアプノフスペクトラムでは一般に //λ//1 から値が大きい順に並べる金子、津田高安。最大値である //λ//1 を、特に最大リアプノフ指数と呼ぶ船越。記事冒頭で述べたように、相空間上の2つの軌道上の時刻 //t// における点の間の距離、すなわちずれを // とする。リアプノフスペクトラム //λi// は以下のように定義される合原高安ベルジュほか。

一般に //λi// は初期値 // に依存する合原。しかし1次元離散力学系の場合と同様に、ほとんどすべての初期位置 //0 から同一の //λi// を得ることができる合原高安。//λi// の定義式にある //αi// は、次式で定義される //k// × //k// 正定値行列 // の固有値である合原高安。

さらに // は、// の解を次の形式で表したときの // に対する乗数として得られる合原高安ベルジュほか。

系が連続力学系の場合、 //k// 個の状態変数 {//x//1, //x//2, ..., //xk//}、常微分方程式 {//f//1, //f//2, ..., //fk//} から成る常微分方程式系が与えられる小室。// が線形近似可能な場合、// のヤコビ行列を用いて、

と表すことができる合原高安。また、系が離散力学系の場合、 //k// 個の状態変数、常差分方程式から成る差分方程式系 // = //)が与えられる小室。同じく、差分方程式系 // が線形近似可能な場合、// のヤコビ行列を用いて、

と表すことができる合原下條。ここに、// は以下に示すようなヤコビ行列による線形写像で、軌道 // に依存し、すなわち初期値 //、時間 //t// に依存して変化する合原高安。

常微分方程式系の場合は、を解いてを得ることで、上記の定義で出てきた正方行列 // を得ることができる高安ベルジュほか。差分方程式系の場合は、// を //t// 回繰り返し適用することで次のような // と // の関係式で書き表すことできるので、// は //Jn// の 0 から //t// − 1 までの総乗として得ることができる合原アリグッドほか。

3. 単位
-----

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

4. 基本特性
-------

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

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

Pesin's theorem によれば、正のリアプノフ指数の総和はコルモゴロフ・シナイ・エントロピーの近似値を与える。

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

5. リアプノフ次元
----------

軌道がカオス的振る舞いをみせるストレンジアトラクターは、フラクタル構造を持つことが多い高安。このようなアトラクターのフラクタル次元とリアプノフスペクトラムの間には関係が存在するアリグッドほか。アトラクターのリアプノフスペクトラムが得られたとして、その各リアプノフ指数 //λi// がそれらの値の大きさで降順に並べられているとしたとき、次のようなフラクタル次元の一種 //DL// が定義される合原。

ここで、 //ξj// は、

であり、//j// は // ξj// が負にならない最大値、すなわち、リアプノフ指数を //λ//1 + //λ//2 + //λ//3 ...と順に足していったときに、その総和が負となる直前における、足し合わせたリアプノフ指数の個数である。//ξj// ≥ 0 を満たす //j// が存在しないときは、//DL// = 0、系の次元数 //k// と //j// が一致する場合は、//DL// = //k// とするアリグッドほか。

このように定義されたフラクタル次元 //DL// は、リアプノフ次元 と呼ばれる。リアプノフ次元は、James L. KaplanとJames A. Yorkeにより1979年に提案された。そのためカプラン・ヨーク次元とも呼ばれ、記号 //DKY// とも記される。

上式で定義されるリアプノフ次元は、フラクタル次元の1つである容量次元の考え方をもとにして、次のように導入される。リアプノフ指数の総和は相空間内の //k// 次元の体積要素の体積変化率を与えるので、同様に部分和である //ξj// から、//j// 次元までの体積要素は拡大するが、//j// + 1 次元以上の体積要素は縮小することになる。そのため、アトラクタを収めるためには少なくとも //j// 次元は必要で、考えられるフラクタル次元の下限を与えている。一方で、//j// + 1 次元は考えられるフラクタル次元の上限といえる。//k// 次元相空間上の一辺が //d// の立方体は、時間発展により各辺 exp, exp, exp,... の直方体に写像されるアリグッドほか。ここで、１辺の長さ //ε// が、

である //j// + 1 次元立方体の箱を考え、容量次元と同じように、ストレンジアトラクタの不変集合をこの箱を何個も当てながら全体を覆うことを考えるアリグッドほか。このとき、覆うのに必要な箱の数 //N// は、

と推論できるアリグッドほか。よって、容量次元と同様の定義から、

となり、//j// の定義より −//λ////j// + 1 = |//λ////j// + 1| なので、リアプノフ次元の定義

を得ることができるアリグッドほか。

他のフラクタル次元と比較したリアプノフ次元の利点は、リアプノフスペクトラムさえ得ることができれば簡単に計算可能なことであるアリグッドほか。また、リアプノフ次元は系の情報量次元の上限を表している。

6. 脚注
-----

6.1. 注釈
-------

6.2. 出典
-------

7. 参考文献
-------

Cvitanovi? P., Artuso R., Mainieri R. , Tanner G. and Vattay G.[1] Niels Bohr Institute, Copenhagen 2005

8. 関連項目
-------

エルゴード理論
リウヴィルの定理 
フロケ理論

9. 外部リンク
--------