Voraussetzungen

  • Runge-Kutta-Verfahren zur numerischen Integration von gewöhnlichen Differentialgleichungen

Lerninhalte

  • Lösen von mehrdimensionalen gewöhnlichen Differentialgleichungen

Systeme von ODEs#

Aufgabe 1: Schwingung eines Pendels#

Die Schwingungsgleichung eines Pendels lernen Sie auch in dem Kapitel Harmonische Schwingung kennen. Da lösen Sie sie mit ode45.

Nun gehen Sie einen Schritt weiter:

  • Stellen Sie die Schwingungsgleichung für sowohl mit als auch ohne Kleinwinkelnäherung auf.

  • Formen Sie die resultierende ODE zweiter Ordnung in ein System von zwei ODEs um.

  • Implementieren Sie nun eines Ihrer Verfahren aus den Kapitel Integration von ODEs mit expliziten Verfahren und Integration von ODEs mit impliziten Verfahren zum Lösen der von Ihnen hergeleiteten Differentialgleichungen. Dazu muss Ihr Code auch mit Systemen von ODEs umgehen können.

%your code here

Aufgabe 2: Die van-der-Pol Gleichung#

Die van-der-Pol Gleichung lautet

\[\ddot{x} - \varepsilon(1-x^2)\dot{x} + x = 0.\]

Analog zu dem Fadenpendel kann \(x\) als ein Winkel beziehungsweise eine Auslenkung interpretiert werden. Die van-der-Pol Gleichung stellt ein besonderes Problem für numerische Verfahren dar. Testen Sie anhand dieser Gleichung Ihre Lösungsverfahren und plotten Sie Ihr Ergebnis. Nehmen Sie auch die Dauern Ihrer Verfahren auf. Ihre Lösung könnte beispielsweise so aussehen:

../../_images/vdP_comparison_1.png
../../_images/vdP_comparison_2.png

Fig. 13 Schwingung und Phasenraum verschiedener Lösungen der van-der-Pol Gleichung. Das explizite Eulerverfahren (rote Punkte) zeigt eine deutliche höhere Amplitude und verschobene Periode. Die Lösung mit einem Runge-Kutta-Verfahren dritter Ordnung (gelbe Striche) liegt sehr nah bei der sehr genauen Lösung von ode45 (blaue Linie).#

%your code here

Aufgabe 3: Die inhomogene van-der-Pol Gleichung#

Ein schwingfähiges System kann periodisch angeregt werden. Das ist besonders interessant für Resonanzfälle, die gewünscht oder unerwünscht sein können. Für eine periodisch angeregte van-der-Pol Schwingung lautet die Gleichung

\[\ddot{x} - \varepsilon(1-x^2)\dot{x} + x = A \cos(\omega t + \varphi_0).\]