自励系

微分方程式論または力学系理論において自励系（じれいけい、英語: Autonomous system）とは、独立変数を陽に含まない常微分方程式である。自律系（じりつけい）とも呼ぶ。逆に独立変数を陽に含む常微分方程式は、非自励系または非自律系と呼ばれる。独立変数を $t$ とし、従属変数を $x$ とすれば、自励系は

{\frac {d{\boldsymbol {x}}}{dt}}=f({\boldsymbol {x}})

で表され、非自励系は

{\frac {d{\boldsymbol {x}}}{dt}}=f(t,\ {\boldsymbol {x}})

で表される。

自励系は、その解の引数（独立変数）を一定値平行移動させたものもまた解になるという一般的性質を持つのに対し、非自励系ではこの性質は一般に成り立たない。自励系の相空間上の軌道は、他の軌道や自身と交わることはない。

定義

独立変数を $t \in ℝ$ とし、 $n$ 個の従属変数 $(x 1, x 2, \dots, x n) \in ℝ n$ に関する一般な1階連立常微分方程式を、正規形で

{\begin{cases}{\dfrac {dx_{1}}{dt}}=f_{1}(t,\ x_{1},\ x_{2},\cdots ,\ x_{n})\\{\dfrac {dx_{2}}{dt}}=f_{2}(t,\ x_{1},\ x_{2},\cdots ,\ x_{n})\\\vdots \\{\dfrac {dx_{n}}{dt}}=f_{n}(t,\ x_{1},\ x_{2},\cdots ,\ x_{n})\\\end{cases}}

と表す^[1]。これをベクトル記号でひとまとめに表すと、次のようになる。

{\frac {d{\boldsymbol {x}}}{dt}}=f(t,\ {\boldsymbol {x}})

ここで、

{\boldsymbol {x}}=(x_{1},\ x_{2},\cdots ,\ x_{n})^{\top }

{\frac {d{\boldsymbol {x}}}{dt}}=\left({\dfrac {dx_{1}}{dt}},\ {\dfrac {dx_{2}}{dt}},\cdots ,{\dfrac {dx_{n}}{dt}}\right)^{\top }

f=(f_{1},\ f_{2},\cdots ,\ f_{n})^{\top }

であり、右肩の ⊤ は転置行列を意味する^[1]^[2]。力学系分野では、独立変数 $t$ を時間とみなす^[3]。

このような微分方程式系において右辺の函数 $f$ が引数として $t$ を含まないとき、すなわち、ある微分方程式系が

{\frac {d{\boldsymbol {x}}}{dt}}=f({\boldsymbol {x}})

で与えられるとき、この微分方程式系を自励系または自律系と呼ぶ^[1]^[2]。あるいは、このような系を自励的であるという^[4]。例えば、ローレンツ方程式

{\begin{cases}{\dfrac {dx}{dt}}=-\sigma (x-y)\\{\dfrac {dy}{dt}}=-y-xz+rx\\{\dfrac {dz}{dt}}=xy-bz\end{cases}}

は $3$ 次元の自励系の例である^[5]。

逆に、微分方程式系において右辺の函数 $f$ が引数として $t$ を含むとき、すなわち、ある微分方程式系が

{\frac {d{\boldsymbol {x}}}{dt}}=f(t,\ {\boldsymbol {x}})

で与えられるとき、この微分方程式系を非自励系または非自律系と呼ぶ^[6]^[2]。あるいは、このような系を非自励的であるという^[4]。例えば、ダフィング方程式

{\begin{cases}{\dfrac {dx}{dt}}=v\\{\dfrac {dv}{dt}}=-\omega x+\epsilon x^{3}-\gamma v+B\sin \Omega t\end{cases}}

は $2$ 次元の非自励系の例である^[7]。

$1$ 階連立微分方程式に限らずに、 $n$ 階の微分方程式

f\left(t,x,{\frac {dx}{dt}},\cdots ,{\frac {d^{n}x}{dt^{n}}}\right)=0

においても、方程式に独立変数 $t$ が陽に含まれないものを自励系、 $t$ が陽に含まれるものを非自励系と呼ぶ^[2]。任意の $1$ 変数の $n$ 階微分方程式は、 $n$ 変数の $1$ 階連立微分方程式に一般的に変換できる^[8]。

自励系という言葉は、振動学の自励振動に由来するが、自励系は自励振動を起こす系をとくに意味するわけではない^[9]。振動学における自励振動は時間（独立変数） $t$ を含む強制項が存在しない方程式で記述されるため、それに関連して上記の種類の微分方程式が自励系と呼ばれるようになった^[9]。

性質

自励系で成立する基本的定理が、独立変数 $t$ を一定値ずらした解もまた解となる点である^[10]。解 $x$ を $t$ の函数として $x (t)$ と表す。自励系の微分方程式系

{\frac {d{\boldsymbol {x}}(t)}{dt}}=f({\boldsymbol {x}}(t))

は、任意の定数 $c \in ℝ$ を $t$ に加えた解 $x (t + c)$ についても、

{\frac {d{\boldsymbol {x}}(t+c)}{dt}}=f({\boldsymbol {x}}(t+c))

が満たされる^[10]。この性質は非自励系では成り立たず、 $x (t)$ がある非自励系の解であったとしても、 $x (t + c)$ は一般に解にならない^[11]。この性質が、自励系と非自励系の本質的な違いといえる^[12]。

従属変数 $(x 1, x 2, \dots, x n)$ の組でつくられる空間を、力学系分野では相空間という^[13]。相空間上に描くことができる、解を表す曲線を軌道という^[13]。自励系の軌道の接ベクトルは、与えられた微分方程式（系）のベクトル $f (x)$ と等しい^[14]。

$x (t)$ に対して $x (t + c)$ も解になるという性質から、自励系の軌道は相空間上で交わらないという性質が導かれる^[15]。別の言い方をすると、自励系の2つの軌道がある点 $x 0$ を共に通るならば、それら2つの軌道は同一の軌道である^[16]。 $x 0$ を通る軌道の形は、 $x 0$ を通る時刻 $t$ の値に無関係に決まる^[17]。また、自励系の軌道が、自身と交わることもない^[18]。非自励系にこのような一般的性質はなく、相空間上で2つの軌道が交わったり、ある軌道が異なる時刻で自身と交わることがありえる^[19]^[18]。また、非自励系の軌道の形は、初期値 $x 0$ だけでなく、初期時刻 $t 0$ の値にも依存して決まる^[20]。

もし自励系がハミルトン系であれば、各軌道に沿ってハミルトニアンは一定となる^[21]。すなわち、 $2 n$ 個の従属変数

{\boldsymbol {q}}=(q_{1},\ q_{2},\cdots ,\ q_{n})

{\boldsymbol {p}}=(p_{1},\ p_{2},\cdots ,\ p_{n})

に対して実数値関数 $H (q, p)$ を定義し、これらが正準方程式

{\dfrac {dq_{i}}{dt}}={\dfrac {\partial H({\boldsymbol {q}},\ {\boldsymbol {p}})}{\partial p_{i}}},\ {\dfrac {dp_{i}}{dt}}=-{\dfrac {\partial H({\boldsymbol {q}},\ {\boldsymbol {p}})}{\partial q_{i}}},\ (i=1,2,\cdots n)

を構成するとき、ハミルトニアン $H (q, p)$ の値は任意の解（軌道）に沿って一定である^[21]。

自励系への変換

任意の $n$ 次元の非自励系は、 $n + 1$ 番目の従属変数として $x n +1 := t$ を導入することで、 $n + 1$ 次元の自励系に機械的に変換できる^[4]^[22]。すなわち、非自励系

{\begin{cases}{\dfrac {dx_{1}}{dt}}=f_{1}(t,\ x_{1},\ x_{2},\cdots ,\ x_{n})\\{\dfrac {dx_{2}}{dt}}=f_{2}(t,\ x_{1},\ x_{2},\cdots ,\ x_{n})\\\vdots \\{\dfrac {dx_{n}}{dt}}=f_{n}(t,\ x_{1},\ x_{2},\cdots ,\ x_{n})\\\end{cases}}

において、 $x n +1 := t$ と置くことで、

{\begin{cases}{\dfrac {dx_{1}}{dt}}=f_{1}(x_{1},\ x_{2},\cdots ,\ x_{n},\ x_{n+1})\\{\dfrac {dx_{2}}{dt}}=f_{2}(x_{1},\ x_{2},\cdots ,\ x_{n},\ x_{n+1})\\\vdots \\{\dfrac {dx_{n}}{dt}}=f_{n}(x_{1},\ x_{2},\cdots ,\ x_{n},\ x_{n+1})\\{\dfrac {dx_{n+1}}{dt}}=1\\\end{cases}}

という自励系を得ることができる^[4]^[22]。独立変数としての $t$ の方を $τ$ と書き換えて、

{\begin{cases}{\dfrac {d{\boldsymbol {x}}}{d\tau }}=f(t,\ {\boldsymbol {x}})\\{\dfrac {dt}{d\tau }}=1\end{cases}}

のように表すこともある^[23]。

$(x 1, x 2, \dots, x n) \in ℝ n$ の相空間に対して、 $(x 1, x 2, \dots, x n, x n +1) \in ℝ n \times ℝ$ （あるいは $(x 1, x 2, \dots, x n, t) \in ℝ n \times ℝ$ ）で張られる $1$ 次元高い空間を特に拡大相空間と呼ぶ^[24]^[25]。非自励系をこのように自励系に変換した方が、軌道の時間依存性が無くなり、解の一意性についても見通しが良い^[22]。

非自励系は上記のように常に自励系の形に書き換え可能なため、自励系の形の方が一般性が高いといえる^[26]。しかし、非自励系には

{\dfrac {dx_{n+1}}{dt}}=1

の存在によって

\left({\dfrac {dx_{1}}{dt}},\ {\dfrac {dx_{2}}{dt}},\cdots ,{\dfrac {dx_{n}}{dt}},\ {\dfrac {dx_{n+1}}{dt}}\right)^{\top }=0

を満たす平衡点が存在しない^[27]。非自励系の軌道は、拡大相空間上で $t$ 軸方向へ常に流れ続ける^[27]。このため、自励系と非自励系では解析のアプローチを変える必要がある^[27]。

出典

^ ^a ^b ^c 小室 2005, p. 17.
^ ^a ^b ^c ^d 今隆助・竹内康博、2018、『常微分方程式とロトカ・ヴォルテラ方程式』初版、共立出版 ISBN 978-4-320-11348-0 p. 3
^ S. ウィギンス、丹羽敏雄（監訳）、今井桂子・田中茂・水谷正大・森真（訳）、2013、『非線形の力学系とカオス』新装版、丸善出版 ISBN 978-4-621-06435-1 p. 1
^ ^a ^b ^c ^d Jackson 1994, p. 16.
^ 小室 2005, pp. 46.
^ 小室 2005, p. 20.
^ 小室 2005, pp. 34–35.
^ ブラウン 2012, p. 2.
^ ^a ^b 齋藤 2004, p. 10.
^ ^a ^b 齋藤 2004, pp. 11–12.
^ 上田 2008, pp. 76–77.
^ 上田 2008, p. 77.
^ ^a ^b 松葉 2011, p. 7.
^ 高橋陽一郎、2004、『力学と微分方程式』、岩波書店〈現代数学への入門〉 ISBN 4-00-006875-X p. 5
^ 松葉 2011, pp. 25–26.
^ ブラウン 2012, p. 190.
^ 齋藤 2004, p. 14.
^ ^a ^b 船越満明、2008、『カオス』初版、朝倉書店〈シリーズ非線形科学入門3〉 ISBN 978-4-254-11613-7 pp. 135–136
^ Jackson 1994, p. 17.
^ 上田 2008, pp. 16, 77.
^ ^a ^b 柴山允瑠、2016、『重点解説ハミルトン力学系』初版、サイエンス社〈SGCライブラリ 130〉 pp. 25–26
^ ^a ^b ^c 松葉 2011, p. 24.
^ 桑村雅隆、2015、『パターン形成と分岐理論―自発的パターン発生の力学系入門』初版、共立出版〈シリーズ・現象を解明する数学〉 ISBN 978-4-320-11004-5 p. 28
^ Jackson 1994, p. 277.
^ 小室 2005, pp. 12, 20.
^ K. T. アリグッド・T. D. サウアー・J. A. ヨーク、津田一郎（監訳）、星野高志・阿部巨仁・黒田拓・松本和宏（訳）、2012、『カオス　第2巻　力学系入門』、丸善出版 ISBN 978-4-621-06279-1 p. 89
^ ^a ^b ^c 川上博 (2005年). “非線形現象入門”. 徳島大学教育・研究者情報データベース. 2020年5月30日閲覧。 pp. 11–12

参照文献

小室元政、2005、『基礎からの力学系―分岐解析からカオス的遍歴へ』新版、サイエンス社 ISBN 4-7819-1118-8
齋藤利弥、2004、『力学系入門』復刊版、朝倉書店 ISBN 4-254-11722-1
松葉育雄、2011、『力学系カオス』第1版、森北出版 ISBN 978-4-627-15451-3
上田睆亮、2008、『カオス現象論』初版、コロナ社〈現代非線形科学シリーズ12〉 ISBN 978-4-339-02611-5
E. Atlee Jackson、田中茂・丹羽敏雄・水谷正大・森真（訳）、1994、『非線形力学の展望I ―カオスとゆらぎ』初版、共立出版 ISBN 4-320-03325-6
M. ブラウン、一樂重雄・河原正治・河原雅子・一樂祥子（訳）、2012、『微分方程式下 ―その数学と応用』、丸善出版 ISBN 978-4-621-06570-9

外部リンク

Hazewinkel, Michiel, ed. (2001), “Autonomous system”, Encyclopedia of Mathematics, Springer, ISBN 978-1-55608-010-4
Weisstein, Eric W. "Autonomous". mathworld.wolfram.com (英語).

[FOOTNOTE小室200517-1] 小室 2005, p. 17.

[今・竹内-2] 今隆助・竹内康博、2018、『常微分方程式とロトカ・ヴォルテラ方程式』初版、共立出版 ISBN 978-4-320-11348-0 p. 3

[3] S. ウィギンス、丹羽敏雄（監訳）、今井桂子・田中茂・水谷正大・森真（訳）、2013、『非線形の力学系とカオス』新装版、丸善出版 ISBN 978-4-621-06435-1 p. 1

[FOOTNOTEJackson199416-4] Jackson 1994, p. 16.

[FOOTNOTE小室200546-5] 小室 2005, pp. 46.

[FOOTNOTE小室200520-6] 小室 2005, p. 20.

[FOOTNOTE小室200534–35-7] 小室 2005, pp. 34–35.

[FOOTNOTEブラウン20122-8] ブラウン 2012, p. 2.

[FOOTNOTE齋藤200410-9] 齋藤 2004, p. 10.

[FOOTNOTE齋藤200411–12-10] 齋藤 2004, pp. 11–12.

[FOOTNOTE上田200876–77-11] 上田 2008, pp. 76–77.

[FOOTNOTE上田200877-12] 上田 2008, p. 77.

[FOOTNOTE松葉20117-13] 松葉 2011, p. 7.

[高橋-14] 高橋陽一郎、2004、『力学と微分方程式』、岩波書店〈現代数学への入門〉 ISBN 4-00-006875-X p. 5

[FOOTNOTE松葉201125–26-15] 松葉 2011, pp. 25–26.

[FOOTNOTEブラウン2012190-16] ブラウン 2012, p. 190.

[FOOTNOTE齋藤200414-17] 齋藤 2004, p. 14.

[船越-18] 船越満明、2008、『カオス』初版、朝倉書店〈シリーズ非線形科学入門3〉 ISBN 978-4-254-11613-7 pp. 135–136

[FOOTNOTEJackson199417-19] Jackson 1994, p. 17.

[FOOTNOTE上田200816,_77-20] 上田 2008, pp. 16, 77.

[柴山-21] 柴山允瑠、2016、『重点解説ハミルトン力学系』初版、サイエンス社〈SGCライブラリ 130〉 pp. 25–26

[FOOTNOTE松葉201124-22] 松葉 2011, p. 24.

[23] 桑村雅隆、2015、『パターン形成と分岐理論―自発的パターン発生の力学系入門』初版、共立出版〈シリーズ・現象を解明する数学〉 ISBN 978-4-320-11004-5 p. 28

[FOOTNOTEJackson1994277-24] Jackson 1994, p. 277.

[FOOTNOTE小室200512,_20-25] 小室 2005, pp. 12, 20.

[26] K. T. アリグッド・T. D. サウアー・J. A. ヨーク、津田一郎（監訳）、星野高志・阿部巨仁・黒田拓・松本和宏（訳）、2012、『カオス　第2巻　力学系入門』、丸善出版 ISBN 978-4-621-06279-1 p. 89

[川上-27] 川上博 (2005年). “非線形現象入門”. 徳島大学教育・研究者情報データベース. 2020年5月30日閲覧。 pp. 11–12

[1]

[2]

[3]

[4]

[5]

[6]

[7]

[8]

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16]

[17]

[18]

[19]

[20]

[21]

[22]

[23]

[24]

[25]

[26]

[27]