🌀 Odwzorowanie Hénona — dziwny atraktor z prostej iteracji

Odwzorowanie Hénona to dyskretna, dwuwymiarowa iteracja, x_n+1 = 1 − a·x_n² + y_n, y_n+1 = b·x_n, opublikowana przez Michela Hénona w 1976 roku. Zaprojektował ją jako uproszczony zamiennik przekroju Poincarégo systemu Lorenza: równania Lorenza wymagają całkowania trzech sprzężonych równań różniczkowych, by zobaczyć dziwny atraktor, podczas gdy odwzorowanie Hénona produkuje go z garści operacji arytmetycznych na krok, co czyniło je znacznie tańszym do badania na komputerach lat 70., przy zachowaniu tego samego jakościowego zachowania chaotycznego.

Człon x_n² zagina płaszczyznę nieliniowo, rozciągając odległości w jednym kierunku, podczas gdy sprzężenie zwrotne b·x_n i ogólne kurczenie się (odwzorowanie kurczy obszar o czynnik b przy każdej iteracji) zaginają rozciągniętą powierzchnię z powrotem na siebie, warstwa po warstwie, jak ugniatanie ciasta przez piekarza. To powtarzane rozciąganie-i-zaginanie jest dokładnie mechanizmem, który daje czułą zależność od warunków początkowych i dodatni wykładnik Lapunowa — dwa początkowe punkty startujące niezwykle blisko siebie oddalają się wykładniczo już po kilkudziesięciu iteracjach.

Przy klasycznych parametrach (a = 1,4, b = 0,3) atraktor Hénona ma wymiar pudełkowy (pojemnościowy) wynoszący około 1,26 — wartość niecałkowita, umieszczająca go ściśle między jednowymiarową krzywą a dwuwymiarowym obszarem. Przekrój prostopadły do głównej nici atraktora wygląda jak zbiór Cantora: nieskończenie drobny pył punktów z lukami na każdym poziomie powiększenia, co jest znakiem rozpoznawczym fraktalnego dziwnego atraktora, a nie zwykłej gładkiej krzywej.

Gdy a maleje od 1,4, atraktor zwykle przechodzi przez sekwencję bifurkacji zmniejszających okres, przeplataną wąskimi oknami periodycznymi, gdzie chaotyczna chmura zapada się na małą liczbę dyskretnych punktów (stabilną orbitę periodyczną) zamiast wypełniać fraktalną chmurę. Znane okno periodyczne znajduje się w okolicy a ≈ 1,22–1,25 przy b = 0,3. Parametr b ustala tempo kurczenia się obszaru (wyznacznik jakobianu odwzorowania jest stały i równy −b); większe wartości b mniej agresywnie ścieńczają nici atraktora, a bardzo małe wartości b upodabniają odwzorowanie do jednowymiarowego (podobnego do odwzorowania logistycznego), a bardzo duże wartości b mogą powodować rozbieganie orbit do nieskończoności zamiast osiadania na atraktorze.

Ponieważ atraktor Hénona jest zbudowany z nieskończenie wielu operacji rozciągania-i-zaginania nałożonych na siebie, każdy mały fragment atraktora po powiększeniu wygląda statystycznie podobnie do całego atraktora — te same łukowate, bananowe kształty i luki pojawiają się na każdej skali, co nazywa się samopodobieństwem. Regulacja powiększenia w tej symulacji przeskalowuje widok wokół środka atraktora bez usuwania zgromadzonych punktów, więc możesz obserwować, jak pojedyncza gruba nić rozdziela się na kilka cieńszych, równoległych nici w miarę zwiększania powiększenia.

Przekrój Poincarégo bierze ciągły trójwymiarowy przepływ, taki jak system Lorenza, i rejestruje punkty, w których trajektoria przecina ustaloną płaszczyznę, zamieniając ciągły w czasie atraktor w dyskretne w czasie odwzorowanie na tej płaszczyźnie. Hénon zbadał jakościowy kształt takich przekrojów atraktora Lorenza, a następnie skonstruował własne odwzorowanie od podstaw, mające podobne zachowanie zaginające, będąc jednak znacznie prostszym w obliczeniach.

Jeśli zaczniesz dwie iteracje od punktów różniących się tylko o niewielką wartość, np. 10⁻¹⁰, nieliniowe zaginanie odwzorowania Hénona wzmacnia tę różnicę w przybliżeniu wykładniczo z każdą iteracją, gdy orbita znajduje się na chaotycznym atraktorze. Po kilkudziesięciu iteracjach obie trajektorie są całkowicie zdekorelowane, mimo że zaczęły niemal identycznie — to samo zjawisko popularnie znane jako efekt motyla.

Odwzorowanie Hénona nie ma globalnego atraktora dla każdego wyboru a i b — dla wystarczająco dużego a, lub dla punktów początkowych dalekich od basenu przyciągania, iteracje mogą rosnąć bez ograniczeń zamiast osiadać na zbiorze fraktalnym. Ta symulacja ogranicza suwaki do zakresów (a między 0,9 a 1,4, b między 0,1 a 0,4), w których, zaczynając od początku układu, odwzorowanie niezawodnie zbiega do ograniczonego atraktora lub orbity periodycznej.

„Klasyczne (a=1,4, b=0,3)” odtwarza oryginalne parametry Hénona z 1976 roku i kanoniczny dziwny atraktor o wymiarze ≈ 1,26. „Okno periodyczne (a=1,25, b=0,3)” znajduje się w jednym z wąskich okien stabilności, gdzie zachowanie chaotyczne tymczasowo znika, a orbita osiada na małym cyklu periodycznym. „a=1,3 (blisko podwojenia okresu)” znajduje się w obszarze przejściowym między oknem periodycznym a w pełni chaotycznym reżimem.