軌道 (力学系)

力学系における軌道（きどう）とは、初期条件に対して時間発展のルールを適用したときに定まる、相空間上の点の集合である。連続的な時間を仮定した系だと、軌道は相空間内で一本の曲線となり、離散的な時間を仮定した系だと、軌道は相空間内で点列となる。 力学系を定める相空間を X、時間を G、時間発展のルールを ϕ: G × X → X とする。ある t ∈ G に固定したときの ϕ を写像 ϕt と表し、X ∋ x ↦ ϕt(x) ∈ X である。G は結合法則 t1 + t2 (t1, t2 ∈ T) で表される群構造を持ち、ϕt は、 ϕ e = i d {\displaystyle \phi ^{e}=\mathrm {id} } ϕ t 1 ∘ ϕ t 2 = ϕ t 1 + t 2 {\displaystyle \phi ^{t_{1}}\circ \phi ^{t_{2}}=\phi ^{t_{1}+t_{2}}} という性質を満たす。ここで e は G の単位元、id は恒等写像、∘ は写像の合成を意味する。 このような力学系 X, T, ϕt において O ( x 0 ) = { x ∈ X ∣ ϕ t ( x 0 ) , t ∈ T } {\displaystyle O(x_{0})=\{x\in X\mid \phi _{t}(x_{0}),t\in T\}} で定義される X の順序部分集合 O(x0) を軌道（英: orbit、英: trajectory）と呼ぶ。ただし、t が取り得る値は ϕt (x0) が定義されている範囲に限られる。O(x0) は「x0 を通る軌道」と呼ばれる。軌道の記号には、O(x0 の他に、𝒪(x0) 、C(x0) 、γ(x0) 、Γ(x0) 、Orb(x0) などの表記がある。 群論の言葉では軌道は次のように定義される。上記を満たす写像 ϕ を、群 G の集合 X への作用という。この作用 ϕ について、X 上の2点 x0, y が適当な t を選びさえすれば ϕt (x0) = y という関係を満たすとき、x0, y は同値関係 x0 ∼ y にあると定義する。この同値関係によって X を分ける同値類が軌道 O(x0) である。 時間 G が整数 ℤ のときの力学系を離散力学系と呼び、G が実数 ℝ のときを連続力学系と呼ぶ。相空間上のどの点も初期値となりうるので、相空間は何かしらの軌道によって完全に埋め尽くされる。力学系理論の主目的は、系の軌道の性質・振る舞いを調べることにある。特に力学系理論の場合、時間が正または負の無限大に発散するときの漸近的振る舞いを問題とする。軌道同士の相互関係や、系に摂動が加わったときに起こる軌道全体の構造の変化なども力学系理論の題目である。 離散力学系は写像の反復によって定義される。相空間上のある点 x0 ∈ M に写像 f: M → M を繰り返し適用することで、x0, f (x0), f 2(x0), … f n(x), … という点列が得られる。点列は x0, x1 = f (x0), x2 = f 2(x0), … xn = f n(x0), … とも表す。この点列が離散力学系の軌道である。多くの力学系で f は連続写像である。 例えば、ℝ 上の正弦関数 f(x) = sin(x) で定義される離散力学系を考える。x0 = 123 とすると、 x 0 = 123 x 1 = − 0.4599 … x 2 = − 0.4438 … ⋮ x 300 = − 0.0975 … x 301 = − 0.0974 … ⋮ {\displaystyle {\begin{aligned}x_{0}&=123\\x_{1}&=-0.4599\ldots \\x_{2}&=-0.4438\ldots \\\vdots \\x_{300}&=-0.0975\ldots \\x_{301}&=-0.0974\ldots \\\vdots \end{aligned}}} というような数列がその軌道である。 細かく分けると、点列 x0, f (x0), f 2(x0), … は特に前方軌道や正の半軌道と呼ばれ、O+(x0) や O+(x0) のように表す。 O + ( x 0 ) = { f n ( x 0 ) ∣ n = 0 ,   1 , 2 , … } {\displaystyle O_{+}(x_{0})=\{f^{n}(x_{0})\mid n=0,\ 1,\,2,\ldots \}} 一方、f が可逆で逆写像 f −1 を持つとき、f 0 は恒等写像だとして、k < 0 についても写像の反復 f k が定義できる。それによって、x0, f −1(x0), f −2(x0), … という点列が定義でき、O−(x0) などのように表す。 O − ( x 0 ) = { f n ( x 0 ) ∣ n = 0 , − 1 , − 2 , … } {\displaystyle O_{-}(x_{0})=\{f^{n}(x_{0})\mid n=0,-1,-2,\ldots \}} 逆写像によって定まる点列 O−(x0) は後方軌道や負の半軌道と呼ばれる。正の半軌道と負の半軌道を足し合わせた集合 O ( x 0 ) = { f n ( x 0 ) ∣ n ∈ Z } = { … ,   f − 2 ( x 0 ) ,   f − 1 ( x 0 ) ,   x 0 ,   f ( x 0 ) ,   f 2 ( x 0 ) ,   … } = { … ,   x − 2 ,   x − 1 ,   x 0 ,   x 1 ,   x 2 ,   … } {\displaystyle {\begin{aligned}O(x_{0})&=\{f^{n}(x_{0})\mid n\in \mathbb {Z} \}\\&=\{\ldots ,\ f^{-2}(x_{0}),\ f^{-1}(x_{0}),\ x_{0},\ f(x_{0}),\ f^{2}(x_{0}),\ \ldots \}\\&=\{\ldots ,\ x_{-2},\ x_{-1},\ x_{0},\ x_{1},\ x_{2},\ \ldots \}\\\end{aligned}}} を軌道や全軌道と呼ぶ。 連続力学系を定義する一番普通の方法は、微分方程式による定義である。相空間 X がユークリッド空間か多様体だとする。未知関数 x(t) ∈ X の常微分方程式系 d x ( t ) d t = V ( x ( t ) ) {\displaystyle {\frac {dx(t)}{dt}}=V(x(t))} を考える。この微分方程式が初期条件 x(0) = x0 を満たす解を x(t, x0) と表す。微分方程式を決めている関数 V(x) は X 上にベクトル場を与える。この解 x(t, x0) は上の節で一般的に定義した写像 ϕt (x0) と等しい。 微分方程式の解が存在する t の領域を I ⊂ ℝ とする。連続力学系の軌道とは、 O ( x 0 ) = { x ( t ,   x 0 ) ∣   t ∈ I } {\displaystyle O(x_{0})=\{x(t,\ x_{0})\mid \ t\in I\}} で定義される集合である。ただし、O(x0) には t が小さい方から大きい方に向かって向きが付いている。 簡単のために I = (−∞, ∞) だと仮定すれば、連続力学系の正の半軌道は、 O + ( x 0 ) = { x ( t ,   x 0 ) ∣ 0 ≤ t < ∞ } {\displaystyle O_{+}(x_{0})=\{x(t,\ x_{0})\mid 0\leq t<\infty \}} で定義され、負の半軌道は、 O − ( x 0 ) = { x ( t ,   x 0 ) ∣ − ∞ < t ≤ 0 } {\displaystyle O_{-}(x_{0})=\{x(t,\ x_{0})\mid -\infty <t\leq 0\}} で定義される。正の半軌道と負の半軌道を足し合わせた集合 O ( x 0 ) = { x ( t ,   x 0 ) ∣ t ∈ R } {\displaystyle O(x_{0})=\{x(t,\ x_{0})\mid t\in \mathbb {R} \}} を離散力学系と同様に軌道や全軌道と呼ぶ。 連続力学系の x0 を通る軌道は、相空間上の x0 を通る一つの曲線に対応する。この曲線を微分方程式の解曲線とも呼ぶ。軌道（解曲線）上の各点 x にはベクトル場のベクトル V(x) が存在し、軌道に接している。解曲線のことを解軌道という風に呼ぶこともある。微分方程式を満たす解 x(t) を指して解軌道や軌道と呼ぶこともある。 もっとも単純な軌道としては、離散力学系の不動点と連続力学系の平衡点がある。これら2つを共に「不動点」と呼んだり、「平衡点」と呼ぶこともある。以下では区別して記す。 離散力学系の不動点とは、写像 f を適用しても動かない点のことで、f (x0) = x0 を満たす点 x0 である。不動点 x0 での軌道は O(x0) = {x0, x0, x0,…} という定数の列となる。連続力学系の平衡点とは、時間が経っても動かない点である。平衡点 x0 での軌道は O(x0) = {x0} となる。微分方程式で定まる系の場合、定常解 x(t) ≡ x0 のことで、微分方程式の右辺（ベクトル場）V(x0) = 0 を満たす点 x0 が平衡点である。 ひとまとめされることがあるように不動点も平衡点も同じ性質のものだといえる。一般化された定義を与えると、時間発展のルール ϕt が任意の t ∈ G について ϕt(x0) = x0 を満たすときの x0 が不動点・平衡点である。不動点・平衡点を調べることは、一般の軌道を調べるよりも総じて容易であり、与えられた力学系を理解するための重要な手がかりとなる。 もう一つの比較的単純な軌道が周期軌道である。 離散力学系で、非零のある自然数 k > 0 について f k(x0) = x0 を満たす x0 を周期点と呼ぶ。条件を満たす最小の k を周期や最小周期と呼ぶ。そして、ある周期点を通る軌道を周期軌道と呼び、軌道は周期的であるという。周期 k の軌道だと O ( x 0 ) = { x 0 ,   f ( x 0 ) ,   f 2 ( x 0 ) , ⋯ ,   f k − 1 ( x 0 ) ,   x 0 ,   f ( x 0 ) ,   f 2 ( x 0 ) , ⋯ } {\displaystyle O(x_{0})=\{x_{0},\ f(x_{0}),\ f^{2}(x_{0}),\cdots ,\ f^{k-1}(x_{0}),\ x_{0},\ f(x_{0}),\ f^{2}(x_{0}),\cdots \}} のようになる。周期軌道の各点は全て同じ周期の周期点である。 連続力学系の場合、非零のある実数 T > 0 と任意の t について微分方程式の解が x(t + T, x0) = x(t, x0) を満たすとき、解を周期解と呼ぶ。条件を満たす最小の T を周期や最小周期と呼ぶ。このような解の軌道、すなわち集合 O ( x 0 ) = { x ( t ,   x 0 ) ∣ 0 ≤ t ≤ T } {\displaystyle O(x_{0})=\{x(t,\ x_{0})\mid 0\leq t\leq T\}} が連続力学系の周期軌道である。連続力学系の周期軌道は相空間上で閉曲線となり、そのため閉軌道とも呼ばれる。 相空間がトーラスになると、準周期軌道という種類の軌道が存在し得る。2次元トーラス 𝕋2 上の d θ 1 d t = ω 1 {\displaystyle {\frac {d\theta _{1}}{dt}}=\omega _{1}} d θ 2 d t = ω 2 {\displaystyle {\frac {d\theta _{2}}{dt}}=\omega _{2}} という微分方程式を考える。トーラスは 2π を法として得られる商集合 𝕋2 = ℝ2/2π ℤ2 と見なし、(θ1, θ2) ∈ 𝕋2 である。 比 ω2/ω1 が有理数のとき、この連続力学系の軌道はトーラス上で周期軌道となる。一方、ω2/ω1 が無理数のとき、任意の解は 𝕋2 上を稠密に埋めつくす。後者のような解を準周期解、軌道を準周期的であるあるいは準周期軌道と呼ぶ。軌道が準周期的なとき、軌道は閉じることも自己交差することもなく、トーラスに永久に巻きつきながら、トーラス上を軌道で埋め尽くす。一般の n 次元トーラス 𝕋n についても同種のことが成り立つ。 𝕋2 上の準周期軌道をポアンカレ写像によって離散力学系に落とし込むと、ポアンカレ断面でトーラスを切り取った格好となるので、準周期軌道は断面上で閉じた曲線として反映される。θ1 = 0 でポアンカレ写像を構成すると、 f ( θ 2 ) = θ 2 + 2 π ( ω 2 / ω 1 ) ( mod 2 π ) {\displaystyle f(\theta _{2})=\theta _{2}+2\pi (\omega _{2}/\omega _{1}){\pmod {2\pi }}} となり、円周上の点を角度 2π(ω2/ω1) ずつ動かす写像になる。ω2/ω1 が無理数のとき、この写像の軌道は円周を稠密に埋め尽くす。離散力学系のこのような軌道も準周期的、準周期軌道と呼ばれることもある。 ^ Kuznetsov 1998, p. 6. ^ 白石 2014, p. 166. ^ a b Kuznetsov 1998, p. 8. ^ a b ウィギンス 2013, p. 2. ^ a b c 國府 2000, p. 7. ^ 青木・白岩 2013, pp. 16–17. ^ a b c 坂井 2015, p. xiv. ^ a b c d e f Devaney 2003, p. 16. ^ a b 青木 1996, p. 6. ^ a b c d e 松葉 2011, p. 28. ^ 荒井 2020, p. 25. ^ 齋藤 2002, p. 9. ^ 郡・森田 2011, p. 42. ^ a b 今・竹内 2018, p. 145. ^ Jackson 1994, p. 44. ^ a b 白石 2014, p. 167. ^ 白石 2014, pp. 166–167. ^ Kuznetsov 1998, p. 5. ^ Strogatz 2015, p. 8. ^ a b 久保・矢野 2018, p. 166. ^ a b c 青木・白岩 2013, p. 17. ^ 荒井 2020, p. 164. ^ 國府 2000, p. 12. ^ 白石 2014, p. 176. ^ a b c デバニー 2007, p. 20. ^ a b Kuznetsov 1998, p. 7. ^ 小室 2005, p. 23. ^ 川上 1990, p. 139. ^ 川上 1990, pp. 139–140. ^ 青木・白岩 2013, p. 6. ^ Kuznetsov 1998, p. 18. ^ 久保・矢野 2018, p. 28. ^ Kuznetsov 1998, p. 21. ^ ウィギンス 2013, pp. 1–2. ^ 今・竹内 2018. ^ a b c 郡・森田 2011, p. 52. ^ a b 柴山 2016, p. 1. ^ 齋藤 2002, p. 8. ^ 小室 2005, p. 7. ^ 丹羽 2004, p. 36. ^ 船越 2008, p. 78. ^ 今・竹内 2018, p. 107. ^ 千葉 2021, p. 100. ^ 高橋 2004, p. 5. ^ 小室 2005, p. 18. ^ 川上 1990, pp. 140–141. ^ a b 柴山 2016, pp. 2–3. ^ a b c Kuznetsov 1998, p. 9. ^ ウィギンス 2013, pp. 6, 9. ^ 千葉 2021, pp. 111, 208. ^ 今・竹内 2018, pp. 108, 216. ^ a b c 松葉 2011, p. 32. ^ a b デバニー 2007, p. 21. ^ 森・水谷 2009, p. 29. ^ 坂井 2015, p. xv. ^ a b 國府 2000, p. 8. ^ a b デバニー 2007, p. 22. ^ 白石 2014, p. 177. ^ ウィギンス 2013, p. 27. ^ デバニー 2007, p. 23. ^ a b 今・竹内 2018, p. 110. ^ a b 伊藤 1998, p. 26. ^ a b 荒井 2020, p. 37. ^ 郡・森田 2011, p. 47. ^ 齋藤 2002, pp. 11–12. ^ a b c Strogatz 2015, p. 303. ^ 伊藤 1998, p. 138. ^ a b 伊藤 1998, p. 139. ^ a b 伊藤 1998, p. 140. ^ 坂井 2015, p. 299. ^ a b c 井上・秦 1999, p. 72. ^ 井上・秦 1999, p. 73. ^ a b Hirsch, Smale & Devaney 2007, p. 120. ^ アリグッド、サウアー、ヨーク 2012, p. 122. ^ ウィギンス 2013, p. 154. 白石 謙一、2014、『力学系の理論』オンデマンド版、岩波書店 ISBN 978-4-00-730152-0 Yuri A. Kuznetsov (1998). Elements of Applied Bifurcation Theory. Applied mathematical sciences Vol. 112 (Second ed.). Springer. ISBN 0-387-98382-1. https://link.springer.com/book/10.1007/b98848  國府 寛司、2000、『力学系の基礎』初版、朝倉書店〈カオス全書2〉 ISBN 4-254-12672-7 荒井 迅、2020、『常微分方程式の解法』初版、共立出版〈共立講座 数学探検 15〉 ISBN 978-4-320-11188-2 青木 統夫・白岩 謙一、2013、『力学系とエントロピー』復刊、共立出版 ISBN 978-4-320-11043-4 坂井 秀隆、2015、『常微分方程式』初版、東京大学出版会〈大学数学の入門 10〉 ISBN 978-4-13-062960-7 齋藤 利弥、2002、『位相力学 ―常微分方程式の定性的理論―』復刊、共立出版 ISBN 4-320-01712-9 郡 宏・森田 善久、2011、『生物リズムと力学系』初版、共立出版〈シリーズ・現象を解明する数学〉 ISBN 978-4-320-11000-7 S. ウィギンス、丹羽 敏雄（監訳）、今井 桂子・田中 茂・水谷 正大・森 真（訳）、2013、『非線形の力学系とカオス』新装版、丸善出版 ISBN 978-4-621-06435-1 Robert L. Devaney、國府 寛司・石井 豊 ・新居 俊作・木坂 正史（新訂版訳）、後藤 憲一（訳）、2003、『カオス力学系入門』新訂版、共立出版 ISBN 4-320-01705-6 青木 統夫、1996、『力学系・カオス ―非線形現象の幾何学的構成―』初版、共立出版 ISBN 4-320-03340-X E. Atlee Jackson、田中 茂・丹羽 敏雄・水谷 正大・森 真（訳）、1994、『非線形力学の展望Ⅰ ―カオスとゆらぎ―』初版、共立出版 ISBN 4-320-03325-6 松葉 育雄、2011、『力学系カオス』第1版、森北出版 ISBN 978-4-627-15451-3 久保 泉・矢野 公一、2018、『力学系』オンデマンド版、岩波書店 ISBN 978-4-00-730742-3 Steven H. Strogatz、田中 久陽・中尾 裕也・千葉 逸人（訳）、2015、『ストロガッツ　非線形ダイナミクスとカオス ―数学的基礎から物理・生物・化学・工学への応用まで―』、丸善出版 ISBN 978-4-621-08580-6 ロバート・L・デバニー、上江洌 達也・重本 和泰・久保 博嗣・田崎 秀一（訳）、2007、『カオス力学系の基礎』新装版、ピアソン・エデュケーション ISBN 978-4-89471-028-3 小室 元政、2005、『基礎からの力学系 ―分岐解析からカオス的遍歴へ―』新版、サイエンス社 ISBN 4-7819-1118-8 川上 博、1990、『カオスCGコレクション』初版、サイエンス社〈Information & Computing 48〉 ISBN 978-4-7819-0591-4 今 隆助・竹内 康博、2018、『常微分方程式とロトカ・ヴォルテラ方程式』初版、共立出版 ISBN 978-4-320-11348-0 柴山 允瑠、2016、「重点解説 ハミルトン力学系 ―可積分系とKAM理論を中心に―」、『臨時別冊・数理科学2016年12月』（SGCライブラリ 130）、サイエンス社、ISSN 0386-8257 船越 満明、2008、『カオス』初版、朝倉書店〈シリーズ 非線形科学入門3〉 ISBN 978-4-254-11613-7 丹羽 敏雄、2004、『微分方程式と力学系の理論入門 ―非線形現象の解析にむけて―』増補版、遊星社 ISBN 4-7952-6900-9 高橋 陽一郎、2004、『力学と微分方程式』初版、岩波書店〈現代数学への入門〉 ISBN 4-00-006875-X 千葉 逸人、2021、『解くための微分方程式と力学系理論』初版、現代数学社 ISBN 978-4-7687-0570-4 森 真・水谷 正大、2009、『入門力学系 ―自然の振舞いを数学で読みとく―』、東京図書 ISBN 978-4-489-02050-6 伊藤 秀一、1998、『常微分方程式と解析力学』初版、共立出版〈共立講座 21世紀の数学 11〉 ISBN 4-320-01563-0 井上 政義・秦 浩起、1999、『カオス科学の基礎と展開 ―複雑系の理解に向けて』初版、共立出版 ISBN 4-320-03323-X Morris W. Hirsch; Stephen Smale; Robert L. Devaney、桐木 紳・三波 篤朗・谷川 清隆・辻井 正人（訳）、2007、『力学系入門　原著第2版 ―微分方程式からカオスまで』初版、共立出版 ISBN 978-4-320-01847-1 K.T.アリグッド；T.D.サウアー；J.A.ヨーク、津田 一郎（監訳）、星野 高志・阿部 巨仁・黒田 拓・松本 和宏（訳）、2012、『カオス　第1巻　力学系入門』、丸善出版 ISBN 978-4-621-06223-4