Doppelpendel: Unterschied zwischen den Versionen

Aktuelle Version vom 4. April 2021, 16:20 Uhr

Schematische Zeichnung eines Doppelpendels.

Trajektorie eines idealisierten Doppelpendels

Das Doppelpendel ist ein beliebtes Modell zur Demonstration von chaotischen Prozessen. Es ist zugleich eines der einfachsten nichtlinearen Dynamischen Systeme, welches chaotisches Verhalten zeigt. An die Masse $m_{1}$ eines Pendels mit der Länge $L_{1}$ wird ein weiteres Pendel der Länge $L_{2}$ mit Masse $m_{2}$ gehängt. Die Herleitung der Bewegungsgleichung zum Berechnen der Bewegung des Doppelpendels lässt sich vereinfachen, wenn man starre, masselose Pendelstangen und Reibungsfreiheit annimmt.

Ein Merkmal eines chaotischen Systems ist, dass es Anfangsbedingungen $x_{i}$ gibt, sodass ein weiteres Experiment mit nahezu identischen Anfangsbedingungen $x_{i} + Δ x$ , die sich nur um eine infinitesimale Störung $Δ x$ unterscheiden, nach kurzer Zeit ein anderes Verhalten zeigt. Diese sensible Abhängigkeit lässt sich durch Berechnen von Ljapunow-Exponenten der Trajektorien charakterisieren.

Herleitung der Bewegungsgleichungen

Wenn $l_{1}$ und $l_{2}$ die Längen der (masselosen) Verbindungsstangen, $m_{1}$ und $m_{2}$ die Pendelmassen, $θ_{1}, θ_{2}$ die Auslenkung vom Lot und $g$ die Erdbeschleunigung bezeichnet, dann findet man für die Positionen von $m_{1}$ und $m_{2}$ :

\begin{aligned} x_{1} & = l_{1} \sin (θ_{1}) \\ y_{1} & = - l_{1} \cos (θ_{1}) \end{aligned}

und

\begin{aligned} x_{2} & = l_{1} \sin (θ_{1}) + l_{2} \sin (θ_{2}) \\ y_{2} & = - l_{1} \cos (θ_{1}) - l_{2} \cos (θ_{2}) . \end{aligned}

Damit lassen sich die Geschwindigkeiten der Massen, welche für den nächsten Schritt notwendig sind, bestimmen, wobei $\dot{θ_{1}}$ die zeitliche Ableitung von $θ_{1}$ ist:

\begin{aligned} u_{1} & = \frac{\partial x_{1}}{\partial t} = {\dot{θ}}_{1} l_{1} \cos (θ_{1}) \\ v_{1} & = \frac{\partial y_{1}}{\partial t} = {\dot{θ}}_{1} l_{1} \sin (θ_{1}) \end{aligned}

und

\begin{aligned} u_{2} & = \frac{\partial x_{2}}{\partial t} = {\dot{θ}}_{2} l_{2} \cos (θ_{2}) + u_{1} \\ v_{2} & = \frac{\partial y_{2}}{\partial t} = {\dot{θ}}_{2} l_{2} \sin (θ_{2}) + v_{1} . \end{aligned}

Unter Verwendung des Lagrange-Formalismus $L = T - V$ wobei $T$ die kinetische Energie der beiden Pendelmassen und $V$ ihre potentielle Energie im konstanten Gravitationsfeld ist, mit

\begin{aligned} T_{1} & = \frac{1}{2} m_{1} (u_{1}^{2} + v_{1}^{2}) \\ T_{2} & = \frac{1}{2} m_{2} (u_{2}^{2} + v_{2}^{2}) \\ V_{1} & = m_{1} g y_{1} \\ V_{2} & = m_{2} g y_{2} \end{aligned}

erhält man

\begin{aligned} T & = T_{1} + T_{2} = \frac{1}{2} m_{1} {\dot{θ}}_{1}^{2} l_{1}^{2} + \frac{1}{2} (m_{2} {\dot{θ}}_{2}^{2} l_{2}^{2} + m_{2} {\dot{θ}}_{1}^{2} l_{1}^{2} + 2 m_{2} {\dot{θ}}_{1} l_{1} {\dot{θ}}_{2} l_{2} \cos (θ_{1} - θ_{2})) \\ V & = V_{1} + V_{2} = - m_{1} g l_{1} \cos (θ_{1}) - m_{2} g l_{1} \cos (θ_{1}) - m_{2} g l_{2} \cos (θ_{2}) . \end{aligned}

Damit ergibt sich für die Lagrange-Funktion dann

L = \frac{1}{2} (m_{1} + m_{2}) {\dot{θ}}_{1}^{2} l_{1}^{2} + \frac{1}{2} m_{2} {\dot{θ}}_{2}^{2} l_{2}^{2} + m_{2} {\dot{θ}}_{1} l_{1} {\dot{θ}}_{2} l_{2} \cos (θ_{1} - θ_{2}) + (m_{1} + m_{2}) g l_{1} \cos (θ_{1}) + m_{2} g l_{2} \cos (θ_{2}) .

Unter Verwendung der Euler-Lagrange-Gleichung

\frac{d}{d t} \frac{\partial L}{\partial \dot{θ_{i}}} - \frac{\partial L}{\partial θ_{i}} = 0

erhält man damit nach einigen Umformungen

\begin{aligned} {\ddot{θ}}_{1} & = - \frac{m_{2}}{m_{1} + m_{2}} \frac{l_{2}}{l_{1}} ({\ddot{θ}}_{2} \cos (θ_{1} - θ_{2}) + {\dot{θ}}_{2}^{2} \sin (θ_{1} - θ_{2})) - \frac{g}{l_{1}} \sin (θ_{1}) \\ {\ddot{θ}}_{2} & = - \frac{l_{1}}{l_{2}} ({\ddot{θ}}_{1} \cos (θ_{1} - θ_{2}) - {\dot{θ}}_{1}^{2} \sin (θ_{1} - θ_{2})) - \frac{g}{l_{2}} \sin (θ_{2}) \end{aligned}

die Winkelbeschleunigungen für $θ_{1}$ und $θ_{2}$ , welche die Evolution des Pendels beschreiben.

In den Bewegungsgleichungen treten Winkelfunktionen $(\sin, \cos)$ der Zustandsgrößen und auch Ableitungen auf. Es handelt sich also um ein nichtlineares System. Im Spezialfall kleiner Auslenkungen als Anfangsbedingungen $θ_{1} (t = 0), θ_{2} (t = 0) \approx 0$ , lassen sich die Bewegungsgleichungen allerdings mittels der Kleinwinkelnäherung vereinfachen. Dann lassen sich beispielsweise weitere Spezialfälle wie $m_{1} ≪ m_{2}$ oder $m_{1} ≫ m_{2}$ mit analytischen Ansätzen betrachten, die eine näherungsweise harmonische Lösung haben; diese kann auch analytisch bestimmt werden.

Lösung der Bewegungsgleichungen

Die Bewegungsgleichungen für die generalisierten Koordinaten $θ_{1}$ und $θ_{2}$ stellen ein nichtlineares System von zwei gekoppelten Differentialgleichungen dar, welches analytisch nicht lösbar ist. Es kann bei vier bekannten Anfangswerten ( $θ_{1}, θ_{2}, {\dot{θ}}_{1}, {\dot{θ}}_{2}$ ) mit numerischen Verfahren gelöst werden. Hierbei werden also die anfänglichen Auslenkungen (z. B. 30° und 30°) und die anfänglichen Geschwindigkeiten (z. B. $0 \frac{rad}{s}$ und $0 \frac{rad}{s}$ ) eingegeben und damit dann die Evolution des Pendels berechnet.

Mittels Trigonometrie können die Winkel $θ_{1}$ und $θ_{2}$ in die kartesischen Koordinaten $(x_{1}, y_{1}, x_{2}, y_{2})$ der Massenpunkte überführt werden.

Anwendungen

Die Blide nutzt den energetischen Austausch zwischen den Massen der Pendel, um ein kleines Gewicht am äußeren Pendel durch die potentielle Energie eines großen Gewichtes am inneren Pendel zu beschleunigen.
Eine Kirchenglocke mit Klöppel bildet ein Doppelpendel, allerdings mit zusätzlicher Beschränkung für den Winkel $θ_{2}$ .

Auswertung des chaotischen Verhaltens

Zur Betrachtung des chaotischen Verhaltens des Doppelpendels gibt es eine Reihe von Möglichkeiten. Oft kann mittels einfachster Berechnungen eine Aussage über chaotisches Verhalten getroffen werden. Beispiele sind der maximale Ljapunow-Exponent (MLE) oder Bifurkationsdiagramme.

Maximaler Ljapunow-Exponent

Der MLE ist der sog. maximale Ljapunow-Exponent (maximum Lyapunov exponent) und beschreibt die „Stärke“ des chaotischen Verhaltens. Er ist Bestandteil des Ljapunow-Spektrums welches alle Ljapunow-Exponenten (je einer pro Freiheitsgrad) beinhaltet. Man geht davon aus, dass das System eine Störung in der Richtung des MLE hat und da er das größte Wachstum zeigt, ist zu erwarten, dass der MLE nach einer gewissen Zeit die Evolution des Systems dominiert. Ein positiver MLE zeigt normalerweise ein chaotisches System an. Er wird berechnet mit:

λ_{\max} (t) = \frac{1}{t} \ln (\frac{‖ δ (t) ‖}{‖ δ (t = 0) ‖})

Bei zwei Experimenten mit einer anfänglichen Separation von $δ (t = 0) ∽ O (10^{- 8})$ in den Anfangsbedingungen oder sogar weniger, verstärkt sich diese Differenz exponentiell und lässt die Trajektorien divergieren.^[1] Die Separation (der natürliche Logarithmus der obigen Gleichung) kann dann in einem halblogarithmischen Diagramm gegen die Zeit aufgetragen werden. Dann wird mittels linearer Regression die Steigung bestimmt und diese gibt dann den approximierten MLE.

Bifurkationsdiagramm

Darstellung der Bifurkation eines Doppelpendels mit gleichen Längen, gleichen Massen,

g = 9, 81 \frac{m}{s^{2}}

für

θ_{1} = θ_{2}

bei

\dot{θ_{2}} = 0

. Die rote Linie zeigt die Bifurkation des Doppelpendels von harmonischen Oszillationen zu chaotischen Oszillationen.

Bifurkationsdiagramme sind eine Möglichkeit komplexe Informationen über den Phasenraum eines dynamischen Systems in einen zweidimensionalen, visualisierbaren Plot zu komprimieren. Üblicherweise wird die qualitative Änderungen des Verhaltens eines Systems mittels der Variation eines geeigneten Parameters untersucht. So können für das Doppelpendel bspw. das Verhältnis der Massen, das Verhältnis der Längen, die Erdbeschleunigung oder die Anfangsbedingungen herangezogen werden. Durch die kontinuierliche Veränderung des gewählten Bifurkationsparameters wird das System auf Stabilität (periodische, quasi-periodische Lösungen) bzw. auf Chaos geprüft.

Wenn man die anfänglichen Winkel $θ_{1} = θ_{2}$ als den Bifurkationsparameter wählt, lässt sich das qualitativ veränderliche Verhalten des Doppelpendels sehr gut veranschaulichen. Dazu werden die beiden Winkel simultan Stück für Stück erhöht und für jedes Inkrement wird das Doppelpendel erneut integriert (berechnet). Mit diesen Daten lässt sich dann veranschaulichen, wie das System schwingt. Man hat also einen vier-dimensionalen Phasenraum, der sich aus $θ_{1}, θ_{2}, \dot{θ_{1}}, \dot{θ_{2}}$ zusammensetzt. Praktischerweise oszillieren die Winkelgeschwindigkeiten $\dot{θ_{1}}, \dot{θ_{2}}$ , obwohl mit unbestimmter Amplitude, um Null. Daher ist zu erwarten, dass beide immer wieder die Null überqueren. Für ein harmonisch schwingendes System (periodische Lösung) sind die Nullüberquerungen von $\dot{θ_{1}}, \dot{θ_{2}}$ an festen Punkten, da das System immer an bestimmten Punkten ( $θ_{1}, θ_{2}$ ) seine Auf- und Abwärtsbewegung beendet und zurückschwingt. Das ist vergleichbar mit einem normalen starren Pendel. Daher ist im Umkehrschluss zu erwarten, dass das chaotisch schwingende System an allen möglichen Punkten ( $θ_{1}, θ_{2}$ ) die Winkelgeschwindigkeit Null zeigt. Wenn man dann eine "Scheibe" aus dem Phasenraum gesondert betrachtet, bspw. Winkelgeschwindigkeit $\dot{θ_{2}} = 0$ , kann man die Bifurkation des Verhaltens zwei-dimensional darstellen, indem man die Winkelgeschwindigkeit $\dot{θ_{1}}$ gegen die veränderliche Anfangsbedingung aufträgt (siehe rechts).

Siehe auch

Magnetisches Pendel
Multipendel
KAM-Theorem
Kirchenglocke

Weblinks

Commons: Double pendulums – Sammlung von Bildern, Videos und Audiodateien

Java-Doppelpendel (englisch)
Doppelpendel-Simulation in Java und Python (deutsch)
Doppelpendel - Grenzen der Simulation zeigt, dass die Bewegung stets nur für eine kurze Zeitspanne simuliert werden kann
Herleitung der Differentialgleichungen zur Beschreibung des Doppelpendels (englisch)

Einzelnachweise

↑ Wolf, A., Swift, J. B., Swinney, H. L., & Vastano, J. A. (1985). Determining lyapunov exponents from a time series. Physica D: Nonlinear Phenomena, Vol. 16 No. 3: 285 - 317

[1] Wolf, A., Swift, J. B., Swinney, H. L., & Vastano, J. A. (1985). Determining lyapunov exponents from a time series. Physica D: Nonlinear Phenomena, Vol. 16 No. 3: 285 - 317

[1]

@@ Zeile 8: / Zeile 8: @@
 == Herleitung der Bewegungsgleichungen ==
-Wenn <math>l_1</math> und <math>l_2</math> die Längen der (masselosen) Verbindungsstangen, <math>m_1</math> und <math>m_2</math> die Pendelmassen, <math>\theta_1, \theta_2 </math> die Auslenkung vom Lot und <math>g</math> die Erdbeschleunigung bezeichnet, dann findet man für die Positionen von <math>m_1 </math> und <math>m_2 </math>:
+Wenn <math>l_1</math> und <math>l_2</math> die Längen der (masselosen) Verbindungsstangen, <math>m_1</math> und <math>m_2</math> die Pendelmassen, <math>\theta_1, \theta_2 </math> die Auslenkung vom Lot und <math>g</math> die [[Erdbeschleunigung]] bezeichnet, dann findet man für die Positionen von <math>m_1 </math> und <math>m_2 </math>:
 :<math>\begin{align}
@@ Zeile 19: / Zeile 19: @@
 :<math>\begin{align}
 x_2 &= l_1 \sin(\theta_1) + l_2 \sin(\theta_2)\
-y_2 &= -l_1 \cos(\theta_1) -l_2 \cos(\theta_2)
+y_2 &= -l_1 \cos(\theta_1) -l_2 \cos(\theta_2).
 \end{align}</math>
@@ Zeile 25: / Zeile 25: @@
 :<math>\begin{align}
-u_1 &= \frac{\partial x_1}{\partial t} = \dot{\theta_1} l_1 \cos(\theta_1)\
+u_1 &= \frac{\partial x_1}{\partial t} = \dot{\theta}_1 l_1 \cos(\theta_1)\
-v_1 &= \frac{\partial y_1}{\partial t} = \dot{\theta_1} l_1 \sin(\theta_1)
+v_1 &= \frac{\partial y_1}{\partial t} = \dot{\theta}_1 l_1 \sin(\theta_1)
 \end{align}</math>
@@ Zeile 32: / Zeile 32: @@
 :<math>\begin{align}
-u_2 &= \frac{\partial x_2}{\partial t} + u_1 = \dot{\theta_2} l_2 \cos(\theta_2) + u_1\
+u_2 &= \frac{\partial x_2}{\partial t}  = \dot{\theta}_2 l_2 \cos(\theta_2) + u_1\
-v_2 &= \frac{\partial y_2}{\partial t} + v_1 = \dot{\theta_2} l_2 \sin(\theta_2) + v_1
+v_2 &= \frac{\partial y_2}{\partial t}  = \dot{\theta}_2 l_2 \sin(\theta_2) + v_1.
 \end{align}</math>
@@ Zeile 48: / Zeile 48: @@
 :<math>\begin{align}
-T &= T_1 + T_2 = \frac{1}{2} m_1\dot{\theta_1^2} l_1^2+ \frac{1}{2} \left(m_2 \dot{\theta_2^2} l_2^2 + m_2 \dot{\theta_1^2} l_1^2 + 2 m_2 \dot{\theta_1} l_1 \dot{\theta_2} l_2 \cos(\theta_1 - \theta_2)\right)\
+T &= T_1 + T_2 = \frac{1}{2} m_1\dot{\theta}_1^2 l_1^2+ \frac{1}{2} \left(m_2 \dot{\theta}_2^2 l_2^2 + m_2 \dot{\theta}_1^2 l_1^2 + 2 m_2 \dot{\theta}_1 l_1 \dot{\theta}_2 l_2 \cos(\theta_1 - \theta_2)\right)\
-V &= V_1 + V_2 = -m_1 g l_1 \cos(\theta_1) -m_2 g l_1 \cos(\theta_1) - m_2 g l_2 \cos(\theta_2)
+V &= V_1 + V_2 = -m_1 g l_1 \cos(\theta_1) -m_2 g l_1 \cos(\theta_1) - m_2 g l_2 \cos(\theta_2).
 \end{align}</math>
 Damit ergibt sich für die Lagrange-Funktion dann
-:<math>{L = \frac{1}{2} (m_1 + m_2) \dot{\theta_1}^2 l_1^2 + \frac{1}{2} m_2 \dot{\theta_2}^2 l_2^2 + m_2 \dot{\theta_1} l_1 \dot{\theta_2} l_2 \cos(\theta_1 - \theta_2) + (m_1 + m_2) g l_1 \cos(\theta_1) + m_2 g l_2 \cos(\theta_2)}</math>
+:<math>{L = \frac{1}{2} (m_1 + m_2) \dot{\theta}_1^2 l_1^2 + \frac{1}{2} m_2 \dot{\theta}_2^2 l_2^2 + m_2 \dot{\theta}_1 l_1 \dot{\theta}_2 l_2 \cos(\theta_1 - \theta_2) + (m_1 + m_2) g l_1 \cos(\theta_1) + m_2 g l_2 \cos(\theta_2)}.</math>
 Unter Verwendung der [[Variationsrechnung|Euler-Lagrange-Gleichung]]
@@ Zeile 63: / Zeile 63: @@
 :<math>\begin{align}
-\ddot{\theta_1} &= -\frac{m_2}{m_1 + m_2} \frac{l_2}{l_1} (\ddot{\theta_2} \cos(\theta_1 - \theta_2) + \dot{\theta_2^2} \sin(\theta_1 - \theta_2)) - \frac{g}{l_1} \sin(\theta_1)\
+\ddot{\theta}_1 &= -\frac{m_2}{m_1 + m_2} \frac{l_2}{l_1} (\ddot{\theta}_2 \cos(\theta_1 - \theta_2) + \dot{\theta}_2^2 \sin(\theta_1 - \theta_2)) - \frac{g}{l_1} \sin(\theta_1)\
-\ddot{\theta_2} &= - \frac{l_1}{l_2} (\ddot{\theta_1} \cos(\theta_1 - \theta_2) - \dot{\theta_1^2} \sin(\theta_1 - \theta_2)) - \frac{g}{l_2} \sin(\theta_2)\end{align}</math>
+\ddot{\theta}_2 &= - \frac{l_1}{l_2} (\ddot{\theta}_1 \cos(\theta_1 - \theta_2) - \dot{\theta}_1^2 \sin(\theta_1 - \theta_2)) - \frac{g}{l_2} \sin(\theta_2)\end{align}</math>
-Das sind die Winkelbeschleunigungen für <math>\theta_1</math> und <math>\theta_2</math> und beschreiben die Evolution des Pendels.
+die Winkelbeschleunigungen für <math>\theta_1</math> und <math>\theta_2</math>, welche die Evolution des Pendels beschreiben.
 In den Bewegungsgleichungen treten Winkelfunktionen <math>(\sin, \cos )</math> der Zustandsgrößen und auch Ableitungen auf. Es handelt sich also um ein nichtlineares System. Im Spezialfall kleiner Auslenkungen als Anfangsbedingungen <math>\theta_1(t=0), \theta_2(t=0) \approx 0 </math>, lassen sich die Bewegungsgleichungen allerdings mittels der [[Kleinwinkelnäherung]] vereinfachen. Dann lassen sich beispielsweise weitere Spezialfälle wie <math>m_1 \ll m_2 </math> oder <math>m_1 \gg m_2</math> mit analytischen Ansätzen betrachten, die eine näherungsweise harmonische Lösung haben; diese kann auch analytisch bestimmt werden.
@@ Zeile 72: / Zeile 72: @@
 == Lösung der Bewegungsgleichungen ==
-Die Bewegungsgleichungen für die [[generalisierte Koordinate|generalisierten Koordinaten]] <math display="inline">{\theta_{1}}</math> und <math display="inline">{\theta_{2}}</math> stellen ein nichtlineares System von zwei gekoppelten [[Differentialgleichungen]] dar, welches analytisch nicht lösbar ist. Es kann bei vier bekannten [[Anfangswertproblem|Anfangswerten]] (<math>\theta_1, \theta_2, \dot{\theta_1}, \dot{\theta_2} </math>) mit [[Liste_numerischer_Verfahren#Numerik_gewöhnlicher_Differentialgleichungen|numerischen Verfahren]] gelöst werden. Hierbei werden also die anfänglichen Auslenkungen (z.&nbsp;B. ''30°'' und ''30°'') und die anfänglichen [[Winkelgeschwindigkeit|Geschwindigkeiten]] (z.&nbsp;B. <math display="inline">0 \mathrm\frac{rad}{s}</math> und <math display="inline">0 \mathrm\frac{rad}{s}</math>) eingegeben und damit dann die Evolution des Pendels berechnet.
+Die Bewegungsgleichungen für die [[generalisierte Koordinate|generalisierten Koordinaten]] <math display="inline">{\theta_{1}}</math> und <math display="inline">{\theta_{2}}</math> stellen ein nichtlineares System von zwei gekoppelten [[Differentialgleichungen]] dar, welches analytisch nicht lösbar ist. Es kann bei vier bekannten [[Anfangswertproblem|Anfangswerten]] (<math>\theta_1, \theta_2, \dot{\theta}_1, \dot{\theta}_2 </math>) mit [[Liste_numerischer_Verfahren#Numerik_gewöhnlicher_Differentialgleichungen|numerischen Verfahren]] gelöst werden. Hierbei werden also die anfänglichen Auslenkungen (z.&nbsp;B. ''30°'' und ''30°'') und die anfänglichen [[Winkelgeschwindigkeit|Geschwindigkeiten]] (z.&nbsp;B. <math display="inline">0 \mathrm\frac{rad}{s}</math> und <math display="inline">0 \mathrm\frac{rad}{s}</math>) eingegeben und damit dann die Evolution des Pendels berechnet.
 Mittels [[Trigonometrie]] können die Winkel <math>\theta_{1}</math> und <math>\theta_{2}</math> in die kartesischen Koordinaten <math>(x_1, y_1, x_2, y_2)</math> der Massenpunkte überführt werden.
 == Anwendungen ==
-Eine [[Kirchenglocke]] mit Klöppel bildet ein Doppelpendel.
+* Die [[Blide]] nutzt den energetischen Austausch zwischen den Massen der Pendel, um ein kleines Gewicht am äußeren Pendel durch die potentielle Energie eines großen Gewichtes am inneren Pendel zu beschleunigen.
+* Eine [[Kirchenglocke]] mit Klöppel bildet ein Doppelpendel, allerdings mit zusätzlicher Beschränkung für den Winkel <math>\theta_2</math>.
 == Auswertung des chaotischen Verhaltens ==
@@ Zeile 83: / Zeile 84: @@
 === Maximaler Ljapunow-Exponent ===
-Der MLE ist der sog. ''maximaler [[Ljapunow-Exponent]]'' (''maximum Lyapunov exponent'') und beschreibt die „Stärke“ des chaotischen Verhaltens. Er ist Bestandteil des Ljapunow-Spektrums welches alle Ljapunow-Exponenten (je einer pro Freiheitsgrad) beinhaltet. Man geht davon aus, dass das System eine Störung in der Richtung des MLE hat und da er das größte Wachstum zeigt, ist zu erwarten, dass der MLE nach einer gewissen Zeit die Evolution des Systems dominiert. Ein positiver MLE zeigt normalerweise ein chaotisches System an. Er wird berechnet mit:
+Der MLE ist der sog. ''maximale [[Ljapunow-Exponent]]'' (''maximum Lyapunov exponent'') und beschreibt die „Stärke“ des chaotischen Verhaltens. Er ist Bestandteil des Ljapunow-Spektrums welches alle Ljapunow-Exponenten (je einer pro Freiheitsgrad) beinhaltet. Man geht davon aus, dass das System eine Störung in der Richtung des MLE hat und da er das größte Wachstum zeigt, ist zu erwarten, dass der MLE nach einer gewissen Zeit die Evolution des Systems dominiert. Ein positiver MLE zeigt normalerweise ein chaotisches System an. Er wird berechnet mit:
 :<math>\lambda_\mathrm{max}(t) = \frac{1}{t} \ln \left(\frac{\|\delta (t)\|}{\| \delta (t=0)\|} \right)</math>
@@ Zeile 89: / Zeile 90: @@
 Bei zwei Experimenten mit einer anfänglichen Separation von <math display="inline">\delta(t=0) \backsim \mathcal{O}(10^{-8} )</math> in den Anfangsbedingungen oder sogar weniger, verstärkt sich diese Differenz exponentiell und lässt die Trajektorien divergieren.<ref>Wolf, A., Swift, J. B., Swinney, H. L., & Vastano, J. A. (1985). Determining lyapunov exponents from a time series. Physica D: Nonlinear Phenomena, Vol. 16 No. 3: 285 - 317</ref> Die Separation (der natürliche Logarithmus der obigen Gleichung) kann dann in einem halblogarithmischen Diagramm gegen die Zeit aufgetragen werden. Dann wird mittels linearer Regression die Steigung bestimmt und diese gibt dann den approximierten MLE.
-=== Bifurkationsdiagram ===
+=== Bifurkationsdiagramm ===
 [[Datei:Bifurkationsdiagram Doppelpendel.png|miniatur|Darstellung der Bifurkation eines Doppelpendels mit gleichen Längen, gleichen Massen, <math display="inline">g=9{,}81\mathrm\frac{m}{s^2}</math> für <math>\theta_1 = \theta_2</math> bei <math>\dot{\theta_2} = 0</math>. Die rote Linie zeigt die Bifurkation des Doppelpendels von harmonischen Oszillationen zu chaotischen Oszillationen.]]
-[[Bifurkation (Mathematik)|Bifurkationsdiagramme]] sind eine Möglichkeit komplexe Informationen über den [[Phasenraum]] eines dynamischen Systems in einen zweidimensionalen, visualisierbaren Plot zu komprimieren. Üblicherweise wird die qualitative Änderungen des Verhaltens eines Systems mittels der Variation eines geeigneten Parameters untersucht. So können für das Doppelpendel bspw. das Verhältnis der Massen, das Verhältnis der Längen, die Erdbeschleunigung oder die Anfangsbedingungen herangezogen werden. Durch die kontinuierliche Veränderung des gewählten Bifurkationsparameters wird das System auf Stabilität (periodische, quasi-periodische Lösungen) bzw. auf Chaos geprüft.
+[[Bifurkation (Mathematik)|Bifurkationsdiagramme]] sind eine Möglichkeit komplexe Informationen über den [[Phasenraum]] eines dynamischen Systems in einen zweidimensionalen, visualisierbaren Plot zu komprimieren. Üblicherweise wird die qualitative Änderungen des Verhaltens eines Systems mittels der Variation eines geeigneten Parameters untersucht. So können für das Doppelpendel bspw. das Verhältnis der Massen, das Verhältnis der Längen, die Erdbeschleunigung oder die Anfangsbedingungen herangezogen werden. Durch die kontinuierliche Veränderung des gewählten Bifurkationsparameters wird das System auf Stabilität (periodische, quasi-periodische Lösungen) bzw. auf [[Chaosforschung|Chaos]] geprüft.
 Wenn man die anfänglichen Winkel <math display="inline">\theta_1 = \theta_2  </math> als den Bifurkationsparameter wählt, lässt sich das qualitativ veränderliche Verhalten des Doppelpendels sehr gut veranschaulichen. Dazu werden die beiden Winkel simultan Stück für Stück erhöht und für jedes Inkrement wird das Doppelpendel erneut integriert (berechnet). Mit diesen Daten lässt sich dann veranschaulichen, wie das System schwingt. Man hat also einen vier-dimensionalen Phasenraum, der sich aus <math display="inline">\theta_1, \theta_2, \dot{\theta_1}, \dot{\theta_2}</math> zusammensetzt. Praktischerweise oszillieren die Winkelgeschwindigkeiten <math display="inline">\dot{\theta_1}, \dot{\theta_2}</math>, obwohl mit unbestimmter Amplitude, um Null. Daher ist zu erwarten, dass beide immer wieder die Null überqueren. Für ein harmonisch schwingendes System (periodische Lösung) sind die Nullüberquerungen von <math>\dot{\theta_1}, \dot{\theta_2}</math> an festen Punkten, da das System immer an bestimmten Punkten (<math display="inline">\theta_1, \theta_2</math>) seine Auf- und Abwärtsbewegung beendet und zurückschwingt. Das ist vergleichbar mit einem [[Pendel|normalen starren Pendel]]. Daher ist im Umkehrschluss zu erwarten, dass das chaotisch schwingende System an allen möglichen Punkten (<math display="inline">\theta_1, \theta_2</math>) die Winkelgeschwindigkeit Null zeigt. Wenn man dann eine "Scheibe" aus dem Phasenraum gesondert betrachtet, bspw. Winkelgeschwindigkeit <math display="inline">\dot{\theta_2} = 0</math>, kann man die Bifurkation des Verhaltens zwei-dimensional darstellen, indem man die Winkelgeschwindigkeit <math display="inline">\dot{\theta_1}</math> gegen die veränderliche Anfangsbedingung aufträgt (siehe rechts).