로렌츠 끌개: 인터랙티브 탐색

이 시각화에서 살펴보는 로렌츠 끌개

로렌츠 방정식은 에드워드 로렌츠가 1963년 대류 모형에서 이끌어 낸 세 상미분방정식입니다[1]. 이 페이지는 $\sigma = 10$ 와 $\beta = 8/3$ 을 고정하고 가열 정도 $\rho$ 를 조절하게 합니다. $\rho = 28$ 에서는 거의 모든 해가 같은 유계의 나비 모양 집합을 끝없이 감돌며 정착하거나 반복하지 않습니다. 이 집합이 로렌츠 끌개입니다.

\dot x = \sigma\,(y - x),\qquad \dot y = x\,(\rho - z) - y,\qquad \dot z = xy - \beta z,\qquad \sigma = 10,\ \beta = \tfrac83

각 그림이 보여 주는 것

큰 그림인 3D 위상 공간에서 각 점은 상태 $(x, y, z)$ 입니다. 옅은 회색은 $(1, 1, 1)$ 에서 출발한 하나의 긴 궤도이고, 움직이는 점과 꼬리는 기준 궤도이며, 파란 고리는 출발 기준점입니다. 드래그해 시점을 회전하고, 꼬리를 눌러 다음 점 구름의 출발 위치를 고르세요.

빨간 출발 버튼은 기준점에서 각기 다른 방향으로 정확히 $\varepsilon$ 만큼 떨어진 시작점 구름을 놓고 모두 따라갑니다. 점의 색은 기준 궤도에서 얼마나 멀어졌는지 나타냅니다. 이 색과 빛, 잔상, 먼 쪽의 어둡기, 출발 전에 흐르는 옅은 입자는 모두 시각적 표현일 뿐입니다.

패널 위쪽 화면의 위에는 기준 궤도와 $\varepsilon$ 만큼 떨어져 시작한 궤도의 $x(t)$ , 즉 두 ‘예측’이 있습니다. 노란 선은 두 $x$ 의 차이가 처음 1을 넘은 때를 표시하며, 0.005 시간 단위마다 확인합니다(‘예측이 갈라졌다’를 나타내기 위해 여기서 정한 기준이며, 계 자체에 내재한 한계는 아닙니다). 아래의 $\frac1m \sum_i \ln\bigl(d_i(t)/d_i(0)\bigr)$ 는 각 점의 거리 증가 배율의 로그를 평균한 값입니다(축은 배율을 표시합니다). 띠는 점들의 10백분위수부터 90백분위수까지이며, 읽음값의 ‘대표 거리’는 기하평균입니다. $\rho = 28$ 에서는 거리가 끌개 크기에 가까워질 때까지 점선으로 장기 평균 기울기 $\lambda \approx 0.906$ 을 보여 줍니다.

다시 출발은 같은 점, 같은 방향, 같은 난수 시드에서 현재 $\varepsilon$ 로 다시 출발합니다. 이전 곡선은 비교할 수 있도록 옅게 남습니다. 이렇게 다시 출발해야 두 $\varepsilon$ 값을 공정하게 비교할 수 있습니다.

왜 이렇게 움직일까요

늘이기와 접기: 끌개 근처에서는 가까운 상태들이 한 방향으로 벌어지고 다른 방향으로 모입니다. 흐름은 늘어난 얇은 면을 다시 접습니다. 점 구름은 둘 다 보여 줍니다. 처음에는 실이 되고, 이어서 그 실이 거듭 접힙니다.

로렌츠 사상( $\rho$ 옆의 창 버튼)은 하나의 긴 궤도에서 $z$ 의 각 극댓값을 다음 극댓값과 짝지어 그립니다. 점들은 가늘고 뾰족한 곡선 위에 놓입니다(로렌츠 1963, 그림 4). 이는 근사적인 일차원 표현이며, 정확한 차원 축소로 증명된 것은 아닙니다. $\rho = 28$ 에서 첨점을 벗어난 같은 가지의 표본 기울기는 절댓값이 1보다 큽니다. 따라서 가까웠던 두 극댓값은 다음 회전에서 더 멀어집니다(첨점의 양쪽에 있는 두 점은 그렇지 않을 수 있습니다).

$\rho > 1$ 에서는 원점 외에 두 고정점 $C_\pm = \bigl(\pm\sqrt{\beta(\rho - 1)},\ \pm\sqrt{\beta(\rho - 1)},\ \rho - 1\bigr)$ 이 있으며, $\rho < 1$ 에서는 원점이 모든 해를 끌어당깁니다. $C_\pm$ 는 $\rho = \sigma(\sigma + \beta + 3)/(\sigma - \beta - 1) = 470/19 \approx 24.74$ 에서 아임계 호프 분기를 통해 안정성을 잃습니다. 다른 거동 변화값은 수치적으로 알려져 있습니다. 약 13.926에서는 원점에서 호모클리닉 폭발이 일어나고, 약 24.06에서는 혼돈 끌개가 나타납니다[4, 5]. 약 24.06(수치 결과)부터 $470/19 \approx 24.74$ 까지는 혼돈 끌개와 두 안정 고정점이 공존합니다. 99.65는 주기 거동의 창 안에 있으며, 160에서는 해가 주기 궤도로 수렴합니다.

수치가 의미하는 것

그림은 궤도의 수치 근사입니다. 오차가 시간에 따라 커지므로, 오랜 시간이 지난 뒤의 정확한 위치를 주어진 초기값에 대한 믿을 만한 예측으로 볼 수는 없습니다. 거리 증가와 장기 통계는 간격, 정밀도, 표집을 바꿔 따로 확인해야 합니다.

모든 궤도는 사용자의 컴퓨터에서 고전적인 사차 룽게–쿠타 방법으로 적분합니다. 프레임 속도와 무관한 고정 간격 0.005와 배정밀도를 사용합니다. 점 구름은 기준 궤도와의 차이로 적분하므로 $\varepsilon = 10^{-12}$ 에서도 유효 자릿수를 유지합니다.

화면의 퍼짐은 하나의 동일한 수치 적분기에서 시작값의 차이가 만들어 낸 결과입니다. 이것만으로 정확한 방정식도 같은 민감성을 가진다고 증명할 수는 없습니다. 이론이 필요합니다(고전적인 매개변수에서는 터커의 증명[2]). 일반적인 그림자 추적 보조정리는 로렌츠 끌개에 적용되지 않습니다[6].

$\lambda \approx 0.906$ 는 통계적 추정값입니다. 고전적인 매개변수에서의 장기 평균 비율로, 스프로트가 독립적으로 계산했습니다[3]. 모든 $\rho$ 에서 일정한 값은 아닙니다. ‘정밀도를 십만 배 높이면 약 $\ln 10^5 / 0.906 \approx 12.7$ 시간 단위를 더 얻는다’는 말은 전형적인 지수적 증가에 대한 추정입니다. 한 번의 출발에서는 크게 달라질 수 있습니다.

헷갈리기 쉬운 점

그림에서 교차하는 듯한 선들은 실제로 교차하지 않습니다. 한 상태는 하나의 궤도에만 속하므로 서로 다른 궤도는 교차할 수 없습니다. 교차하는 것은 화면에 비친 모습입니다.
‘나비 효과’라는 이름과 끌개의 나비 모양은 별개의 이야기입니다. 이름은 로렌츠가 1972년 나비의 날갯짓과 토네이도를 주제로 한 강연에서 나왔습니다.
로렌츠가 0.506127 대신 0.506을 입력해 다른 예측을 얻었다는 일화는 1961년 열두 변수 날씨 모형에서 일어난 일입니다. 이 세 방정식의 이야기는 아닙니다[7].
$\rho = 15$ 에서는( $C_\pm$ 가 안정적입니다) 이 페이지의 시작점에서 출발한 점 구름이 퍼지지 않고, 궤도와 함께 줄어들며 고정점으로 들어갑니다.

역사

로렌츠는 1963년 잘츠만의 모형을 바탕으로 대류를 크게 단순화한 방정식을 발표하고, 그 해가 비주기적이며 시작점에 민감하게 의존함을 보였습니다[1]. 스패로의 1982년 저서는 $\rho$ 에 따라 거동이 어떻게 달라지는지 정리했습니다[4]. 컴퓨터에서 보이는 끌개가 정확한 방정식에서도 정말 존재하는지는 터커의 컴퓨터 보조 증명(1999년 발표, 2002년 출판)까지 미해결 문제였습니다[2].

참고 문헌

E. N. Lorenz, “Deterministic nonperiodic flow”, Journal of the Atmospheric Sciences 20 (1963) 130–141.
W. Tucker, “A rigorous ODE solver and Smale’s 14th problem”, Foundations of Computational Mathematics 2 (2002) 53–117.
J. C. Sprott, Chaos and Time-Series Analysis, Oxford University Press (2003).
C. Sparrow, The Lorenz Equations: Bifurcations, Chaos, and Strange Attractors, Applied Mathematical Sciences 41, Springer (1982).
J. L. Kaplan, J. A. Yorke, “Preturbulence: a regime observed in a fluid flow model of Lorenz”, Communications in Mathematical Physics 67 (1979) 93–108.
M. Komuro, “Lorenz attractors do not have the pseudo-orbit tracing property”, Journal of the Mathematical Society of Japan 37 (1985) 489–514.
E. N. Lorenz, The Essence of Chaos, University of Washington Press (1993).

더 읽어 보기: Wikipedia: Lorenz system (영어); MacTutor 수학사 아카이브: Edward Lorenz (영어).