ローレンツ方程式
ρ = 28 , σ = 10 , β = 8 / 3 の時のローレンツアトラクターの解
ローレンツ方程式 (ローレンツほうていしき)とは、数学者・気象学者であるエドワード・ローレンツ (Edward Lorenz)が最初に研究した非線型常微分方程式 である。特定のパラメータ値と初期条件に対してカオス 的な解を持つことで注目されている。特に、ローレンツ方程式のカオス解の集合はローレンツ・アトラクター と呼ばれる。いわゆるバタフライ効果 の説明に用いられることが多く、決定論的な連立常微分方程式が初期値鋭敏性を持つことは驚きをもって迎えられ、カオス研究の端緒となった。
概要
1963年、エドワード・ローレンツは、数値シミュレーションや数値計算を担当したエレン・フェッター と、ローレンツ方程式の発見に至る初期の数値計算を担当したソフトウェアエンジニア、マーガレット・ハミルトン の協力を得て、大気変動 の簡易数学モデルを開発した[ 1] [ 2] 。このモデルが、現在ローレンツ方程式として知られている以下の3つの常微分方程式の系である:
d
x
d
t
=
σ σ -->
(
y
− − -->
x
)
,
d
y
d
t
=
x
(
ρ ρ -->
− − -->
z
)
− − -->
y
,
d
z
d
t
=
x
y
− − -->
β β -->
z
.
{\displaystyle {\begin{aligned}{\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}}&=xy-\beta z.\end{aligned}}}
この方程式は、下から暖められ、上から冷やされる2次元の流体層の特性に関するもので、3つの量の時間に対する変化率を記述しており、 x は対流速度に、y は水平温度変化に、z は垂直温度変化に比例する。 [ 3] また定数 σ , ρ , β はそれぞれ プラントル数 、レイリー数 に関する不安定度を表すパラメータ、臨界水平波数に関するパラメータである。
ローレンツ方程式はレーザー 、[ 4] 発電機 、[ 5] サーモサイフォン 、 [ 6] ブラシレスDCモーター 、 [ 7] 電気回路 、 [ 8] 化学反応 [ 9] 、正浸透 [ 10] などの簡易モデルで生じうる。また、マルクス水車 のフーリエ空間での支配方程式でもある [ 11] [ 12] 。すなわちマルクス水車はカオス運動を示し、一定速度で一方向に回転するのではなく、その回転が加速したり減速したり停止したり方向転換したり、それらの組み合わせで前後に振動したりと予測不能の動きをする。
解析
通常、パラメータ σ , ρ , β は正であると仮定する。ローレンツはσ =10 、β =8 / 3 、ρ =28 という値を使用し、これらの値(およびその近傍の値)に対して、系がカオス的な振る舞いをすることを示している[ 13] 。
もしρ <1 なら、均衡点は1つだけであり、それは原点である。この点は対流がないことに対応する。すべての軌道は原点に収束し、広域的なアトラクター となる。[ 14] .
ρ = 1 でピッチフォーク分岐 が起こり、ρ > 1 でさらに下記の2つの臨界点が現れる。
(
β β -->
(
ρ ρ -->
− − -->
1
)
,
β β -->
(
ρ ρ -->
− − -->
1
)
,
ρ ρ -->
− − -->
1
)
and
(
− − -->
β β -->
(
ρ ρ -->
− − -->
1
)
,
− − -->
β β -->
(
ρ ρ -->
− − -->
1
)
,
ρ ρ -->
− − -->
1
)
.
{\displaystyle \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).}
これらは定常対流に相当する。この二つの平衡点は
ρ ρ -->
<
σ σ -->
σ σ -->
+
β β -->
+
3
σ σ -->
− − -->
β β -->
− − -->
1
,
{\displaystyle \rho <\sigma {\frac {\sigma +\beta +3}{\sigma -\beta -1}},}
の場合にのみ安定である。これはσ > β + 1 の時のみ、 ρ は正となりうる。また臨界値では、両平衡点はホップ分岐 を経て安定性を失う。[ 15]
ρ =28 、σ =10 、β =8 / 3 のとき、ローレンツ方程式はカオス解を持つ(ただし全ての解がカオスであるとは言えない)。ほぼ全ての初期点は、3つの平衡に関して、不変集合 – ローレンツアトラクター – 、ストレンジアトラクター 、フラクタル 、自己励起アトラクタに傾くことになる。そのハウスドルフ次元 は、上からリアプノフ指数 によって 2.06± 0.01 と見積もられる[ 16]
。また相関次元 は 2.05± 0.01 と推定されている[ 17] 。
グローバルアトラクターの正確なリアプノフ指数の公式は、パラメータの古典的な制限の下で解析的に求めることができ、次に示す[ 18] [ 16] [ 19] 。
3
− − -->
2
(
σ σ -->
+
β β -->
+
1
)
σ σ -->
+
1
+
(
σ σ -->
− − -->
1
)
2
+
4
σ σ -->
ρ ρ -->
{\displaystyle 3-{\frac {2(\sigma +\beta +1)}{\sigma +1+{\sqrt {\left(\sigma -1\right)^{2}+4\sigma \rho }}}}}
ローレンツアトラクターは解析が難しいが、微分方程式のアトラクターへの作用はかなり単純な幾何学モデルで記述され[ 20] 、この証明はスメイルの問題 の14番目の問題であったが、2002年にウォリック・タッカー によって初めて解決された[ 21] 。
ρ の他の値では、系は結び目のある周期的な軌道を示す。例えば、ρ =99.96 ではT (3,2) (トーラス結び目 )となる。
異なるρ に対するローレンツ方程式の解の例
ρ = 14, σ = 10, β = 8 / 3 (拡大)
ρ = 13, σ = 10, β = 8 / 3 (拡大)
ρ = 15, σ = 10, β = 8 / 3 (拡大)
ρ = 28, σ = 10, β = 8 / 3 (拡大)
小さいρ では系は安定し、2つの固定点のうち、いずれかの点アトラクターに進展する。ρ > 24.74 では, 固定点は斥力源となり、軌道はそれらに反発して非常に複雑な形となる。
初期値に対する鋭敏性
Time t = 1 (拡大)
Time t = 2 (拡大)
Time t = 3 (拡大)
ρ = 28 , σ = 10 , β = 8 / 3 の条件で生成されたこれらの画像は二つの軌跡(青と黄色)の時間発展を示している。二つの軌跡の初期値はx 座標のみ10−5 の差がつけられている。初め、二つの軌跡は一致しているように見える(青色の上から黄色が描かれているため黄色の軌跡だけ見える)が、時間経過と共に明らかに分岐していくのがわかる。
テント写像との関連
Mathematica によって作成されたローレンツの結果の再現。赤線より上の点は、ローブが切り替わることに相当する。
ローレンツの論文[ 2] の図4において、 ローレンツは、系が到達したz方向の相対最大値を、z 方向のそれより以前の相対最大値に対してプロットした。この手順は後にローレンツマップとして知られるようになった(軌跡と所定の曲面の交点をプロットするポアンカレプロット と混同しないように)。結果としてこのプロットはテント写像 に非常によく似た形をしており、ローレンツは、z の最大値があるカットオフ値を超えると、系が片方のローブ(軌跡上の片方の円盤)に切り替わることを発見した。これをテント写像で知られているカオスと組み合わせることで、系が2つのローブの間をカオス的に行き来することが判明した。
一般化されたローレンツ方程式
ローレンツ方程式の発見後高次元のローレンツモデルに関する論文が相次ぎ、一般化されたローレンツモデルが作成された[ 22] 。このモデルは、3つの状態変数に対する古典的なローレンツモデル、または5つの状態変数に対する以下の5次元ローレンツモデルに単純化することができる[ 23] 。
d
x
d
t
=
σ σ -->
(
y
− − -->
x
)
,
d
y
d
t
=
x
(
ρ ρ -->
− − -->
z
)
− − -->
y
,
d
z
d
t
=
x
y
− − -->
x
y
1
− − -->
β β -->
z
,
d
y
1
d
t
=
x
z
− − -->
2
x
z
1
− − -->
d
0
y
1
,
d
z
1
d
t
=
2
x
y
1
− − -->
4
β β -->
z
1
.
{\displaystyle {\begin{aligned}{\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}}&=xy-xy_{1}-\beta z,\\[6pt]{\frac {\mathrm {d} y_{1}}{\mathrm {d} t}}&=xz-2xz_{1}-d_{0}y_{1},\\[6pt]{\frac {\mathrm {d} z_{1}}{\mathrm {d} t}}&=2xy_{1}-4\beta z_{1}.\end{aligned}}}
他のパラメータの値によらず、d 0 = 19/3である[ 22] [ 23] 。
シミュレーション
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 };
{ 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 [{{ σ ( y - x ), x ( ρ - z ) - y , x y - β z }, {}}, { x , y , z }, { σ , ρ , β }];
soln [ t_ ] = StateResponse [{ lorenz , { 10 , 10 , 10 }}, { 10 , 28 , 8 / 3 }, { t , 0 , 50 }];
ParametricPlot3D [ soln [ t ], { t , 0 , 50 }]
応用
大気の対流モデル
原論文で言及される通り[ 24] 、ローレンツ系はバリー・サルツマンが以前に研究したより大きな系を縮小したものであり[ 25] 、ローレンツ方程式は下から一様に加熱され、上から一様に冷却される浅い流体層における流体循環を記述する方程式をブシネスク近似 から導いたものである。[ 24] この流体循環はレイリー・ベナール対流 と呼ばれる。流体は2次元(垂直と水平)に循環すると仮定し、矩形の周期的境界条件を設定する。[ 24]
系の流れ関数 と温度をモデル化したこの偏微分方程式を,スペクトルガラーキン近似 を用いて,流体力学場をフーリエ級数で展開し,流れ関数については1次,温度については2次で切り捨てる。これにより、方程式は3つの連立した非線形常微分方程式に縮小される。詳細な導出は, Hilborn (2000) による非線形力学の教科書、Appendix C; Bergé, Pomeau & Vidal (1984) , Appendix D、または Shen (2016)[ 26] を参照。
大気中のカオスと秩序の性質を示すモデル
低次元ローレンツモデルに見られるカオス的な特徴が、地球大気の特徴(気象のカオス性)を表しうることが認められている[ 27] [ 28] [ 29] 。一方で、一般化されたローレンツモデルと初期のローレンツモデルにおいてカオスと予測可能な振る舞いが共存していることから[ 22] [ 30] [ 31] 、Shenとその共著者[ 29] [ 32] は「気象はカオスと、明確な予測可能性を持つ秩序の両方を持っている」という改訂見解を提案した。従来の見解を発展させたこの見解は、「理論的なローレンツモデルに見られるカオスと規則的な特徴は、地球の大気の特徴をよりよく表しうる」と示唆するために用いられている。
スメイルの14番目の問題
スメイルの14番目の問題は「ローレンツアトラクターはストレンジアトラクター の性質を持つか?」というものであった。この問題は2002年にウォリック・タッカーによって肯定的に解決された[ 21] 。証明には 区間演算 、正準系 などの厳密な数値計算が用いられた。初めにタッカーは流れの軌跡によって横方向に切断したものである断面積
Σ Σ -->
⊂ ⊂ -->
{
x
3
=
r
− − -->
1
}
{\displaystyle \Sigma \subset \{x_{3}=r-1\}}
を定義した。ここから、各
x
∈ ∈ -->
Σ Σ -->
{\displaystyle x\in \Sigma }
に対して
x
{\displaystyle x}
の軌跡が初めて
Σ Σ -->
{\displaystyle \Sigma }
に交わる点
P
(
x
)
{\displaystyle P(x)}
とするfirst-return写像
P
{\displaystyle P}
を定義できる。
続く証明は3つのポイントに分かれてなされ、ストレンジアトラクターの存在を示唆する[ 33] 。
the first-return写像で不変、すなわち
P
(
N
)
⊂ ⊂ -->
N
{\displaystyle P(N)\subset N}
なる区域
N
⊂ ⊂ -->
Σ Σ -->
{\displaystyle N\subset \Sigma }
が存在する。
The return写像はforward invariant cone fieldを認める
この不変円錐場内のベクトルは、return写像の微分
D
P
{\displaystyle DP}
によって一様に拡大される。
第一段階では断面積
Σ Σ -->
{\displaystyle \Sigma }
が
P
(
Σ Σ -->
)
{\displaystyle P(\Sigma )}
によって二つの弧に分割される事実を用いる[ 33] 。タッカーはこの二つの弧を小さな長方形
R
i
{\displaystyle R_{i}}
で覆い、これらの長方形の集合が
N
{\displaystyle N}
を与えることを考えた。このことを証明するためには、
N
{\displaystyle N}
の全ての点が
N
{\displaystyle N}
の中の
Σ Σ -->
{\displaystyle \Sigma }
に戻ってくるのをみれば良い。そのために、
Σ Σ -->
{\displaystyle \Sigma }
の下方に
h
{\displaystyle h}
という小さな距離で
Σ Σ -->
′
{\displaystyle \Sigma '}
を想定し、
R
i
{\displaystyle R_{i}}
の中心
c
i
{\displaystyle c_{i}}
とオイラーの積分法を用いて、点
c
i
{\displaystyle c_{i}}
からの流れが
Σ Σ -->
′
{\displaystyle \Sigma '}
内に与える新たな点
c
i
′
{\displaystyle c_{i}'}
を計算できる。そうしてテイラー展開によって
Σ Σ -->
{\displaystyle \Sigma }
内の点が
Σ Σ -->
′
{\displaystyle \Sigma '}
内のどこに写されるかを知ることができる。これによって中心が
c
i
{\displaystyle c_{i}}
で与えられる新たな長方形
R
i
′
{\displaystyle R_{i}'}
が得られる。したがって
R
i
{\displaystyle R_{i}}
内のすべての点は
R
i
′
{\displaystyle R_{i}'}
の中に写される。あとは軌跡の流れが
Σ Σ -->
{\displaystyle \Sigma }
内に戻ってくるまでこの方法を再帰的に実行し、
P
(
R
i
)
⊂ ⊂ -->
R
f
i
{\displaystyle P(R_{i})\subset Rf_{i}}
となる
Σ Σ -->
{\displaystyle \Sigma }
に入る長方形
R
f
i
{\displaystyle Rf_{i}}
を得れば良いのであるが、問題はこの工程を何度か繰り返すうちに推定が不正確になることである。そこでタッカーは
R
i
′
{\displaystyle R_{i}'}
をより小さい長方形
R
i
,
j
{\displaystyle R_{i,j}}
に分割し、処理を再帰的に適応した。もう一つの問題は、このアルゴリズムを適用しているうちに、流れがより「水平」になってしまい[ 33] 、不正確さが飛躍的に増大することで、これを防ぐためにはアルゴリズムでは断面の向きを変え、水平または垂直になるようにする必要があった。
ギャラリー
A solution in the xz 平面で描写された高解像度のローレンツアトラクター
ローレンツ方程式での複数の軌跡のアニメーション
金属製ワイヤーで製作されたローレンツアトラクターの3Dモデル
ローレンツ系の近傍解の発散を示すアニメーション
断続的なサイクルにあるローレンツアトラクターを可視化したもの
ローレンツ系における2つの流線(ρ = 0 から ρ = 28 , σ = 10 , β = 8 / 3 )
ローレンツ系のρ 依存性のアニメーション
Brain Dynamics Toolboxのローレンツアトラクターのアニメーション
[ 34]
関連項目
脚注
^ Lorenz (1960)
^ a b Lorenz (1963)
^ Sparrow (1982)
^ Haken (1975)
^ Knobloch (1981)
^ Gorman, Widmann & Robbins (1986)
^ Hemati (1994)
^ Cuomo & Oppenheim (1993)
^ Poland (1993)
^ Tzenov (2014) [要出典 ]
^ Kolář & Gumbs (1992)
^ Mishra & Sanghi (2006)
^ Hirsch, Smale & Devaney (2003) , pp.303-305
^ Hirsch, Smale & Devaney (2003) , pp.306+307
^ Hirsch, Smale & Devaney (2003) , pp. 307–308
^ a b Kuznetsov, N.V.; Mokaev, T.N.; Kuznetsova, O.A.; Kudryashova, E.V. (2020). “The Lorenz system: hidden boundary of practical stability and the Lyapunov dimension”. Nonlinear Dynamics 102 (2): 713–732. doi :10.1007/s11071-020-05856-4 .
^ Grassberger & Procaccia (1983)
^ Leonov et al. (2016)
^ Kuznetsov, Nikolay; Reitmann, Volker (2021). Attractor Dimension Estimates for Dynamical Systems: Theory and Computation . Cham: Springer. https://www.springer.com/gp/book/9783030509866
^ Guckenheimer, John; Williams, R. F. (1979-12-01). “Structural stability of Lorenz attractors” . Publications Mathématiques de l'Institut des Hautes Études Scientifiques 50 (1): 59–72. doi :10.1007/BF02684769 . ISSN 0073-8301 . http://www.numdam.org/item/PMIHES_1979__50__59_0/ .
^ a b Tucker (2002)
^ a b c Shen, Bo-Wen (2019-03-01). “Aggregated Negative Feedback in a Generalized Lorenz Model” . International Journal of Bifurcation and Chaos 29 (3): 1950037–1950091. Bibcode : 2019IJBC...2950037S . doi :10.1142/S0218127419500378 . ISSN 0218-1274 . https://www.worldscientific.com/doi/abs/10.1142/S0218127419500378 .
^ a b Shen, Bo-Wen (2014-04-28). “Nonlinear Feedback in a Five-Dimensional Lorenz Model” . Journal of the Atmospheric Sciences 71 (5): 1701–1723. Bibcode : 2014JAtS...71.1701S . doi :10.1175/jas-d-13-0223.1 . ISSN 0022-4928 . https://doi.org/10.1175/jas-d-13-0223.1 .
^ a b c Lorenz (1963)
^ Saltzman (1962)
^ 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. Bibcode : 2015NPGeo..22..749S . doi :10.5194/npg-22-749-2015 . ISSN 1607-7946 . https://npg.copernicus.org/articles/22/749/2015/ .
^ Ghil, Michael; Read, Peter; Smith, Leonard (2010-07-23). “Geophysical flows as dynamical systems: the influence of Hide's experiments” . Astronomy & Geophysics 51 (4): 4.28–4.35. Bibcode : 2010A&G....51d..28G . doi :10.1111/j.1468-4004.2010.51428.x . ISSN 1366-8781 . https://doi.org/10.1111/j.1468-4004.2010.51428.x .
^ Read, P. (1993). Application of Chaos to Meteorology and Climate. In The Nature of Chaos; Mullin, T., Ed . Oxford, UK: Oxford Science Publications. pp. 220–260. ISBN 0198539541
^ a b Shen, Bo-Wen; Pielke, Roger; Zeng, Xubin; Cui, Jialin; Faghih-Naini, Sara; Paxson, Wei; Kesarkar, Amit; Zeng, Xiping et al. (2022-11-12). “The Dual Nature of Chaos and Order in the Atmosphere” (英語). Atmosphere 13 (11): 1892. Bibcode : 2022Atmos..13.1892S . doi :10.3390/atmos13111892 . ISSN 2073-4433 .
^ Yorke, James A.; Yorke, Ellen D. (1979-09-01). “Metastable chaos: The transition to sustained chaotic behavior in the Lorenz model” (英語). Journal of Statistical Physics 21 (3): 263–277. Bibcode : 1979JSP....21..263Y . doi :10.1007/BF01011469 . ISSN 1572-9613 . https://doi.org/10.1007/BF01011469 .
^ Shen, Bo-Wen; Pielke, R. A.; Zeng, X.; Baik, J.-J.; Faghih-Naini, S.; Cui, J.; Atlas, R.; Reyes, T. A. L. (2021), Skiadas, Christos H.; Dimotikalis, Yiannis, eds., “Is Weather Chaotic? Coexisting Chaotic and Non-chaotic Attractors Within Lorenz Models” (英語), 13th Chaotic Modeling and Simulation International Conference (Cham: Springer International Publishing): pp. 805–825, doi :10.1007/978-3-030-70795-8_57 , ISBN 978-3-030-70794-1 , https://link.springer.com/10.1007/978-3-030-70795-8_57 2022年12月22日 閲覧。
^ Shen, Bo-Wen; Pielke, Roger A.; Zeng, Xubin; Baik, Jong-Jin; Faghih-Naini, Sara; Cui, Jialin; Atlas, Robert (2021-01-01). “Is Weather Chaotic?: Coexistence of Chaos and Order within a Generalized Lorenz Model” (英語). Bulletin of the American Meteorological Society 102 (1): E148–E158. Bibcode : 2021BAMS..102E.148S . doi :10.1175/BAMS-D-19-0165.1 . ISSN 0003-0007 . https://journals.ametsoc.org/view/journals/bams/102/1/BAMS-D-19-0165.1.xml .
^ a b c Viana (2000)
^ Heitmann, S., Breakspear, M (2017-2022) Brain Dynamics Toolbox. bdtoolbox.org doi.org/10.5281/zenodo.5625923
参考文献
Bergé, Pierre; Pomeau, Yves; Vidal, Christian (1984). Order within Chaos: Towards a Deterministic Approach to Turbulence . New York: John Wiley & Sons . ISBN 978-0-471-84967-4
Cuomo, Kevin M.; Oppenheim, Alan V. (1993). “Circuit implementation of synchronized chaos with applications to communications”. Physical Review Letters 71 (1): 65–68. Bibcode : 1993PhRvL..71...65C . doi :10.1103/PhysRevLett.71.65 . ISSN 0031-9007 . PMID 10054374 .
Gorman, M.; Widmann, P.J.; Robbins, K.A. (1986). “Nonlinear dynamics of a convection loop: A quantitative comparison of experiment with theory”. Physica D 19 (2): 255–267. Bibcode : 1986PhyD...19..255G . doi :10.1016/0167-2789(86)90022-9 .
Grassberger, P.; Procaccia, I. (1983). “Measuring the strangeness of strange attractors”. Physica D 9 (1–2): 189–208. Bibcode : 1983PhyD....9..189G . doi :10.1016/0167-2789(83)90298-1 .
Haken, H. (1975). “Analogy between higher instabilities in fluids and lasers”. Physics Letters A 53 (1): 77–78. Bibcode : 1975PhLA...53...77H . doi :10.1016/0375-9601(75)90353-9 .
Hemati, N. (1994). “Strange attractors in brushless DC motors”. IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications 41 (1): 40–45. doi :10.1109/81.260218 . ISSN 1057-7122 .
Hilborn, Robert C. (2000). Chaos and Nonlinear Dynamics: An Introduction for Scientists and Engineers (second ed.). Oxford University Press . ISBN 978-0-19-850723-9
Hirsch, Morris W. ; Smale, Stephen ; Devaney, Robert (2003). Differential Equations, Dynamical Systems, & An Introduction to Chaos (Second ed.). Boston, MA: Academic Press . ISBN 978-0-12-349703-1
Knobloch, Edgar (1981). “Chaos in the segmented disc dynamo”. Physics Letters A 82 (9): 439–440. Bibcode : 1981PhLA...82..439K . doi :10.1016/0375-9601(81)90274-7 .
Kolář, Miroslav; Gumbs, Godfrey (1992). “Theory for the experimental observation of chaos in a rotating waterwheel”. Physical Review A 45 (2): 626–637. Bibcode : 1992PhRvA..45..626K . doi :10.1103/PhysRevA.45.626 . PMID 9907027 .
Leonov, G.A.; Kuznetsov, N.V.; Korzhemanova, N.A.; Kusakin, D.V. (2016). “Lyapunov dimension formula for the global attractor of the Lorenz system”. Communications in Nonlinear Science and Numerical Simulation 41 : 84–103. arXiv :1508.07498 . Bibcode : 2016CNSNS..41...84L . doi :10.1016/j.cnsns.2016.04.032 .
Lorenz, Edward Norton (1963). “Deterministic nonperiodic flow”. Journal of the Atmospheric Sciences 20 (2): 130–141. Bibcode : 1963JAtS...20..130L . doi :10.1175/1520-0469(1963)020<0130:DNF>2.0.CO;2 .
Mishra, Aashwin; Sanghi, Sanjeev (2006). “A study of the asymmetric Malkus waterwheel: The biased Lorenz equations”. Chaos: An Interdisciplinary Journal of Nonlinear Science 16 (1): 013114. Bibcode : 2006Chaos..16a3114M . doi :10.1063/1.2154792 . PMID 16599745 .
Pchelintsev, A.N. (2014). “Numerical and Physical Modeling of the Dynamics of the Lorenz System”. Numerical Analysis and Applications 7 (2): 159–167. doi :10.1134/S1995423914020098 .
Poland, Douglas (1993). “Cooperative catalysis and chemical chaos: a chemical model for the Lorenz equations”. Physica D 65 (1): 86–99. Bibcode : 1993PhyD...65...86P . doi :10.1016/0167-2789(93)90006-M .
Saltzman, Barry (1962). “Finite Amplitude Free Convection as an Initial Value Problem—I”. Journal of the Atmospheric Sciences 19 (4): 329–341. Bibcode : 1962JAtS...19..329S . doi :10.1175/1520-0469(1962)019<0329:FAFCAA>2.0.CO;2 .
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.
Sparrow, Colin (1982). The Lorenz Equations: Bifurcations, Chaos, and Strange Attractors . Springer
Tucker, Warwick (2002). “A Rigorous ODE Solver and Smale's 14th Problem” . Foundations of Computational Mathematics 2 (1): 53–117. doi :10.1007/s002080010018 . https://link.springer.com/content/pdf/10.1007/s002080010018.pdf .
Tzenov, Stephan (2014). "Strange Attractors Characterizing the Osmotic Instability". arXiv :1406.0979v1 [physics.flu-dyn ]。
Viana, Marcelo (2000). “What's new on Lorenz strange attractors?”. The Mathematical Intelligencer 22 (3): 6–19. doi :10.1007/BF03025276 .
Lorenz, Edward N. (1960). “The statistical prediction of solutions of dynamic equations.” . Symposium on Numerical Weather Prediction in Tokyo . http://eaps4.mit.edu/research/Lorenz/The_Statistical_Prediction_of_Solutions_1962.pdf .
詳細
外部リンク
Hazewinkel, Michiel, ed. (2001), “Lorenz attractor” , Encyclopaedia of Mathematics , Springer , ISBN 978-1-55608-010-4 , http://eom.springer.de/p/l060890.htm
Weisstein, Eric W. "Lorenz attractor" . mathworld.wolfram.com (英語).
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
ローレンツアトラクタを描画する、あるいはそれに類似することをするときのために (英語版 )