Wpisany przez Tomasz Lubiński
piątek, 16 stycznia 2009 21:20
Fraktal Płonący Statek (ang. Burning Ship) został po raz pierwszy opisany w 1992 roku. Jego twórcami są Michael Michelitsch oraz Otto Rössler. By opisać fraktal, zdefiniujemy najpierw dla danego punktu p na płaszczyźnie zespolonej nieskończony ciąg liczb zespolonych z0, z1, z2, ... o wartościach zdefiniowanych następująco:
z0 = 0
zn+1 = ( | Re(zn) | + i | Im(zn) | )2 + p
Widać tutaj podobieństwo do wzoru na ciąg opisujący zbiór Mandelbrota (ang. Mandelbrot Set), różnica polega na tym, że w przypadku Płonącego Statku do potęgi podnoszona jest liczba o wartościach bezwzględnych części rzeczywistej oraz urojonej. Fraktalem jest brzeg zbioru liczb zespolonych p takich, że zdefiniowany powyżej ciąg nie dąży do nieskończoności. W praktyce by narysować fraktale oblicza się kolejne przybliżenia zbioru, które oznacza się różnymi kolorami. I tak kolejne przybliżenia zdefiniujemy jako zbiór liczb zespolonych p takich, że:
z0 = 0
zn+1 = ( | Re(zn) | + i | Im(zn) | )2 + p
Widać tutaj podobieństwo do wzoru na ciąg opisujący zbiór Mandelbrota (ang. Mandelbrot Set), różnica polega na tym, że w przypadku Płonącego Statku do potęgi podnoszona jest liczba o wartościach bezwzględnych części rzeczywistej oraz urojonej. Fraktalem jest brzeg zbioru liczb zespolonych p takich, że zdefiniowany powyżej ciąg nie dąży do nieskończoności. W praktyce by narysować fraktale oblicza się kolejne przybliżenia zbioru, które oznacza się różnymi kolorami. I tak kolejne przybliżenia zdefiniujemy jako zbiór liczb zespolonych p takich, że:
- 1 przybliżenie: wszystkie punkty
- 2 przybliżenie: |z1| < 2
- 3 przybliżenie: |z1| < 2 oraz |z2| < 2
- 4 przybliżenie: |z1| < 2 oraz |z2| < 2 oraz |z3| < 2
- ...
- n-te przybliżenie: |z1| < 2 oraz |z2| < 2, ... |zn-1| < 2
Zatem funkcję obliczającą z jakim maksymalnym przybliżeniem dany punkt p należy do zbioru Płonącego Statku możemy zdefiniować następująco (gdzie maxIter to maksymalne przybliżenie z jakim chcemy wyznaczać zbiór):
przyblizenie(p)
begin
iter := 0;
z := 0;
repeat
iter := iter + 1;
z = (|zr| + i|zi|)^2 + p;
until (|z| < 2) and (iter < maxIter)
przyblizenie = iter;
end;
Przypomnijmy jeszcze działania na liczbach zespolonych jakie będziemy potrzebować podczas obliczeń. Liczba zespolona z składa się z części rzeczywistej zr oraz części urojonej zi, czyli z = zr + i zi.
Potęgowanie definiujemy następująco:
z2 = (zr2 - zi2) + i(2 zr zi)
Dodawanie definiujemy następująco:
a + b = (ar + br) + i(ai + bi)
Moduł z liczby zespolonej definiujemy następująco:
,dlatego też w praktyce warunek |z| < 2 zastępuje się równoważną nierównością (zr2 + zi2) < 4. Pozbywamy się tutaj czasochłonnego obliczania pierwiastka kwadratowego.
W funkcji przedstawionej powyżej nie należy mylić modułu z liczby zespolonej z wartością bezwzględną jej części rzeczywistej: |zr| oraz części urojonej |zi|
Dla kolejnych punktów na płaszczyźnie, obliczamy przybliżenia zgodnie z podanym algorytmem i wzorami. Oś X oznacza wartości rzeczywiste, natomiast os Y wartości urojone. Przedstawiając kolejne przybliżenia na płaszczyźnie (lewy górny róg ma współrzędne -1.85 + -0.1i, dolny prawy róg ma współrzędne -1.7 + 0.05i) i oznaczając je kolorem różnej jasności otrzymujemy wynik - Fraktal Płonący Statek. Na obrazie poniżej barwy kolejnych przybliżeń wyznaczono zgodnie z modelem HSV, zmieniając parametr V (jasność) logarytmicznie w zależności od przybliżenia do którego należy dany punkt. Można też użyć do tego celu odcieni szarości, bądź innego modelu barw.
Przykład w JavaScript:
Zaznaczając obszar uzyskasz jego powiększony obraz. Kliknięcie prawym klawiszem (bądź dotknięcie dwoma palcami na urządzeniach z ekranem dotykowym) spowoduje powtórne pokazanie całego zbioru.
| Autor | Język programowania | Komentarz | Otwórz | Pobierz | Ocena |
| Tomasz Lubiński | C# | MS Visual Studio .net |  | |
  / 1 |
| Tomasz Lubiński | C/C++ | Borland Builder 6 |  | |
/ 1 |
| Tomasz Lubiński | Delphi/Pascal | Borland Delphi 5 |  | |
/ 1 |
| Tomasz Lubiński | Java Script | Firefox 3.0+, Safari 3.0+, Chrome 3.0+, Opera 9.5+, IE 9.0+ |  | |
/ 0 |
Poprawiony: wtorek, 21 czerwca 2011 19:21