미분방정식 · 위상 초상

위상 초상

위상 초상은 연립미분방정식의 각 해를 변수 평면 위의 곡선으로 나타냅니다. 선형 시스템 x′ = Ax에서는 A의 대각합과 행렬식으로 고윳값이 정해집니다. 두 고윳값이 서로 다르면 이 값들로 그림의 유형도 정해집니다.

이 시각화의 위상 초상

위상 초상은 연립미분방정식의 해를 변수 평면 위의 곡선으로 나타냅니다. 각 시작점을 지나는 곡선이 하나씩 있습니다. 이 페이지에서는 2×22 \times 2 행렬 AA로 주어지는 선형 시스템 x′=Axx' = Ax를 보여 줍니다. 대각합–행렬식 지도에서 (τ,Δ)(\tau, \Delta) 위치의 점을 움직이면 안장점, 결절점, 나선점, 중심점이 나타납니다. 옆의 고윳값 평면에는 각 그림을 결정하는 두 고윳값이 표시됩니다.

시스템, 대각합, 행렬식
(xy)′=(abcd)(xy),τ=a+d,Δ=ad−bc\begin{pmatrix} x \\ y \end{pmatrix}' = \begin{pmatrix} a & b \\ c & d \end{pmatrix} \begin{pmatrix} x \\ y \end{pmatrix}, \qquad \tau = a + d, \quad \Delta = ad - bc

모든 곡선은 방향장을 따릅니다. 화살표는 길이가 모두 같아 방향만 나타냅니다. 속력은 같은 시간 간격의 표시로 알 수 있으며, 움직임이 느린 곳에서는 표시가 모여 있습니다. 일부 해는 원점을 지나는 직선 위에서 움직입니다. Av=λvAv = \lambda v이면 x(t)=eλtvx(t) = e^{\lambda t} v는 해입니다. 그 도함수 λeλtv\lambda e^{\lambda t} v가 이 해에 AA를 적용한 것과 정확히 같기 때문입니다. 따라서 이 직선들은 AA의 고유방향이며, 각 방향의 변화율은 λ\lambda입니다. λ>0\lambda > 0이면 바깥으로, λ<0\lambda < 0이면 안으로 움직이고, λ=0\lambda = 0이면 움직이지 않습니다.

고윳값은 λ2−τλ+Δ=0\lambda^2 - \tau\lambda + \Delta = 0의 해입니다. 따라서 τ\tau는 두 고윳값의 합, Δ\Delta는 곱, τ2−4Δ\tau^2 - 4\Delta는 판별식입니다. 지도의 각 경계는 고윳값 평면에서 일어나는 하나의 변화에 대응합니다. Δ=0\Delta = 0 아래에서는 고윳값의 부호가 서로 달라 안장점이 됩니다. 이 선 위에서는 하나의 고윳값이 0이므로 직선 전체의 점이 정지합니다. 포물선 τ2=4Δ\tau^2 = 4\Delta 위에서는 두 고윳값이 만나고, 그 위쪽에서는 복소수 쌍이 됩니다. Δ>0\Delta > 0일 때 τ\tau의 부호가 감쇠와 성장을 구분하며, τ=0\tau = 0 위에서는 두 고윳값이 허수축에 놓입니다. 따라서 τ\tau와 Δ\Delta으로 고윳값과 그림의 유형이 정해집니다. 다만 중복 고윳값에서는 AA가 단위행렬의 상수배인지에 따라서도 유형이 달라집니다.

복소 고윳값이 α±βi\alpha \pm \beta i이고 β>0\beta > 0이면, 영해를 제외한 모든 해는 시간 2π/β2\pi/\beta 동안 원점을 한 바퀴 돌며 그동안 크기가 e2πα/βe^{2\pi\alpha/\beta}배가 됩니다. 그래서 한 바퀴 간격의 표시는 원점에서 뻗는 하나의 반직선 위에 놓이고, 매번 같은 비율로 크기가 변합니다. α=0\alpha = 0이면 배율은 1이며, 영해를 제외한 모든 궤도가 타원으로 닫힙니다. 이것이 중심점입니다. Δ=1\Delta = 1이고 A=(01−1τ)A = \begin{pmatrix} 0 & 1 \\ -1 & \tau \end{pmatrix}이면 이 시스템은 ζ=−τ/2\zeta = -\tau/2인 감쇠 진동자 x′′+2ζx′+x=0x'' + 2\zeta x' + x = 0입니다.

포물선을 가로지르면 두 직선이 가위처럼 오므라들어 하나가 되고, 그림이 나선으로 바뀌면서 사라집니다. 꼭짓점을 제외한 포물선 위의 전형적인 행렬에는 이 직선 하나만 있으며, 나머지 해는 모두 원점에서 이 직선에 접하며 들어오거나 나갑니다. (τ,Δ)=(−2,1)(\tau, \Delta) = (-2, 1)에서 행렬 −I-I는 고윳값이 같지만 모든 방향이 직선입니다. 기저를 A↦PAP−1A \mapsto PAP^{-1}로 바꾸면 그림의 모양은 변해도 τ\tau, Δ\Delta, 고윳값, 유형은 그대로입니다. 단위행렬의 상수배는 바뀌지 않습니다.

눈여겨볼 점

  • 결절점의 두 직선은 대개 서로 수직이 아닙니다. 두 직선 밖의 해는 느린 방향을 따라 원점에 들어오거나 원점을 떠나며, 멀리서는 빠른 방향 쪽으로 휘어집니다.
  • 파란 곡선은 시작점을 지나는 해 전체입니다. 노란 경로는 시작점을 놓은 이후에 해당하는 부분입니다. 따라서 AA가 바뀌면 이전 해의 기록이 아니라 새 해를 보여 줍니다.

조작 방법

  • 위상 평면의 빈 곳을 누르면 시작점이 추가됩니다. 직선 가까이를 누르면 점이 그 직선 위에 놓입니다.

관련 개념

더 읽어 보기: Wikipedia: Phase portrait (영어).