Classic Iterated Function Systems

Larry Riddle, Agnes Scott College

$Z_2$ Heighway Dragon

Description


 

Construction
Animation

Let $K = \{f_1,f_2\}$ be the two functions in the Golden Dragon IFS. Let $G = \{g_1,g_2\}$ be the $Z_2$ cyclic group where $g_1$ is the identity element and $g_2$ is rotation through 180°. We form a new IFS by applying each of the symmetries in G to each of the functions in K to form the set of affine transformations $GK = \{g_1f_1,g_2f_1,g_1f_2,g_2f_2\}$. The first iteration would transform the horizontal line segment into the following figure using the four functions in this set. The red line segments would be the results of the usual Heighway dragon IFS and the blue line segments are the rotations by 180°.

2 red and 2 blue line segments symmetric with respect to the origin.

Recall that \(r = \frac{1}{{\sqrt 2 }}\). Now continue to repeat the construction on each of the four new segments to get the $Z_2$ Heighway dragon.


[Enlarge]


Iterated
Function
System


 
 

IFS
Animation
(with squares)
The IFS for the $Z_2$ Heighway dragon will consist of four functions: the original two functions from the Heighway dragon IFS and the rotation of those two functions through 180°.
\( h_1(\mathbf{x})= g_1f_1(\mathbf{x})= \left[ \begin{array}{cc} 1/2 & -1/2\\ 1/2 & 1/2 \end{array} \right]\mathbf{x} \)
scale by $r$, rotate by 45°
\( h_2(\mathbf{x})= g_2f_1(\mathbf{x})= \left[ \begin{array}{cc} -1/2 & 1/2\\ -1/2 & -1/2 \end{array} \right]\mathbf{x} \)
scale by $r$, rotate by −135°
\( h_3(\mathbf{x})= g_1f_2(\mathbf{x})= \left[ \begin{array}{cc} -1/2 & -1/2\\ 1/2 & -1/2 \end{array} \right]\mathbf{x} + \left[ \begin{array}{c} 1\\ 0 \end{array} \right] \)
scale by $r$, rotate by 135°
\( h_4(\mathbf{x})= g_2f_2(\mathbf{x})= \left[ \begin{array}{cc} 1/2 & 1/2\\ -1/2 & 1/2 \end{array} \right]\mathbf{x} + \left[ \begin{array}{c} -1\\ 0 \end{array} \right] \)
scale by $r$, rotate by −45°

where \(r = \frac{1}{{\sqrt 2 }}\). If H is the unique attractor for this IFS, then

\[H = h_1(H) \cup h_2(H) \cup h_3(H) \cup h_4(H) \]

Rotating H by 180° counterclockwise is the same as applying the function $g_2$ from the cyclic group $Z_2$ to the set H. Because $g_2g_2 = g_1$ (identity) and $g_2g_1 = g_2$, we get

\[\begin{align} {g_2}(H) &= {g_2}\left( {{h_1}(H) \cup {h_2}(H) \cup {h_3}(H) \cup {h_4}(H)} \right)\\ \\ &= {g_2}{g_1}{f_1}(H) \cup {g_2}{g_2}{f_1}(H) \cup {g_2}{g_1}{f_2}(H) \cup {g_2}{g_2}{f_2}(H)\\ \\ &= {g_2}{f_1}(H) \cup {g_1}{f_1}(H) \cup {g_2}{f_2}(H) \cup {g_1}{f_2}(H)\\ \\ &= h_2(H) \cup h_1(H) \cup h_4(H) \cup h_3(H) = H \end{align}\]

This shows that H has 180° rotational symmetry. Also notice that $h_1(H) = h_2(H)$. That is because $h_2$ does the same scaling to H as $h_1$ but rotates the scaled object 135° clockwise rather than 45° counterclockwise. Because of the 180° rotational symmetry of H, however, the result will be the same in both cases. You can see what happens by clicking on each of the buttons to the left to cover H with the four scaled and rotated copies.

Special
Properties

The copy of the $Z_2$ Heighway dragon below lies on top of a grid of size 1/3 by 1/3. It demonstrates that the $Z_2$ Heighway dragon lies within a rectangle with $-7/3 \le x \le 7/3$ and $-5/3 \le x \le 5/3$, i.e. a rectangle of width 4.667 and height 3.333 [Rectangle Details].

The figure below shows the orbit of the point (0,0) under iteration by the fourth function in the IFS. The orbit spirals into the fixed point of $h_4$. That fixed point is located at (−1,1) since

\[{h_4}\left( {\left[ {\begin{array}{c}{ - 1}\\1\end{array}} \right]} \right) = \left[ {\begin{array}{cc}{1/2}&{1/2}\\{ - 1/2}&{1/2}\end{array}} \right]\left[ {\begin{array}{c}{ - 1}\\1\end{array}} \right] + \left[ {\begin{array}{c}{ - 1}\\0\end{array}} \right] = \left[ {\begin{array}{c}0\\1\end{array}} \right] + \left[ {\begin{array}{c}{ - 1}\\0\end{array}} \right] = \left[ {\begin{array}{c}{ - 1}\\1\end{array}} \right]\]

Spiral
Animation
This fixed point is the tip of the spiral of the dragon in the second quadrant. Watch the animation to see how the iteration of the IFS for an initial line segment follows the path of the points in the figure. By symmetry, the tip of the spiral of the dragon in the fourth quadrant is located at (1,−1).

As shown above, the $Z_2$ Heighway dragon has symmetric spirals in the second and the fourth quadrants. This next figure shows a close-up zoom of the spiral in the second quadrant (on a grid of size 1/6 x 1/6).

The white "lakes" inside the red dragon spiral through 45° angles (clockwise), with each lake just touching the one before and the one after as indicated at the dots in the figure. The dot in the upper right where the first lake begins is located at the point (−2/3,1) [Lake Details]. The other dots are the orbit of this point under iteration by the function $h_4$ in the IFS. The first few are (−5/6, 5/6), (−1, 5/6), (−13/12, 11/12), and (−13/12, 1). The orbit converges to the fixed point at (−1,1).

The white lakes in the fourth quadrant exhibit similar behavior with points that are symmetrically opposite from those in the second quadrant.

Variations

The IFS for the Lévy Dragon is very similar to that of the Heighway Dragon, yet the two fractals are quite different. It is perhaps surprising, then, that if the $Z_2$ symmetry group is applied to the Lévy IFS as it was to the Heighway IFS, the resulting figure is the same as the $Z_2$ Heighway dragon, only rotated clockwise by 135° and scaled by a factor of \(1/\sqrt{2}\).
$Z_2$ Heighway dragon
$Z_2$ Lévy dragon

The reason for this is that after composing the transformations in $Z_2$ with the functions in the Heighway IFS and the Lévy IFS, respectively, the new iterated function systems both involve scaling by \(1/\sqrt{2}\) and rotations by 45°, 135°, −45°, and −135°, so they have the same matrices in both systems. The difference is that the translation vectors for the $Z_2$ Lévy dragon are rotated 135° clockwise from those for the $Z_2$ Heighway dragon and scaled by \(1/\sqrt{2}\). This causes the attractor for the $Z_2$ Lévy dragon to be rotated clockwise by 135° and scaled by \(1/\sqrt{2}\) relative to the $Z_2$ Heighway dragon [Rotation Proof].

Red Heighway vectors along x-axis and blue Levy vectors rotated 45 degrees
Red: Heighway
Blue: Lévy


IFS
Animation

Here is the IFS for the $Z_2$ Lévy dragon, where $r = \frac{1}{\sqrt 2}$.

\( k_1(\mathbf{x})= \left[ \begin{array}{cc} 1/2 & -1/2\\ 1/2 & 1/2 \end{array} \right]\mathbf{x} \)
scale by $r$, rotate by 45°
\( k_2(\mathbf{x})= \left[ \begin{array}{cc} -1/2 & 1/2\\ -1/2 & -1/2 \end{array} \right]\mathbf{x} \)
scale by $r$, rotate by −135°
\( k_3(\mathbf{x})= \left[ \begin{array}{cc} 1/2 & 1/2\\ -1/2 & 1/2 \end{array} \right]\mathbf{x} + \left[ \begin{array}{c} 1/2\\ 1/2 \end{array} \right] \)
scale by $r$, rotate by −45°
\( k_4(\mathbf{x})= \left[ \begin{array}{cc} -1/2 & -1/2\\ 1/2 & -1/2 \end{array} \right]\mathbf{x} + \left[ \begin{array}{c} -1/2\\ -1/2 \end{array} \right] \)
scale by $r$, rotate by 135°

Note that the matrix for function $k_3$ corresponds to the matrix for function $h_4$, and that the matrix for function $k_4$ corresponds to the matrix for function $h_3$. Also notice that the translation vectors for these functions satisfy the following calculations for a clockwise rotation of 135° and a scaling by \(1/\sqrt{2}\):

\(\dfrac{1}{{\sqrt 2 }}\left[ {\begin{array}{cc} {\cos ( - {{135}^ \circ })} & { - \sin ( - {{135}^ \circ })} \\ {\sin ( - {{135}^ \circ })} & {\cos ( - {{135}^ \circ })} \\ \end{array}} \right]\left[ {\begin{array}{c} 1 \\ 0 \\ \end{array}} \right] = \left[ {\begin{array}{cc} { - 1/2} & {1/2} \\ { - 1/2} & { - 1/2} \\ \end{array}} \right]\left[ {\begin{array}{c} 1 \\ 0 \\ \end{array}} \right] = \left[ {\begin{array}{c} { - 1/2} \\ { - 1/2} \\ \end{array}} \right]\)

and

\(\dfrac{1}{{\sqrt 2 }}\left[ {\begin{array}{cc} {\cos ( - {{135}^ \circ })} & { - \sin ( - {{135}^ \circ })} \\ {\sin ( - {{135}^ \circ })} & {\cos ( - {{135}^ \circ })} \\ \end{array}} \right]\left[ {\begin{array}{c} -1 \\ 0 \\ \end{array}} \right] = \left[ {\begin{array}{cc} { - 1/2} & {1/2} \\ { - 1/2} & { - 1/2} \\ \end{array}} \right]\left[ {\begin{array}{c} -1 \\ 0 \\ \end{array}} \right] = \left[ {\begin{array}{c} { 1/2} \\ { 1/2} \\ \end{array}} \right]\)

Click for more details on symmetric fractals or see the book by Field and Golubitsky.

 

References

  1. Michael Field and Martin Golubitsky, "Symmetric Fractals" (chapter 7) in Symmetry in Chaos: A Search for Pattern in Mathematics, Art and Nature, Oxford University Press (Edition 1, 1992/1995) and Siam (Edition 2, 2009) [see Preview at Google Books].