동역학계

계산물리 수치해석 물리 마지막 수정: 2026-08-09 04:12:41

1. 개요[편집]

동역학계 Dynamical system
구성	상태공간 + 시간에 따른 진화 규칙
연속 시간	흐름 φ_t (상미분방정식)
이산 시간	사상 f (반복 적용)
주요 대상	불변집합 · 극한집합 · 끌개
동일성 기준	위상적 공액 / 위상적 동등
수치해석과의 관계	시간전진 솔버 = 이산 동역학계

동역학계(dynamical system)는 가능한 상태들의 집합(상태공간)과, 한 상태가 시간이 지나면 어느 상태로 가는지를 정하는 결정론적 규칙의 짝이다. 규칙이 연속 시간이면 흐름(flow) $\varphi_t$ , 이산 시간이면 사상(map) $f$ 라 부르고, 둘 다 “지금 상태만 알면 미래가 정해진다”는 마르코프성 위에 서 있다.

이 문서가 심위키에 있는 이유는 단순하다. 시간을 전진시키는 모든 수치 솔버는 그 자체로 하나의 이산 동역학계이기 때문이다. 우리가 풀고 싶은 미분방정식이 동역학계이고, 그것을 이산화해서 얻은 코드도 동역학계인데, 둘의 동역학이 같다는 보장은 어디에도 없다. 이 간극에서 나오는 것이 가짜 고정점, 뒤집힌 안정성, 존재하지 않는 카오스다. 분기 이론, 카오스 이론, 랴푸노프 지수, 수치 연속법이 전부 이 문서의 자식들이다.

2. 흐름과 사상[편집]

연속 시간 계는 자율 상미분방정식 $\dot{\mathbf x} = \mathbf f(\mathbf x)$ 로 주어지고, 그 해가 정의하는 사상족 $\varphi_t : \mathbb R^n \to \mathbb R^n$ 을 흐름이라 한다. 흐름은 군 성질을 만족한다.

\varphi_0 = \mathrm{id}, \qquad \varphi_{t+s} = \varphi_t \circ \varphi_s

이산 시간 계는 사상 $f$ 를 반복하는 것이다. 즉 $\mathbf x_{n+1} = f(\mathbf x_n)$ 이고 궤도는 $\{f^n(\mathbf x_0)\}$ . 둘은 완전히 다른 세계가 아니라 서로 오간다.

흐름 → 사상: 시간 간격 $h$ 마다 샘플링하면 시간- $h$ 사상 $\varphi_h$ 가 나온다. 이것이 정확히 시간 적분기가 근사하려는 대상이다.
흐름 → 사상 (횡단면): 주기궤도를 가로지르는 초곡면(푸앵카레 단면)을 잡고, 궤도가 그 면으로 되돌아오는 지점을 대응시키면 되돌아옴 사상이 되어 차원이 하나 줄어든다. 흐름의 주기궤도가 사상의 고정점으로 바뀌므로 분기 판정이 훨씬 쉬워진다.

차원 제약도 다르다. 평면 자율 흐름은 푸앵카레-벤딕손 정리 때문에 카오스가 불가능하다 — 유계 궤도의 극한집합은 평형점이거나 주기궤도뿐이라 연속계에서 카오스를 보려면 최소 3차원이 필요하다. 반면 사상은 1차원짜리 로지스틱 사상 하나로도 카오스를 만든다. 궤도가 자기 자신을 가로지를 수 없다는 조르당 곡선 제약이 사상에는 없기 때문이다.¹

3. 궤도, 극한집합, 불변집합[편집]

한 점에서 출발한 궤도가 결국 어디로 가느냐를 담는 것이 극한집합이다.

\omega(\mathbf x) = \{\mathbf y : \varphi_{t_k}(\mathbf x) \to \mathbf y \ \text{인 } t_k \to +\infty \ \text{존재}\}

시간을 거꾸로 돌린 것이 $\alpha$ -극한집합이다. 유계 궤도의 $\omega$ -극한집합은 항상 공집합이 아니고, 연결되어 있으며, 불변이다. 집합 $S$ 가 불변이라는 것은 $\varphi_t(S) = S$ , 즉 한 번 들어가면 영원히 못 나온다는 뜻이다. 불변집합의 종류가 곧 동역학계 도감이다.

불변집합	정체	대표 사례
평형점	$\mathbf f(\mathbf x^*)=0$	정상해, 수렴한 CFD 해
주기궤도	$\varphi_T = \mathrm{id}$ 인 폐곡선	극한 순환, 카르만 와열
불변 토러스	준주기 운동, 무리수 회전수	KAM 정리, 적분가능계
이상 끌개	프랙탈 구조의 끌개	로렌츠 계, 프랙탈 차원
안장형 집합	끌지 않는 카오스 집합	과도 카오스, 말굽

끌개(attractor)는 이 중 “실제로 관측되는” 부류다. 흔한 정의는 (i) 불변이고 (ii) 자기 근방에서 지수적으로 끌어당기며 (iii) 더 작은 불변 부분집합으로 쪼갤 수 없는(추이적) 집합이다. 끌개로 들어가는 초기조건 전체를 흡인 유역이라 하고, 유역 경계는 대체로 안장형 불변집합과 그 불변 다양체로 이루어진다. “같은 해석 조건인데 왜 케이스마다 답이 다르냐”의 절반은 초기값이 다른 유역에 떨어졌다는 뜻이며, 이건 수렴 실패가 아니라 계의 성질이다.

4. 안정성의 언어[편집]

“안정하다”는 말이 이 바닥에서는 세 가지 다른 뜻으로 쓰이므로 구분해 두는 편이 좋다.

랴푸노프 안정. 아무리 작은 허용 오차 $\varepsilon$ 을 줘도, 그 안에 영원히 머물게 하는 초기 오차 $\delta$ 가 존재한다. 진자의 아래 평형점처럼 가까이 있기만 하고 수렴은 안 하는 경우가 여기 해당한다.
점근 안정. 랴푸노프 안정이면서 근방에서 실제로 수렴한다. 수렴이 지수적이면 지수 안정이라 하고, 공학에서 “안정”이라 하면 대개 이 뜻이다.
궤도 안정. 주기궤도에 대해서는 점 대 점 비교가 무의미하다(위상이 조금만 어긋나도 거리가 안 줄어든다). 그래서 집합으로서의 궤도에 가까워지는가를 묻고, 판정은 모노드로미 행렬의 플로케 승수로 한다. 승수 하나는 항상 1이며(궤도를 따라 미는 방향), 나머지 승수가 전부 단위원 안이면 안정이다.

판정 도구는 두 갈래다. 하나는 야코비안 고유값을 보는 선형화 방식으로, 그 정당성을 하트만-그로브만 정리가 대 준다. 다른 하나는 랴푸노프 함수 — 평형점에서 0이고 주위에서 양수이며 궤도를 따라 감소하는 스칼라 함수를 하나 찾으면 선형화 없이 안정성이 증명된다. 선형계에서 이 함수를 이차형식으로 잡는 순간 나오는 것이 리아푸노프 방정식 $A^{\mathsf T}P + PA = -Q$ 이고, 제어 이론에서 이차형 비용과 함께 최적 제어로 이어진다. 문제는 비선형 계에서 랴푸노프 함수를 찾는 일반 방법이 없다는 것 — 물리계라면 에너지에서 출발하는 것이 국룰이다.

5. 보존이냐 소산이냐[편집]

벡터장의 발산이 부피 변화율을 준다. 유역 $V(t)$ 에 대해

\frac{dV}{dt} = \int_{V} \nabla \cdot \mathbf f \; d\mathbf x

이므로 $\nabla\cdot\mathbf f \equiv 0$ 이면 상공간 부피가 보존된다. 해밀턴 계가 정확히 이 경우이고(리우빌 정리), 그래서 해밀토니안 역학에는 끌개가 없다 — 부피가 줄지 않는데 무언가로 수렴할 수는 없다. 반대로 소산계는 부피가 지수적으로 줄고, 그 극한이 측도 0의 끌개다. 로렌츠 계는 발산이 상수 $-(\sigma+1+b)$ 라서 부피가 균일하게 수축하는데, 그러면서도 궤도는 한 점으로 안 모이고 프랙탈 위에서 논다. 이 조합이 “이상(strange)“의 정체다.

이 구분이 적분기 선택으로 직결된다. 보존계를 룽게-쿠타법 같은 일반 적분기로 오래 돌리면 에너지가 서서히 새서 궤도가 나선을 그린다. 그래서 분자동역학과 천체역학은 심플렉틱 적분기와 기하 수치적분을 쓴다. “오차가 작다”가 아니라 “구조를 깨지 않는다”가 기준이 되는 순간이다.

6. 언제 두 계를 같다고 하는가[편집]

동역학계론의 동일성은 좌표가 아니라 궤도의 위상 구조로 정의한다. 두 흐름 $\varphi_t, \psi_t$ 사이에 위상동형 $h$ 가 있어

h \circ \varphi_t = \psi_t \circ h \qquad (\forall t)

이면 위상적 공액이라 하고, 시간 재매개화까지 허용하면 위상적 동등이라 한다. 이 개념 위에서 두 가지 큰 정리가 선다. 하나는 하트만-그로브만 정리 — 쌍곡 평형점 근방에서 계는 자기 선형화와 국소 공액이다. 다른 하나는 구조적 안정성 — 계를 조금 흔들어도 위상 구조가 안 바뀌면 구조적으로 안정하다고 하고, 그렇지 않은 지점이 곧 분기점이다. 쌍곡성이 깨지는 자리에서는 중심 다양체 정리와 정규형 이론이 바통을 받는다.

7. 수치 솔버의 동역학은 다르다[편집]

여기가 이 문서의 결론부다. 단계 크기 $h$ 를 고정한 적분기는 사상 $\Phi_h$ 를 정의하고, 우리는 그것이 $\varphi_h$ 와 닮았기를 바란다. 그런데 다음 일들이 실제로 벌어진다.

가짜 고정점과 가짜 주기해. 로지스틱 방정식 $\dot x = x(1-x)$ 에 전진 오일러를 적용하면

x_{n+1} = x_n + h\,x_n(1-x_n)

이 되는데, 이것은 선형 좌표변환으로 로지스틱 사상 $y_{n+1} = r\,y_n(1-y_n)$ , $r = 1+h$ 와 정확히 공액이다. 원래 미분방정식은 모든 양의 초기값에서 단조롭게 $x=1$ 로 간다. 그런데 수치해는 $h>2$ 에서 주기 2로 갈라지고 $h \approx 2.57$ 부터 카오스에 빠진다.² 미분방정식에 없던 주기배가 캐스케이드가 순전히 이산화 때문에 생긴 것이다. 일반적인 룽게-쿠타 계열도 마찬가지로 원래 방정식에 대응물이 없는 가짜 고정점·가짜 주기 2 해를 갖고, 그중 일부는 안정하기까지 하다(Iserles, 1990).

안정성이 뒤집힌다. $\dot x = \lambda x$ 에 전진 오일러를 쓰면 수치 고정점의 승수는 $1+h\lambda$ 이므로, $\lambda<0$ 인데도 $h|\lambda|>2$ 이면 발산한다. 반대로 과도한 수치 소산은 진짜 불안정을 눌러 없애서 존재하는 진동을 못 보게 만든다. 그래서 강성 방정식에서 음함수 적분기를 쓰고, 폰 노이만 안정성 해석으로 도식의 증폭인자를 미리 따진다.

그럼에도 희망은 있다. 후진 오차 해석은 수치 사상이 원래 방정식이 아니라 약간 다른 방정식(수정 방정식)의 정확한 흐름과 지수적으로 가깝다는 결과를 준다. 심플렉틱 적분기의 수정 방정식이 다시 해밀턴 계라는 사실이 장기 에너지 보존의 진짜 이유다. 카오스 쪽에서는 그림자 정리가 “네 궤적 근처에 진짜 궤적이 하나 있다”는 위안을 주지만, 그건 궤적이 맞다는 보장이 아니라 통계량만 믿으라는 권고에 가깝다.³ 자료동화와 앙상블 예보가 단일 궤적을 포기한 이유가 여기에 있다.

실무 지침으로 압축하면 이렇다. 정상해를 구했으면 야코비안 고유값을 한 번은 봐라. 잔차가 내려갔다는 것은 고정점에 갔다는 뜻이지 그 고정점이 안정하다는 뜻이 아니다.⁴ 그리고 시간 정확도가 필요 없는 문제라 해서 $h$ 를 무한정 키우면, 어느 순간부터 당신이 푸는 것은 원래 방정식이 아니라 그 이산화가 만든 다른 동역학계다.

8. 관련 문서[편집]

그래서 “우리 2차원 정상 해석 코드에서 카오스가 나왔다”는 보고는 십중팔구 카오스가 아니라 수렴 실패다. 물론 비자율(시간 의존 강제항이 있는) 2차원 계는 사실상 3차원이라 얼마든지 카오스가 난다 — 강제 더핑 방정식이 그 표본이다. ↩
$r_\infty \approx 3.569946$ 이므로 $h_\infty = r_\infty - 1 \approx 2.5699$ . 즉 “적당히 큰 시간 간격”이라는 이유로 $h=3$ 을 넣는 순간, 당신의 개체수 모형은 결정론적 카오스를 뱉는다. 물리는 그대로인데 결과만 미쳐 있는 상황이라 디버깅이 제일 괴로운 부류다. ↩
이 바닥의 국룰: 카오스 계에서는 궤적을 자랑하지 말고 평균·분산·스펙트럼을 자랑해라. 궤적을 자랑하는 순간 리뷰어가 “초기조건 끝자리 하나 바꿔서 다시 돌려보라”고 한다. ↩
정상 솔버가 불안정 가지 위에 얌전히 앉아 있는 일은 생각보다 흔하다. 수치 소산이 실제 불안정 모드를 눌러 놓았을 뿐이라, 격자를 조밀하게 하거나 비정상 해석으로 넘어가는 순간 갑자기 진동이 터진다. 그리고 그 시점은 대체로 보고서 제출 전날이다. ↩