Синхронізація хаотичних осциляторів

1. Відкриття 1990 року

У Військово-морській дослідницькій лабораторії США Луїс Пекора та Томас Керрол досліджували, чи можна змусити хаотичні кола узгоджуватися одне з одним. Загальноприйнята думка казала «ні»: хаос визначається чутливою залежністю від початкових умов — позитивний показник Ляпунова означає, що будь-яка мала розбіжність між двома хаотичними траєкторіями зростає експоненційно, точно як в атракторі Лоренца.

Пекора та Керрол знайшли лазівку: якщо взяти хаотичну систему, розбити її на підсистеми та використати сигнал однієї підсистеми, щоб керувати ідентичною копією решти підсистем, керована копія може зійтися до тієї ж самої траєкторії, що й оригінал — навіть попри те, що вся система залишається хаотичною. Їхня стаття 1990 року «Synchronization in Chaotic Systems» (Physical Review Letters) є засадничим результатом синхронізації хаосу.

2. Розклад ведучий-ведений

Розбиваємо вектор стану хаотичної системи на дві частини: ведучу підсистему u та ведену підсистему v, так щоб повна динаміка мала вигляд ẋ = f(u, v). Будуємо другу, незалежну копію веденої підсистеми, v′, і подаємо на неї той самий сигнал ведучої u(t) з оригіналу:

Якщо v і v′ починаються з різних початкових умов, чи буде v′(t) → v(t) при t → ∞? Пекора та Керрол показали, що відповідь повністю залежить від стійкості динаміки похибки e = v′ − v, а не від того, чи є вся система хаотичною.

3. Умовні показники Ляпунова

Лінеаризація динаміки похибки навколо траєкторії ведучої u(t) дає умовні (або підсистемні) показники Ляпунова веденої підсистеми — швидкості зростання малих збурень v за умови, що u зафіксовано ведучим сигналом.

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

4. Приклад: синхронізація систем Лоренца

Оригінальний приклад Пекори та Керрола використовував систему Лоренца. Розбиваємо (x, y, z) на ведучу змінну x і ведену пару (y, z):

Ведена репліка має власні (y′, z′), запущені з зовсім інших початкових умов. Оскільки умовні показники Ляпунова підсистеми (y, z), керованої x, обидва від'ємні для класичних параметрів Лоренца (σ=10, ρ=28, β=8/3), (y′, z′) → (y, z) експоненційно швидко — незалежно від того, наскільки далеко одна від одної починалися дві репліки. Ведена «забуває» власну початкову умову й зливається з хаотичною траєкторією провідної.

5. Поза межами ідентичної синхронізації

• Повна (ідентична) синхронізація — випадок вище: v′(t) − v(t) → 0 для ідентичних систем.
• Узагальнена синхронізація — v′(t) сходиться до деякого функціоналу F(v(t)), а не до самого v(t), корисно, коли ведуча і ведена не є ідентичними системами.
• Фазова синхронізація — лише фази двох хаотичних осциляторів злипаються, тоді як їхні амплітуди залишаються некорельованими й хаотичними (актуально для слабко зв'язаних систем, зустрічається також у моделі Курамото зв'язаних осциляторів).
• Синхронізація із затримкою — v′(t) ≈ v(t − τ) для деякої фіксованої затримки τ, поширена в системах із затримкою передачі.

6. Застосування: захищений зв'язок

У 1990-х та 2000-х синхронізація хаосу надихнула хвилю досліджень хаотичного захищеного зв'язку: передавач запускає хаотичне коло і додає повідомлення до свого хаотичного виходу; приймач з ідентичним (синхронізованим) хаотичним колом віднімає відтворений хаотичний сигнал, відновлюючи повідомлення. Оскільки хаотичний носій виглядає як шум для перехоплювача без відповідного кола, це давало фізичний рівень шифрування.

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

7. Псевдокод

// Провідна система Лоренца: x, y, z
// Ведена репліка, керована x(t) провідної: y2, z2
const σ = 10, ρ = 28, β = 8/3
const dt = 0.005

function masterStep(x, y, z):
  dx = σ*(y - x)
  dy = x*(ρ - z) - y
  dz = x*y - β*z
  return [x + dx*dt, y + dy*dt, z + dz*dt]

// Ведена інтегрує лише y2, z2 — керована x провідної
function slaveStep(xDrive, y2, z2):
  dy2 = xDrive*(ρ - z2) - y2
  dz2 = xDrive*y2 - β*z2
  return [y2 + dy2*dt, z2 + dz2*dt]

[x, y, z]   = [0.1,  0,  0]   // початкова умова провідної
[y2, z2]    = [20, -15]  // ведена — зовсім інша початкова умова

loop (кожен кадр):
  [x, y, z] = masterStep(x, y, z)
  [y2, z2]  = slaveStep(x, y2, z2)
  error = Math.hypot(y2 - y, z2 - z)  // → 0 експоненційно