The story so far
- Classic Mandelbrot and Julia
- More on the Classic Sets
- Other Polynomials
- Trigonometric and Exponential Functions
After a long pause, here is the fifth instalment of this series on fractals with Octave. This time I will look at a fractal called the Burning Ship because of its output that looks like a ship going up in flames.
Burning Ship
The fractal is a Mandelbrot set that has a generating series defined by:
zn+1=(|Re(zn)|+i|Im(zn)|)2+c
where the Re and Im are the functions that return the real and imaginary parts of a complex number. The initial condition is:
z0=0
This fractal can be directly produced with the mandelbrot
function we created in previous articles. In addition, we will take r=200 as a divergence limit.
octave-3.0.1:1> Mbs=mandelbrot(-2.5-2i,1.5+i,320,64, > @(z,c) (abs(real(z))+i*abs(imag(z))).^2.+c,200); octave-3.0.1:2> imagesc(Mbs)
Octave should open the Gnuplot window and display an image similar to this:
Now let's zoom on this fractal centre left by a factor of 30:
octave-3.0.1:1> Mbsz=mandelbrot(-1.8-0.06i,-1.7+0.02i,320,64, > @(z,c) (abs(real(z))+i*abs(imag(z))).^2.+c,200); octave-3.0.1:2> imagesc(Mbsz)
Octave should open the Gnuplot window and display an image similar to this:
As usual with a fractal like this, the zoomed image looks very similar to the original. Now let's have a look at a related fractal called the Bird of Prey.
Bird of Prey
In the same way that you can produce interesting variations on the classic Mandelbrot set by changing the power used in the series, we can do the same with the Burning Ship series. So let's increase the power in the series above from 2 to 3:
zn+1=(|Re(zn)|+i|Im(zn)|)3+c
This cubic series produces a fractal called the Bird of Prey.
octave-3.0.1:1> Mbs3=mandelbrot(-1.5-1.5i,1.5+1.5i,320,64, > @(z,c) (abs(real(z))+i*abs(imag(z))).^3.+c,200); octave-3.0.1:2> imagesc(Mbs3)
You'll note that I chose the boundary values for the complex plane so that I would get a square image centered on the origin. Here is the resulting image:
Next
That's it for the Burning Ship. Next in this series, we will look at the original fractal that Gaston Julia was interested in.
No comments:
Post a Comment