중심 다양체 정리

1. 개요[편집]

중심 다양체 정리 Center manifold theorem
대상	야코비안이 허수축 위에 고유값을 갖는 평형점
결론	중심 고유공간에 접하는 불변다양체 Wc 존재
핵심 도구	불변성 방정식 (테일러 차수별 풀이)
주의	유일하지 않음 · 해석적이지 않을 수 있음
매끄러움	임의 유한 Ck, 단 근방이 k에 따라 줄어듦
짝 이론	정규형 이론

중심 다양체 정리(center manifold theorem)는 평형점에서 자코비안 행렬의 고유값을 안정($\mathrm{Re}\,\lambda<0$)·불안정($>0$)·중심($=0$) 세 무리로 나눌 때, 중심 고유공간 $E^c$ 에 접하는 국소 불변다양체 $W^c$ 가 존재하고, 평형점 근방의 본질적 동역학이 전부 그 위로 환원된다는 정리다. 한 문장으로 줄이면 “허수축에 붙어 있는 방향만 남기고 나머지는 버려도 된다”이며, 100차원 계의 분기 해석이 결국 1~2차원 방정식 한두 줄로 끝나는 이유가 이것이다.

왜 필요한가. 고유값이 전부 허수축 밖에 있으면 하트만-그로브만 정리가 국소 위상적 등가로 선형화를 정당화해 준다. 즉 비선형은 신경 쓸 필요가 없다. 그런데 분기가 일어나는 순간은 정의상 고유값이 허수축 위에 올라온 순간이라 하트만-그로브만이 침묵한다. 이때 “그럼 어느 방향의 비선형을 봐야 하나”에 답하는 것이 중심 다양체 정리이고, “그 방향에서 어떤 비선형 항이 본질인가”에 답하는 것이 정규형 이론이다. 둘은 반드시 짝으로 쓴다.

2. 정리의 진술[편집]

좌표를 잘 잡아 계를 블록으로 나눈다.

$$\dot{\mathbf u} = A\mathbf u + f(\mathbf u,\mathbf v), \qquad \dot{\mathbf v} = B\mathbf v + g(\mathbf u,\mathbf v)$$

여기서 $\mathbf u \in \mathbb R^{c}$ 이고 $A$ 의 고유값은 전부 허수축 위, $\mathbf v \in \mathbb R^{s}$ 이고 $B$ 의 고유값은 전부 실부가 음수, $f, g$ 는 원점에서 값과 1계 미분이 0인 $C^k$ 함수다. 그러면 원점 근방에서

$$W^c = \{(\mathbf u,\mathbf v) : \mathbf v = h(\mathbf u),\ h(0)=0,\ Dh(0)=0\}$$

꼴의 국소 불변다양체가 존재하고, 그 위의 흐름은 $c$ 차원 방정식

$$\dot{\mathbf u} = A\mathbf u + f(\mathbf u, h(\mathbf u))$$

로 주어진다. 이 축약 방정식이 환원된 계(reduced system)다.

환원 원리(reduction principle, 쇼시타이슈빌리)는 여기서 한 걸음 더 나간다. 원래 계는 원점 근방에서

$$\dot{\mathbf u} = A\mathbf u + f(\mathbf u,h(\mathbf u)), \qquad \dot{\mathbf v}_s = -\mathbf v_s, \qquad \dot{\mathbf v}_u = +\mathbf v_u$$

와 국소 위상적으로 동등하다. 즉 안정·불안정 방향은 표준 안장으로 완전히 규격화되고, 계의 개성은 오로지 중심 방향에만 남는다. 실용적 따름정리도 바로 나온다. 불안정 고유값이 없을 때, 환원된 계의 원점이 점근 안정이면 원래 계도 점근 안정이고, 불안정이면 불안정이다.

3. 불변성 방정식과 실제 계산[편집]

$W^c$ 는 존재만 보장될 뿐 닫힌 형태로 주어지지 않는다. 그래서 실무는 늘 테일러 계수를 차수별로 푸는 방식으로 간다. $\mathbf v = h(\mathbf u)$ 를 시간 미분하면 $\dot{\mathbf v} = Dh(\mathbf u)\dot{\mathbf u}$ 이므로, 두 식을 맞붙이면 다음 불변성 방정식(invariance equation)이 나온다.

$$Dh(\mathbf u)\big[A\mathbf u + f(\mathbf u, h(\mathbf u))\big] = B\,h(\mathbf u) + g(\mathbf u, h(\mathbf u))$$

이 함수방정식을 $h(\mathbf u) = \sum_{|\mathbf m|\ge 2} h_{\mathbf m}\mathbf u^{\mathbf m}$ 로 전개해 차수별로 계수를 맞추면, 각 차수마다 유일하게 풀리는 선형계가 나온다. $B$ 의 고유값이 허수축 밖에 있어 관련 연산자가 가역이기 때문이다. 2차 계수를 구하면 환원된 계의 3차 항까지가 확정되고, 그것이 곧 정규형 이론의 제1 랴푸노프 계수 계산에 들어간다.

교과서 예제 하나면 감이 온다.

$$\dot x = x y, \qquad \dot y = -y - x^2$$

원점에서 고유값은 $0, -1$ 이라 $c=1, s=1$ 이다. $h(x) = a_2x^2 + a_3x^3+\cdots$ 로 놓고 불변성 방정식 $h'(x)\,x\,h(x) = -h(x) - x^2$ 에 대입하면, $x^2$ 차수에서 좌변은 0이고 우변은 $-a_2 - 1$ 이므로 $a_2 = -1$. 따라서 $h(x) = -x^2 + O(x^3)$ 이고, 환원된 계는

$$\dot x = x\,h(x) = -x^3 + O(x^4)$$

가 되어 원점은 점근 안정이라는 결론이 나온다. 여기서 $h$ 를 무시하고 $y \equiv 0$ 로 놓았다면 $\dot x = 0$, 즉 “중립”이라는 완전히 틀린 답을 얻었을 것이다. 중심 다양체의 굽음(curvature)이 답의 부호를 만든다는 것이 이 계산의 교훈이다.

3.1. 파라미터를 상태로 승격시키기[편집]

분기 해석에서는 파라미터 $\mu$ 가 함께 움직여야 한다. 표준 요령은 $\dot\mu = 0$ 이라는 자명한 방정식을 계에 추가해 $\mu$ 를 상태변수로 승격시키는 것이다. 그러면 $\mu$ 방향의 고유값이 0이 되어 자동으로 중심 방향에 편입되고, 얻어지는 $W^c$ 는 파라미터에 의존하는 중심 다양체 족이 된다. 이 한 줄짜리 트릭 덕분에 앞의 절차가 파라미터 있는 문제에 그대로 적용된다.¹

4. 유일하지 않고, 해석적이지도 않다[편집]

정리의 진술에 “유일”이라는 단어가 없는 데는 이유가 있다.

유일성 반례. 다음 계를 보자.

$$\dot x = x^2, \qquad \dot y = -y$$

불변성 방정식은 $h'(x)x^2 = -h(x)$ 이고, 변수분리로 풀면 $h(x) = c\,e^{1/x}$ 다. $x<0$ 에서는 $x\to 0^-$ 일 때 $e^{1/x}$ 와 그 모든 도함수가 0으로 가므로, 이것을 $x\ge0$ 에서 $h\equiv0$ 과 이어붙이면 모든 $c$ 에 대해 $C^\infty$ 인 중심 다양체가 만들어진다. 즉 중심 다양체는 무한히 많다. 다행인 것은 이들 전부가 같은 테일러 급수(여기서는 항등적으로 0)를 갖는다는 점이다. 서로 다른 중심 다양체는 평평 접촉(flat contact)으로만 갈리므로, 급수 계산으로 얻는 결론 — 안정성, 분기 종류 — 은 어느 것을 골라도 같다.²

해석성 반례. 더 유명한 것은 이쪽이다.

$$\dot x = -x^3, \qquad \dot y = -y + x^2$$

$h = \sum a_n x^n$ 을 불변성 방정식 $-x^3h'(x) = -h(x) + x^2$ 에 넣고 차수를 맞추면 $a_2 = 1$ 과 점화식 $a_k = (k-2)a_{k-2}$ 가 나온다. 홀수 계수는 전부 0이고 짝수 계수는 $a_2=1,\ a_4=2,\ a_6=8,\ a_8=48,\dots$ 로 계승 속도로 폭발한다. 수렴 반경은 0이다. 즉 이 계의 중심 다양체는 $C^\infty$ 이지만 해석적이지 않으며, 테일러 급수는 형식적 점근 전개일 뿐이다.

여기서 실무 지침이 나온다. 차수를 무작정 올린다고 정확해지지 않는다. 필요한 차수(보통 분기 계수를 확정하는 3차 정도)에서 끊고, 남는 것은 오차로 인정하는 것이 정직한 태도다. 매끄러움도 마찬가지여서, 임의의 유한 $k$ 에 대해 $C^k$ 중심 다양체가 존재하긴 하지만 그 다양체가 정의되는 근방의 크기가 $k$ 가 커질수록 줄어들 수 있다. ”$C^\infty$ 중심 다양체가 존재한다”는 일반 정리는 없다.

5. 무한차원으로의 확장[편집]

같은 논리는 편미분방정식에도 적용된다. 편미분방정식의 선형 작용소가 (i) 허수축 위에 유한개의 고유값만 갖고 (ii) 나머지 스펙트럼이 허수축에서 유한한 간격만큼 떨어져 있으면, 무한차원 상태공간 안에 유한차원 중심 다양체가 존재한다. 결과는 극적이다 — 무한 자유도 PDE의 불안정 개시 근방 거동이 상미분방정식 몇 개로 요약된다.

• 진폭 방정식. 레일리-베나르 대류나 전단류 불안정의 임계점 근방에서 나오는 스튜어트-란다우 방정식, 긴즈부르크-란다우 방정식이 사실상 중심 다양체 축약 + 정규형의 산물이다.
• 관성 다양체. 소산성 PDE에서 모든 궤도가 지수적으로 끌려 들어가는 유한차원 불변다양체다. 중심 다양체가 평형점 근방의 국소 이야기인 것과 달리 관성 다양체는 대역적 끌개 전체를 담으려 한다는 점이 다르고, 존재 증명에는 훨씬 강한 스펙트럴 간격 조건이 필요하다.
• 축소차수모델. POD/DMD 같은 데이터 기반 축약은 “저차원 불변 구조가 있다”는 가정 위에 서 있는데, 그 가정에 이론적 근거를 대 주는 것이 중심·관성 다양체다. 반대로 말하면, 스펙트럴 간격이 없는 문제에서 저차원 ROM이 잘 안 되는 것은 코드 탓이 아니라 구조 탓이다.³

느린-빠른 계의 슬로우 매니폴드(페니첼 이론)와는 사촌 관계다. 저기서는 고유값이 정확히 0이 아니라 작은 파라미터 $\epsilon$ 만큼 느릴 뿐이고, 정상 쌍곡성이 있으면 근처에 진짜 불변다양체가 살아남는다. 형식은 다르지만 “느린 방향으로 환원한다”는 정신과 계산 절차(불변성 방정식 급수 풀이)는 거의 같다.

6. 수치적 근사[편집]

• 급수 전개. 가장 흔하다. 다중선형 형식을 자동 미분으로 뽑고, 차수별 선형계를 푼다. 임계 방향 때문에 계가 특이해지는 지점은 테두리 선형계로 처리한다.
• 그래프 변환. $h$ 를 격자 위 함수로 놓고 흐름을 짧게 전진시킨 뒤 다시 그래프로 사영하는 반복이다. 수축 사상이라 수렴하지만, $c$ 가 2를 넘어가면 격자 비용이 급격히 커진다.
• 매개화 방법(parameterization method). 다양체를 그래프가 아니라 매개화 $K(\theta)$ 로 두고 불변성 방정식 $\dot K = X(K)$ 를 푸는 방식이다. 그래프로 표현 못 하는 접힌 다양체까지 다룰 수 있고, 고차 계수를 안정적으로 뽑는 데 유리하다.
• 연속법과의 결합. 실무에서는 대개 수치 연속법 패키지가 분기점을 잡고, 그 자리에서 중심 다양체 축약과 정규형 계수를 한꺼번에 계산한다. 사용자가 보는 것은 “H(호프), LP(극한점), BT(보그다노프-타켄스)” 같은 라벨과 그 옆의 계수 한 줄뿐이지만, 그 안에서 돌아가는 것이 이 문서의 내용이다.⁴

7. 관련 문서[편집]

• 정규형 이론 · 분기 이론 · 호프 분기
• 자코비안 행렬 · 고유값 문제 · 안장점
• 축소차수모델 · 관성 다양체 · 편미분방정식
• 수치 연속법 · 호장법 · 자동 미분
• 극한 순환 · 카오스 이론 · 테일러 급수

8. Footnotes[편집]

1. $\dot\mu = 0$ 을 추가하는 순간 계의 차원이 늘어나는데도 문제가 쉬워진다는 점이 처음엔 킹받는다. 이유는 간단하다. 파라미터를 상태로 올리면 “파라미터마다 다양체를 새로 구한다”가 “다양체 하나를 한 번에 구한다”로 바뀌기 때문이다. ↩

2. 이 성질이 없었다면 중심 다양체 이론은 실용성이 없었을 것이다. “다양체가 무한히 많은데 어느 걸 쓰죠?”라는 질문의 답이 “아무거나 — 어차피 급수가 같다”라서 살아남은 이론이다. ↩

3. 난류처럼 스펙트럼이 촘촘하게 이어져 있는 문제에서 ROM 모드 수를 늘려도 오차가 잘 안 떨어지는 경험은 이 바닥의 국룰인데, 그 배경에는 “잘라낼 자리가 애초에 없다”는 구조적 사실이 있다. ↩

4. 그래서 초심자가 MatCont 출력에서 l1 = -0.0032 같은 숫자를 보고 “이게 뭐죠?”라고 묻는 일이 반복된다. 정직한 답은 “당신 계를 중심 다양체 위로 내리고 정규형으로 바꿨을 때 남은 3차 계수이고, 부호만 보면 된다”이다. ↩