Iteriertes Funktionensystem

Ein iteriertes Funktionensystem (IFS) dient meist der Konstruktion von Fraktalen, diese werden dementsprechend auch als IFS–Fraktale bezeichnet. Bekannte Vertreter dieser Klasse von Fraktalen sind das Sierpinski-Dreieck und die Koch-Kurve wie auch die Grenzmengen von Lindenmayer-Systeme.

Sierpinski_dreieck.png

Kochkurve.png

Menger_5.PNG]]

Drachen_p.png]]

Diese Art der Fraktalkonstruktion wurde 1981 von Hutchinson erfunden und später von Barnsley mit seinem Buch Fractals Everywhere popularisiert. Dort gab Barnsley auch den Collage-Satz an, welcher die Grundlage der fraktalen Bildkompression bildet. Diese Art, Bilder effizient mittels Datenstrukturen zu kodieren, hat sich jedoch nie richtig durchsetzen können und wird heute im wesentlichen nur noch als Hybridverfahren in Kombination mit einer Wavelet-Transformation untersucht.

Der Begriff iteriertes Funktionensystem bezeichnet eine Menge $\mathcal F$ von Funktionen, die denselben Raum M als Definitions- und Wertebereich haben und unter Verknüpfung abgeschlossen sind. Also

\mathcal F\circ\mathcal F\subset\mathcal F\;

d.h.

\; \forall f,g\in\mathcal F:\; f\circ g\in\mathcal F.

Um für eine solche Funktionenmenge Eigenschaften ableiten zu können, muss sie weitere Voraussetzungen erfüllen. Üblicherweise, wenn von IFS gesprochen wird, werden diese Voraussetzungen stillschweigend als gegeben angenommen. Diese Voraussetzungen sind

dass das IFS endlich erzeugt ist, also endlich viele Funktionen enthält, aus welchen die anderen durch wiederholte (iterierte) Verknüpfung zusammengesetzt werden können,
dass der Raum M ein vollständiger metrischer Raumes mit Metrik d ist, und
dass jede Funktion des IFS kontraktiv ist, mit 1.) reicht es, dies von den Erzeugenden zu verlangen.

Unter diesen Umständen gibt es eine invariante, selbstähnliche Menge X⊆M.

Die Teilmenge X ist invariant, wenn sie von jeder Funktion des IFS wieder in sich abgebildet wird.
Die Teilmenge X ist selbstähnlich, wenn jeder Punkt aus X in der Bildmenge F(X) einer Funktion $F\in\mathcal{F}$ liegt.

Selbstähnliche Mengen werden auch als Fraktal bezeichnet, deshalb die Bezeichnung IFS–Fraktal. Man könnte auch weitergehend den Begriff der Selbstähnlichkeit durch die Forderung der Existenz eines IFS definieren.

Existenz und Eindeutigkeit der invarianten Menge

Mathematisch gesehen handelt es sich bei der Theorie der iterierten Funktionensysteme, wie auch die Begrifflichkeit vermuten läßt, um eine direkte Anwendung des Banachschen Fixpunktsatzes, wobei mehrere Funktionen statt einer betrachtet werden und, statt eines eindeutigen Fixpunktes, sich eine invariante, meist fraktale, Teilmenge des Raumes M ergibt. Zur Illustration wird meist das zweidimensionale Einheitsquadrat M=^*x^* mit dem euklidischen Abstand gewählt.

Wir beginnen also mit einer endlichen Menge von Funktionen eines kompakten metrischen Raumes (M,d) in sich selbst:

\mathcal F_1:=\{\phi_1,\dots,\phi_r:M\to M\},

von denen wir voraussetzen, dass es eine Kontraktionskonstante 0 gibt mit
$\forall x,y\in M,\phi\in\mathcal F_1:\;d(\phi(x),\phi(y))\le c\,d(x,y)$
Durch Iteration setzen wir $\mathcal F_1$ zu einem IFS $\mathcal F$ fort, es sei
$\mathcal F_{n+1}:=\mathcal F_1\circ \mathcal F_n := \{ \phi\circ F:\;\phi\in\mathcal F_1,\;F\in\mathcal F_n\}$
und erhalten schließlich
$\mathcal F:=\bigcup_{n=1}^\infty \mathcal F_n$ .
Satz: Sind alle Funktionen $\phi_1,\dots,\phi_r$ in $\mathcal F_1$ kontraktiv, so gibt es eine invariante Teilmenge X, welche die Fixpunktgleichung
$X=\bigcup_{i=1}^r \phi_i(X)$
erfüllt. Für diese gilt:
Zu jedem $F\in\mathcal F$ gibt es genau einen Fixpunkt. Die invariante Menge X ist der topologische Abschluss der Menge aller Fixpunkte
$\{x\in M|\;\exists F\in\mathcal F:\,x=F(x)\}$ .
Ist y∈ M ein beliebiger Punkt, so gilt für den Abstand dieses Punktes
$d(X,F(y))\le c_F\,d(X,y)$ für jedes $F\in\mathcal F$ .
Es gilt die Abschätzung $c_F\le c^m$ , falls F eine m-fache Verkettung $F\in\mathcal F_m$ der Ausgangsfunktionen ist.
Damit kann X durch Iteration einer beschränkten Ausgangsmenge
$X_0\subset M,\quad X_{n+1}=\bigcup_{i=1}^r \phi_i(X_n)$
beliebig gut angenähert werden.

Approximation der Grenzmenge

Chaosspiel

Die Gestalt der fraktalen Menge X kann durch ein so genanntes Chaosspiel visualisiert werden. Dabei wird zunächst ein Fixpunkt $\vec x$ von $\vec x=\phi_1(\vec x)$ aufgesucht und auf diesen in zufälliger Reihenfolge die definierenden Funktionen angewandt. Als Algorithmus kann dies wie folgt aussehen:
Weise 100 mal hintereinander $\vec x:=\phi_1(\vec x)$ zu
Wiederhole beliebig oft
Wähle zufällig ein i in {1,...,r}
Weise $\vec x:=\phi_i(\vec x)$ zu
Zeichne den Punkt x
Anmerkung: 1) Es ist in den ersten, blinden, Iterationen unwesentlich, welche Funktion gewählt wird, da in jedem Schritt der Abstand zur fraktalen Menge X reduziert wird. Ist z.B. die Kontraktionskonstante c=0.5 und die Grundmenge M das Einheitsquadrat, welches mit 1024x1024 Pixeln dargestellt wird, so ist bereits nach 12 blinden Iterationen der Fehler unter die Pixelgröße gesunken.
2) Es werden im Allgemeinen bessere Darstellungen erzielt, wenn die Wahrscheinlichkeit des Aufrufs jeder der Funktionen $\phi_i$ in etwa proportional zum Volumen von $\phi_i(M)$ ist.
Rekursion

Eine weitere Möglichkeit der Darstellung, vorzugsweise für die affinen Fraktale, ist die rekursive Approximation der Menge X. Dies wird meist anschaulich mittels eines Fotokopierers erklärt: Man macht verschiedene Verkleinerungen eines Ausgangsbildes, fixiert diese nach Vorschrift auf einem neuen Blatt und benutzt dieses dann als Ausgangsbild des nächsten Schrittes.
Auch die Turtle-Grafik, die zur Konstruktion der L-Systeme verwendet wird, folgt einer ähnlichen Idee.
Als Algorithmus braucht man dazu eine rekursiv aufrufbare Funktion, welche die Zuordnung $F(n)=\bigcup_{i=1}^r \phi_i(F(n-1))$ bei einer beliebigen Menge $F(0)$ realisiert. Die Implementierung benötigt einen Stackspeicher, in welchem das jeweils aktuelle Koordinatensystem als affine Koordinatentransformationen festgehalten wird. Damit ergibt sich als Algorithmus
Koch-rekursiv.png Figur(n):
Falls n=0
zeichne die Basisfigur (z.B. eine Strecke, einen Buchstaben, ein schwarzes Rechteck)
sonst:
Für i:=1 bis r
Lege aktuelles Koordinatensystem auf Stack ab
Transformiere das aktuelle Koordinatensystem entsprechend $\phi_i$
Rufe Figur(n-1) auf
Stelle Koordinatensystem vom Stack wieder her
Fraktal:
Rufe Figur(10) auf (10 als Beispiel)
Beispiele für iterierte Funktionensysteme

Affine Funktionen

Die erzeugenden Funktionen des IFS seien affin-lineare Funktionen des zweidimensionalen Einheitsquadrates in sich selbst. Jede Funktion φ_k ist gegeben durch eine 2×2–Matrix A_k und einen Verschiebungsvektor b_k.
Das Koch-Fraktal wird z.B. von folgendem System von 2 Funktionen erzeugt:
$\phi_1\binom{x}{y}= \frac1{\sqrt3}\begin{pmatrix}\cos30^o&\sin30^o\\ \sin30^o&-\cos30^o\end{pmatrix}\, \binom{x}{y}$ ,
$\phi_2\binom{x}{y}= -\frac1{\sqrt3}\begin{pmatrix}-\cos30^o&\sin30^o\\ \sin30^o&\cos30^o\end{pmatrix} \binom{x}{y}+\frac1{\sqrt3}\binom{\cos30^o}{\sin30^o}$
Die klassische Methode zur Erzeugung der Koch-Kurve benutzt 4 Funktionen
$\phi_1(x,y)=\left(\frac{x}3,\frac{y}3 \right)$
$\phi_2(x,y)=\left(\frac{2+x-\sqrt3 y}6,\frac{\sqrt3 x+y}6\right)$
$\phi_3(x,y)=\left(\frac{3+x+\sqrt3 y}6,\frac{\sqrt3-\sqrt3 x+y}6\right)$
$\phi_4(x,y)=\left(\frac{2+x}3 ,\frac{y}3\right)$
Das rechtwinklige Sierpinski-Dreieck wird erzeugt von
$\phi_1(x,y)=\left(\frac{x}2,\frac{y}2\right)$
$\phi_2(x,y)=\left(\frac{x+1}2,\frac{y}2\right)$
$\phi_3(x,y)=\left(\frac{x}{2},\frac{y+1}{2}\right)$
Collagen

Grundlage für die Begeisterung für solche IFS–Fraktale war das Collage–Theorem von Barnsley. Es besagt, dass jede kompakte Menge – jede Gestalt – durch ein IFS–Fraktal beliebig genau angenähert werden kann. Die Grundlage dafür sind folgende Beobachtungen:
Jede endliche Menge ist ein IFS–Fraktal. Die zugehörigen Funktionen sind diejenigen konstanten Funktionen, welche den gesamten Raum auf jeweils einen der endlich vielen Punkte abbilden.
Jede kompakte Menge K hat für jedes ε>0 ein ε-Netz, wird also durch endlich viele Kugeln vom Radius ε überdeckt.
⇒ Das IFS–Fraktal der Kugelmittelpunkte enthält die Ausgangsmenge in einer ε–Umgebung

Anschaulicher: Haben wir in einem 100×100–Pixelbild eine Figur von 500 schwarzen Pixelpunkten, so können wir das Bild um den Faktor 100 auf die Größe eines Pixels verkleinern und mit diesem einen schwarzen Punkt dann wieder die Figur malen, indem wir ihn auf jeden der zugehörigen 500 Pixel abbilden. Diese Vorgehensweise ist bei weitem nicht optimal, hier wäre das einfache Speichern der Positionen der 500 Pixel einfacher. Aber wenn wir für den gleichen Zweck mit nur fünf Abbildungen auskämen, wäre eine Datenreduktion erzielt.
Wir sind auch nicht auf einfache Schwarzweißbilder eingeschränkt. Bei einem Graustufenbild kann der Grad der Schwärzung als dritte Koordinate des Punktes aufgefasst werden, es ergibt sich eine kompakte Fläche im dreidimensionalen Raum, auf welche wieder das Collage–Theorem angewendet werden kann. Mit systematischen Verfahren zur Konstruktion eines IFS–Fraktals mit möglichst wenigen Funktionen befaßt sich die Fraktale Bildkompression.
Weblinks

Eine Sammlung fraktaler Kunst von Paul Bourke
Farbige Erweiterungen des Chaosspiels fraktale Flammen, unter math die Theorie dazu.
Zu Geschichte und Verallgemeinerungen: Cabrelli, Molter: Generalized Self-Similarity ...
1957 Selbstähnliche Funktionen von Bajraktarevic und de Rham untersucht
1981 Selbstähnliche Mengen und Kurven von Hutchinson untersucht
1986 definiert Barnslay fraktale Funktionen
Fraktale Geometrie | Topologie | Computerkunst
Iterated function system | Iteratitaj sistemoj de funkcioj | Système de fonctions itérées | Itererande funktionssystem

Gourt Home::Gourt :: Iteriertes Funktionensystem

Existenz und Eindeutigkeit der invarianten Menge

Approximation der Grenzmenge

Chaosspiel

Rekursion

Beispiele für iterierte Funktionensysteme

Affine Funktionen

Collagen

Weblinks