Біфуркаційні діаграми: подвоєння періоду, Фейгенбаум та шлях до хаосу

1. Логістичне відображення

Логістичне відображення — це найпростіше нелінійне різницеве рівняння, що демонструє повний шлях від порядку до хаосу:

Запроваджене Робертом Меєм (1976) у контексті популяційної динаміки, воно моделює взаємодію між розмноженням (лінійне за x) та конкуренцією за ресурси (член −x²). Попри уявну простоту, воно кодує всю математичну структуру одновимірного хаосу.

Ключові властивості f(x) = rx(1−x):

• Воно унімодальне: один локальний максимум при x* = 1/2, значення f(1/2) = r/4
• Для r ≤ 4, f відображає [0,1] у [0,1] (інваріантність уперед)
• Критична точка x = 1/2 відіграє особливу роль: її орбіта визначає біфуркаційну структуру
• Для r > 4 деякі орбіти втікають у −∞ (хаотичні перехідні процеси)

Біфуркаційну діаграму будують так: (1) обираємо стартове x₀ ≈ 0,5, (2) ітеруємо 1000 разів, щоб перехідні процеси згасли, (3) відкладаємо наступні 200 ітерацій як точки в горизонтальній позиції r. Повторення цього для тисяч значень r дає знамените «дерево подвоєння періоду».

2. Нерухомі точки та аналіз стійкості

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

Нерухома точка задовольняє x* = f(x*):

Стійкість визначається похідною |f'(x*)| = |r(1 − 2x*)|:

• x* = 0: f'(0) = r. Стійка для r < 1, нестійка для r > 1.
• x* = 1 − 1/r: f'(x*) = 2 − r. Стійка для |2 − r| < 1, тобто 1 < r < 3.

Втрата стійкості при r = 3

При r = 3 похідна f'(x*) = −1 точно. Нерухома точка стає перекидним сідлом періоду 1. Для r трохи вище 3 малі збурення зростають (поперемінно вище/нижче x*), і система переходить у цикл періоду 2: пару точок {x₁, x₂}, що задовольняють x₂ = f(x₁) та x₁ = f(x₂), тобто x = f(f(x)) = f²(x).

Це біфуркація подвоєння періоду (перекидна): стійка нерухома точка породжує стійкий 2-цикл при r = 3 у процесі, керованому перетином власним значенням одиничного кола в точці −1.

3. Каскад подвоєння періоду

Сам цикл періоду 2 втрачає стійкість при r₂ ≈ 3,449, породжуючи цикл періоду 4. Цикл періоду 4 біфуркує при r₃ ≈ 3,544 у період 8 і так далі. Послідовність точок біфуркації {rₙ} збігається до точки накопичення r∞ ≈ 3,5699456... — початку хаосу.

Відношення δₙ збігається до сталої Фейгенбаума δ ≈ 4,6692. Проміжки між послідовними точками біфуркації звужуються геометрично: кожен інтервал становить приблизно 1/δ ≈ 21,4% від попереднього. Цикл періоду 2ⁿ має 2ⁿ гілок на біфуркаційній діаграмі, усі створені одним і тим самим перекидним механізмом.

Самоподібність: піддіаграма поблизу кожної точки біфуркації r = rₙ є масштабованою копією всієї діаграми. Масштабування осі x на множник α ≈ −2,5029 (друга стала Фейгенбаума) дає ідентичну копію. Це характерна ознака нерухомої точки ренормгрупи у функціональному просторі.

4. Універсальність Фейгенбаума

Мітчелл Фейгенбаум виявив (1975–1978), що стала δ ≈ 4,6692 є універсальною: вона зʼявляється в кожному гладкому унімодальному відображенні з єдиним квадратичним максимумом, незалежно від конкретної формули. Логістичне відображення, синусоїдальне відображення (xₙ₊₁ = r sin(πxₙ)), намет-відображення та незліченні фізичні системи — усі дають одне й те саме δ.

Пояснення через ренормгрупу

Каскад подвоєння періоду є нерухомою точкою оператора ренормалізації T, що діє на функціональному просторі:

Лінеаризація T у нерухомій функції g має одне нестійке власне значення δ (усі інші < 1). Цей єдиний нестійкий напрям пояснює, чому всі гладкі унімодальні відображення наближаються до того самого каскаду: усі вони притягуються до стійкого многовиду g, а потім рухаються вздовж нестійкого напряму зі швидкістю δ. Клас універсальності аналогічний класам універсальності у статистичній механіці поблизу критичних точок — глибокий звʼязок, формалізований ренормгрупою Вілсона.

Експериментальне підтвердження

δ було виміряно у фізичних експериментах:

• Конвекція Релея–Бенара (Лібшабер і Маурер, 1980): δ виміряно у ртутній комірці з добрим узгодженням
• Нелінійні електричні кола (Лінсей, 1981): подвоєння періоду у керованому RLC-колі
• Акустична турбулентність: каскад Фейгенбаума у керованій кавітації рідини
• Лазерні коливання: модульований CO₂-лазер, що демонструє універсальний каскад

5. Хаос, періодичні вікна та кризи

За межею r∞ ≈ 3,5699 система переважно хаотична — але не рівномірно. Біфуркаційна діаграма після початку хаосу перемежовується періодичними вікнами: вузькими діапазонами r, де стійкий цикл знову виринає з хаосу.

Періодичні вікна

Найпомітніше вікно поблизу r ≈ 3,828 — це вікно періоду 3. Його існування гарантується теоремою Лі–Йорке (1975): «період 3 означає хаос». Точніше, якщо неперервне відображення має орбіту періоду 3, воно має орбіти всіх періодів.

Кожне періодичне вікно саме проходить підкаскад подвоєнь періоду, що веде до власного режиму міні-хаосу. Патерн повторюється на кожному масштабі — фрактальна структура у просторі параметрів. Зʼявляються вікна всіх періодів, упорядковані за порядком Шарковського:

Внутрішні кризи

За певних значень r хаотичний атрактор раптово змінює розмір — внутрішня криза. При знаменитому r ≈ 3,6786 атрактор розширюється з двох неперетинних інтервалів до єдиного інтервалу, що зумовлено зіткненням нестійкої орбіти періоду 3 з межею басейну. Ці переходи (Гребоджі, Отт і Йорке, 1982) є раптовими, але структурно зрозумілими через нестійкі періодичні орбіти.

6. Накладання показника Ляпунова

Показник Ляпунова λ(r) кількісно оцінює середнє експоненційне розходження й відрізняє порядок (λ < 0) від хаосу (λ > 0):

При накладанні на біфуркаційну діаграму:

• λ = 0 точно в кожній точці біфуркації r₁, r₂, r₃, ... — чітка сигнатура
• λ < 0 на всьому протязі кожного періодичного вікна — поглиблюючись до центру кожного вікна
• λ → −∞ на суперстійких орбітах (де f'(x*) = 0, тобто орбіта проходить через x = 1/2)
• λ > 0 у хаотичних смугах — найбільше значення при r = 4, де λ = ln 2 ≈ 0,693

При r = 4 логістичне відображення спряжене з намет-відображенням T(x) = 1 − |2x − 1| через заміну змінних x = sin²(πθ/2). Намет-відображення є кусково-лінійним з нахилом ±2 усюди, що дає λ = ln 2 точно. Це максимально можливий показник Ляпунова для цієї родини.

7. Класифікація типів біфуркацій

Подвоєння періоду логістичного відображення — один із кількох канонічних типів біфуркацій, що зʼявляються в усіх динамічних системах:

Дискретні відображення (1D)

Потоки (ЗДР, 2D+)

8. Інтерактивна біфуркаційна діаграма

Рендерер нижче будує біфуркаційну діаграму логістичного відображення у реальному часі. Перетягуйте повзунки діапазону, щоб наблизитися до будь-якої області — самоподібна структура видима на кожному масштабі.

// Основний рендерер біфуркаційної діаграми
function renderBifurcation(canvas, rMin, rMax, showLyap) {
  const ctx  = canvas.getContext('2d');
  const W    = canvas.width, H = canvas.height;
  const N_r  = W * 2;          // вибірки r (2 на піксель)
  const WARM = 500;            // перехідні ітерації, які відкидаємо
  const PLOT = 300;            // ітерації для відображення на кожне значення r

  ctx.fillStyle = '#0a0a14';
  ctx.fillRect(0, 0, W, H);

  for (let i = 0; i < N_r; i++) {
    const r = rMin + (rMax - rMin) * (i / N_r);
    const px = (i / N_r) * W;

    let x = 0.5;
    // Відпрацювати перехідні процеси
    for (let k = 0; k < WARM; k++) x = r * x * (1 - x);

    // Зібрати орбіту та суму Ляпунова
    let lyapSum = 0;
    for (let k = 0; k < PLOT; k++) {
      const deriv = Math.abs(r * (1 - 2 * x));
      lyapSum += deriv > 1e-12 ? Math.log(deriv) : -30;
      x = r * x * (1 - x);

      // Відкласти точку орбіти
      const py = H - x * H;
      ctx.fillStyle = 'rgba(99,102,241,0.35)';
      ctx.fillRect(px, py, 1.5, 1.5);
    }

    // Необовʼязкове накладання показника Ляпунова
    if (showLyap) {
      const lam = lyapSum / PLOT;
      const lamPy = H / 2 - lam * (H / 4);  // масштаб: -2..+2 → H..0
      ctx.fillStyle = lam > 0 ? 'rgba(239,68,68,0.7)' : 'rgba(34,197,94,0.6)';
      ctx.fillRect(px, Math.min(Math.max(lamPy, 0), H - 1), 1.5, 1.5);
    }
  }

  // Нульова лінія як орієнтир для Ляпунова
  if (showLyap) {
    ctx.strokeStyle = 'rgba(255,255,255,0.15)';
    ctx.beginPath(); ctx.moveTo(0, H / 2); ctx.lineTo(W, H / 2); ctx.stroke();
  }
}