Mathematica 12 - WhenEvent - Warum wird der Vorzeichenwechsel nicht angenommen

• Erstellt am 26.4.2021
• Überarbeitet am 26.8.2024
• Software: Mathematica 12

Fordern Sie hier ein Angebot für die aktuelle Version von Wolfram Mathematica an.

Als ein grundlegendes Beispiel für die Funktion WhenEvent finden Sie auf der Hilfeseite von Wolfram zu dieser Funktion die Simuliation eines springenden Balls, der bei jedem Sprung 95% seiner Geschwindigkeit beibehält. Wenn Sie das Beispiel so abändern, dass der Ball nur 80% seiner Geschwindigkeit beibehält, wird ein Graph angezeigt, in dem der die Höhe 0 dennoch unterschritten wird, obwohl hier die Geschwindigkeit gemäß dem Befehl das Vorzeichen wechseln müsste.

Was ist der Grund hierfür?

Erläuterung

Die Differentialgleichung

$\mathrm{y}\text{'}\text{'}\left(\mathrm{t}\right)=-\mathrm{g}$

mit $\mathrm{g}>0$ hat als allgemeine Lösung

$\mathrm{y}\left(\mathrm{t}\right)=-\frac{\mathrm{g}}{2}*{\mathrm{t}}^2+\mathrm{a}*\mathrm{t}+\mathrm{b}$

mit zwei Konstanten a und b. Die Nullstellen dieser Funktionen werden nach dem Umformen der Gleichung

$-\frac{\mathrm{g}}{2}*{\mathrm{t}}^2+\mathrm{a}*\mathrm{t}+\mathrm{b}=0$

zu

${\mathrm{t}}^2-\frac{2*\mathrm{a}}{\mathrm{g}}*\mathrm{t}-\frac{2*\mathrm{b}}{\mathrm{g}}=0$

gemäß der p-q-Formel berechnet als

${\mathrm{t}}_1=\frac{\mathrm{a}}{\mathrm{g}}-\sqrt{\frac{{\mathrm{a}}^2}{{\mathrm{g}}^2}+\frac{2*\mathrm{b}}{\mathrm{g}}}={\mathrm{t}}_{\mathrm{mitte}}-\frac{1}{2}{\mathrm{t}}_{\mathrm{abstand}}$

${\mathrm{t}}_2=\frac{\mathrm{a}}{\mathrm{g}}+\sqrt{\frac{{\mathrm{a}}^2}{{\mathrm{g}}^2}+\frac{2*\mathrm{b}}{\mathrm{g}}}={\mathrm{t}}_{\mathrm{mitte}}+\frac{1}{2}{\mathrm{t}}_{\mathrm{abstand}}$

Wenn die beiden Nullstellen reelle Zahlen sind, ist ihr Abstand gleich

${\mathrm{t}}_2-{\mathrm{t}}_1=2*\sqrt{\frac{{\mathrm{a}}^2}{{\mathrm{g}}^2}+\frac{2*\mathrm{b}}{\mathrm{g}}}={\mathrm{t}}_{\mathrm{abstand}}$

Wir können die Gleichungen der Funktion und Ableitung auch schreiben als

$\mathrm{y}\left(\mathrm{t}\right)=-\mathrm{g}*(\mathrm{t}-{\mathrm{t}}_1)*(\mathrm{t}-{\mathrm{t}}_2)$

Des Weiteren gilt

und damit

$\mathrm{y}\text{'}\left(0\right)=\mathrm{a}$

Die Nullstelle der ersten Ableitung liegt bei $\frac{\mathrm{a}}{\mathrm{g}}={\mathrm{t}}_{\mathrm{mitte}}$. Die ersten Ableitungen an den Nullstellen sind

$\mathrm{y}\text{'}\left({\mathrm{t}}_1\right)=-\mathrm{g}*({\mathrm{t}}_{\mathrm{mitte}}-\frac{1}{2}*{\mathrm{t}}_{\mathrm{abstand}})+\mathrm{a}=-\mathrm{g}*{\mathrm{t}}_{\mathrm{mitte}}+\mathrm{a}+\mathrm{g}*\frac{1}{2}*{\mathrm{t}}_{\mathrm{abstand}}=\mathrm{y}\text{'}\left({\mathrm{t}}_{\mathrm{mitte}}\right)+\mathrm{g}*\frac{1}{2}*{\mathrm{t}}_{\mathrm{abstand}}=\mathrm{g}*\frac{1}{2}*{\mathrm{t}}_{\mathrm{abstand}}$

$\mathrm{y}\text{'}\left({\mathrm{t}}_2\right)=-\mathrm{g}*({\mathrm{t}}_{\mathrm{mitte}}+\frac{1}{2}*{\mathrm{t}}_{\mathrm{abstand}})+\mathrm{a}=-\mathrm{g}*{\mathrm{t}}_{\mathrm{mitte}}+\mathrm{a}-\mathrm{g}*\frac{1}{2}*{\mathrm{t}}_{\mathrm{abstand}}=\mathrm{y}\text{'}\left({\mathrm{t}}_{\mathrm{mitte}}\right)-\mathrm{g}*\frac{1}{2}*{\mathrm{t}}_{\mathrm{abstand}}=-\mathrm{g}*\frac{1}{2}*{\mathrm{t}}_{\mathrm{abstand}}=-\mathrm{y}\text{'}\left({\mathrm{t}}_1\right)$

Jede der beiden Ableitungen sind proportional zu dem Abstand zwischen den Nullstellen:

$\frac{2}{\mathrm{g}}*\mathrm{y}\text{'}\left({\mathrm{t}}_1\right)={\mathrm{t}}_{\mathrm{abstand}}=-\frac{2}{\mathrm{g}}*\mathrm{y}\text{'}\left({\mathrm{t}}_2\right)$

Für den Abstand zwischen den beiden Nullstellen zweier durch einen Übergang ${\mathrm{y}}_{\mathrm{neu}}\text{'}\left(\mathrm{z}\right)=-\mathrm{q}*{\mathrm{y}}_{\mathrm{alt}}\text{'}\left(\mathrm{z}\right)$ mit $0<q<1$ folgt

${\mathrm{t}}_{\mathrm{abstand},\mathrm{neu}}=\mathrm{q}*{\mathrm{t}}_{\mathrm{abstand},\mathrm{alt}}$

Damit gilt

${\mathrm{t}}_{\mathrm{abstand},2}=\mathrm{q}*{\mathrm{t}}_{\mathrm{abstand},1}$

${\mathrm{t}}_{\mathrm{abstand},3}=\mathrm{q}*{\mathrm{t}}_{\mathrm{abstand},2}={\mathrm{q}}^2*{\mathrm{t}}_{\mathrm{abstand},1}$

${\mathrm{t}}_{\mathrm{abstand},4}=\mathrm{q}*{\mathrm{t}}_{\mathrm{abstand},3}={\mathrm{q}}^3*{\mathrm{t}}_{\mathrm{abstand},1}$

und allgemein

${\mathrm{t}}_{\mathrm{abstand},\mathrm{i}}={\mathrm{q}}^{\mathrm{i}-1}*{\mathrm{t}}_{\mathrm{abstand},1}$

Hier ist ${\mathrm{t}}_{\mathrm{abstand},\mathrm{i}}$ der Abstand zwischen den Nullstellen der $\mathrm{i}$-ten komplett innerhalb des angegebenen Definitionsbereiches liegenden Teilkurve. Die Summe aller dieser Abstände ist

$\sum _{\mathrm{i}=0}^{\infty }{\mathrm{q}}^{\mathrm{i}}*{\mathrm{t}}_{\mathrm{abstand},1}=\frac{1}{1-\mathrm{q}}*{\mathrm{t}}_{\mathrm{abstand},1}$

weil es sich um eine geometrische Reihe handelt.

Wegen ${\mathrm{y}}^{\text{'}}\left(0\right)=0$ liegt die allererste Kurve noch nicht komplett im angegebenen Definitionsbereich, sondern genau zur Hälfte, denn Ihre Mitte liegt bei $\mathrm{t}=0$. Die Gesamtlänge aller Teilkurven ab dem Minimum des angegebenen Definitionsbereichs ist

$\frac{1}{2}*{\mathrm{t}}_{\mathrm{abstand},0}+\frac{1}{1-\mathrm{q}}*{\mathrm{t}}_{\mathrm{abstand},1}=\frac{1}{2}*{\mathrm{t}}_{\mathrm{abstand},0}*{\mathrm{t}}_{\mathrm{abstand},0}+\frac{1}{1-\mathrm{q}}*\mathrm{q}*{\mathrm{t}}_{\mathrm{abstand},0}={\mathrm{t}}_{\mathrm{abstand},0}*(\frac{1}{2}+\frac{\mathrm{q}}{1-\mathrm{q}})$

Für $\mathrm{g}=9.81$ und $\mathrm{q}=0.8$ ist die Gleichung der ersten Lösungskurve

$y\left(\mathrm{t}\right)=-4.905*{\mathrm{t}}^2+5$

mit den Nullstellen

${\mathrm{t}}_1=-1.0096375546923047$

${\mathrm{t}}_2=1.0096375546923047$

und

${\mathrm{t}}_{\mathrm{abstand},0}=2.0192751093846093$

sodass die Gesamtbreite innerhalb des angegebenen Definitionsbereichs durch

$2.0192751093846093*(\frac{1}{2}+\frac{0.8}{1-0.8})=2.0192751093846093*4.5=9.08674$

berechnet wird. Das bedeutet, dass alle der unendlich vielen Ereignisse alle vor dem Wert $9.08674$ passieren, die unendlich vielen Teilkurven, aus denen sich die Lösung der Differentialgleichung zusammensetzt, in der Summe nur eine Breite von $\mathrm{9.08674<10}$ haben. In jeder Annäherung der Lösung durch endlich viele Teilkurven verletzt die letzte der Teilkurven das Ereignis.

Wenn Sie bei Ihre Differentialgleichung mit WhenEvent ein ähnliches Phänomen beobachten, überprüfen Sie, ob es sich analog mit der geometrischen Reihe erklären lässt.