로렌츠 방정식

계산물리 물리 수치해석 마지막 수정: 2026-08-10 04:14:27

1. 개요[편집]

로렌츠 방정식 Lorenz equations (Lorenz-63)
제안	Edward N. Lorenz, 1963
차원	자율 3차원 상미분방정식
고전 파라미터	σ = 10, β = 8/3, ρ = 28
대칭	(x, y, z) → (−x, −y, z)
발산	∇·f = −(σ + 1 + β) = −13.667 (상수)
끌개 차원	D_KY ≈ 2.06
엄밀 증명	W. Tucker, 2002 (스메일의 14번 문제)

로렌츠 방정식(Lorenz equations)은 에드워드 로렌츠가 1963년 논문 Deterministic Nonperiodic Flow에서 제시한 자율 3차원 상미분방정식계로, 결정론적인데도 비주기적으로 영원히 방황하는 해가 존재한다는 것을 처음으로 명확히 보여준 모형이다.

\dot{x} = \sigma(y - x), \qquad \dot{y} = x(\rho - z) - y, \qquad \dot{z} = xy - \beta z

우변에 비선형이라고는 곱 항 $xz$ , $xy$ 딱 두 개뿐이다. 이 초라한 방정식이 카오스 이론의 표준 실험동물이 된 이유는, 카오스가 “복잡한 계에서 생기는 현상”이 아니라 비선형성과 3차원만 있으면 생기는 기본값임을 가장 싸게 증명해 주기 때문이다.¹

이 문서는 방정식 자체의 분기 구조와 끌개의 기하에 집중한다. 예측 가능성 지평이나 랴푸노프 스펙트럼의 계산법은 랴푸노프 지수, 카오스의 일반 정의와 카오스로 가는 경로는 카오스 이론, 상공간 언어 일반은 동역학계에 있다.

2. 어디서 나왔나[편집]

로렌츠는 방정식을 하늘에서 받아온 것이 아니다. 출발점은 아래에서 데우고 위에서 식히는 유체층의 대류(레일리-베나르 문제)를 부시네스크 근사로 쓴 편미분방정식이고, 살츠먼(Barry Saltzman, 1962)이 이를 푸리에 모드로 전개해 놓은 것을 로렌츠가 모드 세 개만 남기고 전부 잘라낸 것이 로렌츠-63이다. 변수의 의미는 각각 이렇다.

$x$ — 대류 롤의 회전 세기(부호는 회전 방향).
$y$ — 상승류와 하강류 사이의 온도차.
$z$ — 연직 온도 분포가 선형에서 벗어난 정도.

파라미터도 물리량이다. $\sigma$ 는 프란틀 수, $\rho$ 는 레일리 수를 임계 레일리 수로 나눈 값, $\beta = 4/(1+a^2)$ 는 대류 셀의 종횡비에서 오는 기하 인자다. 고전값 $\sigma=10$ 은 물에 가까운 프란틀 수이고, $\beta = 8/3$ 은 $a^2 = 1/2$ 일 때 나온다.

여기서 오해를 하나 정리하자. 로렌츠-63은 대류의 모형이 아니다. 모드 세 개짜리 절단은 $\rho \approx 28$ 처럼 임계값의 28배나 되는 영역에서는 물리적으로 전혀 정당화되지 않는다. 로렌츠 본인도 이 점을 알고 있었고, 논문의 주장은 “대류가 이렇게 움직인다”가 아니라 “절단된 결정론적 방정식조차 비주기해를 낳는다”였다. 이 모형의 값어치는 기상 예측이 아니라 예측 가능성의 한계를 논증하는 데 있다.

3. 고정점과 분기 사슬[편집]

$\sigma, \beta$ 를 고정하고 $\rho$ 만 올리면 계의 운명이 순서대로 바뀐다. 먼저 고정점을 구한다. $\dot x = 0$ 에서 $y = x$ , 이를 $\dot y = 0$ 에 넣으면 $x(\rho - 1 - z) = 0$ 이고, $\dot z = 0$ 에서 $z = x^2/\beta$ 다. 따라서

원점 $O = (0,0,0)$ — 항상 존재. 대류가 없는 순수 열전도 상태.
$\rho > 1$ 일 때 추가로 $C^{\pm} = \left(\pm\sqrt{\beta(\rho-1)},\ \pm\sqrt{\beta(\rho-1)},\ \rho-1\right)$ — 시계/반시계 방향으로 도는 정상 대류 롤 한 쌍.

원점에서 선형화하면 $z$ 방향이 분리되어 고유값 $-\beta$ 를 주고, $(x,y)$ 블록은 $\lambda^2 + (\sigma+1)\lambda + \sigma(1-\rho) = 0$ 을 만족한다. $\rho < 1$ 이면 세 고유값이 모두 음수이므로 원점이 전역 안정이고(랴푸노프 함수 $V = x^2 + \sigma y^2 + \sigma z^2$ 로 전역성까지 보인다), $\rho = 1$ 에서 고유값 하나가 0을 지나며 초임계 갈래 분기가 일어나 $C^\pm$ 가 태어난다. $\rho > 1$ 이면 원점은 1차원 불안정 다양체를 가진 안장이 된다. $\rho = 28$ 에서의 고유값은 $\{11.83,\ -2.667,\ -22.83\}$ 이다.

$C^\pm$ 의 안정성이 이 계의 하이라이트다. 특성방정식을 라우스-후르비츠로 따지면 안정 조건이 딱 하나의 임계값으로 정리된다.

\rho_H = \frac{\sigma(\sigma + \beta + 3)}{\sigma - \beta - 1} = \frac{10 \times (10 + 8/3 + 3)}{10 - 8/3 - 1} = 24.7368\ldots

$1 < \rho < \rho_H$ 에서 $C^\pm$ 는 안정한 초점이고, $\rho = \rho_H$ 에서 켤레 고유값 쌍이 허수축을 가로지르며 호프 분기가 일어난다. 그런데 이 호프는 아임계다. 즉 $\rho < \rho_H$ 쪽에 불안정 주기궤도가 존재하다가 $\rho$ 가 올라가면서 $C^\pm$ 쪽으로 수축해 들어와 이들을 삼켜 버린다. 결과적으로 $\rho > \rho_H$ 에서는 세 고정점이 전부 불안정하고, 궤도는 유계 영역 안에 갇혀 있는데 갈 곳이 없다. 남은 선택지는 하나 — 끌개가 점도 순환도 아닌 무언가여야 한다.²

정리하면 $\sigma=10,\ \beta=8/3$ 에서 $\rho$ 에 따른 시나리오는 이렇다.

$\rho$ 구간	상태
$\rho < 1$	원점만 존재, 전역 안정 (대류 없음)
$1 < \rho \lesssim 13.926$	$C^\pm$ 안정, 궤도는 둘 중 하나로 수렴
$\rho \approx 13.926$	원점의 호모클리닉 폭발 — 카오스 안장 탄생
$13.926 \lesssim \rho \lesssim 24.06$	$C^\pm$ 안정 + 과도 카오스(긴 방황 후 수렴)
$24.06 \lesssim \rho < 24.7368$	카오스 끌개와 안정한 $C^\pm$ 가 공존(쌍안정)
$\rho > 24.7368$	아임계 호프 이후 $C^\pm$ 도 불안정, 나비 끌개만 남음

$\rho \approx 13.926$ 에서 원점의 불안정 다양체 두 갈래가 정확히 원점으로 되돌아오는 대칭 호모클리닉 쌍이 생기고, 그 순간 무한히 많은 주기궤도가 한꺼번에 태어난다. 이 사건의 기하는 호모클리닉 궤도 문서가 다룬다. 참고로 $\rho$ 를 더 키운다고 카오스가 영원히 이어지지도 않는다 — $\rho \approx 100$ 근처를 비롯해 여러 곳에 주기 창이 있고, $\rho$ 가 아주 커지면 전역적으로 안정한 주기궤도로 정착한다(Sparrow, 1982).

4. 위상부피는 지수적으로 사라진다[편집]

로렌츠 끌개가 왜 “이상한(strange)” 물건인지를 가장 짧게 보여주는 계산이 벡터장의 발산이다.

\nabla \cdot \mathbf{f} = \frac{\partial}{\partial x}\big[\sigma(y-x)\big] + \frac{\partial}{\partial y}\big[x(\rho-z)-y\big] + \frac{\partial}{\partial z}\big[xy - \beta z\big] = -\sigma - 1 - \beta

$\sigma=10,\ \beta=8/3$ 이면 $-13.6667$ 이다. 중요한 것은 이 값이 상태와 무관한 상수라는 점이다. 리우빌 정리의 소산계 버전에 따라 상공간 영역 $V(t)$ 의 부피는

\frac{dV}{dt} = \int_V \nabla\cdot\mathbf{f}\,dV = -(\sigma+1+\beta)\,V \quad\Longrightarrow\quad V(t) = V(0)\,e^{-(\sigma+1+\beta)t}

로 균일한 비율로 지수 수축한다. 시간 1이 지날 때마다 부피가 $e^{-13.67} \approx 1.2\times10^{-6}$ 배가 된다. 초기 부피가 아무리 커도 $t\to\infty$ 에서 0이므로, 끌개는 3차원 르베그 측도가 0인 집합이다. 동시에 궤도는 유계 영역(부등식 $\rho x^2 + \sigma y^2 + \sigma(z-2\rho)^2 \le C$ 꼴의 흡수 타원체)을 벗어나지 못하므로 한 점으로 도망갈 수도 없다.

여기서 자연스러운 오해가 나온다. “그럼 끌개는 두께 없는 면이겠네?” 아니다. 나비 그림을 보면 날개가 매끄러운 2차원 면처럼 보이지만, 진짜 면이라면 두 날개가 만나는 안쪽에서 궤도가 서로 교차해야 한다 — 결정론적 자율계에서 해의 유일성이 금지하는 일이다. 로렌츠 본인이 1963년 논문에서 이 모순을 지적하며 끌개를 “무한히 많은 면의 복합체(infinite complex of surfaces)“라고 불렀다. 날개는 하나의 면이 아니라, 단면을 확대할 때마다 층이 갈라져 나오는 칸토어 집합 × 면 구조다.

수치로 확인하면 랴푸노프 지수 스펙트럼이 $(\lambda_1,\lambda_2,\lambda_3) \approx (0.906,\ 0,\ -14.572)$ 이고 — 세 값의 합이 위에서 계산한 $-13.667$ 과 정확히 맞아야 한다 —, 카플란-요크 차원은

D_{KY} = 2 + \frac{\lambda_1}{|\lambda_3|} \approx 2 + \frac{0.906}{14.572} \approx 2.062

다. 2보다 크니 면이 아니고, 3보다 작으니 부피도 아니다. 부피 0과 차원 2.06을 동시에 만족하는 물건, 그게 프랙탈 끌개다.

5. 로렌츠가 손으로 만든 1차원 사상[편집]

1963년 논문에서 진짜 놀라운 대목은 그림이 아니라 표다. 로렌츠는 궤도가 $z$ 의 극대값을 찍을 때마다 그 값 $M_n$ 을 기록해 $M_{n+1}$ 대 $M_n$ 을 찍어 보았고, 점들이 뾰족한 봉우리를 가진 1차원 곡선 위에 거의 정확히 얹히는 것을 발견했다. 이 곡선의 기울기 절댓값이 어디서나 1보다 크다는 것이 관찰되면, 안정한 주기궤도가 존재할 수 없다는 결론이 바로 따라 나온다 — 어떤 주기점이든 승수의 절댓값이 1을 넘기 때문이다.

이것이 지금 로렌츠 사상이라 불리는 물건이고, 3차원 흐름을 푸앵카레 단면으로 잘라 이산 사상으로 환원하는 표준 기법의 원형이다. 왜 2차원 단면이 아니라 1차원 곡선으로 뭉개지는가? 위 문단의 수축률 때문이다. 횡방향 수축이 $e^{-14.57}$ 규모로 압도적이라 단면상의 점구름이 사실상 두께 0인 곡선으로 보이는 것이다. 물론 진짜로 두께가 0은 아니고, 확대하면 칸토어 층이 나온다. 로렌츠가 본 “1차원 사상”은 훌륭한 근사였지 정리가 아니었다는 뜻이다.

6. 그래서 그건 진짜로 존재하는가 — 터커의 증명[편집]

여기가 이 문서에서 가장 흥미로운 지점이다. 1963년부터 1990년대까지, 로렌츠 끌개에 대해 인류가 가진 것은 그림과 수치 실험뿐이었다. 우리가 보는 나비가 진짜 카오스 끌개인지, 아니면 주기가 어마어마하게 긴 안정 주기궤도를 유한 정밀도로 잘못 본 것인지 아무도 증명하지 못했다. 카오스계의 수치 궤적은 정의상 지수적으로 오차가 자라므로, “돌려봤더니 나오더라”는 증명이 될 수 없다.

1970년대에 구켄하이머와 윌리엄스는 기하학적 로렌츠 모형을 도입했다. 실제 방정식 대신, 로렌츠 사상이 가진 성질을 공리로 가정한 추상적 흐름을 만들고 그 위에서 카오스를 엄밀히 증명한 것이다. 이 우회로 덕분에 “로렌츠형 끌개”의 이론은 완성되었지만, 원래 방정식이 그 공리를 만족하는가는 여전히 미해결이었다. 스티븐 스메일이 1998년에 꼽은 21세기 수학 문제 18개 중 14번이 바로 이것이다.

2002년 워윅 터커가 이를 해결했다. 전략은 두 갈래를 붙이는 것이다.

원점 근처는 정규형(normal form)으로 해석적으로 처리한다. 안장이 있어서 통과 시간이 무한대로 발산하는 구간이라, 수치 적분으로는 손댈 수 없기 때문이다.
나머지 영역은 구간 해석으로 적분한다. 초기 조건을 점이 아니라 작은 상자로 놓고, 반올림 오차와 절단 오차까지 포함한 상자를 전진시켜 “참 해가 반드시 이 상자 안에 있다”는 포함 관계를 컴퓨터가 유지하게 한다. 상자가 커지면 쪼개고, 되돌아오는 상자들이 출발한 단면을 덮는지를 확인한다.

이 계산으로 로렌츠 사상에 해당하는 되돌아옴 사상이 실제로 팽창적이고 불변 원뿔장을 갖는다는 것, 즉 로렌츠 끌개가 강건한(robust) 특이 쌍곡 끌개이며 SRB 측도를 지지한다는 것이 증명되었다. 요점은 이것이 “컴퓨터가 그림을 예쁘게 그렸다”가 아니라 부동소수점의 모든 오차를 상계로 감싼 논리적 증명이라는 데 있다. 같은 정신의 앞선 성과로 미샤이코프와 므로제크가 1995년 코늘리 지표를 이용해 로렌츠계에 이동사상이 심겨 있음을 컴퓨터 보조로 증명한 것이 있다.

7. 수치적으로 다룰 때[편집]

적분기. RK4에 $\Delta t = 0.005 \sim 0.01$ 이면 충분하다. 이 계는 강성이 아니라서 음함수 적분기를 쓸 이유가 없다. 다만 저차 적분기(전진 오일러 등)는 인공 소산이 커서 $\lambda_1$ 을 실제보다 작게 만들거나 끌개를 부풀린다.
궤적을 믿지 마라. 배정밀도로 적분해도 신뢰 가능한 궤적은 시간 단위 40 남짓이다. 그 이후의 궤적은 내가 넣은 초기값의 답이 아니다. 반대로 통계량(끌개 모양, 차원, $\lambda$ , 시간 평균)은 그림자 정리 덕분에 신뢰할 수 있다.
좋은 단위 테스트. 랴푸노프 스펙트럼 합이 $-(\sigma+1+\beta)$ 와 안 맞으면 야코비안 구현이나 재직교화 간격이 틀린 것이다. 검산이 이렇게 친절한 계는 흔치 않다.
대칭 활용. 계가 $(x,y,z)\to(-x,-y,z)$ 에 불변이므로 코드가 이 대칭을 깨면(예: 비대칭 오차 누적) 버그 신호일 수 있다.

8. 오해와 여담[편집]

나비 효과와 나비 모양은 무관하다. 초기조건 민감성을 뜻하는 나비 효과라는 표현은 1972년 강연 제목에서 왔고, 끌개의 날개 모양은 그 뒤에 붙은 우연한 시각적 일치다. 로렌츠가 원래 쓴 비유는 갈매기였다는 이야기가 있다.
“로렌츠 끌개는 두 개의 소용돌이”가 아니다. 궤도는 한쪽 날개를 몇 바퀴 돌다 반대쪽으로 넘어가는데, 몇 바퀴 도는지의 수열이 사실상 난수다. 이 수열이 앞 절의 1차원 사상의 기호 동역학이다.
로렌츠-96은 다른 물건이다. 1996년에 로렌츠가 따로 만든 $N$ 변수 순환 모형으로, 자료동화·앙상블 예보 알고리즘의 시험대로 쓰인다. 이름이 비슷해서 자주 혼동된다.
회로로도 돈다. 연산증폭기와 아날로그 곱셈기로 로렌츠계를 그대로 배선하면 오실로스코프에 나비가 뜬다(쿠오모-오펜하임, 1993). 카오스가 소프트웨어 버그가 아니라는 가장 저렴한 증거다.³

9. 관련 문서[편집]

그래서 이 방정식은 “3줄짜리 코드로 재현 가능한 유일한 유명 물리”라는 지위를 누린다. 학부 수치해석 과제로 나오고, 논문 그림으로 나오고, 티셔츠로도 나온다. 로렌츠 본인은 1963년 논문이 20년 가까이 거의 인용되지 않는 것을 지켜봐야 했다. 기상학 저널에 실린 순수 수학 논문의 숙명이었다. ↩
아임계 호프라는 사실이 실무적으로 잔인한 이유는 호프 분기 문서에 정리되어 있다. 요약하면 임계점을 넘는 순간 진폭이 0에서 살살 자라는 것이 아니라 궤도가 멀리 있는 다른 끌개로 점프한다. 로렌츠계에서는 그 “다른 끌개”가 하필 카오스였다. ↩
반대로 말하면, 아날로그 회로 두 대를 똑같이 만들어도 두 오실로스코프의 파형은 몇 초 만에 갈라진다. 부품 공차가 초기조건 오차 역할을 하기 때문이다. 이걸 역이용해 카오스 동기화로 비밀통신을 하겠다는 연구가 1990년대에 유행했고, 지금은 조용하다. ↩