Das Wachstum des Fehlers bei der Approximation einer Differentialgleichung

Beim Lösen der Differentialgleichung θ ¨ + ω 2 Sünde θ = 0 , mein Lehrbuch nähert sich dem an Sünde θ = θ und löst die Näherungsgleichung θ ¨ + ω 2 θ = 0 .

Aber woher weiß ich, wie schnell der Fehler zwischen der realen Lösung und der angenäherten Lösung wächst?

Der Fehler hängt mit \theta zusammen. Wenn also in Ihrer Näherungsgleichung die Lösung beschränkt ist, ist dies auch der Fehler. In diesem Fall wird es im Laufe der Zeit nicht explodieren. Bei einigen anderen Differentialgleichungen ist es jedoch viel schwieriger zu sagen, ob der Fehler bei einer Annäherung explodiert oder nicht.
Dies ist eine sehr weit gefasste Frage (meiner Meinung nach eher für Mathematik SE geeignet), hängt jedoch stark von der Gleichung, den Randbedingungen und der numerischen Implementierung des Algorithmus ab.

Antworten (1)

Die exakte Lösung kann durch elliptische Integrale ausgedrückt werden , und man kann in diesem Fall tatsächlich das Wachstum des Fehlers berechnen. Als Formel allein wird es nicht besonders aufschlussreich sein. Vielleicht kann eine bessere Schätzung gefunden werden, indem die Differenz in der Schätzung des Zeitraums zwischen der tatsächlichen und der angenäherten Lösung aufgetragen wird . Wir können sehen, dass der Periodenfehler ungefähr quadratisch mit dem Anfangswinkel wächst.

Dies sagt uns nicht direkt die Wachstumsrate des Fehlers im Laufe der Zeit, aber man könnte argumentieren, dass, wenn der obige Fehler ist E ( θ 0 ) = | θ ^ ( T ( θ 0 ) ) θ ( T ( θ 0 ) ) | für eine wahre Zeit T ( θ 0 ) dann sollte die Wachstumsrate ungefähr sein E ' ( T ) E ( θ 0 ) / T ( θ 0 ) . Dies wird eine ungefähr quadratische Funktion sein ( E ( θ 0 ) ) geteilt durch eine weitere ungefähr lineare Funktion für kleine Winkel ( T ( θ 0 ) ), so dass die Wachstumsrate im Laufe der Zeit für kleine Winkel ungefähr linear mit dem Anfangswinkel zunimmt.

Mit zunehmendem Winkel wird der Fehler immer schlimmer, zumal die tatsächliche Flugbahn bei Annäherung anders aussieht als eine Sinuskurvenschwingung θ π und die Periode geht gegen unendlich.