Fraktal Nowa: Gdy Poszukiwanie Pierwiastków Rysuje Fraktal

Korzenie problemu stają się widoczne

Metoda Newtona znajduje pierwiastek funkcji, zaczynając od przybliżenia z₀ i wielokrotnie iterując w kierunku, gdzie styczna do wykresu funkcji przecina oś zerową:

z_{n+1} = z_n − f(z_n) / f'(z_n) Dla funkcji rzeczywistej jest to historia jednowymiarowa: od dobrego początkowego przybliżenia sekwencja zbiega się do najbliższego pierwiastka, podwajając liczbę poprawnych cyfr w każdym kroku (zbieżność kwadratowa). Fragmety Newtona powstaje z zadania bardziej ryzykownego: zastosowania dokładnej samej iteracji z liczbami zespolonymi na całej płaszczyźnie zespolonej do wielomianu takiego jak f(z) = zⁿ − 1, którego n pierwiastków są to n-te korzenie jedności rozmieszczone równomiernie po okręgu jednostkowym. Kolorujemy każdy punkt początkowy z₀, do którego pierwiastek ostatecznie zbiega, i cieniujemy kolor w zależności od liczby iteracji, które zajęło to zbiegnięcie. Rezultatem jest jeden z najbardziej skomplikowanych fraktali w matematyce, a wszystko to prosto wyłania się z rachunku różniczkowego wysokiej szkoły.

z_{n+1} = z_n − f(z_n) / f'(z_n)

Baseny przyciągania i dlaczego granica jest fraktalna

W pobliżu każdego pierwiastka metoda Newtona zachowuje się dokładnie tak, jak można by się spodziewać: cała okolica punktów startowych szybko zbiega się do tego jednego pierwiastka, tworząc duży, dobrze zdefiniowany basen przyciągania. W odległości dużych od korzeni historia ulega zmianie. Ponieważ f'(z) może być małe lub korekta kroku może przebiec w potęgę, punkty, które są ze sobą bardzo blisko, mogą zostać skierowane do różnych pierwiastków lub potrzebują znacznie różnej liczby iteracji, aby się zbiegnąć. Zbiór punktów znajdujących się na granicy między dwoma lub więcej basenami – gdzie nawet niewielkie zakłócenia powodują zmianę pierwiastka, na którym lądujesz – to dokładnie ten zestaw o ekstremalnej wrażliwości, a dla metody Newtona zastosowanej do wielomianu stopnia trzeciego lub wyższego, ta graniczna struktura jest fraktalem z samopodobną strukturą na każdym poziomie: powiększając dowolny punkt graniczny, znajdziesz mniejsze kopie całego wzoru, nieskończony regres splątanych basenów.

Dlaczego kubiki i wyższe zachowują się bardziej skomplikowanie niż kwadraty

Dla kwadratu, takiego jak z² − 1, dwa obszary przyciągania są po prostu pół płaszczyznami lewymi i prawymi, rozdzielone prostą linią — w ogóle nie występuje żadne zachowanie fraktalne. Problem pojawia się na trzecim stopniu. Dla z³ − 1 trzy obszary przyciągania splatają się tak skomplikowanie, że każdy punkt graniczny między dwoma obszarami przyciągania jest również punktem granicznym trzeciego obszaru — to zaskakujący fakt topologiczny zwany właściwością Jezior Wady: trzy (lub więcej) obszarów mogą dzielić wspólną granicę wszędzie, nie tylko w izolowanych punktach. To część tego, dlaczego obraz wygląda w ten sposób: kolory z trzech pierwiastków splatają się ze sobą niezależnie od tego, jak bardzo zbliżysz się do dowolnej strefy granicznej.

Połączenie z zbiorami Julia

Metoda Newtona dla wielomianu f jest samą sobą daną funkcją racyjną z ↦ N(z) = z − f(z)/f'(z), a iteracja dowolnej funkcji racyjnej na płaszczyźnie zespolonej to precyzyjnie ustawienie dynamiki złożonej, które generuje zbiory Julia i podział Fatou/Julia. Pierwiastki wielomianu są superprzyciągającymi punktami stałymi N (pochodna N tam zanika, co dokładnie wyjaśnia kwadratową konwergencję), zbiorami siatki są zbiory Fatou N, a fraktalny próg między zbiorami to zbiór Julia funkcji N. Newtonowski fraktal nie jest jedynie podobny do zbioru Julia — dosłownie jest to zbiór Julia iteracji funkcji Newtona, co wyjaśnia jego tę samą samo-podobieństwo i patologie graniczne jak te produkowane przez iterację zⁿ + c w pobliżu zbioru Mandelbrota.

Generowanie go

Obraz jest generowany piksel po pikselu: przyporządkuj każdy piksel do liczby zespolonej z₀, wykonaj krok Newtona przez określoną liczbę iteracji (lub do momentu, gdy |f(z)| spada poniżej tolerancji), a zarejestruj, który pierwiastek sekwencja zbliżyła się najbliżej oraz ile iteracji zajęło. Pierwiastek tożsamości wybiera odcienie; liczba iteracji wybiera jasność, która wyżyna finezę paski widoczne wewnątrz każdej miski, ponieważ punkty dalej od pierwiastka wymagają więcej kroków do zbiegu. Punkty dokładnie na granicy zbioru Julia nie zbiegają się nigdy, więc rzeczywiste implementacje ograniczają liczbę iteracji i kolorują wszystko, co wtedy jeszcze nie uległo zbiegnięciu, charakterystycznym kolorem "niezbieżny", który jest ciemną siatką widoczną w najdrobniejszej skali obrazu.

Frequently asked questions

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Newton's Fractal i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.