Mandelbrot with GNU Octave
In preparation for another project idea I have (to be posted in the future), I wrote a Mandelbrot fractal generator in Octave. Octave is great for just trying things out and prototyping your algorithms. It is very slow, however. Here is the code,
function mandel_img = mandel () ## Parameters w = [-2.5 1.5]; # Domain h = [-1.5 1.5]; # Range s = 0.005; # Step size it = 64; # Iteration depth ## Prepare the complex plane [wa ha] = meshgrid (w(1):s:w(2), h(1):s:h(2)); complex_plane = wa + ha * i; ## Preallocate image mandel_img = zeros( length(h(1):s:h(2)), length(w(1):s:w(2))); ## Generate mandelbrot for wi = 1:size(mandel_img, 2) for hi = 1:size(mandel_img, 1) z = 0; k = 0; while k < it && abs(z) < 2 z = z^2 + complex_plane (hi, wi); k = k + 1; end mandel_img (hi, wi) = k - 1; end ## Display progress waitbar (wi/size(mandel_img, 2)); end end
You may need to comment out the
(This information is out of date.)
waitbar line if you do not have
Octave-Forge installed properly (as is the case with Octave
2.9 on Debian as of this writing) or at all. You will also need
Octave-Forge if you want to use the image functions described below.
You can find the same code all over the Internet for many different
languages. The advantage with Octave is that it knows about complex
numbers so that this can be expressed directly with
z = z^2 + c and
Now, this code just generates a matrix of the escape iteration numbers for each pixel. To visualize this, you will need to use the image functions. The simplest thing to do is view the data as a boring greyscale image.
octave> m = mandel(); # Generate the data octave> imshow(m);
You should see something like this,
You can save this as an image with
octave> imwrite("mandel.png", m*4)
*4 part is because the iteration depth was set to 64. The image
being written will have values between 0 and 255. This allows the data
to use the full dynamic range of the image.
If you want more interesting images, you can apply a
colormap. Octave-Forge has two handy color maps,
ocean (cool). To make the inside of the fractal
black, which are the points that are part of the set and never
escaped, stick black on the end of the colormap. This can be done like
this (viewing and saving),
octave> cmap = [hot(63); 0 0 0]; # The colormap octave> imshow(m + 1, cmap); octave> imwrite("mandel.png", m + 1, cmap);
m is between 0 and 63. We add one to it to put it between 1 and 64.
Then we take the colormap of length 63 and stick black on the end. If
hot, you will get a nice blue version.
You can modify the code above to try to get different fractals. For
z = z^4 + c instead,
More on fractals another time.blog comments powered by Disqus