Numerisches Lösen von Differentialgleichungen 6. Ordnung

Ein sich in drei Dimensionen bewegendes Objekt wird mit 6 Variablen beschrieben: $x$, $y$, $z$, $v_x$, $v_y$, and $v_z$. Ist die Kraft aus das Objekt bekannt, kann die Bewegung mithilfe von 6 Differentialgleichungen erster Ordnung beschrieben werden:

$\large \frac{dx}{dt}=v_x$   $\large \frac{dv_x}{dt}=F_x(x,y,z,v_x,v_y,v_z,t)/m$

$\large \frac{dy}{dt}=v_y$   $\large \frac{dv_y}{dt}=F_y(x,y,z,v_x,v_y,v_z,t)/m$

$\large \frac{dz}{dt}=v_z$   $\large \frac{dv_z}{dt}=F_z(x,y,z,v_x,v_y,v_z,t)/m$

Dabei sind $F_x$, $F_y$, and $F_z$ die drei Komponenten der Kraft, $m$ die Masse und $t$ die Zeit. Das folgende Formular kann zum numerischen Integrieren dieser Gleichungen für eine Gesamtzahl von $N_{steps}$ Schritten der Schrittweite $\Delta t$ benutzt werden. Die Kraft kann als Funktion von $x$, $y$, $z$, $v_x$, $v_y$, $v_z$ und $t$ angegeben werden.

 Numerisches Lösen von Differentialgleichungen 6. Ordnung 

$ \large \frac{dx}{dt}=$

$v_x$

$ \large \frac{dv_x}{dt}=$

$ \large \frac{dy}{dt}=$

$v_y$

$ \large \frac{dv_y}{dt}=$

$ \large \frac{dz}{dt}=$

$v_z$

$ \large \frac{dv_z}{dt}=$

Anfangsbedingungen: 

$t_0=$

 

$\Delta t=$

$x(t_0)=$

 

$N_{steps}$

$v_x(t_0)=$

 

Graphische Darstellung:

vs.

$y(t_0)=$

 

$v_y(t_0)=$

 

$z(t_0)=$

 

$v_z(t_0)=$

 

 

 $t$   $x$   $v_x$   $y$   $v_y$   $z$   $v_z$

Die Idee der numerischen Integration ist, dass für bekannte Anfangsbedingungen $x(t_0)$ und $v_x(t_0)$ eine gute Schätzung für $x$ und $v_x$ eine kurze Zeit $\Delta t$ später ist:

$\large x(t_0+\Delta t) \approx \frac{dv_x}{dt}|_{t_0}\Delta t$ and $\large v_x(t_0+\Delta t) \approx F(x(t_0),v_x(t_0),t_0)\Delta t/m.$

Sobald Schätzungen für Position und Geschwindigkeit des Objekts zur Zeit $t_0 + \Delta t$ berechnet wurden, können diese benutzt werden, um die Position und Geschwindigkeit zur Zeit $t_0 + 2\Delta t$ zu bestimmen.

Im obigen Formular kann die Kraft mithilfe der mathematischen Standardfunktionen abs(x), acos(x), asin(x), atan(x), cos(x), exp(x), pi = 3.141592653589793, pow(x,y) = xy, round(x), sin(x), sqrt(x), and tan(x) angegeben werden. Beachten Sie, daß die Multiplikation mit einem '*' Symbol angegeben werden muss, also 3*cos(x) statt 3cos(x). Potenzen werden mit der 'pow' Funktion spezifiert: x² ist pow(x,2) statt x^2.

Die numerische Integration wird instabil, wenn der Zeitschritt zu lang wird. Der Zeitschritt sollte einen Faktor 10 bis 100 mal kleiner sein als jede berechnete, charakteristische Zeit der Bewegung. Wird der Orbit der Erde um die Sonne berechnet, dann ist ein Zeitschritt von 3 Tagen (ca. 4E6 Sekunden) angemessen. Wird die Bewegung eines Elektrons vorbei an einem Ion berechnet, ist ein Zeitschritt von 1 ps sinnvoll. Wird die Routine instabil, wird nichts graphisch dargestellt. In diesem Fall sollte die Berechnung mit einem kürzeren Zeitschritt versucht werden.

Es gibt viele Routinen, um eine numerische Integration vorzunehmen. Die oben beschriebene Methode nennt sich die Euler Methode. Konkret jedoch benutzt diese APP die genauere Runge-Kutta Methode vierter Ordnung mit fixer Schrittweite. Eine Routinen zur numerischen Integration berechnen den optimalen Zeitschritt automatisch.