Multibrot set

Multibrot 3 at the bottom-left of the main part.
Detail of Multijulia 8.
Multibrot 4.
Multibrot exponent 0 - 8

In mathematics, a Multibrot set is the set of values in the complex plane whose absolute value remains below some finite value throughout iterations by a member of the general monic univariate polynomial family of recursions.[1][2][3] The name is a portmanteau of multiple and Mandelbrot set. The same can be applied to the Julia set, this being called Multijulia set.

where d ≥ 2. The exponent d may be further generalized to negative and fractional values.[4]

Examples[5][6]

[edit]

The case of

is the classic Mandelbrot set from which the name is derived.

The sets for other values of d also show fractal images[7] when they are plotted on the complex plane.

Each of the examples of various powers d shown below is plotted to the same scale. Values of c belonging to the set are black. Values of c that have unbounded value under recursion, and thus do not belong in the set, are plotted in different colours, that show as contours, depending on the number of recursions that caused a value to exceed a fixed magnitude in the Escape Time algorithm.

Positive powers

[edit]

The example d = 2 is the original Mandelbrot set. The examples for d > 2 are often called multibrot sets. These sets include the origin and have fractal perimeters, with (d − 1)-fold rotational symmetry.

zz2 + c
zz3 + c
zz4 + c
zz5 + c
zz6 + c
zz96 + c
zz96 + c detail x40

Negative powers

[edit]

When d is negative the set appears to surround but does not include the origin, However this is just an artifact of the fixed maximum radius allowed by the Escape Time algorithm, and is not a limit of the sets that actually have a shape in the middle with an no hole (You can see this by using the Lyapunov exponent [No hole because the origin diverges to undefined not infinity because the origin {0 or 0+0i} taken to a negative power becomes undefined]). There is interesting complex behaviour in the contours between the set and the origin, in a star-shaped area with (1 − d)-fold rotational symmetry. The sets appear to have a circular perimeter, however this is an artifact of the fixed maximum radius allowed by the Escape Time algorithm, and is not a limit of the sets that actually extend in all directions to infinity.

zz−2 + c
zz−3 + c
zz−4 + c
zz−5 + c
zz−6 + c

Fractional powers

[edit]

Rendering along the exponent

[edit]

An alternative method is to render the exponent along the vertical axis. This requires either fixing the real or the imaginary value, and rendering the remaining value along the horizontal axis. The resulting set rises vertically from the origin in a narrow column to infinity. Magnification reveals increasing complexity. The first prominent bump or spike is seen at an exponent of 2, the location of the traditional Mandelbrot set at its cross-section. The third image here renders on a plane that is fixed at a 45-degree angle between the real and imaginary axes. [8]

Multibrot rendered with real value along horizontal axis and exponent along vertical axis, imaginary value fixed at zero
Multibrot rendered with imaginary value on horizontal axis and exponent on vertical axis, real value fixed at zero
Multibrot rendered with exponent on vertical axis along a plane angled 45-degrees between the real and imaginary axes.

Rendering images

[edit]

All the above images are rendered using an Escape Time algorithm that identifies points outside the set in a simple way. Much greater fractal detail is revealed by plotting the Lyapunov exponent,[9] as shown by the example below. The Lyapunov exponent is the error growth-rate of a given sequence. First calculate the iteration sequence with N iterations, then calculate the exponent as

and if the exponent is negative the sequence is stable. The white pixels in the picture are the parameters c for which the exponent is positive aka unstable. The colours show the periods of the cycles which the orbits are attracted to. All points colored dark-blue (outside) are attracted by a fixed point, all points in the middle (lighter blue) are attracted by a period 2 cycle and so on.

Enlarged first quadrant of the multibrot set for the iteration zz−2 + c rendered with the Escape Time algorithm.
Enlarged first quadrant of the multibrot set for the iteration zz−2 + c rendered using the Lyapunov exponent of the sequence as a stability criterion rather than using the Escape Time algorithm. Periodicity checking was used to colour the set according to the period of the cycles of the orbits.

Pseudocode

[edit]
ESCAPE TIME ALGORITHM  for each pixel on the screen do     x = x0 = x co-ordinate of pixel     y = y0 = y co-ordinate of pixel        iteration := 0     max_iteration := 1000        while (x*x + y*y ≤ (2*2) and iteration < max_iteration do         /* INSERT CODE(S)FOR Z^d FROM TABLE BELOW */         iteration := iteration + 1        if iteration = max_iteration then         colour := black     else         colour := iteration        plot(x0, y0, colour) 

The complex value z has coordinates (x,y) on the complex plane and is raised to various powers inside the iteration loop by codes shown in this table. Powers not shown in the table can be obtained by concatenating the codes shown.

z−2 z−1 z2
(for Mandelbrot set)
z3 z5 zn

d=x^4+2*x^2*y^2+y^4
assert d != 0
xtmp = (x^2-y^2)/d+a
y = -2*x*y/d+b
x = xtmp

d=x^2+y^2
assert d != 0
x = x/d + a
y= -y/d + b

xtmp=x^2-y^2 + a
y=2*x*y + b
x=xtmp

xtmp=x^3-3*x*y^2 + a
y=3*x^2*y-y^3 + b
x=xtmp

xtmp=x^5-10*x^3*y^2+5*x*y^4 + a
y=5*x^4*y-10*x^2*y^3+y^5 + b
x=xtmp

xtmp=(x*x+y*y)^(n/2)*cos(n*atan2(y,x)) + a
y=(x*x+y*y)^(n/2)*sin(n*atan2(y,x)) + b
x=xtmp

References

[edit]
  1. ^ "Definition of multibrots". Retrieved 2008-09-28.
  2. ^ "Multibrots". Retrieved 2008-09-28.
  3. ^ Wolf Jung. "Homeomorphisms on Edges of the Mandelbrot Set" (PDF). p. 23. The Multibrot set Md is the connectedness locus of the family of unicritical polynomials zd + c, d ≥ 2
  4. ^ "WolframAlpha Computation Knowledge Engine".
  5. ^ "23 pretty JavaScript fractals". 23 October 2008. Archived from the original on 2014-08-11.
  6. ^ "Javascript Fractals". Archived from the original on 2014-08-19.
  7. ^ "Animated morph of multibrots d = −7 to 7". Retrieved 2008-09-28.
  8. ^ Fractal Generator, "Multibrot Slice"
  9. ^ Ken Shirriff (Sep 1993). "An Investigation of Fractals Generated by z → 1/zn + c". Computers & Graphics. 17 (5): 603–607. doi:10.1016/0097-8493(93)90012-x. Retrieved 2008-09-28.