🔮 Eksplorator fraktali Julii — interaktywna dynamika zespolona

Dla ustalonej liczby zespolonej c zbiór Julii to granica między punktami startowymi z₀, których orbity pod działaniem z ↦ z² + c pozostają ograniczone, a tymi, które uciekają do nieskończoności. Symulacja koloruje obszar ucieczki zależnie od tego, jak szybko punkt rozbiega, pozostawiając wypełniony zbiór ograniczony na czarno.

Każdy piksel to liczba zespolona z, a shader na karcie graficznej iteruje z = z² + c aż do limitu Max iterations. Jeśli moduł z pozostaje poniżej promienia ucieczki równego 256, piksel należy do zbioru i jest czarny; w przeciwnym razie jest kolorowany zależnie od płynnej wartości czasu ucieczki. Wykonywanie tego dla każdego piksela na karcie graficznej pozwala renderować obraz w czasie rzeczywistym.

Ustalają część rzeczywistą i urojoną stałej c, w zakresie od -1 do 1 z krokiem 0,001. Każdy zbiór Julii jest zdefiniowany pojedynczą wartością c, więc przeciąganie dowolnego suwaka w czasie rzeczywistym przekształca fraktal w zupełnie inny kształt.

Odwzorowanie ma postać z(n+1) = z(n)² + c, stosowane wielokrotnie od punktu startowego każdego piksela. Aby uniknąć pasmowania kolorów, symulacja używa płynnej liczby iteracji: mu = n + 1 - log(log(modułu z)) / log(2), co interpoluje między całkowitymi wartościami ucieczki, dając ciągłe cieniowanie.

Przesuwa stałą c po okręgu o promieniu 0,7885 na płaszczyźnie zespolonej — to klasyczna animacja zbiorów Julii. Ponieważ kształt zależy wyłącznie od c, fraktal nieustannie rozkwita i się fałduje, przechodząc przez sekwencję powiązanych form, gdy punkt przemierza okrąg.

Zbiór Mandelbrota jest katalogiem wszystkich wartości c, dla których odpowiadający zbiór Julii jest spójny, a nie rozproszony w niespójny pył. Zgodnie z twierdzeniem Douady–Hubbarda zbiór Julii J(c) jest spójny wtedy i tylko wtedy, gdy c leży wewnątrz zbioru Mandelbrota, więc oba zbiory są ze sobą głęboko powiązane.

Każdy preset przenosi do słynnej wartości c. Dendrite przy -0,8, 0,156 to rozgałęziające się drzewo, Siegel disc przy -0,7269, 0,1889 zawiera dysk rotacyjny, Rabbit przy -0,4, 0,6 to królik Douady'ego, Quadratic przy 0,285, 0,01 daje pogrubioną spiralę, a Cantor dust przy 0,45, 0,1428 leży poza zbiorem Mandelbrota, więc rozpada się na nieskończenie wiele drobin.

Liczba iteracji mieści się w zakresie od 50 do 1000. Punkty leżące blisko granicy zbioru potrzebują wielu iteracji, by ujawnić, czy uciekają, więc niski limit rozmywa cienkie nitki, a wysoki wyostrza je kosztem szybkości renderowania. Głębsze przybliżenia zwykle wymagają większej liczby iteracji, by odsłonić detale.

Tak, w granicach skończonej precyzji. Shader wykonuje dokładną iterację zespoloną z = z² + c z promieniem ucieczki 256 i standardową matematyką płynnego kolorowania. Jedynymi przybliżeniami są skończony limit iteracji oraz pojedyncza precyzja zmiennoprzecinkowa, które ograniczają, jak głęboko można powiększać obraz, zanim pojawią się artefakty.

Na granicy wypełnionego zbioru dwa dowolnie bliskie sobie punkty startowe mogą mieć zupełnie różne losy podczas iteracji: jeden pozostaje ograniczony, a drugi ucieka. Ta czuła zależność od warunków początkowych, połączona z w pełni deterministyczną regułą, jest charakterystycznym sygnaturą chaosu deterministycznego.

Poza sztuką generatywną teoria iteracji leżąca u podstaw zbiorów Julii i Mandelbrota znajduje zastosowanie w projektowaniu anten, kompresji obrazu i sygnału oraz w modelowaniu układów dynamicznych w fizyce i biologii. Są też standardowym przykładem dydaktycznym dla geometrii fraktalnej, analizy zespolonej i matematyki chaosu.