지수 이분성

1. 개요[편집]

지수 이분성 Exponential dichotomy
대상	선형 비자율계 ẋ = A(t)x (또는 xn+1 = Anxn)
자료	사영족 P(t) + 상수 K, α > 0
상수계에서는	쌍곡성(Re λ ≠ 0)과 동치
핵심 정리	거칠기(roughness) — 작은 섭동에 살아남음
동치 조건	유계 비제차항 ↔ 유일한 유계해 (페론 조건)
주요 문헌	페론(1930) · 마세라-셰퍼(1966) · 코펠(1978)
수치적 의미	절단 경계값 문제의 조건수 = K/α

지수 이분성(exponential dichotomy)은 선형 비자율계 $\dot{\mathbf x} = A(t)\mathbf x$ 의 해공간이 매 시각마다 지수적으로 수축하는 부분공간과 지수적으로 팽창하는 부분공간으로 갈라지고, 그 갈라짐이 시간에 대해 균등한 상수로 통제되는 성질이다. 상수계에서 “고유값의 실수부가 0이 아니다”라고 말하던 쌍곡성을, 고유값이라는 개념 자체가 없는 시간의존 계로 옮긴 것이다.

왜 이게 필요한지는 안정 다양체 정리를 보면 바로 나온다. 평형점 근방의 랴푸노프-페롱 구성은 $e^{At}P^s$, $e^{At}P^u$ 의 지수 감쇠 부등식 두 줄만 쓰고, $A$ 가 상수라는 사실은 그 두 줄 이외의 어디에도 등장하지 않는다. 그러니 두 줄을 공리로 승격시키면 증명이 통째로 옮겨간다. 스턴버그 정리나 위상 공액의 비자율 판본, 그림자 보조정리, 연결 궤도의 수치 추적이 전부 같은 방식으로 이 개념을 경유한다.

이 문서는 정의와 그 동치 조건, 그리고 이 개념이 실제로 값어치를 하는 지점인 거칠기 정리와 수치 조건수에 집중한다. 평형점 주위의 불변 다양체 구성 자체는 안정 다양체 정리가, 사영 경계조건을 쓴 호모클리닉 연속화는 호모클리닉 궤도가 이미 다루므로 여기서 반복하지 않는다.

2. 정의[편집]

$A(t)$ 가 구간 $J$($\mathbb R$, $\mathbb R_+$, $\mathbb R_-$ 중 하나) 위에서 연속·유계이고 $\Phi(t)$ 가 기본해행렬($\Phi(0)=I$)이라 하자. 사영행렬 $P$($P^2=P$)와 상수 $K\ge1$, $\alpha>0$ 이 존재해 모든 $t,s\in J$ 에 대해

$$\bigl\lVert \Phi(t)P\,\Phi(s)^{-1} \bigr\rVert \le K e^{-\alpha(t-s)}\quad (t\ge s), \qquad \bigl\lVert \Phi(t)(I-P)\,\Phi(s)^{-1} \bigr\rVert \le K e^{-\alpha(s-t)}\quad (t\le s)$$

가 성립하면 이 계가 $J$ 위에서 지수 이분적이라 한다. $P(t) := \Phi(t)P\Phi(t)^{-1}$ 를 두면 이것이 흐름과 교환하는 사영족이 되고, $\mathcal S(t) = \operatorname{im}P(t)$ 는 미래로 지수 감쇠하는 해들의 초기값, $\mathcal U(t) = \ker P(t)$ 는 과거로 지수 감쇠하는 해들의 초기값이다. 매 시각 $\mathbb R^n = \mathcal S(t)\oplus\mathcal U(t)$ 이고, 이 분해는 흐름을 따라 옮겨간다.

읽을 때 놓치기 쉬운 것 셋.

• $K$ 와 $\alpha$ 가 $t,s$ 에 무관하다. 이 균등성이 정의의 전부다. 각 해가 개별적으로 지수 감쇠한다는 것과는 전혀 다른 요구다(아래 §랴푸노프 지수와 다른 점).
• 부등식이 두 방향으로 쓰여 있다. 안정 쪽은 미래로, 불안정 쪽은 과거로 감쇠한다. 그래서 $\Phi(s)^{-1}$ 가 등장하는데, 불안정 성분을 과거로 되돌려도 폭발하지 않는다는 뜻이지 계를 역방향으로 풀라는 뜻이 아니다.
• $P$ 의 계수(rank)는 구간에 따라 달라질 수 있다. 반직선 $\mathbb R_+$ 위에서는 $\mathcal S$ 만 유일하게 결정되고 $\mathcal U$ 는 임의의 보충공간으로 골라도 되지만, $\mathbb R$ 전체에서는 양쪽이 모두 유일하다. 연결 궤도 계산에서 두 반직선의 이분성을 따로 잡은 다음 그 교집합의 차원을 세는 논법이 여기서 나온다.

이산 판본도 문자 그대로 같다. $x_{n+1}=A_n x_n$ 에서 전이행렬 $\Phi(n,m)$ 을 놓고 $\lVert\Phi(n,m)P_m\rVert\le K\lambda^{n-m}$, $\lVert\Phi(n,m)(I-P_m)\rVert\le K\lambda^{m-n}$ ($0<\lambda<1$)로 쓴다. 수치 궤도는 어차피 이산이므로 실무에서 실제로 쓰이는 것은 이쪽이다.

3. 상수계·주기계로 되돌리면[편집]

$A(t)\equiv A$ 이면 $\Phi(t)=e^{At}$ 이고, 위 두 부등식은 정확히 ”$A$ 의 스펙트럼이 허수축과 만나지 않는다”와 동치다. 즉 쌍곡성의 완전한 일반화다. $P$ 는 안정 스펙트럼 사영이고 $\alpha$ 는 스펙트럼 간극이며, $K$ 는 고유벡터 기저의 조건수를 반영한다. 비정규 행렬에서 $K$ 가 커진다는 사실이 바로 뒤에서 볼 “거칠기 문턱이 $K^2$ 로 나빠진다”의 정체다.

주기계 $A(t+T)=A(t)$ 이면 플로케 이론이 답을 준다. 모노드로미 행렬 $M=\Phi(T)$ 의 고유값(플로케 승수)이 단위원 위에 하나도 없을 때, 그리고 오직 그때만 지수 이분적이다. 극한주기해 위의 변분방정식은 접선 방향 승수가 항상 1이라 이분적일 수 없고, 그래서 그쪽은 이분성 대신 삼분성(중심 방향을 따로 떼는 exponential trichotomy)으로 다룬다.

4. 그린 함수와 유계해 — 동치 조건[편집]

정의의 부등식 두 줄은 사실 훨씬 쓰기 좋은 형태로 압축된다. $\mathbb R$ 위의 이분성이 있으면

$$G(t,s) = \begin{cases} \ \ \Phi(t)P\Phi(s)^{-1}, & t \ge s\\[2pt] -\Phi(t)(I-P)\Phi(s)^{-1}, & t < s \end{cases}$$

가 $\lVert G(t,s)\rVert \le Ke^{-\alpha\lvert t-s\rvert}$ 를 만족하는 **그린 함수**가 되고, 비제차 방정식 $\dot{\mathbf x}=A(t)\mathbf x+\mathbf f(t)$ 의 유계해가

$$\mathbf x(t) = \int_{-\infty}^{\infty} G(t,s)\,\mathbf f(s)\,ds$$

로 한 방에 나온다. 적분이 수렴하는 이유는 커널이 대각선에서 지수적으로 죽기 때문이고, $\lVert \mathbf x\rVert_\infty \le (2K/\alpha)\lVert\mathbf f\rVert_\infty$ 라는 사전 평가가 공짜로 딸려 온다.

놀라운 것은 역도 성립한다는 점이다(페론 조건, 또는 admissibility).

$\mathbb R$ 위에서: 임의의 유계 연속 $\mathbf f$ 에 대해 $\dot{\mathbf x}=A(t)\mathbf x+\mathbf f(t)$ 가 유계해를 정확히 하나 가진다 ⟺ 지수 이분적이다.

즉 지수 이분성은 미분작용소 $L=\frac{d}{dt}-A(t)$ 가 유계함수 공간 사이의 동형사상이라는 말과 같은 뜻이다. 이 재해석이 결정적인 이유는, 부등식 두 줄을 확인하는 대신 선형 작용소 하나가 가역인지만 보면 되고, 가역성은 바나흐 공간 위의 섭동론이 통째로 적용되는 성질이기 때문이다. 아래의 거칠기 정리도, 멜니코프 함수가 왜 하필 그 형태인지도 여기서 나온다. $L$ 이 전사지만 단사가 아니면 유계해 공간의 차원이 $W^u$ 와 $W^s$ 의 초과 교차 차원이 되고, 그 여핵이 멜니코프 방법의 적분에 등장하는 수반해다.

5. 거칠기 — 이 개념의 존재 이유[편집]

$$\dot{\mathbf x} = \bigl(A(t)+B(t)\bigr)\mathbf x, \qquad \sup_t \lVert B(t)\rVert \le \delta$$

거칠기 정리(roughness, 코펠 1978)는 $\delta$ 가 충분히 작으면 섭동계도 지수 이분적이고, 지수는 $\alpha-O(K\delta)$ 로만 나빠지며, 사영족도 $O(K^2\delta/\alpha)$ 만큼만 움직인다고 말한다. 문턱은 판본마다 상수가 다르지만 구조는 항상 같다 — $\delta \lesssim \alpha/K^2$. 여기서 $\alpha$ 는 한 번, $K$ 는 두 번 등장한다는 점이 중요하다. 지수 간극이 좁거나 안정·불안정 부분공간의 사잇각이 작으면($K$ 가 크면) 이분성은 아주 쉽게 깨진다.¹

이 정리가 없으면 수치해석에서 할 말이 사라진다.

• 이산화 오차. 수치 적분기가 푸는 것은 $A(t)$ 가 아니라 $A(t)+O(h^p)$ 다. 거칠기가 있어야 “계산된 궤도의 안정 다양체”라는 표현이 의미를 가진다.
• 그림자 정리. 유사궤도(pseudo-orbit)의 변분방정식이 이산 지수 이분성을 가지면 근처에 진짜 궤도가 있고, 그림자 거리는 $K/\alpha$ 곱하기 국소 오차로 평가된다. 팔머의 판본이 정확히 이 형태다.
• 모델 오차. 물리 모형을 조금 고쳐도 쌍곡 구조가 유지된다는 구조적 안정성의 선형 부분이 곧 거칠기다.
• 비자율 하트만-그로브만. 팔머(1973)는 $\dot{\mathbf x}=A(t)\mathbf x+\mathbf g(t,\mathbf x)$ 에서 $A(t)$ 가 이분적이고 $\mathbf g$ 가 작으면 선형계와 위상 공액임을 보였다. 상수계용 하트만-그로브만 정리의 증명에서 “쌍곡”을 “이분”으로 바꿔 쓴 것이 전부다.

6. 랴푸노프 지수와 다른 점 — 균등성이 전부다[편집]

여기서 사고가 가장 많이 난다. 모든 랴푸노프 지수가 0이 아닌 것과 지수 이분성은 같은 말이 아니다. 전자가 훨씬 약하다.

랴푸노프 지수는 $\lim_{t\to\infty}\frac1t\ln\lVert\Phi(t)\mathbf v\rVert$ 라는 점근적 양이라, 유한 시간 구간에서 무슨 짓을 하든 상관하지 않는다. 반면 이분성은 모든 $t\ge s$ 쌍에 대해 같은 $K$ 로 통제할 것을 요구한다. 페론이 만든 고전적 반례가 이 간극을 보여준다 — 모든 지수가 음수인데도 비선형 항을 조금만 얹으면 해가 발산하는 계를 만들 수 있다. 이유는 뻔하다. 지수가 음수라도 도중에 $\lVert\Phi(t)\Phi(s)^{-1}\rVert$ 가 $e^{100}$ 까지 치솟았다 내려오면, 그 순간 비선형 항이 증폭돼 판이 뒤집힌다.

정리하면 이렇다.

개념	통제하는 것	섭동에 강한가
랴푸노프 지수	점근 성장률(점)	아니다(불연속적으로 변할 수 있다)
새커-셀 이분성 스펙트럼	이분성이 깨지는 λ 집합(구간)	그렇다(상반연속)
지수 이분성	균등 지수 부등식	그렇다(거칠기)

가운데 줄이 둘을 잇는다. **새커-셀 스펙트럼**은 $\dot{\mathbf x}=(A(t)-\lambda I)\mathbf x$ 가 지수 이분적이지 않은 $\lambda$ 들의 집합으로 정의되고, 최대 $n$ 개의 서로 떨어진 닫힌 구간의 합집합이 된다는 것이 정리다. 이때 원래 계가 이분적이라는 것은 $0\notin$ 스펙트럼이라는 것과 같다. 랴푸노프 지수들은 이 구간들 안에 들어가지만, 구간이 점으로 찌그러지지는 않는다 — 시간 평균이 수렴하지 않는 계에서 구간의 폭이 그대로 남는다. 페신 이론이 다루는 비균등 쌍곡성은 이 구간이 0을 포함하는데도 거의 모든 궤도의 지수가 0이 아닌 상황이며, 다양체는 여전히 존재하지만 그 크기가 궤도를 따라 지수적으로 쪼그라들어 균등한 상수가 없다. 수치적으로 다양체를 잡으려는 순간 이 차이가 곧바로 고통으로 돌아온다.

7. 수치적으로[편집]

1. 조건수가 곧 $K/\alpha$ 다. 무한 구간 문제를 $[-T,T]$ 로 자르고 양 끝에 사영 경계조건을 걸면 유한차원 선형계가 나오는데, 그 계의 조건수 상계가 그린 함수 노름, 즉 $K/\alpha$ 로 표현된다. 그래서 지수 이분성은 “안정 다양체가 존재한다”는 정성적 진술이 아니라 미리 계산 가능한 오차 예산이다. 반대로 $\alpha$ 가 작은 문제(느린 방향이 섞인 계)에서 다중 슈팅법이 왜 갑자기 발산하는지도 같은 식이 설명한다. 절단 오차 자체는 $T$ 에 대해 $O(e^{-\alpha T})$ 로 지수적으로 작으므로, 구간을 조금만 늘려도 이득이 크다.

2. 사영을 실제로 어떻게 구하나. 정공법은 $\Phi(t)$ 를 직접 적분한 뒤 분해하는 것인데, $\Phi$ 자체가 $e^{\alpha T}$ 규모로 폭발하므로 유한정밀도에서 즉사한다. 표준 처방은 기본해행렬을 통째로 들고 다니지 않고 **연속 QR 분해**로 가는 것이다. $\Phi(t)=Q(t)R(t)$ 로 두고 $Q$ 에 대한 미분방정식(직교군 위의 흐름)과 $R$ 의 대각성분에 대한 방정식을 따로 적분하면, 지수적 스케일 차이가 $R$ 의 대각성분에 로그로 흡수돼 $Q$ 는 얌전히 유계로 남는다. 디에치와 판블렉의 연속·이산 QR 알고리즘이 이 방식으로 새커-셀 스펙트럼 구간의 상하한을 직접 추정한다. 실무 코드가 그람-슈미트 재직교화를 주기적으로 끼워 넣는 것도 같은 이유다.

3. 이분성 판정은 스펙트럼 구간 문제다. “이 궤도가 쌍곡인가”를 수치로 묻는 것은 고유값 하나를 보는 일이 아니라 구간이 0을 덮는지 보는 일이다. 유한 시간 QR 추정치는 진짜 구간의 상하한을 주므로, 상한이 음수이거나 하한이 양수인 방향으로 딱 갈리면 안전하게 이분적이라 말할 수 있고, 구간이 0을 물고 있으면 판정 보류가 정답이다. 여기서 “랴푸노프 지수가 $-0.003$ 이니까 안정하다”고 선언하는 것이 왜 위험한지가 드러난다.

4. 어디서 쓰이나. 절단 구간 경계값 문제로 호모클리닉·헤테로클리닉 궤도를 잡는 수치 연속법 코드(AUTO, MatCont의 HomCont), 이동파 해의 스펙트럼 안정성 판정에 쓰는 에번스 함수(그 정의 자체가 양 반직선 이분성 부분공간의 행렬식이다), 유한시간 국소 다양체(LCS) 추출, 그리고 데이터 기반 동역학에서 유한시간 리아푸노프 지수를 해석할 때의 경고문까지 — 전부 같은 부등식 두 줄 위에 서 있다.²

8. 무엇을 포기하고 무엇을 얻는가[편집]

지수 이분성은 공짜가 아니다. 균등성을 요구하는 대가로 적용 범위가 좁아진다. 대부분의 카오스 끌개는 균등 쌍곡이 아니고(로렌츠 끌개조차 고전적 의미로는 아니다), 그래서 “수치 궤도가 진짜 궤도를 그림자한다”는 주장을 로렌츠계에 무조건 적용할 수는 없다. 실제로 성립하는 것은 유한 시간·유한 정밀도 판본이며, 그 조건을 구간마다 확인하는 것이 정직한 태도다.

대신 얻는 것이 크다. 부등식 두 줄이면 평형점도, 주기궤도도, 아무 기준궤도도, 심지어 이산 수치궤도도 같은 문장으로 처리된다. 상수계의 행렬 지수함수와 리아푸노프 방정식이 하던 일을, 스펙트럼이라는 단어 없이 해내는 것. 동역학계 이론에서 “비자율로 확장한다”는 말이 나오면 열에 아홉은 이 개념을 끼워 넣었다는 뜻이다.³

9. 관련 문서[편집]

• 안정 다양체 정리 · 불변 다양체 · 중심 다양체 정리
• 하트만-그로브만 정리 · 스턴버그 정리 · 위상 공액
• 플로케 이론 · 랴푸노프 지수 · 구조적 안정성
• 호모클리닉 궤도 · 멜니코프 방법 · 수치 연속법
• 그린 함수 · 바나흐 공간 · QR 분해 · 그람-슈미트
• 조건수 · 행렬 지수함수 · 리아푸노프 방정식 · 동역학계

10. Footnotes[편집]

1. $K$ 가 제곱으로 들어오는 게 억울해 보이지만 이유는 단순하다. 섭동된 사영을 만들 때 안정 쪽 추정과 불안정 쪽 추정을 한 번씩 곱해 쓰기 때문이다. 실무적 교훈은 하나다 — 스펙트럼 간극만 보지 말고 부분공간 사잇각도 같이 보라. 간극이 넉넉해도 사잇각이 1도면 $K$ 가 60을 넘고, 그러면 문턱이 3600배 나빠진다. ↩

2. 유한시간 리아푸노프 지수(FTLE) 그림을 예쁘게 뽑아 놓고 “여기가 쌍곡 구조입니다”라고 말하기 전에, 그 값이 적분 시간에 얼마나 민감한지 한 번은 찍어 봐야 한다. 적분 시간을 두 배로 늘렸는데 그림이 달라진다면 그건 구조가 아니라 유한시간 효과다. 이 검사를 통과하지 못한 FTLE 그림이 학회 포스터의 절반쯤 된다는 것이 이 바닥의 공공연한 비밀. ↩

3. 페론이 1930년에 이 아이디어의 씨앗을 심고, 마세라와 셰퍼가 1950~60년대에 함수공간 언어로 다시 쓰고, 코펠이 1978년 강의노트에서 “Dichotomies in Stability Theory”라는 제목으로 정리해 지금의 표준형이 됐다. 90쪽짜리 얇은 노트인데 인용 수만 놓고 보면 웬만한 두꺼운 교과서를 다 이긴다. 얇을수록 많이 인용되는 것은 이 바닥의 국룰이다. ↩