Притягач Лорренца: Як метелик змінює погоду

Три ОДУ, одна метелик

Едвард Лорнз, сотрудник MIT, спростив модель з водночас 12 рівнянь атмосферної конвекції до трьох спарних звичайних диференціальних рівнянь, де x — швидкість перекидання конвективних потоків, y — різниця температур між підіймаючимися та опускаючимися течіями, а z — відхилення вертикального температурного профілю від лінійної траєкторії:

dx/dt = σ(y − x) dy/dt = x(ρ − z) − y z/dt = xy − βz

Класичні хаотичні параметри: σ = 10, ρ = 28, β = 8/3

З цими параметрами траєкторія ніколи не повторюється і ніколи не виходить до нескінченності — вона спіралевидним чином обертається навколо однієї лобової частини, непередбачувано перетинаючи іншу, спіралевидно там і повертаючись назад, відтворюючи форму «метелика» з двома крилами, яка стала найвідомішим символом хаосу. Лорнз опублікував свої результати у 1963 році; на лекції 1972 року він поставив знамените питання — чи можуть крильця метелика в Бразилії спровокувати торнадо в Техасі?

dx/dt = σ(y − x)
dy/dt = x(ρ − z) − y
dz/dt = xy − βz
classic chaotic values: σ = 10, ρ = 28, β = 8/3

жива демонстрація · пов'язана симуляція● LIVE

Ліapunovівські експоненти та межа передбачуваності

Як швидко розходяться дві сусідні траєкторії? Відповідь - найбільший ліapunovівський експонент λ₁ ≈ +0.906 для класичних параметрів — невелика початкова помилка δ₀ зростає як δ(t) ≈ δ₀·e^(λ₁t). Оскільки це експоненціально, межа передбачуваності T_p ≈ (1/λ₁)·ln(Δ/δ₀) росте лише логарифмічно з точністю вимірювання: зменшення початкової помилки на 0.76 одиниць часу.

Це пояснює, чому глобальні моделі погоди з 10⁹ градусами вільності все ще втрачають передбачуваність через приблизно 7–10 днів — це винуватий хаос сам по собі, а не недосконале обчислювальне потенціал.

λ₁ ≈ +0.906 λ₂ ≈ 0.000 λ₃ ≈ −14.572 D_KY = 2 + (λ₁+λ₂)/|λ₃| ≈ 2.062 ← Каплановський-Йоркський фрактальний вимір Фрактальна розмірність аттрактора ледь перевищує 2 — тонка, як поверхня, але з нескінченною тонкошаровою структурою. Операційні прогнористи реагують на цей фундаментальний ліміт із використанням ансамблевого прогнозування: запуском десятків симуляцій з незначно відмінних початкових умов і читанням невизначеності, як далеко вони розходяться.

λ₁ ≈ +0.906   λ₂ ≈ 0.000   λ₃ ≈ −14.572
D_KY = 2 + (λ₁+λ₂)/|λ₃| ≈ 2.062   ← Kaplan-Yorke fractal dimension

Frequently asked questions

Які три рівняння Лорренца?

dx/dt = σ(y−x), dy/dt = x(ρ−z)−y, dz/dt = xy−βz. З класичними параметрами σ=10, ρ=28, β=8/3 система ніколи не стабілізується в періодичній орбіті та ніколи не розходяться — траєкторія назавжди утримується на фрактальній «дивній привабливій» формі, що нагадує метелика.

Що таке ефект метелика, точніше?

Це твердження про те, що близькі траєкторії в хаотичній системі експоненціально швидко відокремлюються, швидкість якого визначається найбільшим лінійним показником (λ₁ ≈ 0.906 для класичних параметрів Лорренца). Навіть невелике відхилення менше за діаметр атома зрештою робить довгострокове прогнозування неможливим, тому Лорренц запитав, чи можуть крильця метелика в Бразилії викликати торнадо в Техасі.

Чому не можна прогнозувати погоду більше ніж приблизно на 10 днів?

Атмосфера є хаотичною системою з позитивним лінійним показником, тому будь-яка похибка вимірювань початкових умов експоненціально зростає. Оскільки межа передбачуваності росте логарифмічно зі збільшенням точності вимірювань, то подвоєння точності інструментів купує лише кілька додаткових днів майстерності, а не місяці.

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте the simulation і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.

▶ Відкрити симуляцію the simulation