안정 다양체 정리

1. 개요[편집]

안정 다양체 정리 Stable manifold theorem
다른 이름	아다마르-페롱 정리
전제	쌍곡 고정점 (Re λ ≠ 0, 사상은 |λ| ≠ 1)
결론	Wsloc, Wuloc 존재 · Es, Eu에 접함
매끄러움	벡터장과 같은 Cr (해석적이면 해석적)
유일성	국소 조각은 유일 (중심 다양체와 대비)
증명 노선	그래프 변환 · 적분방정식 (둘 다 축소사상)
일반화	지수 이분성 → 비자율계 · 쌍곡 불변집합

안정 다양체 정리는 쌍곡 고정점 근방에서 그 점으로 지수적으로 수렴하는 점들의 집합이 안정 고유공간 $E^s$ 에 접하는 매끄러운 다양체를 이루고, 대칭적으로 과거 방향에 대해서도 같은 일이 성립한다는 정리다. 아다마르(1901)와 페롱(1928)의 이름을 따 아다마르-페롱 정리라고도 부른다.

하트만-그로브만 정리와 짝을 이루지만 주는 물건이 다르다. 하트만-그로브만은 “궤도 그림이 선형화와 같다”는 위상적 사전을 주고, 그 사전은 일반적으로 미분 가능하지 않다. 안정 다양체 정리는 그림 전체가 아니라 골격 두 개만 떼어 주는 대신 그것이 벡터장과 같은 $C^r$ 급임을 보장한다. 그래서 실무에서 실제로 계산되는 것은 언제나 후자다 — 아무도 하트만-그로브만의 $h$ 를 수치적으로 구하지 않지만, $W^s$ 와 $W^u$ 는 매일 계산된다.

불변 다양체 일반의 지형(중심 다양체, 법선 쌍곡성, 느린 다양체, 유역 경계로서의 의미)은 불변 다양체가 다루므로 여기서 반복하지 않는다. 이 문서는 정리 자체의 진술과 두 증명 노선, 그리고 그 증명이 그대로 알고리즘이 되는 방식에 집중한다.

2. 정리의 진술[편집]

$\dot{\mathbf x} = \mathbf f(\mathbf x)$ 가 $C^r$($r\ge1$)이고 $\mathbf x^*$ 가 쌍곡 평형점이라 하자. $A = D\mathbf f(\mathbf x^*)$ 의 스펙트럼을 $\mathrm{Re}\,\lambda<0$ 쪽과 $>0$ 쪽으로 갈라 $\mathbb R^n = E^s\oplus E^u$ 로 두면, $\mathbf x^*$ 의 근방 $U$ 가 존재해

$$W^s_{\text{loc}} = \{\mathbf x\in U:\ \varphi_t(\mathbf x)\in U\ \ \forall t\ge0,\ \ \varphi_t(\mathbf x)\to\mathbf x^*\}$$

가 $\dim E^s$ 차원의 $C^r$ 매장 원판이고, $\mathbf x^*$ 에서 $E^s$ 에 접한다. 시간을 뒤집으면 $W^u_{\text{loc}}$ 에 대해 같은 결론이 나온다. 수렴은 지수적이며, $0<\alpha<\min|\mathrm{Re}\,\lambda_s|$ 인 임의의 $\alpha$ 에 대해 $\|\varphi_t(\mathbf x)-\mathbf x^*\|\le Ce^{-\alpha t}$ 가 성립한다.

세 가지를 눈여겨봐야 한다.

• 국소 조각은 유일하다. 위 집합은 “근방에 머물면서 수렴하는 점 전부”로 정의되므로 정의상 유일하다. 무한히 많은 중심 다양체가 존재했던 중심 다양체 정리와 정면으로 대비되며, 계산 결과를 검증할 수 있느냐 없느냐를 가르는 차이다.
• 매끄러움을 잃지 않는다. $\mathbf f$ 가 $C^r$ 이면 다양체도 $C^r$, 해석적이면 해석적이다. 중심 다양체는 $C^k$ 마다 근방이 줄어들고 $C^\infty$ 보장이 없었는데, 여기서는 그런 잔소리가 없다. 이유도 뒤에서 보듯 명확하다 — 관련 선형계가 비특이하다.
• 표현은 그래프다. $E^s$ 좌표를 $\mathbf u$, $E^u$ 좌표를 $\mathbf v$ 라 하면 $W^s_{\text{loc}} = \{\mathbf v = h^s(\mathbf u)\}$, $h^s(0)=0$, $Dh^s(0)=0$ 이다. 접한다는 말이 곧 $Dh^s(0)=0$ 이다.

사상 버전도 같다. $F$ 가 $C^r$ 국소 미분동형이고 $|\lambda|\ne1$ 이면 같은 결론이 성립하며, 더 나아가 쌍곡 불변집합(말굽 같은 것)의 각 점마다 안정·불안정 다양체가 정의되어 국소 곱 구조를 만든다. 스메일 말굽의 칸토어 집합 위에서 기호열이 좌표 노릇을 할 수 있는 것이 이 확장 덕분이다.

3. 두 개의 증명 노선[편집]

두 노선 모두 원하는 다양체를 어떤 축소사상의 유일한 고정점으로 구성한다는 점이 같다. 다른 것은 무엇을 미지수로 두느냐다.

3.1. 아다마르의 그래프 변환[편집]

미지수를 다양체 자체(그래프 함수)로 둔다. $E^u$ 위에서 정의되고 립시츠 상수가 1 이하인 함수 $g:E^u\to E^s$ 들의 집합 $\mathcal G$ 를 잡자. 그래프 $\mathrm{graph}(g)$ 를 시간 $T$ 만큼 흘려보내면 — 컷오프로 비선형 항을 작게 만들어 둔 덕분에 — 그 상도 다시 $E^u$ 위의 그래프가 된다. 이렇게 정의된 그래프 변환 $\Gamma: g\mapsto g'$ 가 $\mathcal G$ 위의 축소사상임을 보이고, 바나흐 고정점 정리로 유일한 $g^*$ 를 얻으면 $\mathrm{graph}(g^*) = W^u_{\text{loc}}$ 다.

수축이 일어나는 이유는 기하적으로 뻔하다. 불안정 방향은 늘어나고 안정 방향은 줄어드니, 그래프의 기울기가 반복마다 $E^u$ 쪽으로 눌린다. 립시츠 상수 1짜리 콘(cone) 안에 들어온 것은 다시 나가지 못한다는 콘 조건이 이 논증의 심장이고, 같은 조건이 나중에 법선 쌍곡성 이론으로 그대로 확장된다.

3.2. 랴푸노프-페롱 적분방정식[편집]

미지수를 궤도로 둔다. $\dot{\mathbf x} = A\mathbf x + \mathbf g(\mathbf x)$ 에서 $t\ge0$ 동안 유계인 해를 찾는다고 하고, 상수변화법을 안정·불안정 사영 $P^s, P^u$ 로 쪼개 쓰면 다음 적분방정식이 나온다.

$$\mathbf x(t) = e^{At}P^s\boldsymbol\xi + \int_0^t e^{A(t-s)}P^s\,\mathbf g(\mathbf x(s))\,ds - \int_t^{\infty} e^{A(t-s)}P^u\,\mathbf g(\mathbf x(s))\,ds$$

핵심은 적분 구간의 방향을 성분마다 다르게 잡은 것이다. 안정 성분은 과거에서 현재로, 불안정 성분은 미래에서 현재로 적분해야 각 지수함수가 감쇠 쪽으로 작동한다. 이렇게 짜 놓으면 우변이 지수 감쇠 함수공간 위의 축소사상이 되고, 유일한 해 $\mathbf x(\cdot;\boldsymbol\xi)$ 를 얻는다. $h^s(\boldsymbol\xi) = P^u\mathbf x(0;\boldsymbol\xi)$ 가 찾던 그래프다.

이쪽 노선이 실용적으로 더 중요하다. 적분방정식의 피카르 반복이 그대로 수치 알고리즘이고, 위 식의 사영-적분 구조가 뒤에서 볼 사영 경계조건의 정체이기 때문이다. $\boldsymbol\xi$ 에 대한 미분가능성을 반복 안에서 같이 추적하면 $C^r$ 매끄러움도 따라 나온다.¹

4. 지수 이분성 — 비자율계로의 확장[편집]

평형점이 없는 문제, 예컨대 어떤 기준 궤도 $\bar{\mathbf x}(t)$ 주위의 변분방정식 $\dot{\mathbf y} = A(t)\mathbf y$ 를 다루려면 “고유값”이라는 말을 버려야 한다. 대체물이 지수 이분성(exponential dichotomy)이다. 기본해행렬 $\Phi(t)$ 와 사영족 $P(t)$ 가 있어

$$\|\Phi(t)P(s)\Phi(s)^{-1}\| \le Ke^{-\alpha(t-s)}\ (t\ge s),\qquad \|\Phi(t)(I-P(s))\Phi(s)^{-1}\| \le Ke^{-\alpha(s-t)}\ (t\le s)$$

가 성립하면 이 선형계가 지수 이분적이라 한다. 상수계에서 쌍곡성과 정확히 같은 조건이고, 랴푸노프-페롱 구성은 문자 그대로 이 부등식만 쓴다. 그래서 증명이 통째로 비자율계로 옮겨간다.

실무적으로 더 중요한 것은 거칠기(roughness)다. $A(t)$ 를 조금 흔들어도 지수 이분성은 살아남고 사영족도 조금만 움직인다. 이 정리가 없으면 “이산화 오차를 먹은 수치 궤도의 안정 다양체”라는 말 자체가 성립하지 않는다. 절단 구간 위의 경계값 문제가 잘 세워진 문제(well-posed)가 되는 근거, 그림자 정리의 조건, 호모클리닉 궤도 추적 코드의 수렴 증명이 전부 지수 이분성의 거칠기 위에 서 있다.

5. 람다 보조정리[편집]

람다 보조정리(inclination lemma, 팔리스)는 대역 다양체가 어떻게 서로를 축적하는지를 말한다. $p$ 가 쌍곡 고정점이고 $D$ 가 $W^s(p)$ 를 어떤 점에서 횡단하는 $\dim E^u$ 차원 원판이면, 반복상 $f^n(D)$ 안에 $W^u_{\text{loc}}(p)$ 에 $C^1$ 으로 임의로 가까운 원판이 들어 있다.

결론이 강하다. 불안정 다양체는 자기 근처를 지나는 모든 것을 자기 쪽으로 정렬시킨다. 그래서 아무 초기조건이나 전진 적분하면 궤도가 곧 $W^u$ 에 얹히고, 반대로 $W^s$ 는 전진 적분으로 절대 보이지 않는다. 스메일-버코프 정리에서 횡단 호모클리닉 점 하나로부터 말굽을 뽑아낼 때 실제로 삽질을 하는 것도 이 보조정리다 — $W^u$ 조각이 $W^s$ 를 가로지르며 되돌아오는 것을 반복하면 필요한 띠 배치가 자동으로 생긴다.

대역 다양체 자체는 국소 조각을 흘려 정의한다. $W^s = \bigcup_{t\le0}\varphi_t(W^s_{\text{loc}})$ 이며, 이때 매장성이 깨질 수 있다는 점(자기 자신에게 감기는 얽힘)은 불변 다양체 쪽에 정리돼 있다.

6. 수치적으로 성장시키기[편집]

정리는 존재만 말하고 위치는 말하지 않는다. 실제 계산은 다음 순서로 간다.

1. 씨앗 놓기. $W^u$ 의 1차원 가지는 불안정 고유벡터 $\mathbf v$ 방향으로 $\mathbf x_0 = \mathbf x^* + \delta\mathbf v$ 를 잡고 전진 적분하면 얻어진다. 그런데 다양체는 $E^u$ 에 접할 뿐 일치하지 않으므로 이 씨앗은 진짜 다양체에서 $O(\delta^2)$ 만큼 벗어나 있고, 상대오차는 $O(\delta)$ 다. $\delta$ 를 줄이면 정확해지지만 평형점에서 빠져나오는 데 걸리는 시간이 $\ln(1/\delta)/\lambda_u$ 로 늘어 적분 오차가 쌓인다. 양쪽으로 손해 보는 트레이드오프라, 진지한 계산은 씨앗을 고차로 놓는다.

2. 고차 씨앗 — 매개화. 다양체를 그래프가 아니라 매개화 $K(s) = \mathbf x^* + s\mathbf v + \sum_{k\ge2}K_k s^k$ 로 두고, 그 위의 동역학이 $\dot s = \lambda s$ 가 되도록 요구하면 불변성 방정식 $\mathbf f(K(s)) = K'(s)\lambda s$ 가 나온다. 차수별로 정리하면 매 차수마다

$$(A - k\lambda I)\,K_k = (\text{낮은 차수로 이미 결정된 항})$$

를 푸는 일이 된다. 이 선형계가 비특이한 조건이 $k\lambda \notin \sigma(A)$, 즉 비공명이다. 안정 다양체 정리가 중심 다양체와 달리 잔소리 없이 매끄러운 진짜 이유가 이것이고, 공명이 걸릴 때 무슨 일이 생기는지는 스턴버그 정리와 정규형 이론이 다룬다. 계수를 $m$ 차까지 구해 씨앗을 놓으면 오차가 $O(\delta^{m+1})$ 로 떨어진다. 이 방식을 매개화 방법이라 부른다.

3. 기본 영역으로 덮기. 씨앗 $\mathbf x_0$ 와 그 상 $\varphi_T(\mathbf x_0)$ 를 잇는 호를 기본 영역이라 한다. 이 호 하나를 촘촘히 채우고 반복상을 취하면 $W^u$ 가지 전체가 중복 없이 덮인다. 무작정 오래 적분하는 것보다 점 분포가 훨씬 고르다.

4. 안정 다양체는 역방향으로. 시간을 뒤집으면 $W^s$ 는 $W^u$ 가 되므로 같은 절차를 쓴다. 다만 뒤집는 순간 원래 안정하던 방향이 가장 강한 팽창 방향이 되어, 강한 감쇠가 섞인 계에서는 역방향 적분이 곧바로 발산한다. 차원이 2 이상이면 궤도가 지배적 방향으로 몰려 다양체가 한쪽만 그려지는 문제도 생긴다. 이 지점부터는 적분을 버리고 경계값 문제로 가는 것이 정석이며, 측지 준위집합 방법 같은 구체적 처방은 불변 다양체에 정리돼 있다.

5. 사영 경계조건. 무한 구간을 $[-T, T]$ 로 자를 때 끝점에 무엇을 요구할 것인가. 순진하게 $\mathbf x(-T)=\mathbf x^*$ 라고 쓰면 평형점은 유한 시간에 도달하지 않으므로 해가 없다. 바른 처방은 끝점을 고유공간에 사영하는 것이다.

$$L^s\big(\mathbf x(-T)-\mathbf x^*\big) = 0, \qquad L^u\big(\mathbf x(T)-\mathbf x^*\big) = 0$$

여기서 $L^s, L^u$ 는 각각 안정·불안정 고유공간을 소거하는 행렬이다. 즉 왼쪽 끝은 $W^u_{\text{loc}}$ 의 선형 근사 위에, 오른쪽 끝은 $W^s_{\text{loc}}$ 의 선형 근사 위에 놓으라는 뜻이며, 이것이 랴푸노프-페롱 적분방정식의 사영 구조가 이산 세계에서 취하는 모습이다. 바인(1990)이 이 절단의 오차가 $T$ 에 대해 지수적으로 작다는 것을 보였기 때문에, $T$ 를 조금만 키워도 정확도가 확 좋아진다.

6. 연결 궤도 찾기. 호모클리닉·헤테로클리닉 연결은 결국 $W^u$ 와 $W^s$ 의 교차다. 실무 표준은 린(Lin)의 방법이다. 단면 $\Sigma$ 를 하나 잡고, 불안정 쪽에서 출발한 궤도 조각과 안정 쪽에서 거슬러 온 궤도 조각을 각각 $\Sigma$ 까지 적분한 뒤, 두 도착점의 차이를 미리 고른 방향 $Z$ 로만 허용한다. 이 린 간극 $\eta(\mu)$ 는 매개변수의 매끄러운 함수이고, $\eta(\mu)=0$ 이 곧 연결의 존재 조건이다. 다양체를 통째로 그려 눈으로 교차를 찾는 대신 스칼라 하나의 영점을 수치 연속법으로 추적하는 문제로 바뀌는 것이 이 방법의 값어치이며, MatCont와 AUTO의 호모클리닉 연속화가 이 위에 서 있다.²

7. 역할 분담 정리[편집]

같은 쌍곡 평형점을 놓고 네 정리가 서로 다른 질문에 답한다. 헷갈리면 이 표만 기억하면 된다.

정리	주는 것	매끄러움
하트만-그로브만	선형화와의 국소 위상 공액	C⁰ (뇔더)
안정 다양체	두 불변 다양체와 접선 방향	벡터장과 동급
스턴버그	비공명 조건 아래 매끄러운 선형화	Cᵏ / C∞
중심 다양체	비쌍곡일 때 남는 방향	유한 Cᵏ, 비유일

정리의 힘이 셀수록 가정도 세다. 안정 다양체 정리가 유독 널리 쓰이는 이유는 가정은 하트만-그로브만과 똑같으면서 결론은 매끄럽기 때문이다. 공짜로 얻는 매끄러움이 아니라, 그림 전체가 아니라 골격만 요구하는 대가로 얻는 매끄러움이다.³

8. 관련 문서[편집]

• 불변 다양체 · 하트만-그로브만 정리 · 중심 다양체 정리
• 스턴버그 정리 · 위상 공액 · 정규형 이론
• 호모클리닉 궤도 · 스메일 말굽 · 안장점 · 흡인 유역
• 수치 연속법 · MatCont · 뉴턴-랩슨법
• 동역학계 · 구조적 안정성 · 분기 이론
• 자코비안 행렬 · 고유값 문제 · 테일러 급수

9. Footnotes[편집]

1. 두 노선의 성격 차이가 재미있다. 그래프 변환은 “다양체를 통째로 하나의 점으로 보는” 기하적 증명이라 그림이 잘 그려지고, 적분방정식 쪽은 계산은 지저분하지만 가정을 부등식 몇 줄로 압축해서 일반화가 잘 된다. 교과서가 둘 다 싣는 데는 이유가 있다. ↩

2. 다양체를 그려서 교차를 찾겠다는 접근은 3차원까지는 예쁘고 4차원부터는 재앙이다. 그림이 필요한 것은 논문의 그림 4이지 연속화 코드가 아니다. 간극 하나를 뉴턴으로 0으로 보내는 쪽이 훨씬 빠르고 정확한데도, 처음 배우는 사람은 대개 다양체를 먼저 그려 보려 한다. ↩

3. 그래서 “선형화가 정당한가”라는 질문에는 하트만-그로브만을, “유역 경계가 어디인가”라는 질문에는 안정 다양체 정리를 인용해야 한다. 둘을 섞어 쓰면 심사위원에게 정확히 그 지적을 받는다. ↩