next up previous contents
Next: Contraction mappings and hyperbolic Up: Iterated Function Systems Previous: Contractions and Fixed Points

Linear Contractions

For iterated function systems, the mappings are affine linear functions. Therefore, it is important to determine when these are contractions and what there contractivity factors are. Choose an affine linear map:

displaymath8076

The first thing to observe is that if we compose w with an isometry T, the resulting map tex2html_wrap_inline8082 changes distances in the same way:

displaymath8084

Isometries consist of translations and reflections and rotations. We have already seen how to choose a translation so that

displaymath8086

So we shall assume that w already satisfies this equation. That comes down to e=f=0. That is,

displaymath8092

The easiest case is when b=c=0. Then w(x,y)=(ax,dy) corresponds to a diagonal matrix. We see that w(1,0)=a(1,0) and w(0,1)=d(0,1). The map w has a separate scaling factor along the x-axis and along the y-axis. Then

align2062

So the contractivity factor is tex2html_wrap_inline8108 if both diagonal coefficients have absolute value less than 1.

In linear algebra courses, we spend a great deal of time learning how to change coordinates so that general linear transformations become diagonal ones. This is the theory of eigenvalues (the eventual entries on the diagonal) and eigenvectors. This is related to the problem of determining whether or not a linear transformation is a contraction; however, we will approach the problem from the point of view of finding the maximum scaling factor. Writing points as column vectors, suppose we have tex2html_wrap_inline8110 and tex2html_wrap_inline8112 . Suppose the matrix of our transformation w is tex2html_wrap_inline8116 . Then

displaymath8118

using matrix multiplication. The two distances we must compare are

align2073

using the usual notation for transpose of a matrix, and

align2079

remembering that transpose switches the order of matrix multiplication. We want the maximum of

displaymath8120

tex2html_wrap_inline8122 is a symmetric matrix with the property that tex2html_wrap_inline8124 for all vectors v. Let's suppose that tex2html_wrap_inline8128 and that tex2html_wrap_inline8130 . Then we want the maximum of

displaymath8132

subject to tex2html_wrap_inline8130 . Multiplying out, the function we must maximize is

displaymath8136

This problem is best handled by the theory of Lagrange multipliers, which says that at any extreme point of F subject to tex2html_wrap_inline8130 we must have

align2088

This means

align2094

In matrix form, this may be written

displaymath8142

This is exactly what it means for tex2html_wrap_inline8144 to be an eigenvector of tex2html_wrap_inline8122 , and the number tex2html_wrap_inline6624 is the associated eigenvalue.

The scaling factor at the solution tex2html_wrap_inline8144 is actually tex2html_wrap_inline8152 . So we can determine whether or not w is a contraction by computing the eigenvalues of tex2html_wrap_inline8122 . We end this section with an example. Consider

displaymath8158

(this actually arises in the fern IFS). Then

displaymath8160

The eigenvalue equation is

displaymath8162

The roots are

displaymath8164

The square root of the larger value is the contractivity factor; this is roughly 0.34.

In summary, an affine map with tex2html_wrap_inline6172 coefficient matrix A is a contraction provided both the eigenvalues of tex2html_wrap_inline8122 are strictly less than 1. The best choice of contractivity factor is the square root of the larger eigenvalue.


next up previous contents
Next: Contraction mappings and hyperbolic Up: Iterated Function Systems Previous: Contractions and Fixed Points

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