a = 1.4 、b = 0.3 の場合のエノンアトラクターa = 1.4 、b = 0.3 の場合のエノンアトラクター数学 において、エノン写像は 離散時間 力学系 です。これは、カオス的挙動 を示す力学系の中で最もよく研究されている例の一つです。エノン写像は、平面上の点( x n , y n ) を新しい点に写像します。{ x n + 1 = 1 − 1 x n 2 + y n y n + 1 = b x n {\displaystyle {\begin{cases}x_{n+1}=1-ax_{n}^{2}+y_{n}\\y_{n+1}=bx_{n}\end{cases}}} このマップは、2 つのパラメータa とbに依存しており、 古典的なエノンマップ では、a = 1.4 、b = 0.3 の値が取られます。[ 1 ] 古典的な値の場合、エノンマップはカオス的です。a と b の他の値の場合、 マップ はカオス的、断続的、または 周期軌道 に収束する可能性があります。さまざまなパラメータ値でのマップの挙動の概要は、軌道図 で確認できます。
このマップは、ミシェル・エノンによって ローレンツ系 のポアンカレ部分 の簡略化されたモデルとして導入されました。[ 1 ] 古典的なマップの場合、平面上の初期点は、エノンのストレンジアトラクター として知られる点の集合に近づくか、無限に発散します。エノンアトラクターはフラクタル であり、一方向には滑らかで、別の方向にはカントール集合です 。[ 1 ] 古典的なマップのストレンジアトラクターのフラクタル次元 の数値推定では、相関次元 が 1.21 ± 0.01 [ 2 ] 、ボックスカウント次元 が 1.261 ± 0.003 [ 3 ]となります。
ダイナミクス
アトラクター エノン写像は、定数ヤコビ行列式 を持つ2次元微分同相写像 である。写像のヤコビ行列は次のとおりである。J = [ − 2 1 x 1 b 0 ] {\displaystyle J={\begin{bmatrix}-2ax&1\\b&0\end{bmatrix}}} この行列の行列式はd e t ( J ) = − b {\displaystyle det(J)=-b} 写像は散逸的であるため(つまり、反復によって体積が縮小するため)、行列式は -1 から 1 の間になければなりません。Hénon 写像は1 の 場合に散逸的です。[ 4 ] 古典的なパラメータの場合1 = 1.4 、 b = 0.3 {\displaystyle a=1.4,b=0.3} 行列式は-0.3なので、マップは一定の割合で領域を縮小します。反復ごとに領域は0.3倍ずつ縮小します。
この収縮は、伸縮と折り畳み動作と組み合わさって、エノンアトラクターの特徴的なフラクタル構造を生み出します。古典的なパラメータの場合、ほとんどの初期条件は、このブーメランのような形状を描き出す軌道につながります。アトラクターには、その構造の根幹となる無数の不安定な周期軌道が含まれています。[ 5 ]
固定点 この地図には、写像によって変化しない2つの固定点があります。これらは、 x = 1 - ax² + y とy = bxを解くことによって求められます。2 番目 の式を最初の式に代入すると、次の二次方程式が得られます。1 x 2 + ( 1 − b ) x − 1 = 0 {\displaystyle ax^{2}+(1-b)x-1=0} 解(固定点のx座標)は以下のとおりです。x = − ( 1 − b ) ± ( 1 − b ) 2 + 4 1 2 1 {\displaystyle x={\frac {-(1-b)\pm {\sqrt {(1-b)^{2}+4a}}}{2a}}} 古典的なパラメータa = 1.4 およびb = 0.3 の場合、2 つの固定点は次のようになります。
x 1 ≈ 0.631 、 y 1 ≈ 0.189 {\displaystyle x_{1}\approx 0.631,\quad y_{1}\approx 0.189}
x 2 ≈ − 1.131 、 y 2 ≈ − 0.339 {\displaystyle x_{2}\approx -1.131,\quad y_{2}\approx -0.339}
これらの点の安定性は、固定点におけるヤコビ行列Jの 固有値によって決定されます。古典的な写像では、最初の固定点は 鞍点 (不安定)であり、2 番目の固定点は反発点(これも不安定)です。[ 6 ] 最初の固定点の不安定多様体 は、ストレンジ アトラクター自体を生成する重要な構成要素です。[ 6 ]
エノン写像の軌道図 。b = 0.3 に固定し、a を変化させた場合を示す。この図は、周期倍分岐を経て カオスに至る古典的な経路を示している。
分岐図 エノン写像は、そのパラメータが変化すると複雑な挙動を示します。これを視覚化する一般的な方法は、分岐図 を使用することです。bを 一定 (例えば 0.3) に保ち、a を変化させると、写像は規則的 (周期的) な挙動からカオス的な挙動へと遷移します。この遷移は、ロジスティック写像 と同様に、周期倍分岐カスケード によって発生します。[ 4 ]
分岐図の3Dビュー。パラメータbを 新しい軸に沿って変化させている。これにより、アトラクターが存在するパラメータ空間領域が「ブーメラン」形状をしていることがわかる。 a の値が小さい場合、システムは単一の安定な固定点に収束します。a が 増加すると、この点は不安定になり、安定な 2 サイクルに分裂します。このサイクルは不安定になり、4 サイクル、8 サイクルなどに分裂し、システムが完全カオス状態になるaの臨界値に達するまで続きます。カオス領域内には、特定の a の範囲で安定な軌道が再び現れる周期性の「ウィンドウ」も存在します。[ 6 ]
クープマン演算子解析 エノン写像の近似クープマンモード。このモードのレベルセットは、力学系の複雑な幾何学的構造を明らかにします。濃い青色の領域は、ストレンジアトラクターの吸引領域 を示し、その内部の複雑な構造は、アトラクターの安定多様体 に対応します。 エノン写像のような力学系を解析する別の方法として、クープマン演算子 法があります。このアプローチは、非線形力学に線形的な視点を提供します。位相空間における個々の点の進化を研究する代わりに、システムの作用を「観測可能な」関数g ( x , y ) の空間上で考えます。クープマン演算子Uは、観測可能な関数 g を次の時間ステップにおける値に写像する線形演算子です。( U g ) ( x n ) = g ( x n + 1 ) = g ( 1 − 1 x n 2 + y n 、 b x n ) {\displaystyle (Ug)(\mathbf {x} _{n})=g(\mathbf {x} _{n+1})=g(1-ax_{n}^{2}+y_{n},bx_{n})} 演算子U は線形ですが、無限次元の関数空間に作用します。解析の鍵は、 U φ k = λ k φ k を満たすこの演算子の固有関数 φ k と固有値 λ k を見つけることです。これらの固有関数はクープマンモード とも呼ばれ、対応する固有値にはシステムのダイナミクスに関する重要な情報が含まれています。[ 7 ]
エノン写像のようなカオス系の場合、固有関数は一般的に複雑なフラクタル関数です。これらは解析的に求めることはできず、動的モード分解(DMD) などの方法を用いて数値的に計算する必要があります。[ 8 ] クープマンモードのレベルセットは、安定多様体や不安定多様体、吸引領域などのシステムの不変構造を明らかにし、ダイナミクスの全体像を提供します。[ 9 ]
分解 15回の反復後の古典的なエノンマップ。分解によって記述される伸縮と折り畳みの動作を示している。 エノン写像は、3つのより単純な幾何学的変換のシーケンスに分解できます。これは、写像が位相空間をどのように伸縮、圧縮、折り畳むかを理解するのに役立ちます。[ 1 ] 写像T ( x , y ) = (1 - ax² + y , bx )は 、3つの関数の合成T = R ∘ C ∘ B と見なすことができます。
曲げ: y 方向の面積保存型非線形曲げ: ( x 1 、 y 1 ) = ( x 、 1 − 1 x 2 + y ) {\displaystyle (x_{1},y_{1})=(x,1-ax^{2}+y)} 収縮: x 方向の収縮: ( x 2 、 y 2 ) = ( b x 1 、 y 1 ) {\displaystyle (x_{2},y_{2})=(bx_{1},y_{1})} 反射 :直線y = x に関する反射: ( x 3 、 y 3 ) = ( y 2 、 x 2 ) {\displaystyle (x_{3},y_{3})=(y_{2},x_{2})} 最終点は( x n +1 , y n +1 ) = ( x 3 , y 3 ) です。この分解により、面積保存折り畳み動作 (ステップ 1) と散逸収縮 (ステップ 2) が分離されます。[ 4 ]
一般化
エノンの3DマップHitzlとZeleは、Hénonマップの3D一般化を提案した。[ 13 ]
s ( n + 1 ) = [ s 1 ( n + 1 ) s 2 ( n + 1 ) s 3 ( n + 1 ) ] = [ 1 − α s 1 2 ( n ) + s 3 ( n ) − β s 1 ( n ) β s 1 ( n ) + s 2 ( n ) ] {\displaystyle \mathbf {s} (n+1)={\begin{bmatrix}s_{1}(n+1)\\s_{2}(n+1)\\s_{3}(n+1)\end{bmatrix}}={\begin{bmatrix}1-\alpha s_{1}^{2}(n)+s_{3}(n)\\-\beta s_{1}(n)\\\beta s_{1}(n)+s_{2}(n)\end{bmatrix}}} 特定のパラメータ (例:α = 1.07 {\displaystyle \alpha =1.07} そしてβ = 0.3 {\displaystyle \beta =0.3} )、このマップはカオスアトラクターを生成する。[ 13 ]
四次元拡張 4次元におけるエノン地図。bの 範囲は-1.5から0.5、a の範囲は-2.3から1.0です。この動画では、点が無限遠に発散しない3次元断面図を示しています。エノンマップは、そのパラメータa とbを追加の軸として扱うことで、 4 次元空間 にプロットできます。これにより、パラメータ空間全体にわたるマップの挙動を視覚化できます。この 4 次元構造を視覚化する 1 つの方法は、一連の 3 次元スライスをレンダリングすることです。各スライスは、1 つのパラメータ (例えばa ) の固定値を表し、他の 3 つのパラメータ ( x 、y 、b ) が表示されます。次に、4 番目のパラメータを時間変数として変化させることで、変化する 3 次元構造の動画を作成できます。
参考文献 1 2 3 4 5 6 M. Hénon (1976). "奇妙なアトラクターを持つ2次元マッピング" . Communications in Mathematical Physics . 50 (1): 69– 77. Bibcode : 1976CMaPh..50...69H . doi : 10.1007/BF01608556 . S2CID 12772992 . ↑ P. Grassberger; I. Procaccia (1983). "ストレンジアトラクターの奇妙さの測定". Physica D: Nonlinear Phenomena . 9 ( 1– 2): 189– 208. Bibcode : 1983PhyD....9..189G . doi : 10.1016/0167-2789(83)90298-1 . ↑ DA Russell; JD Hanson; E. Ott (1980). "Dimension of strange attractors". Physical Review Letters . 45 (14): 1175. Bibcode : 1980PhRvL..45.1175R . doi : 10.1103/PhysRevLett.45.1175 . 1 2 3 Alligood, KT; Sauer, T.; Yorke, JA (1996). Chaos: An Introduction to Dynamical Systems . Springer. pp. 129–134 . ISBN 978-0-387-94677-1 。↑ Predrag Cvitanović; Gemunu Gunaratne; Itamar Procaccia (1988). "Hénon型ストレンジアトラクターの位相的および計量的性質". Physical Review A . 38 (3): 1503– 1520. Bibcode : 1988PhRvA..38.1503C . doi : 10.1103/PhysRevA.38.1503 . PMID 9900529 . 1 2 3 ストロガッツ、スティーブン H. (2015). 非線形力学とカオス:物理学、生物学、化学、工学への応用 (第 2 版). ウェストビュー プレス. pp. 498–501 . ISBN 978-0-8133-4910-7 。↑ Mezić, Igor (2005). "力学系のスペクトル特性、モデル縮約と分解". Nonlinear Dynamics . 41 ( 1– 3): 309– 325. doi : 10.1007/s11071-005-2824-x . S2CID 121958639 . ↑ Brunton, Steven L. ; Proctor, Joshua L.; Kutz, J. Nathan (2022). "Dynamic mode decomposition for control". Annual Review of Control, Robotics, and Autonomous Systems . 5 : 435– 462. doi : 10.1146/annurev-control-062921-090210 . ↑ コン・チャン;李海鵬。ユエヘン・ラン (2022)。 「クープマン解析による位相空間分割」。 カオス 。 32 (6): 063132. 土井 : 10.1063/5.0079812 。 PMID 35778118 。 ↑ Hénon, M.; Pomeau, Y. (1976). "単純な構造を持つ2つのストレンジアトラクター".Temam, Roger (編). Turbulence and Navier Stokes Equations . Lecture Notes in Mathematics. Vol. 565. Berlin, Heidelberg: Springer, Berlin, Heidelberg. pp. 29–68 . doi : 10.1007/BFb0091446 . ISBN 978-3-540-37516-6 。↑ Pomeau, Y.; Manneville, P. (1980). "散逸力学系における乱流への断続的遷移". Communications in Mathematical Physics . 74 (2): 189–197 . Bibcode : 1980CMaPh..74..189P . doi : 10.1007/BF01197757 . S2CID 123753342 . ↑ プティジラール、ロイク (2004)。 「混沌の誕生 (1970-1985)」 [ 混沌の誕生 (1970-1985) ] 。 theses.univ-lyon2.fr (フランス語) 。 2025年 8月6日 取得 。 1 2 Hitzl, Donald L.; Zele, Frank (1985 年 3 月). 「エノン二次写像の探求」. Physica D: Nonlinear Phenomena . 14 (3): 305– 326. doi : 10.1016/0167-2789(85)90092-2 . ↑ ボルヘス、ヴィニシウス S.。アイゼンクラフト、マルシオ(2022年12月)。 「フィルタリングされたエノンマップ」。 カオス、ソリトン、フラクタル 。 165 112865.arXiv : 2211.16964 。 土井 : 10.1016/j.chaos.2022.112865 。 S2CID 254095983 。 ↑ ボルヘス、ヴィニシウス S.。シルバ、マグノ TM;アイゼンクラフト、マルシオ(2024年4月)。 「FIR フィルターをかけたヘノン マップのカオス プロパティ」。 非線形科学と数値シミュレーションにおけるコミュニケーション 。 131 107845.arXiv : 2401.10281 。 土井 : 10.1016/j.cnsns.2024.107845 。
さらに読む Kuznetsov, Nikolay; Reitmann, Volker (2020).動的システムのアトラクター次元推定: 理論と計算 . Cham: Springer. M.ミケリッチュ。 OE レスラー (1989)。 「ヘノンの地図の新しい特徴」。コンピュータとグラフィックス 。13 (2): 263–265 。土井 : 10.1016/0097-8493(89)90070-8。 。
外部リンク ibiblio.org提供のインタラクティブなエノン地図 Wolfram Demonstrations Projectによる、エノン地図の軌道図。 CNRSによるJavaScriptを用いたエノン地図のシミュレーション。