Die Mandelbrotmenge

Geschichte

Popularisierung der Mandelbrotmenge durch Benoît Mandelbrot

Wie geht die Mandelbrotmenge?

  1. Man nimmt sich einen Startwert cCc\in\mathbb{C}.
  2. Danach bastelt man sich die Folge cc, c2+cc^2+c, (c2+c)2+c(c^2+c)^2+c, \dots.
  3. Man macht das für viele Startwerte cc.

Also hat man im Grunde genommen eine Iteration

zk+1=zk2+cz_{k+1} = z_k^2+c

mit z0=cz_0=c.

Wir schreiben präziser zc,kz_{c,k}.

Wir rechnen ein paar Orbits von Hand:

Also gibt es unterschiedliche Arten von Folgen (konstant, divergent, alternierend)

Theorem 1

Ist zc,k>2|z_{c,k}| > 2 für ein kk, dann ist zc,kk\langle z_{c,k}\rangle_k unbeschränkt.

Proof

Graphischer Beweis mit Fallunterscheidung c>2|c| > 2 und c2|c| \leq 2. Dann guckt man sich den Absolutbetrag als Abstand an und sieht, dass vektoriell die Folge nicht in den r=2r=2-Kreis zurück hinein kommen kann. Für den zweiten Fall gilt dann nach Voraussetzung für ein kk zc,k>2|z_{c,k}| > 2 voraus fast unmittelbar die Behauptung analog zum ersten Fall folgt.

Definition 1: Mandelbrotmenge

Man definiert die Mandelbrotmenge als

M:={cCzc,kk ist beschra¨nkt}\mathbb{M} := \{c\in\mathbb{C} \mid \langle z_{c,k}\rangle_k\text{ ist beschränkt}\}

Man überlege grob den Bereich möglicher Punkte (der 5 bereits bekannten). Es ist aber nicht so offensichtlich, welche Punkte innerhalb des Kreises mit Radius 22 zur Mandelbrotmenge gehören und welche nicht. Ein Computer kann das für relativ viele Startwerte relativ gut abschätzen.

Wir betrachten ein ganz grobes Mandelbrotraster, so ähnlich wie es damals von Brooks realisiert wurde.

Man kann jetzt also die Feinheiten farblich abstufen, je nachdem wie lange die Punkte gebraucht haben, um aus dem Kreis mit Radius 2 zu entfliehen. Also je mehr Iterationen ich wähle, umso detaillierter kann ich die Farbabstufungen machen.

Mandelbrot Plots

Jetzt wollen wir die Mandelbrotmenge plotten...

Exercise 1

Man guckt nach, an welcher Stelle hier jeweils "reingezoomt" wird.

Die Mandelbrotmenge wird im Volksmund auch "Apfelmännchen" genannt. Sieht halt irgendwie nach Apfel aus...

Wie lässt man solche Bilder berechnen?

Grundsätzlich guckt man, ob für einen Wert cCc\in\mathbb{C} der Orbit von f(z)=z2+cf(z) = z^2+c mit Startwert z=0z=0 divergiert. Falls nein, dann wird der Startpunkt schwarz eingefärbt. Falls aber ja, dann wird der Punkt mit einer Farbe eingefärbt, die vom Index des Iterationsschrittes abhängt, an dem die Folge als divergent erklärt wird. Stellt sich also die Frage, wann dass eine Folge als divergent eingestuft werden darf?

Wer oben das mandelbrotProgramm überflogen hat, dem ist aufgefallen, dass als Abbruchbedingung für die Iteration eines Startwerts die Bedingung zk>2|z_k| > 2 verwendet wurde. Man kann nämlich zeigen, dass, falls für einen Wert zkz_k in der Itartionsfolge zk>2|z_k| > 2 gilt, der Orbit sicher divergiert. Diesen Fact macht man sich in der Programmierung zu Nutze.

Theorem 2: Nicht-Mandelbrot-Kriterium

Ist für ein kNk\in\mathbb{N} der Betrag zk>2|z_k| > 2, so divergiert die Folge zk\langle z_k\rangle definiert durch zk+1=zk2+cz_{k+1} = z_k^2+c mit fixem cCc\in\mathbb{C} und Startwert z0=0z_0 = 0.

Proof

Für f(z)=z2+cf(z)=z^2+c mit z0=0z_0=0 haben wir z0=0,z1=c,z2=c2+c,z_0=0, z_1=c, z_2=c^2+c,\dots Ferner ist f(z)z=z2+czz2cz=zcz>z1>1\frac{|f(z)|}{|z|}=\frac{|z^2+c|}{|z|}\geq\frac{|z|^2-|c|}{|z|}=|z|-\frac{|c|}{|z|}>|z|-1>1 wobei die Dreiecksungleichung verwendet wurde. Das heisst wir haben ein Wachstum pro Zeitschritt mit Faktor grösser 11.

Ist nun c<2|c|<2 und für ein kk zk>2|z_k|>2, dann erfüllt z=zkz=z_k die Voraussetzung und die Folge wird divergieren. Ist c>2|c|>2, dann ist c2+cc2c=cc1>c>2|c^2+c|\geq |c|^2-|c|=|c|\cdot|c-1|>|c|>2. Also erfüllt z=z2z=z_2 die Voraussetzung.

Der kommentierte Beweis zur Divergenz im Apfelmännchen von gym math gibts unter eben dem Link.

Das Apfelmännchen ist mit den Walen verwandt!

Betrachte die Iteration zk+1=zk2+cz_{k+1} = z_{k}^2+c und setze eine affine Transformation zk=axk+bz_k = ax_k+b mit a,b,xkRa,b,x_k\in\mathbb{R}. Es ist

axk+1+b=(axk+b)2+c=a2xk2+2abxk+b2+c.ax_{k+1}+b = (ax_k+b)^2+c = a^2x_k^2+2abx_k+b^2+c.

Es folgt

xk+1=2bxk(1+a2bxk)+b2b+ca.x_{k+1} = 2bx_k(1+\frac{a}{2b}x_k)+\frac{b^2-b+c}{a}.
Exercise 2

Zeige obige Folgerung.

Solution

Falls es nicht klappt, dann kannst du dir den Zusammenhang zwischen Mandelbrot und Wale anschauen.

Exercise 3: Parameter-Vergleich

Begründe, dass ein Vergleich mit der logistischen Funktion fr(x)f_r(x) heisst: b=r2b=\frac{r}{2} und a=ra=-r. Ferner zeige man, dass der interessante Bereich 1<r<41<r<4 den Werten 2<c<14-2<c<\frac{1}{4} entspricht.

Exercise 4: Bestimme die c-Bereiche der Fixpunktregionen 1<r<3 und 3<r<1+\sqrt{6}.

Mandelbrotbereiche für Fixpunkte der Ordnung 1 und 2 von frf_r

Abschliessend begreife man folgendes Bild !

Enjoy! Don't Freak out!

Exercise 5

Jetzt kommt eine gemütliche Pflichtaufgabe. Höre zum Video Mandelbrot Zoom, in dem eine Stelle in der komplexen Ebene vergrössert wird, den sinnigen Titel "Functionality" dazu. Mindestens bis Minute 7.5 den Zoom angucken. Dem Video voraus schicken möchte ich, dass die Dimension des Unversiums im Verhältnis zu einem Atomkern in etwa 1040÷110^{40}\div1 ist. In Längen gedacht heisst dies also für Zahlen, welche sich erst an der 4040-sten Stelle nach dem Komma unterscheiden, dass das Universum um einen Atomdurchmesser verschoben wurde. Erinnere dich an diesen Fact, wenn jeweils die Grössenordnung im Video eingeblendet wird.

So, los gehts! Zuerst den Sound einrichten...

  1. Functionality 2 von DJ Dimsa

  2. Mandelbrot Zoom