next up previous contents
Next: Dynamical Systems and Fractal Up: Tilings as Dynamical Systems Previous: Self-similar tilings arising from

Complex numbers and similarity constants

 

The sphinx pattern (see this figure) has a twist to its similarity. For the type X sphinx, there is exactly one smaller type X sphinx inside, rotated at tex2html_wrap_inline5566 with respect to the larger one. The similarity that carries the smaller one to the larger one has this rotation as well as a scale factor of 2. We can express similarities more easily with complex numbers

displaymath6736

corresponding to the point tex2html_wrap_inline6738 . Here i stands for a choice of tex2html_wrap_inline6742 , i.e. a number satisfying tex2html_wrap_inline6744 . We use tex2html_wrap_inline5806 to denote the set of real numbers, and tex2html_wrap_inline6748 to denote the set of complex numbers. Multiplication by complex numbers captures both scaling (multiplication by positive real numbers) and rotation. For a first example, consider

displaymath6750

The point (-y,x) is the rotation of (x,y) by tex2html_wrap_inline5502 counterclockwise.

Rather than include pictures of this sort of rotation here, we offer a MAPLE worksheet instead that will allow the readers to try out these sorts of calculations and plots themselves.

Complex Rotations and Spirals worksheet
Once you download this file and initiate MAPLE, you may load the worksheet by selecting File-Open in the MAPLE menus. Then read the text included in the worksheet. You may modify the calculations in the worksheet at will. A good introduction to MAPLE in the framework of the physical sciences is [Bay94]. This book is accompanied by a fine set of MAPLE worksheets (in a library called tmlib) for each chapter. The worksheet for Chapter 8 covers complex number calculations.

More generally, every nonzero complex number can be written in polar coordinates

displaymath6758

where r>0 is the length tex2html_wrap_inline6762 of z and tex2html_wrap_inline6058 is the angle that the line from 0 to z makes with the positive x-axis. Euler discovered a remarkable formula relating the trigonometric functions and the exponential. One way to approach Euler's formula is to study the Taylor series expansions of the three functions:

align1203

The series for tex2html_wrap_inline6772 and tex2html_wrap_inline6774 resemble the odd and even terms, respectively, of tex2html_wrap_inline6776 . To account for the signs, we consider the powers tex2html_wrap_inline6778 . These powers repeat every fourth term in the pattern

equation1259

Then

align1267

The beautiful aspect of this extension of tex2html_wrap_inline6776 to complex numbers is that the usual properties of tex2html_wrap_inline6776 :

align1324

can be viewed as consequences of just the algebra of the Taylor series for tex2html_wrap_inline6776 . Therefore, these properties also hold for complex values of z and w. This allows us to see the geometric effect of multiplication by complex numbers. Choose two complex numbers tex2html_wrap_inline6798 and tex2html_wrap_inline6800 in polar coordinates. Then

displaymath6802

We are scaling the length tex2html_wrap_inline6804 by the factor r, and we are adding angle tex2html_wrap_inline6058 to tex2html_wrap_inline6810 , that is, we rotate w by an angle tex2html_wrap_inline6058 counterclockwise.

For our sphinx, we must rotate by tex2html_wrap_inline6816 (or tex2html_wrap_inline5566 clockwise), which corresponds to

displaymath6820

Therefore, the complex scaling ratio between the small type X tile and the larger type X tile is tex2html_wrap_inline6826 .

However, simply multiplying by a will not cause the smaller tile to precisely overlap the larger tile. We may have to translate the result. The added bonus of using complex numbers is that translations are also easy to express. Adding fixed amounts to the real and imaginary parts amounts to adding a complex number. That is, a translation of the complex plane takes the form:

displaymath6830

where b is a complex number. Combining multiplications and additions, the general similarity of the complex plane is

displaymath6834

where tex2html_wrap_inline6836 and b are both complex numbers.

We can find the complex similarity of the sphinx if we introduce complex coordinates for the corners of the sphinx. We choose the origin 0 to be the lower left corner, and we choose the bottom side to run along the positive real axis. Finally, suppose that the length of a short side of the sphinx is exactly 1. Let's walk around the edge of the sphinx in the clockwise direction starting at 0. The first vertex we come to is at angle tex2html_wrap_inline5564 or tex2html_wrap_inline6842 radians exactly 1 unit away from 0. This point is

displaymath6844

The remaining vertices are obtained by adding 1, tex2html_wrap_inline5070 , or the conjugate tex2html_wrap_inline6848 , depending on the direction chosen. Thus, the vertices in clockwise order are

displaymath6850

The smaller X tile has vertices at tex2html_wrap_inline6854 and tex2html_wrap_inline6856 as well. The sphinx with the vertices so labelled is shown in this figure. The complex similarity that maps the smaller sphinx to the larger one must satisfy:

align1350

Subtracting the two equations, we have

displaymath6858

Dividing by tex2html_wrap_inline5070 , we find

displaymath6862

as we already knew to be the case. We can also determine

align1357

   figure1361
Figure 34: Labelling the vertices of the sphinx by complex coordinates

We may now exploit this similarity and the deflation of the type X and Y sphinxes to obtain a tiling of the whole plane. The first few generations of this method are shown in this figure. We see a spiral pattern emerging. In the last frame, we draw a spiral passing through the images of the vertex at 3 under iteration of the complex similarity T.

   figure1372
Figure 35: Generations of the sphinx tiling under the complex similarity T.

The spiral converges onto a single point which is the fixed point of the transformation:

displaymath6872

Solving, we obtain the solution

displaymath6874

This is approximately 1.86 + .988 i. This kind of behavior is common for complex similarities; if tex2html_wrap_inline6878 , the sequence of iterates tex2html_wrap_inline6880 , tex2html_wrap_inline6882 , tex2html_wrap_inline6884 etc., generally lies on an infinite spiral curve, known as a ``loxodrome.''

If we change coordinates so that the origin becomes the fixed point tex2html_wrap_inline6886 , the transformation T becomes simply a complex multiplication. We do this by defining

displaymath6890

Then

displaymath6892

Iterates of multiplications are easy to write down. We denote the composition of U with itself n times by tex2html_wrap_inline6898 . Then tex2html_wrap_inline6900 . The spirals that result are parametrized by tex2html_wrap_inline6902 , tex2html_wrap_inline6904 . We can calculate the powers of a complex number by using polar coordinates tex2html_wrap_inline6906 . Then the spiral curve has the form tex2html_wrap_inline6908 for some constants c,r>0 and real angles tex2html_wrap_inline6912 . These are also called exponential spirals.

Our figures of the sphinx tiling in this section were generated with the aid of a MAPLE worksheet. We encourage the reader to examine and experiment with this worksheet which is located at

MAPLE worksheet sphinx.ms

To review, similarities of the plane are the same as complex linear maps T(z)=az+b. Any structure (tiling or fractal) with self-similarities that involve both scaling and rotations generally exhibits spiralling behavior.


next up previous contents
Next: Dynamical Systems and Fractal Up: Tilings as Dynamical Systems Previous: Self-similar tilings arising from

David J. Wright
Mon Aug 19 17:21:15 CDT 1996