ローレンツ方程式

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

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

:  \begin{align}
\frac{\mathrm{d}x}{\mathrm{d}t} &= \sigma (y - x), \\[6pt]
\frac{\mathrm{d}y}{\mathrm{d}t} &= x (\rho - z) - y, \\[6pt]
\frac{\mathrm{d}z}{\mathrm{d}t} &= x y - \beta z.
\end{align} 

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

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

## 解析
通常、パラメータ , ,  は正であると仮定する。ローレンツは、、という値を使用し、これらの値（およびその近傍の値）に対して、系がカオス的な振る舞いをすることを示している。

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

でピッチフォーク分岐が起こり、 でさらに下記の2つの臨界点が現れる。
\left( \sqrt{\beta(\rho-1)}, \sqrt{\beta(\rho-1)}, \rho-1 \right) \quad\text{and}\quad \left( -\sqrt{\beta(\rho-1)}, -\sqrt{\beta(\rho-1)}, \rho-1 \right).  
これらは定常対流に相当する。この二つの平衡点は
\rho 
の場合にのみ安定である。これはの時のみ、 は正となりうる。また臨界値では、両平衡点はホップ分岐を経て安定性を失う。

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

 3 - \frac{2 (\sigma + \beta + 1)}{\sigma + 1 + \sqrt{\left(\sigma-1\right)^2 + 4 \sigma \rho}} 

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

の他の値では、系は結び目のある周期的な軌道を示す。例えば、では（トーラス結び目）となる。

## テント写像との関連

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

## 一般化されたローレンツ方程式
ローレンツ方程式の発見後高次元のローレンツモデルに関する論文が相次ぎ、一般化されたローレンツモデルが作成された。このモデルは、3つの状態変数に対する古典的なローレンツモデル、または5つの状態変数に対する以下の5次元ローレンツモデルに単純化することができる。
 \begin{align}
\frac{\mathrm{d}x}{\mathrm{d}t} &= \sigma (y - x), \\[6pt]
\frac{\mathrm{d}y}{\mathrm{d}t} &= x (\rho - z) - y, \\[6pt]
\frac{\mathrm{d}z}{\mathrm{d}t} &= x y - x y_1 - \beta z, \\[6pt]
\frac{\mathrm{d}y_1}{\mathrm{d}t} &= x z - 2 x z_1 - d_0 y_1, \\[6pt]
\frac{\mathrm{d}z_1}{\mathrm{d}t} &= 2 x y_1 - 4\beta z_1.
\end{align} 

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

## シミュレーション

### MATLAB

```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
標準的な記法：


```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};
{xs, ys, zs} = 
  NDSolveValue[{eq /. pars, init}, {x, y, z}, {t, 0, tend}];
ParametricPlot3D[{xs[t], ys[t], zs[t]}, {t, 0, tend}]
```


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


```mathematica
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年にウォリック・タッカーによって肯定的に解決された。証明には 区間演算、正準系などの厳密な数値計算が用いられた。初めにタッカーは流れの軌跡によって横方向に切断したものである断面積\Sigma\subset \{x_3 = r - 1 \}を定義した。ここから、各x\in\Sigmaに対してx の軌跡が初めて\Sigmaに交わる点P(x)とするfirst-return写像 Pを定義できる。

続く証明は3つのポイントに分かれてなされ、ストレンジアトラクターの存在を示唆する。
*  the first-return写像で不変、すなわちP(N)\subset Nなる区域N\subset\Sigmaが存在する。
* The return写像はforward invariant cone fieldを認める
* この不変円錐場内のベクトルは、return写像の微分DPによって一様に拡大される。

第一段階では断面積\SigmaがP(\Sigma)によって二つの弧に分割される事実を用いる。タッカーはこの二つの弧を小さな長方形R_iで覆い、これらの長方形の集合がNを与えることを考えた。このことを証明するためには、Nの全ての点がNの中の\Sigmaに戻ってくるのをみれば良い。そのために、\Sigmaの下方にhという小さな距離で\Sigma'を想定し、R_iの中心c_iとオイラーの積分法を用いて、点c_iからの流れが\Sigma'内に与える新たな点c_i'を計算できる。そうしてテイラー展開によって\Sigma内の点が\Sigma'内のどこに写されるかを知ることができる。これによって中心がc_iで与えられる新たな長方形R_i'が得られる。したがってR_i内のすべての点はR_i'の中に写される。あとは軌跡の流れが\Sigma内に戻ってくるまでこの方法を再帰的に実行し、P(R_i)\subset Rf_iとなる\Sigmaに入る長方形Rf_iを得れば良いのであるが、問題はこの工程を何度か繰り返すうちに推定が不正確になることである。そこでタッカーはR_i'をより小さい長方形R_{i,j}に分割し、処理を再帰的に適応した。もう一つの問題は、このアルゴリズムを適用しているうちに、流れがより「水平」になってしまい、不正確さが飛躍的に増大することで、これを防ぐためにはアルゴリズムでは断面の向きを変え、水平または垂直になるようにする必要があった。

## ギャラリー

File:Lorenz system r28 s10 b2-6666.png|A solution in the 平面で描写された高解像度のローレンツアトラクター
File:Lorenz attractor.svg|SVGでレンダリングされたローレンツアトラクター
File:A Lorenz system.ogv|ローレンツ方程式での複数の軌跡のアニメーション
File:Lorenzstill-rubel.png|金属製ワイヤーで製作されたローレンツアトラクターの3Dモデル
File:Lorenz.ogv|ローレンツ系の近傍解の発散を示すアニメーション
File:Intermittent Lorenz Attractor - Chaoscope.jpg|断続的なサイクルにあるローレンツアトラクターを可視化したもの
File:Lorenz apparition small.gif|ローレンツ系における2つの流線（ から , , ）
File:Lorenz(rho).gif|ローレンツ系の依存性のアニメーション
File:Lorenz Attractor Brain Dynamics Toolbox.gif|Brain Dynamics Toolboxのローレンツアトラクターのアニメーション


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

## 脚注

## 参考文献
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* 
* Shen, B.-W. (2015-12-21). "Nonlinear feedback in a six-dimensional Lorenz model: impact of an additional heating term". *Nonlinear Processes in Geophysics*. **22** (6): 749–764. doi:10.5194/npg-22-749-2015. ISSN 1607-7946.
* 
* 
* 
* 
* 

## 詳細
* 
* 

## 外部リンク

* 
* 
* Lorenz attractor by Rob Morris, Wolfram Demonstrations Project.
* Lorenz equation on planetmath.org
* Synchronized Chaos and Private Communications, with Kevin Cuomo. The implementation of Lorenz attractor in an electronic circuit.
* Lorenz attractor interactive animation (you need the Adobe Shockwave plugin)
* 3D Attractors: Mac program to visualize and explore the Lorenz attractor in 3 dimensions
* Lorenz Attractor implemented in analog electronic
* Lorenz Attractor interactive animation (implemented in Ada with GTK+. Sources & executable)
* Web based Lorenz Attractor (implemented in JavaScript/HTML/CSS)
* Interactive web based Lorenz Attractor made with Iodide
* ローレンツアトラクタを描画する、あるいはそれに類似することをするときのために （英語版）