
Apfelmännchen algebraisch
Von: shadowking
Datum: Fr. 16. Oktober 2015 19:43:06 Thema: Mathematik
|
Wie die Mandelbrotmenge aussieht, brauche ich jemandem, der auf diesen Seiten regelmäßig aktiv ist, nicht mehr zu erklären. Diese äußerst komplizierte fraktale Menge ist zu einem Aushängeschild für die moderne rechnergestützte Mathematik und Algorithmik geworden. Ob ein wissenschaftlicher Themenbereich "populär" wird, ist weniger eine Frage seines Inhalts oder seiner Bedeutung, sondern eine Frage der Reklame. Und dafür hatten Herr Mandelbrot und seine Mitstreiter definitiv die schöneren Bilder produziert - und das in einer Zeit, da Rechenmaschinen noch groß wie Schränke waren und weniger konnten als etwa heute ein Smartphone. Da die numerische Rechenpower heute leicht zu bekommen und kostengünstig ist, siegt allzu oft die schiere Rechnerkraft über eine Betrachtungsweise, die den Phänomenen mit klassischer Mathematik auf den Grund geht.
Inhalt:
1. Einleitung
Ich verfolge mit diesem Artikel die Absicht, mich diesem Objekt einmal nicht von der numerischen, sondern der algebraischen Sichtweise her zu nähern, womit gemeint ist, die Bestimmung von Punkten und Kurven soweit zu treiben, wie dies mit algebraischer Exaktheit (und nicht nur numerischer Präzision) möglich ist. Die ersten Mathematiker, die sich dem Thema "quadratischer Iteration" näherten, waren die Franzosen Gaston Julia und Pierre Fatou um 1910, und die hatten schließlich auch noch keinen Rechner, um sich die fraktalen Gebilde veranschaulichen zu können, mit denen sie sich beschäftigten. Sehr weit wird man mit der reinen Kraft der Algebra in diesem Gebiet allerdings dennoch nicht kommen, denn die zugrundeliegende quadratische Dynamik impliziert ein exponentielles Anwachsen der Grade der Polynome, mit denen man es zu tun bekommt und die schon für kleine Zyklenlängen nicht nur nicht mehr praktisch, sondern wegen prinzipieller algebraischer Gründe (man denke an den Satz von Abel-Ruffini) nicht mehr handhabbar sind. Wer die Berechnungen, die in diesem Artikel vorgenommen werden, nachvollziehen möchte, der sollte also über ein leistungsfähiges Computeralgebrasystem verfügen. Ich habe alles mit wxMaxima durchgerechnet, eine freie und dennoch sehr brauchbare Software. Die Plots sind mittels Gnuplot erstellt, und so gute habe ich selten gesehen.
2. Notationen und Definitionen
Gegeben sei ein komplexer Parameter aus sowie die quadratische Abbildung
.
Wir betrachten die Iteration dieser Abbildung für verschiedene Startwerte aus .
Es gibt mehrere Möglichkeiten:
- Es tut sich gar nichts; es handelt sich um einen Fixpunkt.
- Die Iterierten streben einem Grenzwert zu.
- Die Iterierten streben einem Grenzwert nicht nur für zu, sondern erreichen ihn exakt nach einer endlichen Anzahl Schritte.
- Die Iterierten drehen sich im Kreis; verhalten sich zyklisch.
- Die Iterierten werden nach einer endlichen Anzahl Schritte periodisch.
- Die Iterierten nähern sich einem Grenzzyklus asymptotisch an.
- Die Iterierten wachsen betragsmäßig über jede Grenze.
- Keiner der bisherigen Fälle tritt ein; die Iterierten verhalten sich scheinbar regellos (deterministisch chaotisch).
Die Julia-Menge zum komplexen Parameter ist die Menge aller Punkte der komplexen Ebene, für die der vorletzte Fall nicht eintritt, die Iteriertenfolge also nicht betragsmäßig gegen Unendlich divergiert.
Das Verhalten hängt kritisch von dem Parameter ab. Wenn wir die Betrachtung vom Raum der Folgenglieder auf den Raum der Parameter richten, können wir setzen und die Menge derjenigen betrachten, für die die spezielle Folge , also nicht beliebig anwächst. Diese Punktmenge ist als Mandelbrotmenge definiert.
Definition: Die Menge aller Parameter, für die die Iteration für nicht gegen Unendlich divergiert, heißt die Mandelbrotmenge.
Bei genauerer Betrachtung fällt auf, daß unzählige kleine (näherungsweise) Kopien von sich selbst enthält, also eine selbstähnliche fraktale Menge ist. Wir werden dies im folgenden streng beweisen und auch die Gründe für dieses Phänomen näher beleuchten. Zunächst sei hier nur der sprachlichen Bequemlichkeit halber festgelegt, daß diese "kleinen Kopien von " innerhalb von Satelliten genannt werden.
Die psychedelisch-bunten Bilder, die man von überall zu sehen bekommt, ergeben sich dadurch, daß man die Bereiche, für die die Iteriertenfolge divergiert, unterschiedlich einfärbt, je nachdem, wie schnell diese Divergenz erfolgt. Doch in aller algebraischen Strenge muß darauf hingewiesen werden, daß ein Schwarzweißbild die Situation klarer wiedergibt, da ein Parameter nur entweder zur Mandelbrotmenge dazugehört oder nicht.
3. Mandelbrot und Feigenbaum
Beschränkt man sich auf die reellen Zahlen, so ähneln sich die Iterationen und die der logistischen Abbildung so sehr, daß man die Ergebnisse der Untersuchung ihrer Parameterräume durch eine Umrechnung bzw. von der einen auf die andere übertragen kann. Das bekannte Feigenbaumdiagramm mit seinen Periodenverdopplungskaskaden und sich in "chaotischen" Bereichen unvermittelt öffnenden "Fenstern der Ordnung" überträgt sich darum quasi unverändert auf die Mandelbrotmenge. Die Bereiche stabiler 2-, 4-, 8- etc.-Zyklen entspricht der Ansammlung von "Dutts" auf dem "Kopf" von , während der größte Satellit in der "Antenne", die bis reicht, dem größten Fenster der Ordnung im Feigenbaumdiagramm entspricht.
Folgende markante Punkte seien hier mitsamt der exakt bestimmten Parameterwerte gegenübergestellt:
Man kann die Umrechnungsfunktion entweder empirisch ermitteln und mit den einander entsprechenden Parameterwerten einen quadratischen Fit veranstalten, wie ich das zunächst gemacht habe - oder wieder algebraisch vorgehen. Dazu wähle man das dritte Parameterpaar . Das Polynom annulliert , also setze man eine quadratische Ansatzfunktion für ein und vergleiche das Ergebnis koeffizientenweise mit dem Polynom , welches annulliert. Hierbei wird man praktischerweise das eine Polynom mit multiplizieren, dadurch werden die konstanten Glieder gleich und man hat sofort . Dann ergeben sich die anderen beiden Parameter leicht als .
Bild 1: Die Vermittlung zwischen den Parameterwerten des Feigenbaumdiagramms und der reellen Achse der Mandelbrotmenge erfolgt über die Funktion
4. Die Mandelbrotmenge als Bifurkationsdiagramm im Parameterraum der Julia-Mengen
In gewisser Weise ist es also angebracht zu sagen, daß die Mandelbrotmenge die Erweiterung des Feigenbaum-Diagramms ins Komplexe, oder, daß das Feigenbaum-Diagramm der "reelle Schatten" der komplexen Mandelbrotmenge ist. stellt eine Landkarte für Julia-Mengen dar, und Bifurkationslinien, an denen sich das Verhalten der Iterierten qualitativ ändert, treten auch in der Mandelbrotmenge auf. Während aber das Feigenbaum-Diagramm nicht nur den Ort dieser Stellen, sondern auch zugleich die Entwicklung des Attraktors darstellt, sieht man in nur die Örter, an denen sich das Systemverhalten ändert, aber nicht, wie es sich ändert. Wie es sich ändert, verraten die Gestaltänderungen der entsprechenden Julia-Mengen bei Überschreiten der Bifurkation.
So hat der im Reellen einzig mögliche Bifurkationstyp mit Kodimension 1, die Heugabelbifurkation, im Komplexen mannigfache Entsprechungen. Ein bislang attraktiver Fixpunkt oder zyklischer Punkt kann im zweidimensionalen nicht nur in zwei, sondern in beliebig viele zyklische Punkte auseinanderspringen. Und zwar geschieht dies, wenn die Randkurve der hyperbolischen Bereiche an Stellen überschritten wird, die bei einer sinnvoll zu wählenden Parametrisierung rationalen Vielfachen von entsprechen, wie in Abschnitt 9 dieses Artikels näher erläutert werden soll.
Das eröffnet die Frage, was mit einem attraktiven Zyklus geschieht, wenn der Parameter die Hauptkardioide von etwa dort verlassen wird, wo dieser Parameter z.B. beträgt. Da er nicht in zyklisches Verhalten umspringen kann, sondern nur "in die Nähe einer Periodizität", könnte das Systemverhalten "fastperiodisch" oder deterministisch chaotisch werden, oder es könnte sogar langfristig Divergenz eintreten, also schon außerhalb von liegen.
5. Zentren und Misiurewicz-Punkte
5.1 Zentren
Jede hyperbolische Komponente in hat genau einen "Mittelpunkt". Dieser ist dadurch gekennzeichnet, daß die Folge der Iterierten zu nach Schritten wieder exakt in 0 zurückläuft, wobei die der Komponente zugeordnete Zyklenlänge ist. Es kann nur einen derartigen Punkt je Komponente geben, da sich bei Variation des Parameters auch die Lage der Punkte des stabilen Zyklus innerhalb der Julia-Menge etwas verschiebt und so 0 nicht mehr im Zyklus liegt. Als das Zentrum der hyperbolischen Menge wird dieser ausgezeichnete Punkt definiert.
Die Bestimmung dieses Zentrums erfolgt mit der Gleichung . Zudem lassen sich für sämtliche Teiler der Zyklenlänge Faktoren ausdividieren, und ist ohnehin für jede dieser Gleichungen eine Lösung.
Beispiel: Die Zentren der hyperbolischen Mengen, die zu attraktiven Zyklen der Länge 4 führen, ergeben sich aus
Algebraisch können sie nicht näher bestimmt werden. Da wir wissen, daß eine hyperbolische Menge zur Zyklenlänge 4 den ersten Dutt auf der Kreisscheibe bildet und eine andere einen kleinen Satelliten weit links in der Antenne, gehören die beiden reellen Lösungen gerade zu diesen. Die weiteren Lösungen sind komplex und liegen daher paarweise vor: gehört zu einem Knospenpaar auf der Hauptkardioide dort, wo der entsprechende Fixpunkt repulsiv wird und in einen attraktiven 4-Zyklus übergeht. gehört dagegen zu einem Paar weiterer Satelliten, die über ein Verzweigungssystem mit den Knospen zu 3-Zyklen in verbunden sind.
Da in jeder hyperbolischen Menge nur genau ein Zentrum existiert, läßt sich die Anzahl der hyperbolischen Mengen in zur Zyklenlänge bequem mit angeben, wobei die Teiler von durchläuft und die zahlentheoretische Möbiusfunktion bezeichnet.
Bild 2: Zentren aller hyperbolischer Mengen zu den Zyklenlängen 4,5 und 6. Zu den hyperbolischen Mengen und sind deren Randkurven rot, hellgrün bzw. dunkelblau eingezeichnet.
5. 2 Echte Vorfixpunkte und präzyklische Punkte (Misiurewicz-Punkte)
Die Mandelbrotmenge setzt sich nicht nur aus den hyperbolischen Mengen zusammen, sondern sie ist darüberhinaus stark verästelt. Die auffälligste Verästelung ist die "Antenne" auf der reellen Achse, die dem rechten Rand des Feigenbaumdiagramms entspricht und sich nach links bis hinzieht. Innerhalb dieser "Antenne" fällt ein Satellit ins Auge, eine Miniaturausgabe der gesamten Menge . Diese liegt exakt da, wo an der entsprechenden Stelle im Feigenbaumdiagramm das "Dreierfenster" liegt. Wir werden im nächsten Abschnitt dessen Parametrisierung auffinden und dabei auch erkennen, daß er keine exakte Kopie darstellt. Wenn von der "Selbstähnlichkeit" von die Rede ist, so ist dies also nicht so exakt zu verstehen wie bei den ausdrücklich als selbstähnlich konstruierten Fraktalen Drachenkurve, Sierpiński-Dreieck, Kochscher Schneeflockenkurve etc., sondern nur asymptotisch gültig.
Das Kennzeichen der hyperbolischen Mengen ist, daß die dazugehörigen Julia-Mengen "gutartig" und zusammenhängend sind und ein Inneres besitzen. Alle Iteriertenfolgen nähern sich asymptotisch einem Zyklus oder erreichen ihn sogar exakt. Außerhalb der hyperbolischen Bereiche finden wir nur noch Parameter vor, die zu repulsiven Zyklen in den entsprechenden Julia-Mengen führen, d.h. ein Punkt der komplexen Ebene, der nicht exakt auf einem zyklischen Orbit oder auf dessen Vorperiode liegt, kann noch so nahe daran liegen, doch wird er sich diesem Zyklus nicht nähern, sondern sich immer weiter davon entfernen. Er kann sich statt dessen einem anderen Zyklus annähern oder auch völlig azyklisch und scheinbar chaotisch hin- und herspringen. Algebraisch faßbar sind allerdings nur noch die Punkte, die entweder sofort zyklisches Verhalten zeigen oder aber nach einer endlichen Vorperiode in ein solches übergehen.
Definition: Ein Punkt , für den existieren, so daß , aber für alle wird, heißt Misiurewicz-Punkt. Das minimale , so daß dies gilt, heißt Periodenlänge von , das minimale heißt Vorperiodenlänge. Die Iterierten teilen sich auf in strikt vorperiodische Punkte und periodische Punkte. Ist , so besteht die Periode nur aus einem Fixpunkt und die vorperiodischen Punkte heißen auch Vorfixpunkte.
Misiurewicz-Punkte sind damit genau die Parameter, für die 0 ein strikt vorperiodischer Punkt ist.
Rein zyklische Punkte hatten wir als die Zentren hyperbolischer Mengen identifiziert; diese Zyklen sind also immer attraktiv für ihre nahe Umgebung. Ferner ist, da die Iteration mit dem Startwert 0 beginnt, 0 für rein zyklische Punkte nicht präperiodisch. Gestattet man aber die Einschränkung, daß ein Parameter erst nach Durchlaufen einer gewissen Vorperiode in einen Zyklus einmündet, so finden wir abseits der hyperbolischen Bereiche Punkte, die entscheidbar zu gehören.
Zunächst ist festzustellen, daß es Punkte mit der Vorperiode 1 nicht geben kann. Gäbe es einen solchen, so wäre der Beginn der Periode. Die Iteration besteht aber gerade darin, den Ausgabewert zu quadrieren und dann den Parameter zu addieren. Man landet also nach einem Schritt nur dann bei , wenn man im vorangegangenen Schritt bei 0 stand, also ist 0 nicht strikt vorperiodisch. Minimale Vorperiodenlänge ist also .
Zur praktischen Berechnung wird also der Ansatz bzw. gemacht. Die Einschränkungen, daß die Iteration genau nach dem -ten Schritt in die Periode übergeht und nicht bereits vorher, erlaubt das Ausdividieren des Faktors . Ist weiterhin , so kann, da die Periodenlänge minimal definiert ist, zu den Koteilern jedes Primteilers von ein weiterer Faktor ausdividiert werden, und zwar , wobei man darauf zu achten hat, daß man eventuell mehrfach ausdividierte Teiler nach der Methode des Einschluß-Ausschluß-Verfahrens wieder hinzumultipliziert. Insgesamt erreicht man durch das Ausdividieren von Faktoren, daß die zu erwartenden hohen Polynomgrade weitestgehend reduziert werden und man die Nullstellen sicher als Misiurewicz-Punkte des gewünschten Typs identifizieren kann, ohne in Listen nachsehen zu müssen, ob sie evtl. schon bei kleineren Vorperioden- bzw. Periodenlängen aufgetreten sind.
Beispiel: Vorperiode 2, Periode 2
und sind also die einzigen Misiurewicz-Punkte zu .
Ein anderer bemerkenswerter Misiurewicz-Punkt ist die Antennenspitze . Hierbei handelt es sich um einen echten Vorfixpunkt, denn ihre Iteriertenfolge lautet
Die folgende Bilderserie zeigt für einige kleine die dazugehörigen Misiurewicz-Punkte. Vorfixpunkte sind in allen Bildern dunkelblau dargestellt.
Bild 3: Einige Misiurewicz-Punkte mit der Vorperiodenlänge 2. Enthalten sind auch –2, i und –i
Bild 4: Einige Misiurewicz-Punkte mit der Vorperiodenlänge 3
Bild 5: Einige Misiurewicz-Punkte mit Vorperiodenlängen 4, 5 und 6
Soweit ich die berechneten Misiurewicz-Punkte mit Fraktalprogrammen untersucht habe, waren sie alle vom Typ "Spitze" oder "Ende einer verzweigten Struktur". Wikipedia meint, es gäbe sie auch als diskrete Punkte auf der Antenne oder – was noch interessanter aufzufinden wäre – als Verzweigungspunkte in den Dendriten jenseits der Vervielfachungskaskaden der hyperbolischen Mengen. Wenn man letztere nicht abzählen mag, um festzustellen, welche Zyklenlänge zu den Punkten in ihr gehört, so kann man nämlich praktischerweise auch die Äste, die von dem über ihr liegenden Verzweigungspunkt ausstrahlen, abzählen. Der Ast, über den dieser mit dem Hauptkörper von verbunden ist, ist dabei mitzuzählen. Diese Verästelungszahl ist gleich der Zyklenlänge. Die Regel, nach der diese in den hyperbolischen Mengen entlang des Randes von anwächst, kann man sich auf diese Weise nochmals bestätigen.
6. Hyperbolische Mengen. Bestimmung von und
Wir können außer zwischen Divergenz und Nichtdivergenz auch nach dem Vorliegen der anderen unter 2. genannten Möglichkeiten fragen. Das führt auf den Begriff der hyperbolischen Menge. Dies ist ein zusammenhängender Bereich der Mandelbrotmenge mit nichtleerem Inneren, dessen Elemente zu einem ähnlichen Verhalten der Iterierten in der Julia-Menge führen. Wenn wir schauen, bei welchen Parametern nur ein Fixpunkt vorliegt, den die Iterierten entweder irgendwann erreichen oder dem sie asymptotisch zustreben, so müssen folgende zwei Bedingungen erfüllt sein:
- Es existiert ein aus , für das gilt:
- Der Betrag der Ableitung ist betragsmäßig höchstens 1, so daß eine Fixpunktiteration gemäß dem Satz von Banach kontrahierend ist und einer der Lösungen von zustrebt
Während die erste Bedingung zur exakten Gleichung führt, ergibt sich aus der zweiten eine Betragsungleichung bzw. . Um die zu bestimmen, die beides erfüllen, müssen wir die Variable eliminieren, was hier der einfachen Lösung einer quadratischen Gleichung in entspricht, die auf führt. Die Randkurve, auf der diese Ungleichung mit Gleichheit erfüllt wird, ist hierbei das Interessante, denn hier wird, wenn man den Rand von innen nach außen überschreitet, gerade der bis dahin attraktive Fixpunkt repulsiv.
Parametrisiert man den Einheitskreis als , so läßt sich dies in überführen (siehe Bildfolge). Dies ist die Parametrisierung einer Kardioide, der sogenannten Hauptkardioide .
Bild 6: Transformation des Einheitskreises durch Skalieren und Verschieben - Quadrieren - Spiegeln und erneutes Verschieben - in die Hauptkardioide
Grundsätzlich werden wir bei dieser Vorgehensweise bleiben, auch wenn wir noch erheblich schwereres algebraisches Geschütz brauchen werden, um die Randkurve von zu bestimmen.
Doch zunächst sehen wir uns noch den Bereich an, der die Parameter beinhaltet, für die es in den entsprechenden Julia-Mengen einen attraktiven 2-Zyklus gibt. Das heißt nichts anderes, als daß die zweite Iterierte, , einen Fixpunkt besitzt, für den die Ableitung betragsmäßig höchstens 1 wird. ist zwar schon ein Polynom 4. Grades, doch werden mit der Fixpunktgleichung auch bereits die Fixpunkte der Abbildung selbst erfaßt, so daß nach dem Ausdividieren von nur noch ein Polynom 2. Grades übrigbleibt, und zwar . Die Betragsbedingung lautet . Da hier der Grad des Polynoms in Betragsstrichen größer ist als der dessen, welches exakt zu erfüllen ist, können wir eine Polynomdivision durchführen und nur deren Rest weiterbetrachten, was uns sofort auf bringt. Der Bereich der Parameter, denen in der Julia-Menge attraktive 2-Zyklen entsprechen, ist also exakt eine Kreisscheibe mit Radius um . Wir werden noch sehen, daß dies der einzige exakte Kreis im Rand von ist.
Bild 7: So weit haben wir die Mandelbrotmenge nun schon erzeugt. Man sieht: Es wird...
7. Die hyperbolische Menge und darüber hinaus
Nun werden wir mit dieser Methodik noch den Bereich der stabilen 3-Zyklen untersuchen. Auch wenn wir es jetzt nicht mehr mit leichten Einsichten zu tun haben, sondern wir das ganze Instrumentarium von Lösung linearer Systeme über die Auflösung polynomialer Gleichungen bis hin zur Körpertheorie in Anschlag bringen müssen, bleibt uns hierbei doch das Glück hold.
Exakt zu erfüllen:
bzw.
Betragsmäßig zu erfüllen:
bzw.
Nach der Polynomdivision durch
bleibt in Betragsstrichen
.
Wenn wir hierin zur Randkurve übergehen, können wir dies zu gegebener Zeit als
umschreiben.
Als nächstes haben wir die Variable zu eliminieren. Dazu suchen wir ein Polynom, über , welches unter der Voraussetzung, daß die Gleichung 6. Grades erfüllt, die algebraische Zahl annulliert. Wir wissen aus der Körpertheorie, daß dies ebenfalls eine ganz-algebraische Zahl 6. Grades über ist, und machen also den Ansatz
,
in welches wir den zu annullierenden Term einsetzen. Dies ergibt zwar ein Polynom 24. Grades, dessen Koeffizienten von abhängen, doch nachdem wir wiederum ausdividiert haben, bleibt als Rest formal ein Polynom 5. Grades (konkret in diesem Fall sogar 4. Grades), dessen Koeffizienten Linearkombinationen der Koeffizienten und über sind. Diese Linearkombinationen setzen wir alle gleich Null, so daß sich ein lineares Gleichungssystem mit 6 Gleichungen in 6 Unbekannten ergibt, welches wir lösen. Dabei stellt sich heraus, daß 4 dieser 6 Gleichungen redundant sind, d.h. von den 6 Koeffizienten sind 4 frei wählbar und nur der des quadratischen und der des konstanten Gliedes sind von diesen abhängig. Wenn wir schon die freie Wahl haben, machen wir es uns so einfach wie möglich und wählen alle freien Koeffizienten des linearen Systems gleich Null.
Damit ergibt sich für :
Und wir haben noch ein zweites Mal Glück. Dieses Polynom faktorisiert praktischerweise in drei in quadratische Polynome, und zwar ist . Nur der letzte dieser Faktoren hat tatsächlich eine Bedeutung für die Randkurve; die beiden übrigen Faktoren haben wir uns durch die Methode eingehandelt, die der Anwendung einer Polynomfunktion 6. Grades auf beiden Seiten entspricht.
Da wir so gewählt hatten, daß , wandelt sich die Betragsungleichung in , was die Parametrisierung von als gestattet. Wenn wir dies in für einsetzen, definiert dies implizit die Randkurve von , wobei von bis läuft. Lösen wir diese Gleichung nach auf – hierbei haben wie das Cardanosche Verfahren anzuwenden –, so gewinnen wir eine Parametrisierung der Randkurven sämtlicher Zusammenhangskomponenten von .
Hier ist der komplexe Term für die "Hände" von und die kardioidenähnliche Kurve bei . Ein scharfer Vergleich im nächsten Abschnitt wird zeigen, daß es sich um keine exakten Kreise und keine exakte Kardioide handelt.
nimmt wie üblich die Werte 0, 1, 2 an. Einsetzen von für liefert eine explizite Parametrisierung der blauen Bereiche im folgenden Plot.
Bild 8: Wir haben die "Hände" des Apfelmännchens exakt parametrisiert, und auch der Hauptkörper des größten Satelliten in der Antenne zeigt sich bereits. Kenneth Falconer nennt in [3] diese Strukturen "kreisförmig (circular)" - hat er damit recht?
Mit dem hier angewendeten Verfahren, ein nichtlineares System aus zwei Polynomgleichungen in zwei Variablen zu lösen, ließe sich prinzipiell auch die Bestimmung sämtlicher hyperbolischer Mengen angehen. Man hat allerdings zu beachten, daß durch die Anwendung einer höhergradigen Polynomfunktion die Lösungsmenge vergrößert wird, d.h. es wird so sein, daß das aufgefundene Polynom zwar stets faktorisiert, doch daß dessen Faktoren nicht alle zu einer Randkurve einer hyperbolischen Menge in gehören. So findet man für den Rand der hyperbolischen Menge ein Polynom 12. Grades , doch nur dessen Faktor entsprechen tatsächlich die Ränder sämtlicher hyperbolischer Bereiche, deren Elementen in den ihnen entsprechenden Julia-Mengen attraktive Zyklen der Länge 4 zuzuordnen sind. Das rasche Anwachsen der Polynomgrade mit der Zyklenlänge verdeutlicht aber, wie rasch die Grenzen der Bestimmbarkeit exakter Grenzkurven erreicht ist. Die explizite Darstellung des Randes von ist schon aufgrund des Satzes von Abel-Ruffini nicht mehr möglich.
Die Anzahl der Zusammenhangskomponenten dieses Bereichs ist übrigens, wie man aufgrund des Beispiels vermuten könnte, durch den Grad dieses Polynoms in festgelegt. hat Grad 6 in , und es existieren 6 hyperbolische Zusammenhangskomponenten: zwei sind die 4-Knospen an der Hauptkardioide rechts von den "Händen" von , eine bildet den "Dutt" am "Kopf" des Apfelmännchens, eine vierte bildet den kardioidenartigen Hauptkörper eines kleinen Satelliten in der "Antenne" von bei (diese Zahl ist nebenbei für mich das schönste Ergebnis dieser Befassung mit der Materie überhaupt), und die übrigen beiden sind kardioidenförmige hyperbolische Bereiche, die im Komplexen innerhalb von Verästelungen liegen, die von den beiden "Händen" von ausgehen.
Wurzeln
Doch auch wenn es nicht möglich ist, die gesamte Randkurve eines hyperbolischen Bereichs zu bestimmen, um sie sich z.B. plotten zu lassen, so läßt sich zumindest der zu gehörige Punkt bestimmen. Für diesen nimmt den Wert an, wobei die Zyklenlänge ist. Das sich so ergebende Polynom in läßt sich stets in Faktoren zerlegen, denn es gibt für jedes hyperbolische Bereiche auf den Hauptkörpern der Mandelbrotmenge, die an ausgezeichneten Punkten der Randkurven als "Knospen" sitzen. Deren Aufsitzpunkte erfüllen also stets Polynomgleichungen niedrigeren Grades. (Welche das sind, wird in Kapitel 9 näher erörtert.) Für die übrigen hyperbolischen Bereiche - man denke an kleine Satelliten irgendwo in den Verästelungen von - hat man aber weitere Faktoren, die mit den Standardverfahren im Allgemeinen nicht mehr faktorisierbar sind. Doch kann zumindest numerisch im Komplexen dieser Anfangspunkt gefunden werden. Als dieser ist die Wurzel der hyperbolischen Menge definiert.
Zwei Fälle sind also möglich: Die Wurzel sitzt entweder auf der Randkurve eines größeren hyperbolischen Bereichs, dann ist die hyperbolische Menge kreisähnlich, oder sie hängt in den Verästelungen des Randes von und ähnelt dann einer Kardioide, an der sich sämtliche Substrukturen der gesamten Mandelbrotmenge unter leichter Verzerrung wiederentdecken lassen. Die Zyklenlänge, denen die Elemente dieser Substrukturen jeweils entsprechen, ergibt sich dabei aus der Zyklenlänge, die der mit dieser verbundenen Kardioide zugeordnet ist, multipliziert mit der Zyklenlänge der gleichartigen Knospe an der Hauptkardioide .
8. Die Form der Grenzkurven
Noch einmal zurück zu den Knospen von . Deren gesamte Randkurve gliedert sich in drei Zusammenhangskomponenten: eine liegt symmetrisch auf der reellen Achse bei , innerhalb der reellen "Antenne" der Mandelbrotmenge und ähnelt einer verkleinerten Kopie der Hauptkardioide. Die anderen beiden liegen symmetrisch bezüglich der reellen Achse und berühren die Hauptkardioide genau dort, wo diese waagerechte Tangenten hat. Kenneth Falconer behauptet in [3], diese Kurven seien "kreisförmig (circular)". Ich habe mich gefragt: Sind dies wirklich exakte Kreise und, falls ja, warum lassen sich ihre Gleichungen nicht einfacher auffinden, warum muß man dafür so tief in die algebraische Werkzeugkiste greifen?
Um die Frage zu beantworten, verlassen wir kurz die strengen Pfade der Algebra und bedienen uns der komplexen Analysis, um den Schmiegekreis der Zusammenhangskomponente im Kontaktpunkt mit der Hauptkardioide zu bestimmen. Dazu setzen wir an und wählen die drei Parameter , , so, daß für am Kreis die reellen und imaginären Aufpunktkoordinaten, die Richtung des Tangentenvektors sowie die Krümmung ( und sind Real- bzw. Imaginärteil der Parametrisierung der Randkurve) mit denjenigen der Randkurve identisch werden.
Wir erhalten:
Wäre die Zusammenhangskomponente der Randkurve tatsächlich ein Kreis, so müßte sie nicht nur lokal im Kontaktpunkt, sondern global mit ihrem eigenen Schmiegekreis zusammenfallen. Schaut man sich den entsprechenden Plot sehr genau an, so erkennt man aber, daß die Kurven in dem Bereich, der zwischen "halb elf" und "halb zwei" auf dem Schmiegekreis liegt, auseinanderlaufen. Die Ausschnittsvergrößerung bestätigt die Diskrepanz, und ebenso täte dies ein Nachrechnen – welches wir uns aber an dieser Stelle ersparen, da der Augenschein deutlich genug ist. Diese Teile der Randkurve von sind also tatsächlich keine Kreise - wenn auch die Abweichung unter 1% liegt.
Bild 9: Vergleich der oberen "Hand" des Apfelmännchens, eines Teils des Randes von , mit ihrem Schmiegekreis im Kontaktpunkt mit der Hauptkardioide. Bei exakter Kreisform müßten die Kurven überall übereinstimmen; es erweist sich jedoch, daß sie dies nicht tun
Bild 10: Ausschnittvergrößerung des linken oberen Bildbereichs
Man kann sich dies auch anderweitig überlegen, indem man die Entstehung der Kurve aus dem Kreis nachvollzieht, mit dem wir das Argument der Betragsungleichung parametrisiert haben:
.
Dieser wird einer polynomialen Abbildung 2. Grades unterworfen, und solche Abbildungen überführen im Allgemeinen einen Kreis in eine Epizykloide, wie wir am Beispiel der Hauptkardioide bereits sahen. Doch in diesem Fall wird noch eine inverse polynomiale Abbildung nachgeschaltet, und zwar eine, die ein Polynom 3. Grades invertiert – ein Term, der dritte Wurzeln enthält. Die Mehrdeutigkeit dieses Wurzelausdrucks ist auch dafür verantwortlich, daß die Randkurven der höheren hyperbolischen Mengen nicht mehr zusammenhängend sind, sondern in so viele Zusammenhangskomponenten zerfällt, wie der Grad des Polynoms in angibt. Es gilt jedoch: Je kleiner der Radius des Ausgangskreises und je weiter seine Peripherie vom Koordinatenursprung entfernt bleibt, desto geringer fallen die Abweichungen von einem exakten Kreis nach den beiden Transformationen aus. Hieraus können wir den Schluß (zumindest aber die begründete Vermutung) ziehen, daß die beiden spiegelsymmetrischen Knospen von diejenigen unter den nicht-kardioidenförmigen hyperbolischen Mengen in der ganzen Mandelbrotmenge sind, die am stärksten von der exakten Kreisform abweichen.
Bild 11: Vergleich des kardioidenartigen Hauptkörpers des Satelliten innerhalb von mit einer exakten Kardioide, die in der Symmetrieachse die gleiche Breite hat.
9. Zur Selbstähnlichkeit von
Wir werden uns zum Abschluß mit der Frage beschäftigen, welchen Regeln die Knospenbildung auf der Randkurve einer hyperbolischen Menge gehorcht. Es wird dazu hilfreich sein, sich den verwandten Vorgang der Periodenverdopplung im Feigenbaum-Szenario vor Augen zu führen.
Der Fixpunkt bzgl. wird repulsiv, wenn bei stetiger Variation von den Wert 1 überschreitet, wobei rechts und links von ein neues Paar Fixpunkte der zweiten Iterierten von entsteht, die ab dann einen attraktiven 2-Zyklus bilden. Dieser Typ der Bifurkation heißt aufgrund seiner Form Heugabelbifurkation. Es ist dies der einzig mögliche Bifurkationstyp im Reellen mit der Kodimension 1, weshalb sich im Feigenbaum-Szenario Zyklen mit ungerader Länge nur durch spontanes Aufreißen eines Fensters im Chaos ergeben können.
Da in den Julia-Mengen wesentlich die gleiche Situation, allerdings im Komplexen, vorliegt, sind hier auch Bifurkationen mit Kodimension 1 möglich, bei denen sich beim Repulsivwerden eines Fixpunktes oder eines Zyklus mehrere Fixpunkte einer höheren als der zweiten Iterierten von ergeben. Damit z.B. ein 3-Zyklus entsteht, muß die dritte Iterierte von lokal um den Fixpunkt wie aussehen. Den entsprechenden Bifurkationstyp könnte man z.B. "Dreifußbifurkation" nennen.
Sehen wir uns an, was beim Iterieren von in der Nähe von dessen Fixpunkt geschieht, wenn den kritischen Wert annimmt: nicht nur stimmt mit überein, sondern auch wird Eins. -faches Iterieren von ändert an dieser Tatsache wenig, da sich die Funktion sehr eng an die komplexe Hauptdiagonale anschmiegt.
Liegt auf dem Rand der Hauptkardioide , so entspricht ihm ein kritischer Fixpunkt auf dem Kreis . Es ist . Unter welchen Bedingungen geht aus ein attraktiver -Zyklus hervor? Die Ableitung der -ten Iterierten, die sich nach mehrfacher Anwendung der Kettenregel als Produkt gleicher Faktoren entpuppt,
ist gleich 1, denn die Iterierte wird lokal, ebenso wie die Ausgangsfunktion, durch approximiert; erst der -te Taylorkoeffizient der Entwicklung von um ist von Null verschieden, was das Auftreten des neuen -Zyklus bedingt. Folglich ist
mit . Auf die Kardioide umgerechnet:
.
Für ergibt dies gerade die Spitze der Hauptkardioide; diese Lösung kann also fallengelassen werden. Für existieren am Hauptkörper des Apfelmännchens also Knospen ( ist die Eulersche Phi-Funktion, die zu teilerfremde Zahlen zählt), in welchen sich diejenigen Parameter finden, in deren Julia-Mengen attraktive -Zyklen auftreten.
In analoger Weise kann man auch von beliebigen hyperbolischen Mengen über deren Rand in weitere Knospen von vorstoßen. Es gilt dabei, daß die der Knospe entsprechende Zyklenlänge sich als das Produkt aus derjenigen der ihr entsprechenden Knospe, die mit über gemeinsame Randpunkte verbunden ist, und derjenigen des kardioidenförmigen Hauptkörpers, mit dem sie selbst über gemeinsame Randpunkte verbunden ist, ergibt.
Wir betrachten zum Beweis die -te Iteration von am Ort eines -zyklischen Punktes , dessen Zyklus kritisch ist, d.h. der unitäre Ableitung hat. Diese entspricht lokal bis in -te Näherung der Identität; insbesondere ist ihre Ableitung gerade 1. Für deren Ableitung gilt
,
und der Betrag dieser Ableitung ist, so wie wir die Randkurve parametrisiert haben, gleich 1, also ist für ein . Demnach ist und mit (bei liegt der Kontaktpunkt mit dem Hauptkörper oder, falls die Knospe kardioidenähnlich ist, deren Spitze, und bei mit befinden sich Knospen mit geringerer Zyklenlänge), d.h. überall dort, wo der Parameter ein rationales Verhältnis zu aufweist, dessen Nenner ist, befindet sich eine Unterknospe, die solche Parameter enthält, in deren Julia-Mengen Zyklen auftreten, die -mal so lang sind wie diejenigen, die den Punkten in der Oberknospe zugeordnet sind. Die übrigen hyperbolischen Bereiche gleicher Zyklenlänge liegen dann nicht als Knospen vor, sondern als Hauptkörper von Satelliten.
Man kann die hyperbolischen Mengen am Rand eines Hauptkörpers also mit einem Indexpaar indizieren: |
|