🪻 Fraktal metody Newtona — baseny przyciągania
Każdy piksel iteruje metodę Newtona w ℂ ku pierwiastkowi, a jego kolor zależy od tego, który pierwiastek osiągnie. Przełączaj wielomiany zⁿ−1 i z³−2z+2, strój relaksację i przybliżaj brzeg fraktala.
Podobne symulacje
O fraktalu Newtona
Metoda Newtona zastosowana do wielomianów zespolonych tworzy jedną z najpiękniejszych struktur fraktalnych matematyki. Zaczynając od dowolnej liczby zespolonej z₀, iteracja z ← z − a·f(z)/f′(z) (gdzie a to parametr relaksacji, domyślnie 1) zbiega do pierwiastka f(z) = 0. Kolorowanie każdego punktu startowego według tego, do którego pierwiastka dociera — i cieniowanie go według szybkości zbieżności — ujawnia diagram basenów przyciągania. Dla stopnia n ≥ 3 granice między basenami są fraktalnymi zbiorami Julii odwzorowania wymiernego N(z) = z − f(z)/f′(z), nieskończenie złożonymi i samopodobnymi w każdej skali.
Wybierz spośród kilku wielomianów (w tym „chaotycznego” z³ − 2z + 2, który ma przyciągające cykle 2-okresowe), dostrój parametr relaksacji i liczbę iteracji, a następnie kliknij w dowolnym miejscu, aby przybliżyć. Renderer wierszowy aktualizuje płótno w miarę kończenia każdego pasa, dzięki czemu widzisz, jak obraz powstaje w czasie rzeczywistym. Zapisz swój ulubiony widok fraktala jako PNG za pomocą przycisku w panelu sterowania.
Najczęściej zadawane pytania
Jak działa metoda Newtona na płaszczyźnie zespolonej?
Iteracja Newtona z ← z − f(z)/f′(z) została zaprojektowana dla liczb rzeczywistych, ale działa równie dobrze dla liczb zespolonych — algebra jest identyczna, tylko z arytmetyką zespoloną. Zaczynając od zespolonego punktu startowego z₀, ciąg zwykle spiralnie zmierza do jednego z pierwiastków wielomianu. Niezwykłe odkrycie polega na tym, że zbiór punktów startowych zbiegających do danego pierwiastka (basen przyciągania) ma fraktalną granicę, gdy wielomian ma stopień 3 lub więcej.
Co sprawia, że z³ − 2z + 2 jest „chaotyczny”?
Wielomian z³ − 2z + 2 ma trzy pierwiastki w pobliżu z ≈ −1,769, 0,885 ± 0,590i, ale odwzorowanie Newtona N(z) = z − f(z)/f′(z) dla tego wielomianu ma przyciągające cykle 2-okresowe — punkty, w których iteracja odbija się okresowo między dwiema wartościami, nigdy nie zbiegając do żadnego pierwiastka. Obszary płaszczyzny zespolonej uwięzione w tych cyklach pojawiają się jako ciemne, niezbiegające strefy, a granice między obszarami zbieżnymi i niezbieżnymi są szczególnie zawiłe.
Co robi parametr relaksacji a?
Zrelaksowana metoda Newtona wykorzystuje z ← z − a·f(z)/f′(z). Przy a = 1 jest to standardowa metoda Newtona ze zbieżnością kwadratową. Wartości a < 1 spowalniają zbieżność, ale mogą okiełznać rozbieżność. Wartości a > 1 przyspieszają zbieżność dla niektórych punktów startowych, ale powiększają chaotyczną granicę, w której nie osiąga się żadnego pierwiastka. Ustawienie a dalekiego od 1 może tworzyć zupełnie nowe struktury fraktalne, nieobecne w klasycznej metodzie.
Czym są pierwiastki z jedynki i dlaczego pojawiają się tak często?
N-te pierwiastki z jedynki to n rozwiązań równania zⁿ = 1, rozmieszczonych równomiernie na okręgu jednostkowym pod kątami 2πk/n dla k = 0, …, n−1. Są wygodnymi przypadkami testowymi, ponieważ mają proste, dokładne wartości, pochodna f′(z) = n·z^{n−1} łatwo się oblicza, a n-krotna symetria obrotowa sprawia, że fraktal również jest symetryczny. Zmiana suwaka stopnia z 3 na 9 dodaje więcej pierwiastków i basenów, zwiększając złożoność fraktala.
Dlaczego obraz renderuje się stopniowo, wiersz po wierszu?
Każdy piksel wymaga do maxIter kroków metody Newtona (do 80 w tej symulacji), a płótno 480×360 oznacza ~172 800 pikseli — potencjalnie miliony operacji zmiennoprzecinkowych. Aby uniknąć blokowania wątku interfejsu przeglądarki, renderer przetwarza 12 wierszy na klatkę animacji i oddaje kontrolę między partiami. Widzisz, jak fraktal pojawia się od góry do dołu; zmiana dowolnego ustawienia anuluje bieżące renderowanie i rozpoczyna nowe.
Jak działa przybliżanie po kliknięciu?
Każde kliknięcie ponownie wyśrodkowuje widok płaszczyzny zespolonej w klikniętym punkcie i zmniejsza skalę o połowę (skala × 0,5), podwajając powiększenie. Widok jest zdefiniowany przez środek (centerX + centerY·i) oraz szerokość (skala w jednostkach zespolonych). Na każdym poziomie powiększenia działa ta sama iteracja Newtona przy wyższej rozdzielczości — struktura fraktala powtarza się w każdej skali, ponieważ leżąca u jej podstaw matematyka jest niezmiennicza względem skali na granicy fraktala.
Jakie palety kolorów są dostępne i jak działają?
Dostępne są cztery palety. „Vivid” wykorzystuje nasycone kolory HSL z odcieniami rozłożonymi po całym kole barw, po jednym odcieniu na basen. „Pastel” używa niższego nasycenia dla łagodniejszego wyglądu. „Neon” podnosi nasycenie do 100% dla efektu świecenia. „Mono” mapuje baseny na tony skali szarości z lekkimi przesunięciami koloru. We wszystkich przypadkach szybciej zbiegające piksele są jaśniejsze (wyższa luminancja), więc granica fraktala — gdzie zbieżność jest najwolniejsza — pojawia się jako ciemny, misternie wzorzysty obszar.
Jaki jest związek między fraktalami Newtona a zbiorami Julii?
Odwzorowanie Newtona N(z) = z − f(z)/f′(z) jest funkcją wymierną stopnia 2n−2 dla wielomianu stopnia n. Jego zbiór Julii — granica między basenami przyciągania — to dokładnie fraktalna granica widoczna na obrazie fraktala Newtona. Dla z³ − 1 ten zbiór Julii ma wymiar Hausdorffa około 1,29. Zbiór Mandelbrota można traktować jako parametryzujący, które zbiory Julii są „gładkie” (spójne), a które „pyłowe” (podobne do zbioru Cantora).
Czy mogę zapisać mój obraz fraktala?
Tak — przycisk „Zapisz PNG” eksportuje bieżącą zawartość płótna w wewnętrznej rozdzielczości renderowania (do 560×420 pikseli). Obraz zawiera dokładne kolory pikseli obliczone przez iterację Newtona przy bieżącym widoku, wielomianie, palecie i ustawieniach relaksacji. Aby uzyskać eksport w wyższej rozdzielczości, trzeba by zwiększyć RW i RH w kodzie źródłowym i zaakceptować dłuższy czas renderowania.
Czym jest teoria Douady-Hubbarda istotna dla tych fraktali?
Adrien Douady i John Hubbard opracowali rygorystyczne ramy matematyczne dla dynamiki wielomianów zespolonych w latach 80. XX wieku, dowodząc, że zbiór Mandelbrota jest spójny, i klasyfikując topologię zbiorów Julii. Ich praca ma bezpośrednie zastosowanie do fraktali Newtona: każdy basen przyciągania odpowiada natychmiastowemu basenowi zawierającemu jeden pierwiastek, a granica fraktala jest okręgiem topologicznym (dla wielomianów, których wszystkie pierwiastki leżą w tym samym komponencie Fatou) lub zbiorem Cantora. Ta teoria wyjaśnia, dlaczego każdy punkt graniczny jest jednocześnie dostępny ze wszystkich basenów.
Każdy piksel iteruje metodę Newtona w ℂ w kierunku pierwiastka, kolorowany według tego, do którego pierwiastka dociera. Przełączaj wielomiany zⁿ−1 i z³−2z+2, dostrajaj relaksację i przybliżaj granicę fraktala.
3D · Renderer Three.js / WebGL · docelowo 60 FPS · działa w całości po stronie klienta, bez instalacji