logo

自律システムにおける固定点の分類 📂動力学

自律システムにおける固定点の分類

定義

空間 XX と関数 fC1(X,X)f \in C^{1}(X,X) に対して、以下のようなベクトル場微分方程式として与えられているとしよう。 x˙=f(x) \dot{x} = f(x) x\overline{x} がこの自律系の一つの固定点であり、Df(x)D f \left( \overline{x} \right)固有値λ1,,λm\lambda_{1} , \cdots , \lambda_{m} とする。

ハイパーボリック:双曲固定点1

  1. ハイパーボリック:Df(x)D f \left( \overline{x} \right) の全ての固有値の実部が 00 でなければ、x\overline{x}ハイパーボリックだと言われる。 Re(λ1)0,,Re(λm)0 \operatorname{Re} \left( \lambda_{1} \right) \ne 0 , \cdots , \operatorname{Re} \left( \lambda_{m} \right) \ne 0
    1. サドル:x\overline{x} がハイパーボリックであり、Df(x)D f \left( \overline{x} \right) が実部が正の固有値と負の固有値を少なくとも一つずつ持つ場合、x\overline{x}サドルだという。 i,j[1,m]:Re(λi)>0Re(λj)<0 \exists i, j \in [1,m] : \operatorname{Re} \left( \lambda_{i} \right) > 0 \land \operatorname{Re} \left( \lambda_{j} \right) < 0
    2. シンク:Df(x)D f \left( \overline{x} \right) の全ての固有値の実部が負ならば、x\overline{x}安定stableだと言われ、シンクという。 Re(λ1)<0,,Re(λm)<0 \operatorname{Re} \left( \lambda_{1} \right) < 0 , \cdots , \operatorname{Re} \left( \lambda_{m} \right) < 0
    3. ソース:Df(x)D f \left( \overline{x} \right) の全ての固有値の実部が正ならば、x\overline{x}不安定untableだと言われ、ソースという。 Re(λ1)>0,,Re(λm)>0 \operatorname{Re} \left( \lambda_{1} \right) > 0 , \cdots , \operatorname{Re} \left( \lambda_{m} \right) > 0

エリプティック:楕円固定点2

  1. エリプティック、センター:Df(x)D f \left( \overline{x} \right) の全ての固有値が純虚数ならば、x\overline{x}エリプティックだと言われ、センターという。 Im(λ1)=λ1,,Im(λm)=λm \operatorname{Im} \left( \lambda_{1} \right) = \lambda_{1} , \cdots , \operatorname{Im} \left( \lambda_{m} \right) = \lambda_{m}

  • Re\reIm\im はそれぞれ複素数から実部と虚部を取り出す関数である。

説明

定義ではそう言われていないけれども、ハイパーボリックかどうかは、端的に言って、システムが単純かどうかとほぼ一致する。システムを動力学的に分析する場合、主な問題はほとんどの場合、固有値が 00 である場合に生じる。ハイパーボリックであれば、そのような厄介ごとを心配せずに済むため、分析も単純になる。

例として、ダフィング振り子を考えてみよう: x˙=yy˙=xx3δy,δ0 \begin{align*} \dot{x} =& y \\ \dot{y} =& x - x^{3} - \delta y \qquad , \delta \ge 0 \end{align*} ダフィング振り子の固定点は (x,y)=(0,0),(±1,0)(x,y) = (0,0) , (\pm 1 , 0) であり、ヤコビアンDf=[0113x2δ] D \mathbf{f} = \begin{bmatrix} 0 & 1 \\ 1 - 3 x^{2} & - \delta \end{bmatrix} であるため、固定点でのヤコビアンは Df(0,0)=[011δ]Df(±1,0)=[012δ] D \mathbf{f} (0,0) = \begin{bmatrix} 0 & 1 \\ 1 & - \delta \end{bmatrix} \\ D \mathbf{f} (\pm1,0) = \begin{bmatrix} 0 & 1 \\ -2 & - \delta \end{bmatrix} である。これらの固有値を計算すると、(0,0)(0,0) の時 det(Df(0,0)λE)=det[λ11λδ]=λ2+δλ1 \det ( D \mathbf{f} (0,0) - \lambda E ) = \det \begin{bmatrix} -\lambda & 1 \\ 1 & -\lambda - \delta \end{bmatrix} = \lambda^{2} + \delta \lambda - 1 であり、二次方程式の解によって λ1,2=δ±δ2+42 \lambda_{1,2} = {{ -\delta \pm \sqrt{\delta^{2} + 4 } } \over { 2 }} である。Df(0,0)D \mathbf{f} (0,0)固有値は常に正数と負数の一つずつを持つため、固定点 (0,0)(0,0) はサドルである。同様に、Df(±1,0)D \mathbf{f} (\pm1,0) の固有値を計算すると λ1,2=δ±δ282 \lambda_{1,2} = {{ -\delta \pm \sqrt{\delta^{2} - 8 } } \over { 2 }} であるため、Df(±1,0)D \mathbf{f} (\pm1,0)固有値δ>0\delta > 0 の時、全部負であり、固定点 (±1,0)(\pm1,0) はシンクである。しかし、δ=0\delta = 0 の時は純虚数 λ1,2=±2i\lambda_{1,2} = \pm \sqrt{2} i であるため、固定点 (±1,0)(\pm1,0) はセンターである。

上記の例では、各固定点でヤコビアンを求め、パラメーター δ\delta の設定によって安定性がどのように変化するかを見た。このような分析は、システムがベクトル場で表されていれば、動力学のどの論文でも同じように使われる方法である。少なくとも一度は、必ず自分で試してみよう。


  1. Wiggins. (2003). Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition(2nd Edition): p12. ↩︎

  2. Wiggins. (2003). Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition(2nd Edition): p12. ↩︎