Неперервну траєкторію у фазовому просторі важко візуалізувати та класифікувати. Прозріння Анрі Пуанкаре полягало в тому, щоб стробоскопічно семплювати траєкторію — щоразу, коли вона перетинає вибрану поверхню, фіксувати точку перетину. Утворений розсип точок, переріз Пуанкаре, згортає багатовимірний плин у 2-вимірну картину, що одразу розкриває, чи є динаміка періодичною, квазіперіодичною або хаотичною.

1. Означення та побудова

Нехай система з неперервним часом еволюціонує у фазовому просторі ℝⁿ. Виберіть поверхню перерізу Пуанкаре Σ — многовид корозмірності 1, трансверсальний до плину — та угоду про знак напрямку перетину.

Σ = { x ∈ ℝⁿ : g(x) = 0, ṡ(x) · n̂ > 0 } Відображення Пуанкаре: P : Σ → Σ, xₖ₊₁ = P(xₖ) де xₖ — послідовні точки трансверсального проколу

Відображення P кодує одне повне повернення. Вивчення його нерухомих точок, стійкості та закономірностей ітерацій розкриває повну геометрію неперервного плину.

Практичне обчислення

// RK4-інтегратор з виявленням події (перетин нуля)
function integrateUntilCrossing(state, f, gFn, dt) {
  let prev = gFn(state), s = { ...state };
  for (let iter = 0; iter < 1e6; iter++) {
    s = rk4Step(s, f, dt);
    const curr = gFn(s);
    if (prev < 0 && curr >= 0) {
      // бісекція для уточнення часу перетину
      return refine(s_prev, s, f, gFn, dt);
    }
    prev = curr;
  }
}

2. Відображення повернення Пуанкаре та нерухомі точки

Орбіта плину з періодом T відповідає нерухомій точці x* = P(x*) відображення. Орбіта періоду n (n-цикл плину, поділений на T) відображається в орбіту періоду n для P. Стійкість визначається матрицею монодромії M = DP(x*):

Тип орбіти Власні числа M ───────────────────────────────────────────────── Стійкий фокус |λ₁,₂| < 1, комплексно-спряжені Нестійкий фокус |λ₁,₂| > 1, комплексно-спряжені Центр (інтегровний)|λ₁,₂| = 1, комплексні (KAM-тори) Сідло |λ₁| < 1, |λ₂| > 1 (гіперболічне) Еліптичний острів власні числа на одиничному колі

3. KAM-тори та інваріантні криві

Для майже інтегровних гамільтонових систем теорема Колмогорова-Арнольда-Мозера (KAM) гарантує, що більшість інваріантних торів переживають малі збурення. На перерізі Пуанкаре непошкоджені KAM-тори постають як гладкі замкнені криві. Число обертання ρ = ω₁/ω₂ характеризує кожен тор:

Ірраціональне ρ (діофантове) → гладка замкнена крива (KAM-тор виживає) Раціональне ρ = p/q → ланцюг p-періодичних островів + p нестійких нерухомих точок

На межі між стійкістю та хаосом тори розпадаються на канторі — рештки канторової множини, що діють як часткові бар'єри для перенесення.

4. Стандартне відображення (Чирикова-Тейлора)

Стандартне відображення — це канонічне 2-вимірне відображення зі збереженням площі, що моделює штовхані ротатори, утримання плазми та динаміку комет:

pₙ₊₁ = pₙ + K sin(θₙ) (mod 2π) θₙ₊₁ = θₙ + pₙ₊₁ (mod 2π) K = параметр стохастичності K = 0: чисте обертання (інтегровне) K = 0.97163: розпадається останній KAM-тор (критичне значення, Chirikov 1979) K >> 1: глобальний хаос, дифузія ⟨p²⟩ ≈ K²t/2
Значення K Динаміка Вигляд перерізу Пуанкаре
0 Інтегровне обертання Горизонтальні лінії (стале p)
0.5 Острови + KAM-тори Гладкі криві + ланцюги островів
0.97163 Критичне (останній KAM-тор) Фрактальна межа між хаосом та порядком
1.5 Поширений хаос Розсіяні точки із залишковими островами
5.0 Глобальна стохастичність Рівномірна хмара, видима дифузія Арнольда

5. Острови подвоєння періоду та ланцюг Біркгофа

Поблизу раціонального тора ρ = p/q теорія збурень передбачає ланцюг островів Біркгофа: p еліптичних (стійких) островів, оточених p гіперболічними (нестійкими) сідловими точками. Кожен острів сам може містити підострови за раціональнішими числами обертання — це самоподібна ієрархія:

Дерево чисел обертання: 1/2 → 1/3, 2/3 → 1/4, 3/4, … (послідовність Фарея) Кожне раціональне ρ = p/q має резонансний ланцюг островів періоду p Острови біфуркують зі зростанням K: каскад подвоєння періоду → хаос

Ширина резонансної зони масштабується як ΔI ~ K^{1/2}, тоді як проміжок між сусідніми резонансами зменшується, спричиняючи перекриття резонансів — критерій Чирикова для настання глобального хаосу: K_c ≈ 1, коли півширини резонансів перекриваються.

6. Перерізи Пуанкаре системи Енона-Гайлеса

Гамільтоніан Енона-Гайлеса — стандартний тест на хаос у консервативних системах із 2 ступенями вільності:

H = ½(ẋ² + ẏ²) + ½(x² + y²) + x²y − y³/3 = E Переріз: y = 0, ẏ > 0 Змінні на перерізі: (x, ẋ) E = 1/12: переважно регулярна (гладкі KAM-криві) E = 1/8 : змішаний фазовий простір E = 1/6 : настання глобального хаосу (енергія сепаратриси)
Властивість Інтегровна границя (E→0) Хаотичний режим (E≈1/6)
Структура перерізу Вкладені інваріантні криві Розсіяне хаотичне море + острови
Найбільший показник Ляпунова λ₁ ≈ 0 ≈ 0.08
Короткочасний Ляпунов Коливається біля нуля Додатне експоненційне зростання
Час повторення Передбачуваний (теорема Пуанкаре про повернення) Статистика польотів Леві через «липкість»

7. Підкова Біркгофа-Смейла та гомоклінічні плетива

Походження хаосу поблизу гомоклінічної точки — де нестійкий многовид W^u сідла перетинає власний стійкий многовид W^s — зрозумів сам Пуанкаре та формалізував Смейл. Кожен трансверсальний перетин зумовлює нескінченно багато інших, породжуючи геометрію відображення-підкови:

Трансверсальний гомоклінічний перетин ⟹ підкова Смейла ⟹ зліченна нескінченність періодичних орбіт ⟹ незліченна множина обмежених неперіодичних орбіт ⟹ топологічна ентропія h > 0 (додатний показник Ляпунова)

«Гомоклінічне плетиво» — це фрактальна мережа многовидів, що перетинаються, видима поблизу сепаратрис у будь-якому перерізі Пуанкаре хаотичної гамільтонової системи.

8. Інтерактив: дослідник стандартного відображення

Клацніть по полотну, щоб додати початкову умову. Кожне клацання запускає траєкторію, що ітерується протягом 5 000 кроків, наносячи точки свого перерізу Пуанкаре. Налаштуйте стохастичність K, щоб розкрити KAM-тори (мале K), ланцюги островів (помірне K) або глобальний хаос (велике K). Використовуйте Очистити, щоб скинути.

Точки: 0

Кожен колір відповідає іншій початковій умові. Гладкі замкнені криві → KAM-тори; ланцюги точок → раціональні резонанси; розсіяні точки → хаос.