Атрактор Ресслера: спіральний хаос, біфуркації та дивні атрактори

1. Система Ресслера

Атрактор Ресслера визначається автономним тривимірним ОДУ:

Канонічний набір параметрів для добре сформованого спірального атрактора — a = 0,2, b = 0,2, c = 5,7. Єдина нелінійність — це одинарний добуток z·x у третьому рівнянні. Ресслер свідомо обрав цю мінімальну форму, щоб продемонструвати: хаос не вимагає внутрішньої складності — достатньо одного квадратичного члена.

Порівняйте з системою Лоренца, яка має два нелінійних члени (xz та xy). Простота системи Ресслера робить її механізм прозорим: площина x–y містить лінійну спіраль назовні (зумовлену членом a·y), а напрямок z забезпечує механізм згортання та повторного введення. Щоразу, коли орбіта достатньо далеко відхиляється назовні, вона повертається до початку динамікою z, створюючи характерну смугову структуру атрактора.

Нерухомі точки

Встановлюючи всі похідні рівними нулю:

Для канонічних параметрів: z ≈ 0,0394 (нестійка спіраль) та z ≈ 28,46 (сідло-фокус). Атрактор обгортається навколо нестійкої нерухомої точки в площині xy.

2. Геометрія фазового простору

Атрактор Ресслера має форму скрученої стрічки — топологічно подібної до стрічки Мебіуса структури у трьох вимірах. На відміну від двокрилого «метелика» атрактора Лоренца, у Ресслера одна пелюстка, що спірально розширюється у z ≈ 0 та періодично складається назад.

Механізм можна розкласти на три фази орбіти:

1. Спіраль: траєкторія обертається навколо осі z у площині xy, розширюючись назовні через член a·y.
2. Підйом: коли x стає великим, член z·(x−c) стає додатнім, z швидко зростає — орбіта «злітає» з площини xy.
3. Повернення: великий z штовхає dx/dt = −y−z від'ємним, x та y колапсують до початку координат, a z повертається до нуля.

Цей механізм розтягування та складання — фундаментальний інгредієнт усіх дивних атракторів. Версія Ресслера особливо чиста, оскільки складання відбувається лише в одному напрямі: стрічка складається один раз за оберт, формуючи топологічну скрутку π (напівоберт).

Переріз Пуанкаре

Переріз Пуанкаре при y = 0 (з dy/dt > 0) виявляє гладку одновимірну карту атрактора — характерну ознаку хаосу типу Ресслера. Будуючи графік n-го перетину x у порівнянні з (n+1)-м, отримуємо майже параболічну карту xₙ₊₁ = f(xₙ), що нагадує логістичне відображення. Ця структурна карта є рушієм каскаду подвоєнь до хаосу.

3. Шлях до хаосу через біфуркації

Фіксуючи a = 0,2, b = 0,2 та змінюючи c, система Ресслера зазнає канонічного каскаду подвоєнь (сценарій Фейгенбаума) до хаосу:

Відношення послідовних інтервалів подвоєнь сходиться до константи Фейгенбаума δ ≈ 4,669 — тієї ж універсальної константи, що з'являється в логістичному відображенні, відображенні Енона та незліченних інших системах з подвоєнням. Ця універсальність стала одним із найвражаючих відкриттів теорії хаосу (Фейгенбаум, 1978).

Біфуркаційна діаграма

Для побудови діаграми: перебирайте c від 2 до 8, інтегруйте довгий перехідний процес (відкидіть перші 500 оборотів), потім будуйте графік координат x у точках перерізу Пуанкаре. Результуюча діаграма показує чисте дерево подвоєнь, що зливається у хаотичну смугу, пронизану вікнами стабільності (зокрема, при c ≈ 5,3: вікно 3-го порядку).

4. Спіральний хаос проти гвинтового хаосу

У системі Ресслера зі збільшенням c виникають два якісно відмінних типи хаосу:

Спіральний хаос (c ≈ 5,7)

Атрактор складається з однієї когерентної смуги, що спірально розширюється та чисто складається назад в одномодальний спосіб. Карта Пуанкаре є унімодальною (однокупольна). Це топологічно еквівалентно логістичному відображенню на інтервалі — «найпростіший» тип хаосу.

• Один позитивний показник Ляпунова (λ₁ > 0)
• Фрактальна розмірність: трохи вище 2 (D_KY ≈ 2,01 для канонічних параметрів)
• По суті унімодальна карта повернення → сильний зв'язок з 1D хаосом

Гвинтовий хаос (c ≈ 6,5)

При більших c атрактор розростається та смуги зливаються — орбіта тепер кілька разів намотується навколо нерухомої точки перед повверненням. Карта повернення стає мультимодальною (кілька горбів). Це топологічно складніше та дає багатшу структуру атрактора:

• При погляді вздовж осі x атрактор нагадує гвинт або завиток
• Стабільні вікна густіші та вужчі
• Більший позитивний показник Ляпунова → швидше перемішування

Перехід від спірального до гвинтового хаосу проходить через кризу злиття смуг — раптову топологічну зміну структури атрактора при критичному значенні c, де дві раніше роз'єднані смуги орбіт зливаються в одну. Це приклад внутрішньої кризи (Грьобогі, Отт, Йорк, 1982).

5. Показники Ляпунова

Показники Ляпунова вимірюють середню експоненційну швидкість розбіжності (або збіжності) близьких траєкторій вздовж головних напрямків у фазовому просторі. Тривимірна система має три показники λ₁ ≥ λ₂ ≥ λ₃.

Для атрактора Ресслера при канонічних параметрах:

Сума λ₁ + λ₂ + λ₃ = div(F) = a − c ≈ −5,5 (дивергенція векторного поля, що дорівнює усередненій швидкості стиснення фазового простору). Ця від'ємна дивергенція підтверджує, що система дисипує об'єм і підтримує існування дивного атрактора.

Розмірність Каплана–Йорке

Ляпунівська (Каплан–Йорке) розмірність оцінює фрактальну розмірність атрактора:

Ця надзвичайно мала фрактальна розмірність (ледь перевищує 2) відповідає тонкій, стрічкоподібній структурі атрактора — він майже двовимірний, але з нескінченно тонкою структурою.

Чисельне обчислення показників Ляпунова

Стандартний алгоритм (Бенеттін та ін., 1980) інтегрує систему разом із варіаційними рівняннями (якобіан, застосований до набору ортонормованих векторів) з періодичною ортонормалізацією Грама–Шмідта:

6. Порівняння з системою Лоренца

Система Лоренца перемішується швидше (більший λ₁) та хаотично перемикається між двома листами. Система Ресслера спірально розвивається повільніше та передбачуваніше, роблячи її більш м'яким вступом до хаосу. Ресслер також простіший для симуляції: помірне стиснення (λ₃ ≈ −5,4 проти −14,6 у Лоренца) означає, що його повільний многовид менш жорсткий чисельно.

7. Реалізація RK4 на JavaScript

// Атрактор Ресслера — інтегратор RK4 із візуалізацією на Canvas 2D
const canvas = document.getElementById('canvas');
const ctx = canvas.getContext('2d');

// Параметри
const a = 0.2, b = 0.2, c = 5.7;
const dt = 0.01;
const STEPS_PER_FRAME = 50;

// Стан
let x = 0.1, y = 0.0, z = 0.0;

// Похідні Ресслера
function deriv(x, y, z) {
  return {
    dx: -y - z,
    dy:  x + a * y,
    dz:  b + z * (x - c)
  };
}

// Один крок RK4
function rk4Step(x, y, z, h) {
  const k1 = deriv(x, y, z);
  const k2 = deriv(x + h/2*k1.dx, y + h/2*k1.dy, z + h/2*k1.dz);
  const k3 = deriv(x + h/2*k2.dx, y + h/2*k2.dy, z + h/2*k2.dz);
  const k4 = deriv(x + h*k3.dx,   y + h*k3.dy,   z + h*k3.dz);
  return {
    x: x + h/6 * (k1.dx + 2*k2.dx + 2*k3.dx + k4.dx),
    y: y + h/6 * (k1.dy + 2*k2.dy + 2*k3.dy + k4.dy),
    z: z + h/6 * (k1.dz + 2*k2.dz + 2*k3.dz + k4.dz)
  };
}

// Проєкція 3D → 2D (вид у площині xy)
function project(x, y) {
  const cx = canvas.width  / 2;
  const cy = canvas.height / 2;
  const scale = 10;
  return { px: cx + x * scale, py: cy - y * scale };
}

function render() {
  // Затухання сліду
  ctx.fillStyle = 'rgba(10, 10, 20, 0.04)';
  ctx.fillRect(0, 0, canvas.width, canvas.height);

  for (let i = 0; i < STEPS_PER_FRAME; i++) {
    const prev = project(x, y);
    const next = rk4Step(x, y, z, dt);
    x = next.x; y = next.y; z = next.z;
    const curr = project(x, y);

    // Колір за висотою z
    const t  = Math.min(Math.max(z / 25, 0), 1);
    const r  = Math.round(100 + t * 155);
    const g  = Math.round(50  + t * 50);
    const bl = Math.round(200 - t * 100);
    ctx.strokeStyle = `rgb(${r},${g},${bl})`;
    ctx.lineWidth = 0.8;
    ctx.beginPath();
    ctx.moveTo(prev.px, prev.py);
    ctx.lineTo(curr.px, curr.py);
    ctx.stroke();
  }
  requestAnimationFrame(render);
}

ctx.fillStyle = '#0a0a14';
ctx.fillRect(0, 0, canvas.width, canvas.height);
render();
        

Кути огляду

Атрактор Ресслера виглядає по-різному з кожної осі:

• Площина xy (вид по осі z): чітко видно спіральні рукави — канонічна «спіраль Ресслера»
• Площина xz (вид по осі y): видно вертикальну структуру складання — ширина та висота смуги
• Площина yz (вид по осі x): «гвинтова» структура у режимі гвинтового хаосу; структура циклів у періодичному режимі
• 3D обертання: найкраще для порівняння з метеликом Лоренца

8. Застосування та розширення

Хімічні осцилятори

Ресслер розробив свою систему як мінімальну модель хімічного осцилятора з механізмом повторного введення — аналог реакції Бєлоусова–Жаботинського. Абстрактна модель відображає якісну динаміку без вказівки конкретного хімічного механізму. Подальші дослідження (Спарроу, 1982; Леткельє та ін., 1995) пов'язали конкретні режими параметрів Ресслера з реальними осцилюючими хімічними системами.

Пов'язані осцилятори Ресслера та синхронізація

Дві системи Ресслера, пов'язані дифузійним членом:

При малому зв'язку k орбіти незалежні. Вище критичного k_c системи синхронізуються за фазою (фази блокуються, тоді як амплітуди залишаються незалежними) — одна з відмінних рис синхронізації Пекора–Керрола (1990). Це явище спостерігається в парах лазерів, парах нейронів і осцилюючих хімічних реакторах.

Гіперхаос Ресслера

Додавання четвертої змінної створює гіперхаотичну систему Ресслера з двома позитивними показниками Ляпунова:

Гіперхаос (d = 0,05, канонічний c = 0,5) має λ₁ > 0 та λ₂ > 0 — траєкторії розходяться одночасно у двох незалежних напрямках. Це суттєво ускладнює керування, синхронізацію або прогнозування порівняно зі звичайним хаосом.

Аналіз часових рядів і теорема Такенса

Один виміряний часовий ряд (наприклад, x(t), відібраний через рівні інтервали) можна використовувати для реконструкції повного атрактора через занурення із затримкою (Такенс, 1981): будуємо {x(t), x(t+τ), x(t+2τ)} для підходящого τ. Отримана тривимірна хмара апроксимує топологію справжнього атрактора Ресслера, дозволяючи оцінити показники Ляпунова та розмірність атрактора лише з експериментальних даних — без знання базових рівнянь.