자율 시스템에서 고정점의 분류

자율 시스템에서 고정점의 분류

Classification of fixed point

정의

공간 $X$ 와 함수 $f \in C^{1}(X,X)$ 에 대해 다음과 같은 벡터 필드미분 방정식으로 주어져 있다고 하자. $$ x' = f(x) $$ $\overline{x}$ 가 이 자율 시스템의 한 고정점이라 하고 $D f \left( \overline{x} \right)$ 의 아이겐 밸류들을 $\lambda_{1} , \cdots , \lambda_{m}$ 이라 하자.

하이퍼볼릭: 쌍곡 고정점1

  1. Hyperbolic: $D f \left( \overline{x} \right)$ 의 모든 아이겐 밸류들의 실수부가 $0$ 이 아니면 $\overline{x}$ 가 하이퍼볼릭하다고 말한다. $$ \text{Re} \left( \lambda_{1} \right) \ne 0 , \cdots , \text{Re} \left( \lambda_{m} \right) \ne 0 $$
    1. Saddle: $\overline{x}$ 가 하이퍼볼릭하고 $D f \left( \overline{x} \right)$ 가 실수부가 양수인 아이겐 밸류와 음수인 아이겐 밸류를 적어도 하나씩 가지면 $\overline{x}$ 가 새들이라고 한다. $$ \exists i, j \in [1,m] : \text{Re} \left( \lambda_{i} \right) > 0 \land \text{Re} \left( \lambda_{j} \right) < 0 $$
    2. Sink : $D f \left( \overline{x} \right)$ 의 모든 아이겐 밸류들의 실수부가 음수면 $\overline{x}$ 가 스테이블Stable하다고 말하고 싱크라고 한다. $$ \text{Re} \left( \lambda_{1} \right) > 0 , \cdots , \text{Re} \left( \lambda_{m} \right) > 0 $$
    3. Source : $D f \left( \overline{x} \right)$ 의 모든 아이겐 밸류들의 실수부가 양수면 $\overline{x}$ 가 언스테이블Untable하고 말하고 소스라고 한다. $$ \text{Re} \left( \lambda_{1} \right) < 0 , \cdots , \text{Re} \left( \lambda_{m} \right) < 0 $$

일립틱: 타원 고정점2

  1. Elliptic, Center : $D f \left( \overline{x} \right)$ 의 모든 아이겐 밸류들이 순허수면 $\overline{x}$ 가 일립틱하다고 말하고 센터라고 한다. $$ \text{Im} \left( \lambda_{1} \right) = \lambda_{1} , \cdots , \text{Im} \left( \lambda_{m} \right) = \lambda_{m} $$

  • $\Re$ 와 $\Im$ 는 각각 복소수에서 실수부와 허수부만을 취하는 함수들이다.

설명

정의에서는 그렇게 말하지 않지만, 하이퍼볼릭이냐 아니냐는 쉽게 말해 시스템이 간단하냐 간단하지 않냐와 거의 상통한다. 시스템을 동역학적으로 분석함에 있어서 문제가 되는 경우는 아이겐 밸류가 $0$ 인 경우가 대부분으로, 하이퍼볼릭이면 그런 골칫거리를 생각하지 않아도 되기 때문에 분석도 간단해진다.

예시

예로써 더핑 오실레이터를 생각해보자: $$ \begin{align*} x’ =& y \\ y’ =& x - x^{3} - \delta y \qquad , \delta \ge 0 \end{align*} $$ 더핑 오실레이터의 고정점은 $(x,y) = (0,0) , (\pm 1 , 0)$ 이고 자코비안은 $$ D \mathbb{f} = \begin{bmatrix} 0 & 1 \\ 1 - 3 x^{2} & - \delta \end{bmatrix} $$ 이므로 고정점의 자코비안은 $$ D \mathbb{f} (0,0) = \begin{bmatrix} 0 & 1 \\ 1 & - \delta \end{bmatrix} \\ D \mathbb{f} (\pm1,0) = \begin{bmatrix} 0 & 1 \\ -2 & - \delta \end{bmatrix} $$ 이다. 이들의 아이겐 밸류들을 계산해보면 $(0,0)$ 일 때 $$ \det ( D \mathbb{f} (0,0) - \lambda E ) = \det \begin{bmatrix} -\lambda & 1 \\ 1 & -\lambda - \delta \end{bmatrix} = \lambda^{2} + \delta \lambda - 1 $$ 이므로 근의 공식에 따라 $$ \lambda_{1,2} = {{ -\delta \pm \sqrt{\delta^{2} + 4 } } \over { 2 }} $$ 이다. $D \mathbb{f} (0,0)$ 의 아이겐 밸류는 $\delta \ge 0$ 에 대해 항상 양수 하나와 음수 하나이므로 고정점 $(0,0)$ 은 새들이다. 비슷한 방법으로 $D \mathbb{f} (\pm1,0)$ 의 아이겐 밸류를 계산해보면 $$ \lambda_{1,2} = {{ -\delta \pm \sqrt{\delta^{2} - 8 } } \over { 2 }} $$ 이므로 $D \mathbb{f} (\pm1,0)$ 의 아이겐 밸류는 $\delta > 0$ 일 때 모두 음수고, 고정점 $(\pm1,0)$ 은 싱크다. 그러나 $\delta = 0$ 일 때는 순허수 $\lambda_{1,2} = \pm \sqrt{2} i$ 이므로 고정점 $(\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. ↩︎

댓글