ローレンツ方程式

ローレンツ方程式（ろーれんつほうていしき）とは、数学者・気象学者であるエドワード・ローレンツ（Edward Norton Lorenz|Edward Lorenz）が最初に研究した非線型常微分方程式である。特定のパラメータ値と初期条件に対してカオス的な解を持つことで注目されている。特に、ローレンツ方程式のカオス解の集合はローレンツ・アトラクターと呼ばれる。いわゆるバタフライ効果の説明に用いられることが多く、決定論的な連立常微分方程式が初期値鋭敏性を持つことは驚きをもって迎えられ、カオス研究の端緒となった。
概要.
1963年、エドワード・ローレンツは、数値シミュレーションや数値計算を担当したエレン・フェッターと、ローレンツ方程式の発見に至る初期の数値計算を担当したソフトウェアエンジニア、マーガレット・ハミルトンの協力を得て、大気変動の簡易数学モデルを開発した。このモデルが、現在ローレンツ方程式として知られている以下の3つの常微分方程式の系である：
この方程式は、下から暖められ、上から冷やされる2次元の流体層の特性に関するもので、3つの量の時間に対する変化率を記述しており、 は対流速度に、は水平温度変化に、は垂直温度変化に比例する。 また定数 , , はそれぞれ プラントル数、レイリー数に関する不安程度を表すパラメータ、臨界水平波数に関するパラメータである。
ローレンツ方程式はレーザー、 発電機、サーモサイフォン、 ブラシレスDCモーター、 電気回路、 化学反応、正浸透 などの簡易モデルで生じうる。また、マルクス水車のフーリエ空間での支配方程式でもある 。すなわちマルクス水車はカオス運動を示し、一定速度で一方向に回転するのではなく、その回転が加速したり減速したり停止したり方向転換したり、それらの組み合わせで前後に振動したりと予測不能の動きをする。
解析.
通常、パラメータ , , は正であると仮定する。ローレンツは、、という値を使用し、これらの値（およびその近傍の値）に対して、系がカオス的な振る舞いをすることを示している。
もしなら、均衡点は1つだけであり、それは原点である。この点は対流がないことに対応する。すべての軌道は原点に収束し、広域的なアトラクターとなる。.
でピッチフォーク分岐が起こり、 でさらに下記の2つの臨界点が現れる。
formula_2 
これらは定常対流に相当する。この二つの平衡点は
formula_3
の場合にのみ安定である。これはの時のみ、 は正となりうる。また臨界値では、両平衡点はホップ分岐を経て安定性を失う。
、、のとき、ローレンツ方程式はカオス解を持つ（ただし全ての解がカオスであるとは言えない）。ほぼ全ての初期点は、3つの平衡に関して、不変集合ローレンツアトラクター、ストレンジアトラクター、フラクタル、隠れアトラクタ#自己励起アトラクタ#自己励起アトラクタに傾くことになる。そのハウスドルフ次元は、上からリアプノフ指数 によってと見積もられる 
。また相関次元はと推定されている。
グローバルアトラクターの正確なリアプノフ指数の公式は、パラメータの古典的な制限の下で解析的に求めることができ、次に示す。
formula_4
ローレンツアトラクターは解析が難しいが、微分方程式のアトラクターへの作用はかなり単純な幾何学モデルで記述され、この証明はスメイルの問題の14番目の問題であったが、2002年にウォリック・タッカーによって初めて解決された。
の他の値では、系は結び目のある周期的な軌道を示す。例えば、では（トーラス結び目）となる。
テント写像との関連.
ローレンツの論文の図4において、 ローレンツは、系が到達したz方向の相対最大値を、方向のそれより以前の相対最大値に対してプロットした。この手順は後にローレンツマップとして知られるようになった（軌跡と所定の曲面の交点をプロットするポアンカレプロットと混同しないように）。結果としてこのプロットはテント写像に非常によく似た形をしており、ローレンツは、 の最大値があるカットオフ値を超えると、系が片方のローブ（軌跡上の片方の円盤）に切り替わることを発見した。これをテント写像で知られているカオスと組み合わせることで、系が2つのローブの間をカオス的に行き来することが判明した。
一般化されたローレンツ方程式.
ローレンツ方程式の発見後高次元のローレンツモデルに関する論文が相次ぎ、一般化されたローレンツモデルが作成された。このモデルは、3つの状態変数に対する古典的なローレンツモデル、または5つの状態変数に対する以下の5次元ローレンツモデルに単純化することができる。
formula_5
他のパラメータの値によらず、"d"0 = 19/3である。
シミュレーション.
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 = @(t,a) [-sigma*a(1) + sigma*a(2); rho*a(1) - a(2) - a(1)*a(3); -beta*a(3) + a(1)*a(2)];
[t,a] = ode45(f,[0 100],[1 1 1]); % Runge-Kutta 4th/5th order ODE solver
plot3(a(:,1),a(:,2),a(:,3))
Mathematica.
標準的な記法：
tend = 50;
eq = {x'[t] == σ (y[t] - x[t]), 
 y'[t] == x[t] (ρ - z[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};
 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}]
応用.
大気の対流モデル.
原論文で言及される通り、ローレンツ系はバリー・サルツマンが以前に研究したより大きな系を縮小したものであり、ローレンツ方程式は下から一様に加熱され、上から一様に冷却される浅い流体層における流体循環を記述する方程式をブシネスク近似から導いたものである。 この流体循環はレイリー・ベナール対流と呼ばれる。流体は2次元（垂直と水平）に循環すると仮定し、矩形の周期的境界条件を設定する。
系の流れ関数と温度をモデル化したこの偏微分方程式を，スペクトルガラーキン近似を用いて，流体力学場をフーリエ級数で展開し，流れ関数については1次，温度については2次で切り捨てる。これにより、方程式は3つの連立した非線形常微分方程式に縮小される。詳細な導出は， による非線形力学の教科書、Appendix C; , Appendix D、または Shen (2016)を参照。
大気中のカオスと秩序の性質を示すモデル.
低次元ローレンツモデルに見られるカオス的な特徴が、地球大気の特徴（気象のカオス性）を表しうることが認められている。一方で、一般化されたローレンツモデルと初期のローレンツモデルにおいてカオスと予測可能な振る舞いが共存していることから、Shenとその共著者は「気象はカオスと、明確な予測可能性を持つ秩序の両方を持っている」という改訂見解を提案した。従来の見解を発展させたこの見解は、「理論的なローレンツモデルに見られるカオスと規則的な特徴は、地球の大気の特徴をよりよく表しうる」と示唆するために用いられている。
スメイルの14番目の問題.
スメイルの14番目の問題は「ローレンツアトラクターはストレンジアトラクターの性質を持つか？」というものであった。この問題は2002年にウォリック・タッカーによって肯定的に解決された。証明には 区間演算、正準系などの厳密な数値計算が用いられた。初めにタッカーは流れの軌跡によって横方向に切断したものである断面積formula_6を定義した。ここから、各formula_7に対してformula_8 の軌跡が初めてformula_9に交わる点formula_10とするfirst-return写像 formula_11を定義できる。
続く証明は3つのポイントに分かれてなされ、ストレンジアトラクターの存在を示唆する。
第一段階では断面積formula_9がformula_16によって二つの弧に分割される事実を用いる。タッカーはこの二つの弧を小さな長方形formula_17で覆い、これらの長方形の集合がformula_18を与えることを考えた。このことを証明するためには、formula_18の全ての点がformula_18の中のformula_9に戻ってくるのをみれば良い。そのために、formula_9の下方にformula_23という小さな距離でformula_24を想定し、formula_17の中心formula_26とオイラーの積分法を用いて、点formula_26からの流れがformula_24内に与える新たな点formula_29を計算できる。そうしてテイラー展開によってformula_9内の点がformula_24内のどこに写されるかを知ることができる。これによって中心がformula_26で与えられる新たな長方形formula_33が得られる。したがってformula_17内のすべての点はformula_33の中に写される。あとは軌跡の流れがformula_9内に戻ってくるまでこの方法を再帰的に実行し、formula_37となるformula_9に入る長方形formula_39を得れば良いのであるが、問題はこの工程を何度か繰り返すうちに推定が不正確になることである。そこでタッカーはformula_33をより小さい長方形formula_41に分割し、処理を再帰的に適応した。もう一つの問題は、このアルゴリズムを適用しているうちに、流れがより「水平」になってしまい、不正確さが飛躍的に増大することで、これを防ぐためにはアルゴリズムでは断面の向きを変え、水平または垂直になるようにする必要があった。