Logistische Gleichung

Die logistische Gleichung wurde ursprünglich 1837 von Pierre François Verhulst^[1] als demographisches mathematisches Modell eingeführt. Die Gleichung ist ein Beispiel dafür, wie komplexes, chaotisches Verhalten aus einfachen nichtlinearen Gleichungen entstehen kann. Infolge einer richtungsweisenden Arbeit des theoretischen Biologen Robert May aus dem Jahr 1976^[2] fand sie weite Verbreitung. Bereits 1825 stellte Benjamin Gompertz in einem verwandten Zusammenhang eine ähnliche Gleichung vor.^[3]

Die zugehörige Dynamik kann anhand eines sogenannten Feigenbaumdiagramms (siehe unten) veranschaulicht werden. Eine wichtige Rolle spielt dabei die schon 1975 von Mitchell Feigenbaum gefundene Feigenbaum-Konstante.

Das demographische Modell

Für den stetigen Fall siehe logistische Funktion.

Es werden mathematische Gesetzmäßigkeiten gesucht, die die Entwicklung einer Population modellhaft darstellen. Aus der Größe $X_{n}$ der Population zu einem gewissen Zeitpunkt $$ n $$ soll auf die Größe $X_{n+1}$ nach einer Fortpflanzungsperiode (z. B. nach einem Jahr) geschlossen werden.

Das logistische Modell berücksichtigt zwei Einflüsse:

Durch Fortpflanzung vermehrt sich die Population geometrisch. Die Individuenzahl ist im Folgejahr um einen Wachstumsfaktor $q_{\mathrm {f} }$ größer als die aktuelle Population.
Durch Verhungern verringert sich die Population. Die Individuenzahl vermindert sich in Abhängigkeit von der Differenz zwischen ihrer aktuellen Größe und einer theoretischen Maximalgröße Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): G mit der Proportionalitätskonstante Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): q_\mathrm v . Der Faktor, um den sich die Population vermindert, hat also die Gestalt Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): q_\mathrm h = (G-X_n) \, q_\mathrm v .

Um bei der Berechnung der Population im Folgejahr beide Prozesse zu berücksichtigen, multipliziert man die aktuelle Population Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): X_n sowohl mit dem Vermehrungsfaktor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): q_\mathrm f als auch mit dem Hungerfaktor Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): q_\mathrm h . Man erhält damit die logistische Gleichung

X_{n+1}=q_{\mathrm {f} }\,q_{\mathrm {v} }\,X_{n}\,(G-X_{n})

.

Um die folgenden mathematischen Untersuchungen zu vereinfachen, wird die Populationsgröße Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): X_n oft als Bruchteil Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): x_n der Maximalgröße Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): G angegeben:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): x_n = \frac{X_n}{G} .

Außerdem werden Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): G , Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): q_\mathrm f und $q_{\mathrm {v} }$ zusammengefasst zum Parameter Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r :

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r =G\, q_\mathrm f \, q_\mathrm v .

Damit ergibt sich die folgende Schreibweise für die logistische Gleichung:

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): x_{n+1} = r\, x_n\,(1 - x_n / K) .

Hierbei ist Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): K die Kapazität des Biotops. Das heißt, es ist die Population, die bei geeigneter Wahl von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r dem Fixpunkt der Dynamik entspricht.

Das mathematische Modell

Man kann ohne Beschränkung der Allgemeinheit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): K=1 setzen. Dann ergibt sich

x_{n+1}=r\cdot x_{n}\cdot (1-x_{n})

.

Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): x_n ist dabei eine Zahl zwischen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): 0 und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): 1 . Sie repräsentiert die relative Größe der Population im Jahr Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): n . Die Zahl Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): x_0 steht also für die Startpopulation (im Jahr 0). Der Parameter Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r ist immer positiv, er gibt die kombinierte Auswirkung von Vermehrung und Verhungern wieder.

Verhalten in Abhängigkeit von r

Die Animation unten zeigt Zeitreihenentwicklungen der Logistischen Gleichung im Zeit- und Frequenzbereich (Fourier-Analysis), die sich für wachsende Parameter $$ 2<r<4 $$ ergeben. Startwert ist jedes Mal Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): x_0=\tfrac12 .

Datei:Logistic map animation.gif

Gut sichtbar sind die Zonen der Intermittenz innerhalb des deterministischen Chaos.

Bei verschiedenen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r können die folgenden Verhaltensweisen für große Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): n beobachtet werden. Dabei hängt dieses Verhalten nicht vom Anfangswert ab, sondern nur von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r :

Mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r von 0 bis 1 stirbt die Population in jedem Fall aus.
Mit $$ r $$ zwischen 1 und 2 nähert sich die Population monoton dem Grenzwert Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \tfrac{r-1}{r} an.
Mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r zwischen 2 und 3 nähert sich die Population dem Grenzwert Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): \tfrac{r-1}{r} alternierend, d. h. die Werte liegen ab einem bestimmten $$ n $$ abwechselnd über und unter dem Grenzwert.
Mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r zwischen 3 und $1+{\sqrt {6}}$ (etwa 3,45) wechselt die Folge bei fast allen Startwerten (ausgenommen 0, 1 und Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): 1-\tfrac{1}{r} ) zwischen den beiden Umgebungen zweier Häufungspunkte.
Mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r zwischen Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): 1+\sqrt{6} und ungefähr 3,54 wechselt die Folge bei fast allen Startwerten zwischen den Umgebungen von vier Häufungspunkten.
Wird Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r größer als 3,54, stellen sich erst 8, dann 16, 32 usw. Häufungspunkte ein. Die Intervalle mit gleicher Anzahl von Häufungspunkten (Bifurkationsintervalle) werden immer kleiner; das Längenverhältnis zweier aufeinanderfolgender Bifurkationsintervalle nähert sich der Feigenbaumkonstanten. Diese Konstante ist auch in anderen mathematischen Zusammenhängen von Bedeutung. (Zahlenwert: δ ≈ 4,6692016091029906718532038204662016172581…).
Bei $$ r $$ annähernd 3,57 beginnt das Chaos: Die Folge springt zunächst periodisch zwischen den Umgebungen der nun instabilen Häufungspunkte umher. Mit weiter wachsendem Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r verschmelzen diese Intervalle so, dass sich deren Anzahl im Rhythmus der Feigenbaumkonstante halbiert, bis es nur noch ein Intervall gibt, in dem die Folge chaotisch ist. Perioden sind dann nicht mehr erkennbar. Winzige Änderungen des Anfangswertes resultieren in unterschiedlichsten Folgewerten – eine Eigenschaft des Chaos.
Bei vielen Koeffizienten zwischen 3,57 und 4 kommt es zu chaotischem Verhalten, obwohl für bestimmte Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r wieder Häufungspunkte (d. h. stabile periodische Orbits, gegen die fast jeder Anfangswert konvergiert) vorhanden sind. Beispielsweise existieren in der Nähe von Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r = 3{,}82 bei steigendem Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r erst 3, dann 6, 12 usw. Häufungspunkte. Ebenso gibt es r-Werte mit 5 oder mehr Häufungspunkten – alle Periodendauern tauchen auf.
Drei tiefliegende mathematische Sätze besagen folgendes: (1) jedes noch so kleine Intervall von Koeffizienten enthält Parameter, für die es stabile periodische Orbits gibt (so dass die Dynamik eben nicht chaotisch ist): also nicht-chaotische Parameter sind "dicht" im Intervall der Koeffizienten. Chaotische Parameter enthalten also keine Intervalle. Aber (2) die chaotischen Parameter haben positives Maß: also mit echt positiver Wahrscheinlichkeit liefert ein zufälliger Parameter chaotische Dynamik. Schließlich (3) hat fast jeder reelle Koeffizient $$ r $$ (im Sinne voller Wahrscheinlichkeit) entweder einen stabilen periodischen Orbit (gegen den fast jeder Anfangswert konvergiert) oder ist in strengem Sinne "chaotisch". (Weitere dynamische Möglichkeiten gibt es, haben aber Wahrscheinlichkeit null.)
Für Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r größer 4 divergiert die Folge für fast alle Anfangswerte und verlässt das Intervall Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): [0;1] .

Dieser Übergang von konvergentem Verhalten über Periodenverdopplungen zu chaotischem Verhalten ist generell für nichtlineare Systeme typisch, die in Abhängigkeit von einem Parameter chaotisches oder nicht-chaotisches Verhalten zeigen.

Eine Erweiterung des Wertebereiches auf die komplexen Zahlen führt nach einer Koordinatentransformation zur Mandelbrotmenge.

Beispiel

Die „logistische Kurve“ mit einer Wachstumsrate Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r = 1{,}4 verläuft S-förmig. Ab einem Wert um 3,6 bricht Chaos aus, wie die Abbildung mit Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r = 3{,}81 illustriert.^[4]

Datei:Logistische Gleichung kleines r.JPG

„Logistische Kurve“ mit einer Wachstumsrate r = 1,4

Datei:Logistische Gleichung großes r.JPG

„Logistische Kurve“ mit einer Wachstumsrate r = 3,81

Graphische Darstellung

Das folgende Bifurkationsdiagramm, bekannt als Feigenbaum-Diagramm, fasst diese Beobachtungen zusammen. Die horizontale Achse gibt den Wert des Parameters Fehler beim Parsen (MathML mit SVG- oder PNG-Rückgriff (empfohlen für moderne Browser und Barrierefreiheitswerkzeuge): Ungültige Antwort („Math extension cannot connect to Restbase.“) von Server „https://wikimedia.org/api/rest_v1/“:): r an und die vertikale Achse die Häufungspunkte für die Folge $x_{n}$ .

LogisticMap BifurcationDiagram.png

Bifurkationsdiagramm der logistischen Gleichung
Bifurkation wiki.gif

Hochauflösende Version ohne Skala
Hochauflösender Ausschnitt des Bifurkationsdiagramms der logistischen Gleichung

Zusammenhang mit der Mandelbrotmenge (nach Koordinatentransformation)

Analytische Lösung

Für den Parameter $$ r=2 $$ existiert eine analytische Lösung:

x(n)={\frac {1}{2}}-{\frac {1}{2}}(1-2x_{0})^{\left(2^{n}\right)}

.

Für die Parameter $$ r=-2 $$ und $$ r=4 $$ können ebenfalls analytische Lösungen angegeben werden.

Einzelnachweise

↑ Pierre-François Verhulst: Notice sur la loi que la population suit dans son accroissement. In: Correspondance Mathématique et Physique. Band 10, 1838, ZDB-ID 428605-4, S. 113–121.
↑ Robert May: Simple mathematical models with very complicated dynamics Nature V. 261, S. 459–467 (10 Juni 1976)
↑ Benjamin Gompertz: On the Nature of the Function Expressive of the Law of Human Mortality, and on a New Mode of Determining the Value of Life Contingencies. In: Philosophical Transactions of the Royal Society of London. Vol. 115, 1825, ISSN 0260-7085, S. 513–585.
↑ Jürgen Beetz: 1 + 1 = 10. Mathematik für Höhlenmenschen. Springer, Heidelberg 2012, ISBN 978-3-8274-2927-8, S. 313 f.

Weblinks

[1] Pierre-François Verhulst: Notice sur la loi que la population suit dans son accroissement. In: Correspondance Mathématique et Physique. Band 10, 1838, ZDB-ID 428605-4, S. 113–121.

[2] Robert May: Simple mathematical models with very complicated dynamics Nature V. 261, S. 459–467 (10 Juni 1976)

[3] Benjamin Gompertz: On the Nature of the Function Expressive of the Law of Human Mortality, and on a New Mode of Determining the Value of Life Contingencies. In: Philosophical Transactions of the Royal Society of London. Vol. 115, 1825, ISSN 0260-7085, S. 513–585.

[4] Jürgen Beetz: 1 + 1 = 10. Mathematik für Höhlenmenschen. Springer, Heidelberg 2012, ISBN 978-3-8274-2927-8, S. 313 f.

[1]

[2]

[3]

[4]