The Mandelbrot and Julia sets
| The Mandelbrot and Julia sets |
 |
 |
 |
| Wednesday, 14 May 2008 00:00 | ||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Gaston Julia studied the iteration of polynomials and rational functions
in the early twentieth century. If f(x) is a function,
various behaviors can arise when f is iterated. Let's take,
for example, the function
We will iterate this function when initially applied to an initial value of x, say x = a0. Leta1 denote the first iterate f(a0), let a2 denote the second iteratef(a1), which equalsf(f(a0)), and so forth. Then we'll consider the infinite sequence of iterates It may happen that these values stay small or perhaps they don't, depending on the initial value a0. For instance, if we iterate our sample functionf(x) = x2 – 0.75 starting with the initial value a0 = 1.0, we'll get the following sequence of iterates (easily computed with a handheld calculator)
If you extend this table far enough, you'll see the iterates slowly approach the number –0.5. The iterates are above or below –0.5, but they get closer and closer to –0.5. In summary, when the initial value isa0 = 1.0, the iterates stay small, and, in particular, they approach –0.5.  It helps to see what's going on graphically. In the diagram above, the graphy = x2 – 0.75 of our function is drawn in magenta. Also, to help us see the iterates, the line y = x is drawn in green. Then the values a0, a1,a2, a3, anda4 are shown grapically, starting with our first value of a0, namely, 1.0. To find an+1 from an, follow the vertical blue line to the graph of the function (the magenta curve), then follow the horizontal red line across to the diagonal green liney = x, then follow another vertical blue line back to the x-axis. Notice how close a4 is to a2, but just a little to the right. The iteratea5 isn't shown on the graph, but it would be a little to the left of a3. If the rest of the iterates were shown, you would see a rectangular spiral in blue and red closing in on the point where the magenta curve and the green line intersect, that is, at the point (–0.5,–0.5). Thus, the iterates approach –0.5. With different initial values, different things can happen to the iterates. Let's take this time a0 = 2.0. Here's the table of iterates.
Here's the graph of the iterates starting with the initial valuea0 = 2.0. The graph is colored in the same way, but it's drawn to a different scale. This time, the iterates run off to infinity; they don't remain small.  Thus, what happens in the long run depends on the initial valuea0. It's easy to show that if the initial valuea0 is greater than 1.5, then the iterates run off to infinity, like they did when a0 was 2.0. Also, when the initial value a0 is less than –1.5, then, since after one iteration a1 is greater than 1.5, the iterates go off to infinity. With more work, you can show that if –1.5 < a0 < 1.5, then the iterates approach –0.5, like they did when a0 was 1.0. The remaining two initial values, namely 1.5 and –1.5, have all their iterates equal to 1.5. So the set of real numbers is partitioned into two parts, the closed interval [–1.5,1.5] of initial values whose iterates remain bounded, and the rest of the real numbers whose iterates do not remain bounded but approach infinity. This is shown graphically below. The real line is partitioned into the two sets with the interval [–1.5,1.5] drawn in light yellow, the rest in pink. We'll call the set of initial values whose iterates remain bounded the filled-in Julia set, or simply the Julia set for short. Thus, the light yellow set is the real Julia set for the functionf(x) = x2 – 0.75.  Complex numbers. We'll continue studying this same function, but extend our comprehension to complex numbers. For background on complex numbers, see Dave's Short Course on Complex Numbers. When we go to complex functions, we'll use z as the independent variable to emphasize that we're working with complex numbers. So, we write our sample function asf(z) = z2 – 0.75. What happens when the initial value a0 is a complex number instead of a real number? Let's take a0 to be the imaginary number i for an example. The table of iterates of this function looks like this.
Ah, let's stop there. We have the real number –1.75, and we know from the analysis above that it's iterates go to infinity since –1.75 is not in the range from –1.5 to 1.5. Let's try another imaginary number, say 0.5i.
Again, we can stop because we found that the iterate a1, which equals –1, does lie in the range from –1.5 to 1.5. Therefore, it's iterates converge. We could try other pure imaginary numbers ci where c is a real number. We would find that the iterates are bounded whenever –c2 – 0.75 lies in the range from –1.5 to 1.5. A little algebra shows that occurs whenc2 is less than or equal to 0.75, that is, whenci is between 0.866i and –0.766i. Therefore, on the complex axis, we should color the imaginary numbers between o.866i and –0.866i light yellow, the others pink.  We have begun to explore the complex Julia set for this function, but only just begun. Let's try another complex number for a0. Let's take a0 to be the complex number 1/3 + 5i/9. We'll get this table of iterates.
Well, it's unclear what's going to happen in the long run from looking at just these first five iterates. It looks like it's time to turn the computation over to a computer. What we need a computer to do is try lots of different initial complex valuesa0. For each value, iterate the function until it's clear whether the iterates approach infinity. It isn't hard to find how far away from 0 an value has to be in order that further iterates run away to infinity, but it's hard to conclude that all the iterates remain bounded. So what's done here is look at a limited number of iterates. If one gets far enough from 0, then we know the iterates will approach infinity. But if all remain sufficiently close to 0, then we won't draw any certain conclusion. Just below the result of that computation for our functionf(z) = z2 – 0.75. The image shows the part of the complex plane including numbersz = x + iy where both x and y are between ±1.5 is shown.) The black part of the plane indicates that the iterates did not get far from 0, so for an initial valuea0 in that region, all the iterates (probably) remain bounded. The black region includes the parts of the axes in the previous image that were drawn in light yellow. Most of the rest of the plane is colored red shading to some other colors. In this colored region of the plane, the iterates soon left toward infinity. Different colors indicate how soon, with red being quite quickly, yellow taking longer, and green, blue, and magenta taking much longer. This colored region includes the parts of the axes in the previous image that were drawn in pink.  That finishes our example of iterating a function. This was a quadratic function gave an interesting Julia set. Linear functions likef(z) = 3z + 5 don't yield interesting partitions of the complex plane, but quadratic and higher degree polynomials do as do other nonlinear functions.
A natural source of iterated functions comes from the approximation of
roots of functions by Mandelbrot Sets
Consider a whole family of functions parameterized by a variable. Although any family of functions
can be studied, we'll look at the most studied family, that being the family of quadratic polynomialsf(x) = x2 - µ, where µ is a complex parameter. As µ
varies, the Julia set will vary on the complex plane. Some of these Julia sets will be connected,
and some will be disconnected, and so this character of the Julia sets will partition the
µ-parameter plane into two parts. Those values of µ for which the Julia set is connected
is called the Mandelbrot set in the parameter plane. The boundary between the Mandelbrot set and
its complement is often called the Mandelbrot separator curve. The Mandelbrot set is the black shape
in the picture. This is the portion of the plane where x varies from -1 to 2 and y
varies between -1.5 and 1.5. There are some surprising details in this image, and it's well worth exploring.
The bulk of the Mandelbrot set is the black cardioid. (A cardioid is a heart-shaped figure). It's
studded with circles all around its boundary.
If you look close, you'll see the strands of dark blue above the circle under discussion. So, what's going on up there? Here's a blowup of that portion of the figure.
Aha! There's something black up near the top of the picture. What's that?
Each point C in the parameter space specifies the geometric structure of the corresponding Julia set J(C) in the dynamical plane z. If C1 is in the Mandelbrot set, the Julia set is connected. If C2 is not in the Mandelbrot set, the Julia set is a Cantor dust. You see connected 'Douady's rabbit' and 'Cantor dust' below.
More fractals
Alfeld, P. "The Mandelbrot Set." http://www.math.utah.edu/~alfeld/math/mandelbrot/mandelbrot.html.
Bowen, J. P. "Virtual Museum of Computing Mandelbrot Exhibition." http://archive.comlab.ox.ac.uk/other/museums/computing/mandelbrot.html.
Branner, B. "The Mandelbrot Set." In Chaos and Fractals: The Mathematics Behind the Computer Graphics,
Proc. Sympos. Appl. Math., Vol. 39 (Ed. R. L. Devaney and
L. Keen). Providence, RI: Amer. Math. Soc., 75-105, 1989.
Devaney, R. "The Mandelbrot Set and the Farey Tree, and the Fibonacci Sequence."
Amer. Math. Monthly 106, 289-302, 1999.
Dickau, R. M. "Mandelbrot (and Similar) Sets." http://mathforum.org/advanced/robertd/mandelbrot.html.
Douady, A. "Julia Sets and the Mandelbrot Set." In The Beauty of Fractals: Images of Complex Dynamical Systems
(Ed. H.-O. Peitgen and D. H. Richter). Berlin: Springer-Verlag, p. 161,
1986.
Elenbogen, B. and Kaeding, T. "A Weak Estimate of the Fractal Dimension of the
Mandelbrot Boundary." Phys. Lett. A 136, 358-362, 1989.
Eppstein, D. "Area of the Mandelbrot Set." http://www.ics.uci.edu/~eppstein/junkyard/mand-area.html.
Ewing, J. H. and Schober, G. "The Area of the Mandelbrot Set." Numer.
Math. 61, 59-72, 1992.
Fisher, Y. and Hill, J. "Bounding the Area of the Mandelbrot Set." Submitted.
Gleick, J. Chaos: Making a New Science. New York: Penguin Books, center
plate (following p. 114), 1988.
Hill, J. R. "Fractals and the Grand Internet Parallel Processing Project." Ch. 15 in Fractal Horizons: The Future Use of Fractals. New York:
St. Martin's Press, pp. 299-323, 1996.
Lauwerier, H. Fractals: Endlessly Repeated Geometric Figures. Princeton,
NJ: Princeton University Press, pp. 148-151 and 179-180, 1991.
Lei, T. (Ed.) The Mandelbrot Set, Theme and Variations. Cambridge, England:
Cambridge University Press, 2000.
Lesmoir-Gordon, N.; Rood, W.; and Edney, R. Introducing Fractal Geometry. Cambridge, England: Icon
Books, p. 97, 2000.
Mitchell, K. "A Statistical Investigation of the Area of the Mandelbrot Set." 2001. http://www.fractalus.com/kerry/articles/area/mandelbrot-area.html.
Munafo, R. "Mu-Ency--The Encyclopedia of the Mandelbrot Set." http://www.mrob.com/pub/muency.html.
Munafo, R. "Area of the Mandelbrot Set." http://www.mrob.com/pub/muency/areaofthemandelbrotset.html.
Munafo, R. "Pixel Counting." http://www.mrob.com/pub/muency/pixelcounting.html.
Peitgen, H.-O. and Saupe, D. (Eds.). The Science of Fractal Images. New York: Springer-Verlag,
pp. 178-179, 1988.
Shishikura, M. "The Boundary of the Mandelbrot Set has Hausdorff Dimension Two."
Astérisque, No. 222, 7, 389-405, 1994.
Sloane, N. J. A. Sequences A054670, A054671, A098403, and A114448 in "The On-Line Encyclopedia of Integer Sequences."
Wells, D. The Penguin Dictionary of Curious and Interesting Geometry.
London: Penguin, pp. 146-148, 1991.
Now we can go to a Julia orbit trip.
600x600 Mandelbrot set.
See also the
Mandelbrot Set
in Eric W. Weisstein's "World of |
Newton
Subscribe to this comment's feed
