ローレンツ方程式

ローレンツ方程式とは、数学者・気象学者であるエドワード・ローレンツが最初に研究した非線型常微分方程式である。特定のパラメータ値と初期条件に対してカオス的な解を持つことで注目されている。特に、ローレンツ方程式のカオス解の集合はローレンツ・アトラクターと呼ばれる。いわゆるバタフライ効果の説明に用いられることが多く、決定論的な連立常微分方程式が初期値鋭敏性を持つことは驚きをもって迎えられ、カオス研究の端緒となった。

1. 概要
-----

1963年、エドワード・ローレンツは、数値シミュレーションや数値計算を担当したエレン・フェッターと、ローレンツ方程式の発見に至る初期の数値計算を担当したソフトウェアエンジニア、マーガレット・ハミルトンの協力を得て、大気変動の簡易数学モデルを開発した。このモデルが、現在ローレンツ方程式として知られている以下の3つの常微分方程式の系である：

この方程式は、下から暖められ、上から冷やされる2次元の流体層の特性に関するもので、3つの量の時間に対する変化率を記述しており、 x}}は対流速度に、y}}は水平温度変化に、z}}は垂直温度変化に比例する。 また定数 σ}}, ρ}}, β}} はそれぞれ プラントル数、レイリー数に関する不安定度を表すパラメータ、臨界水平波数に関するパラメータである。

ローレンツ方程式はレーザー、 発電機、サーモサイフォン、 ブラシレスDCモーター、 電気回路、 化学反応、正浸透 などの簡易モデルで生じうる。また、マルクス水車のフーリエ空間での支配方程式でもある 。すなわちマルクス水車はカオス運動を示し、一定速度で一方向に回転するのではなく、その回転が加速したり減速したり停止したり方向転換したり、それらの組み合わせで前後に振動したりと予測不能の動きをする。

2. 解析
-----

通常、パラメータ σ}}, ρ}}, β}} は正であると仮定する。ローレンツは1=//σ//=10}}、1=//β//=sfrac、1=//ρ//=28}}という値を使用し、これらの値に対して、系がカオス的な振る舞いをすることを示している。

もし//ρ//<1}}なら、均衡点は1つだけであり、それは原点である。この点は対流がないことに対応する。すべての軌道は原点に収束し、広域的なアトラクターとなる。.

1=//ρ// = 1}}でピッチフォーク分岐が起こり、//ρ// > 1}} でさらに下記の2つの臨界点が現れる。 これらは定常対流に相当する。この二つの平衡点は の場合にのみ安定である。これは//σ// > //β// + 1}}の時のみ、 ρ}}は正となりうる。また臨界値では、両平衡点はホップ分岐を経て安定性を失う。

1=//ρ//=28}}、1=//σ//=10}}、1=//β//=sfracのとき、ローレンツ方程式はカオス解を持つ。ほぼ全ての初期点は、3つの平衡に関して、不変集合ローレンツアトラクター、ストレンジアトラクター、フラクタル、隠れアトラクタ#自己励起アトラクタ#自己励起アトラクタに傾くことになる。そのハウスドルフ次元は、上からリアプノフ指数 によって2.06と見積もられる 。また相関次元は2.05と推定されている。 グローバルアトラクターの正確なリアプノフ指数の公式は、パラメータの古典的な制限の下で解析的に求めることができ、次に示す。

ローレンツアトラクターは解析が難しいが、微分方程式のアトラクターへの作用はかなり単純な幾何学モデルで記述され、この証明はスメイルの問題の14番目の問題であったが、2002年にウォリック・タッカーによって初めて解決された。

ρ}}の他の値では、系は結び目のある周期的な軌道を示す。例えば、1=//ρ//=99.96}}では//T//}}となる。



3. テント写像との関連
------------

ローレンツの論文の図4において、 ローレンツは、系が到達したz方向の相対最大値を、z}}方向のそれより以前の相対最大値に対してプロットした。この手順は後にローレンツマップとして知られるようになった。結果としてこのプロットはテント写像に非常によく似た形をしており、ローレンツは、z}} の最大値があるカットオフ値を超えると、系が片方のローブに切り替わることを発見した。これをテント写像で知られているカオスと組み合わせることで、系が2つのローブの間をカオス的に行き来することが判明した。

4. 一般化されたローレンツ方程式
-----------------

ローレンツ方程式の発見後高次元のローレンツモデルに関する論文が相次ぎ、一般化されたローレンツモデルが作成された。このモデルは、3つの状態変数に対する古典的なローレンツモデル、または5つの状態変数に対する以下の5次元ローレンツモデルに単純化することができる。

他のパラメータの値によらず、//d//0 = 19/3である。

5. シミュレーション
-----------

5.1. MATLAB
-----------

% Solve over time interval [0,100] with initial conditions [1,1,1] % //f// is set of differential equations % //a// is array containing x, y, and z variables % //t// is time variable

sigma = 10; beta = 8/3; rho = 28; f = @ [-sigma*a + sigma*a; rho*a - a - a*a; -beta*a + a*a]; [t,a] = ode45; % Runge-Kutta 4th/5th order ODE solver plot3,a,a)<WtXmlEndTag />

5.2. Mathematica
----------------

標準的な記法：

tend = 50; eq = {x'[t] == σ , y'[t] == x[t]  - y[t], z'[t] == x[t] y[t] - β z[t]}; init = {x[0] == 10, y[0] == 10, z[0] == 10}; pars = {σ->10, ρ->28, β->8/3}; {xs, ys, zs} = NDSolveValue[{eq /. pars, init}, {x, y, z}, {t, 0, tend}]; ParametricPlot3D[{xs[t], ys[t], zs[t]}, {t, 0, tend}]

冗長性を抑えたバージョン：

lorenz = NonlinearStateSpaceModel[, {x, y, z}, {σ, ρ, β}]; soln[t_] = StateResponse[{lorenz, {10, 10, 10}}, {10, 28, 8/3}, {t, 0, 50}]; ParametricPlot3D[soln[t], {t, 0, 50}]

6. 応用
-----

6.1. 大気の対流モデル
-------------

原論文で言及される通り、ローレンツ系はバリー・サルツマンが以前に研究したより大きな系を縮小したものであり、ローレンツ方程式は下から一様に加熱され、上から一様に冷却される浅い流体層における流体循環を記述する方程式をブシネスク近似から導いたものである。 この流体循環はレイリー・ベナール対流と呼ばれる。流体は2次元に循環すると仮定し、矩形の周期的境界条件を設定する。

系の流れ関数と温度をモデル化したこの偏微分方程式を，スペクトルガラーキン近似を用いて，流体力学場をフーリエ級数で展開し，流れ関数については1次，温度については2次で切り捨てる。これにより、方程式は3つの連立した非線形常微分方程式に縮小される。詳細な導出は， Hilbornによる非線形力学の教科書、Appendix C; Bergé, Appendix D、または Shen を参照。

6.2. 大気中のカオスと秩序の性質を示すモデル
------------------------

低次元ローレンツモデルに見られるカオス的な特徴が、地球大気の特徴を表しうることが認められている。一方で、一般化されたローレンツモデルと初期のローレンツモデルにおいてカオスと予測可能な振る舞いが共存していることから、Shenとその共著者は「気象はカオスと、明確な予測可能性を持つ秩序の両方を持っている」という改訂見解を提案した。従来の見解を発展させたこの見解は、「理論的なローレンツモデルに見られるカオスと規則的な特徴は、地球の大気の特徴をよりよく表しうる」と示唆するために用いられている。

6.3. スメイルの14番目の問題
-----------------

スメイルの14番目の問題は「ローレンツアトラクターはストレンジアトラクターの性質を持つか？」というものであった。この問題は2002年にウォリック・タッカーによって肯定的に解決された。証明には 区間演算、正準系などの厳密な数値計算が用いられた。初めにタッカーは流れの軌跡によって横方向に切断したものである断面積を定義した。ここから、各に対して の軌跡が初めてに交わる点とするfirst-return写像 を定義できる。

続く証明は3つのポイントに分かれてなされ、ストレンジアトラクターの存在を示唆する。

the first-return写像で不変、すなわちなる区域が存在する。
The return写像はforward invariant cone fieldを認める
この不変円錐場内のベクトルは、return写像の微分によって一様に拡大される。第一段階では断面積がによって二つの弧に分割される事実を用いる。タッカーはこの二つの弧を小さな長方形で覆い、これらの長方形の集合がを与えることを考えた。このことを証明するためには、の全ての点がの中のに戻ってくるのをみれば良い。そのために、の下方にという小さな距離でを想定し、の中心とオイラーの積分法を用いて、点からの流れが内に与える新たな点を計算できる。そうしてテイラー展開によって内の点が内のどこに写されるかを知ることができる。これによって中心がで与えられる新たな長方形が得られる。したがって内のすべての点はの中に写される。あとは軌跡の流れが内に戻ってくるまでこの方法を再帰的に実行し、となるに入る長方形を得れば良いのであるが、問題はこの工程を何度か繰り返すうちに推定が不正確になることである。そこでタッカーはをより小さい長方形に分割し、処理を再帰的に適応した。もう一つの問題は、このアルゴリズムを適用しているうちに、流れがより「水平」になってしまい、不正確さが飛躍的に増大することで、これを防ぐためにはアルゴリズムでは断面の向きを変え、水平または垂直になるようにする必要があった。

7. ギャラリー
--------

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

バタフライ効果
カオス理論
ロジスティック方程式
レスラー方程式
常微分方程式
Eden's conjecture on the Lyapunov dimension
Lorenz 96 model
List of chaotic maps
Takens' theorem

9. 脚注
-----

10. 参考文献
--------

Shen, B.-W. . "Nonlinear feedback in a six-dimensional Lorenz model: impact of an additional heating term". //Nonlinear Processes in Geophysics//. 22 : 749–764. doi:10.5194/npg-22-749-2015. ISSN 1607-7946.

11. 詳細
------

12. 外部リンク
---------

[1] by Rob Morris, Wolfram Demonstrations Project.
[2] on planetmath.org
[3]. The implementation of Lorenz attractor in an electronic circuit.
[4] 
[5]
[6]
[7] 
[8] 
[9] made with Iodide
[10]