Fraktal Metody Newtona: Baseny Przyciągania w Trakcie Złożonym

Metoda Newtona, uruchamiana na każdym pikselu

Metoda Newtona znajduje pierwiastek funkcji f poprzez iteracyjne śledzenie stycznej w dół do punktu przecięcia z zerem: zaczynając od szacunku z, następny szacąg to z - f(z)/f'(z). Dla funkcji rzeczywistej konwergencja jest szybka w pobliżu pierwiastka i może zachować się nieprawidłowo daleko od niego. Uruchom tę samą iterację z początkowymi wartościami złożonymi z = x + iy, a każdy piksel początkowy zostanie skolorowany przez pierwiastek, do którego jego iteracja zbiega się, co daje Newtonowską fraktalną - obraz zbiorów przyciągania dla każdego pierwiastka, w którym wszystkie zwykłe gwarancje konwergencji zawiodą spektakularnie na granicach zbioru.

Iteracja i konkretny wielomian

Dla najprostszej interesującej się przypadki, f(z) = z^n - 1, którego n pierwiastków to n-te potęgi jedności równomiernie rozmieszczone na okręgu jednostkowym, metoda Newtona upraszcza się do:

z_(k+1) = z_k - (z_k^n - 1) / (n * z_k^n + 1) = ((n - 1) * z_k^n + 1) / (n * z_k^n + 1) Każdy piksel płaszczyzny jest niezależnym wartością początkową z0. Iteruj przez określoną liczbę razy (zwykle 30-60, wystarczającą do prawie wszędzie konwergencji, poza granicami), sprawdź, do którego z n pierwiastków ostateczna wartość jest najbliżej, i pomal piksel tym kolorem pierwiastka. Cieniowanie w zależności od liczby iteracji do konwergencji, zamiast jednolitego koloru, przekształca płaskie baseny w pasmowane pierścienie konturowych kręgów, które ujawniają strukturę fraktalną ukrytą w obszarze wyglądającym na jednolity kolor.

z_(k+1) = z_k - (z_kⁿ - 1) / (n·z_kⁿ⁻¹)
        = ((n-1)·z_kⁿ + 1) / (n·z_kⁿ⁻¹)

Dlaczego granice są fraktalne

W pobliżu pierwiastka metoda Newtona zbiega się kwadratowo – liczba poprawnej cyfr rośnie w przybliżeniu podwójnie na każdym kroku – więc głęboko wewnątrz basenu obraz jest gładki i nudny. Najbardziej interesującą strukturę tworzą całkowicie granice między basenami. Ponieważ funkcja iteracyjna z - f(z)/f'(z) jest funkcją racjonalną z sfery Riemanna, jej zbiór Julia (granica, gdzie bliskie punkty początkowe rozbiegają się do różnych atrakktorów) ma samopodobną, nieskończenie szczegółową strukturę typową dla dynamiki złożonej: powiększając dowolny punkt graniczny, znajdujemy mniejsze kopie tego samego splątania trzech lub więcej basenów, na zawsze.

Jest to specjalny przypadek starszej i zaskakującej rzeczy z dynamiki złożonej: granica między obszarami przyciągania dla dowolnej funkcji racjonalnej z trzema lub więcej odrębnymi stałymi punktami przyciągającymi jest albo połączona, albo ma nieskończenie wiele komponentów, a w przypadku metody Newtona dla ogólnego sześcianu trzy baseny słynnie dzielą tę samą granicę – każdy punkt graniczny znajduje się arbitralnie blisko wszystkich trzech basenów jednocześnie, nie tylko dwóch z nich, to naprawdę poczucie trójwymiarowej patologii żyjącej w obrazie dwuwymiarowym.

Gdzie metoda Newtona w ogóle nie zbiega się

Nie każda wartość początkowa ugruntowuje się. Jeśli pochodna f'(z) jest dokładnie równa zeru podczas iteracji, aktualizacja eksploduje; częściej wartość początkowa wpada w cykliczny układ, który nigdy nie osiąga pierwiastka – dla z^19 - 2z + 2, znana wartość początkowa z0 = 0 słynnie ulega cyklowi pomiędzy 1 a 0 bezbytnemu, klasyczny przykład podważający naiwne twierdzenia o zbieżności. Obraz fraktali oznaczają takie punkty (lub punkty, które nie zbiegają w ramach budżetu iteracji) charakterystycznym kolorem 'brak zbieżności', zwykle czarnym, i to właśnie tam granice basenów są najbardziej gęste i skomplikowane.

Rozluźnione Newtonowskie i zmienne wielomiany

Współczynnik tłumienia (relaksacji) lambda wpływa na aktualizację, zmieniając ją z postaci z(k+1) = z_k - lambda * f(z)/f'(z). Przy lambda=1 otrzymujemy standardowy Newtona; mniejsze lambda spowalniają zbieżność, ale mogą skurczyć lub zmienić kształt granicy fraktali, a złożone lambda powodują pełne obroty zbioru barwnego, ponieważ zmieniają kierunek, w jakim każdy krok jest dopuszczalny względem linii stycznej. Zmiana wielomianu z z-1 na asymetryczny taki jak z-2z+2 usuwa symetrię rotacyjną przypadku pierwiastków jedności i generuje widocznie inny, mniej symetryczny węzeł zbiorów barwnych - demo pozwala na przełączanie się między obiema opcjami i ponowne renderowanie w czasie rzeczywistym.

Frequently asked questions

Wypróbuj na żywo

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