Iterated function system
From Wikipedia, the free encyclopedia
In mathematics, iterated function systems or IFSs are a method of constructing fractals; the resulting constructions are always self-similar.
IFS fractals as they are normally called can be of any number of dimensions, but are commonly computed and drawn in 2D. The fractal is made up of the union of several copies of itself, each copy being transformed by a function (hence "function system"). The canonical example is the Sierpinski gasket. The functions are normally contractive which means they bring points closer together and make shapes smaller. Hence the shape of an IFS fractal is made up of several possibly-overlapping smaller copies of itself, each of which is also made up of copies of itself, ad infinitum. This is the source of its self-similar fractal nature.
Contents |
[edit] Definition
Formally,
where
and
are the functions to be iterated. S is the fixed point of a Hutchinson operator, which is the union of the functions fi.
[edit] Properties
The collection of functions fi together with composition form a monoid. If there are only two such functions, the monoid is the dyadic monoid. Composition can be visualized as an infinite binary tree, where, at each node of the tree, one may compose with the one or the other function (i.e take the left or the right branch). In general, if there are p functions, then one may visualize the composition as a p-adic tree.
Because function composition takes this form, the elements of the monoid can be seen to be isomorphic with the p-adic numbers; that is, each digit of the p-adic number indicates which function is to be composed with.
The automorphism group of the dyadic monoid is the modular group; this can be seen to describe the fractal self-similarity of many fractals, including the Cantor set and the de Rham curves.
[edit] Examples
Sometimes each function fi is required to be linear, or more accurately an affine transformation and hence can be represented by a matrix. However, IFSs may also be built from non-linear functions, including projective transformations and Möbius transformations. The Fractal flame is an example of an IFS with nonlinear functions
The most common algorithm to compute IFS fractals is called the chaos game. It consists of picking a random point in the plane, then iteratively applying one of the functions chosen at random from the function system and drawing the point. An alternative algorithm is to generate each possible sequence of functions up to a given maximum length, and then to plot the results of applying each of these sequences of functions to an initial point or shape.
Each of these algorithms provides a global construction which generates points distributed across the whole fractal. If a small area of the fractal is being drawn, many of these points will fall outside of the screen boundaries. This makes zooming into an IFS construction normally impractical. Where a high degree of detail is required in a small area of the fractal, local construction methods based on calculating forward orbits and the fate of individual points may be more efficient (although no software for solving IFS does this).
[edit] Example: a fractal "fern"
Here is an example of a fern-like image computed using an iterated function system.
The first point drawn is at the origin (x0 = 0, y0 = 0) and then the new points are iteratively computed by randomly applying one of the following four coordinate transformations:
- xn + 1 = 0
- yn + 1 = 0.16 yn.
This coordinate transformation is chosen 1% of the time and maps any point to a point in the line segment shown in green in the figure.
- xn + 1 = 0.2 xn − 0.26 yn
- yn + 1 = 0.23 xn + 0.22 yn + 1.6.
This coordinate transformation is chosen 7% of the time and maps any point inside the black rectangle to a point inside the red rectangle in the figure.
- xn + 1 = −0.15 xn + 0.28 yn
- yn + 1 = 0.26 xn + 0.24 yn + 0.44.
This coordinate transformation is chosen 7% of the time and maps any point inside the black rectangle to a point inside the dark blue rectangle in the figure.
- xn + 1 = 0.85 xn + 0.04 yn
- yn + 1 = −0.04 xn + 0.85 yn + 1.6.
This coordinate transformation is chosen 85% of the time and maps any point inside the black rectangle to a point inside the light blue rectangle in the figure.
The first coordinate transformation draws the stem. The second draws the bottom frond on the left. The third draws the bottom frond on the right. The fourth generates successive copies of the stem and bottom fronds to make the complete fern. The recursive nature of the IFS guarantees that the whole is a larger replica of each frond. Note: The fern is within the range -5 <= x <= 5 and 0 <= y <= 10.
[edit] History
Iterated function systems are a generalisation of the Cantor set, first described in 1884, and of de Rham curves, a type of self-similar curve described by Georges de Rham in 1957. IFS were conceived in their present form by John Hutchinson in 1981 and popularized by Michael Barnsley's book Fractals Everywhere. In 1992 Scott Draves developed the Fractal flame algorithm.