Generating Fractals With Complex Numbers

Learning Outcomes

  • Identify the difference between an imaginary number and a complex number
  • Identify the real and imaginary parts of a complex number
  • Plot a complex number on the complex plane
  • Perform arithmetic operations on complex numbers
  • Graph physical representations of arithmetic operations on complex numbers as scaling or rotation
  • Generate several terms of a recursive relation
  • Determine whether a complex number is part of the set of numbers that make up the Mandelbrot set

Complex Recursive Sequences

Some fractals are generated with complex numbers.  The Mandlebrot set, which we introduced briefly at the beginning of this module, is generated using complex numbers with a recursive sequence. Before we can see how to generate the Mandelbrot set, we need to understand what a recursive sequence is.

Recursive Sequence

A recursive relationship is a formula which relates the next value, [latex]{{z}_{n+1}}[/latex], in a sequence to the previous value, [latex]{{z}_{n}}[/latex]. In addition to the formula, we need an initial value, [latex]{{z}_{0}}[/latex].

The sequence of values produced is the recursive sequence.

Example

Given the recursive relationship [latex]{{z}_{n+1}}={{z}_{n}}+2,\quad{{z}_{0}}=4[/latex], generate several terms of the recursive sequence.

We are given the starting value, [latex]{{z}_{0}}=4[/latex]. The recursive formula holds for any value of n, so if [latex]n = 0[/latex], then [latex]{{z}_{n+1}}={{z}_{n}}+2[/latex] would tell us [latex]{{z}_{0+1}}={{z}_{0}}+2[/latex], or more simply, [latex]{{z}_{1}}={{z}_{0}}+2[/latex].

Notice this defines [latex]{{z}_{1}}[/latex] in terms of the known [latex]{{z}_{0}}[/latex], so we can compute the value:

[latex]{{z}_{1}}={{z}_{0}}+2=4+2=6[/latex].

Now letting [latex]n = 1[/latex], the formula tells us [latex]{{z}_{1+1}}={{z}_{1}}+2[/latex], or [latex]{{z}_{2}}={{z}_{1}}+2[/latex]. Again, the formula gives the next value in the sequence in terms of the previous value.

[latex]{{z}_{2}}={{z}_{1}}+2=6+2=8[/latex]

Continuing,

[latex]{{z}_{3}}={{z}_{2}}+2=8+2=10\\{{z}_{4}}={{z}_{3}}+2=10+2=12[/latex]

Try It

 

The previous example generated a basic linear sequence of real numbers. The same process can be used with complex numbers.

Example

Given the recursive relationship [latex]{{z}_{n+1}}={{z}_{n}}\cdot{i}+(1-i),\quad{{z}_{0}}=4[/latex], generate several terms of the recursive sequence.

We are given [latex]{{z}_{0}}=4[/latex]. Using the recursive formula:

[latex]{{z}_{1}}={{z}_{0}}\cdot{i}+(1-i)=4\cdot{i}+(1-i)=1+3i\\{{z}_{2}}={{z}_{1}}\cdot{i}+(1-i)=(1+3i)\cdot{i}+(1-i)=i+3{{i}^{2}}+(1-i)=i-3+(1-i)=-2\\{{z}_{3}}={{z}_{2}}\cdot{i}+(1-i)=(-2)\cdot{i}+(1-i)=-2i+(1-i)=1-3i\\{{z}_{4}}={{z}_{3}}\cdot{i}+(1-i)=(1-3i)\cdot{i}+(1-i)=i-3{{i}^{2}}+(1-i)=i+3+(1-i)=4\\{{z}_{5}}={{z}_{4}}\cdot{i}+(1-i)=4\cdot{i}+(1-i)=1+3i[/latex]

Notice this sequence is exhibiting an interesting pattern—it began to repeat itself.

Try It

In the following video we show more worked examples of how to generate the terms of a recursive, complex sequence.

Mandelbrot Set

The Mandelbrot Set is a set of numbers defined based on recursive sequences.

Mandelbrot Set

For any complex number c, define the sequence [latex]{{z}_{n+1}}={{z}_{n}}^{2}+c,\quad{{z}_{0}}=0[/latex]

If this sequence always stays close to the origin (within 2 units), then the number c is part of the Mandelbrot Set. If the sequence gets far from the origin, then the number c is not part of the set.

Example

Determine if [latex]c=1+i[/latex] is part of the Mandelbrot set.

We start with [latex]{{z}_{0}}=0[/latex]. We continue, omitting some detail of the calculations

[latex]{{z}_{1}}={{z}_{0}}^{2}+1+i=0+1+i=1+i\\{{z}_{2}}={{z}_{1}}^{2}+1+i={{(1+i)}^{2}}+1+i=1+3i\\{{z}_{3}}={{z}_{2}}^{2}+1+i={{(1+3i)}^{2}}+1+i=-7+7i\\{{z}_{4}}={{z}_{3}}^{2}+1+i={{(-7+7i)}^{2}}+1+i=1-97i[/latex]

We can already see that these values are getting quite large. It does not appear that [latex]c=1+i[/latex] is part of the Mandelbrot set.

Example

Determine if [latex]c=0.5i[/latex] is part of the Mandelbrot set.

We start with [latex]{{z}_{0}}=0[/latex]. We continue, omitting some detail of the calculations

[latex]{{z}_{1}}={{z}_{0}}^{2}+0.5i=0+0.5i=0.5i\\{{z}_{2}}={{z}_{1}}^{2}+0.5i={{(0.5i)}^{2}}+0.5i=-0.25+0.5i\\{{z}_{3}}={{z}_{2}}^{2}+0.5i={{(-0.25+0.5i)}^{2}}+0.5i=-0.1875+0.25i\\{{z}_{4}}={{z}_{3}}^{2}+0.5i={{(-0.1875+0.25i)}^{2}}+0.5i=-0.02734+0.40625i[/latex]

While not definitive with this few iterations, it does appear that this value is remaining small, suggesting that 0.5i is part of the Mandelbrot set.

Try It

Determine if [latex]c=0.4+0.3i[/latex] is part of the Mandelbrot set.

A shaded fractal shown on a graph.If all complex numbers are tested, and we plot each number that is in the Mandelbrot set on the complex plane, we obtain the shape to the right.[1]

The boundary of this shape exhibits quasi-self-similarity, in that portions look very similar to the whole.

Watch the following video for more examples of how to determine whether a complex number is a member of the Mandelbrot set.

In addition to coloring the Mandelbrot set itself black, it is common to the color the points in the complex plane surrounding the set. To create a meaningful coloring, often people count the number of iterations of the recursive sequence that are required for a point to get further than 2 units away from the origin. For example, using [latex]c=1+i [/latex] above, the sequence was distance 2 from the origin after only two recursions.

For some other numbers, it may take tens or hundreds of iterations for the sequence to get far from the origin. Numbers that get big fast are colored one shade, while colors that are slow to grow are colored another shade. For example, in the image below[2], light blue is used for numbers that get large quickly, while darker shades are used for numbers that grow more slowly. Greens, reds, and purples can be seen when we zoom in—those are used for numbers that grow very slowly.

The Mandelbrot set, for having such a simple definition, exhibits immense complexity. Zooming in on other portions of the set yields fascinating swirling shapes.

 

 

Additional Resources

A much more extensive coverage of fractals can be found on the Fractal Geometry site. This site includes links to several Java software programs for exploring fractals.

If you are impressed with the Mandelbrot set, check out this TED talk from 2010 given by Benoit Mandelbrot on fractals and the art of roughness.


  1. http://en.wikipedia.org/wiki/File:Mandelset_hires.png
  2. This series was generated using Scott’s Mandelbrot Set Explorer