Mo Logo [Home] [Lexikon] [Aufgaben] [Tests] [Kurse] [Begleitmaterial] [Hinweise] [Mitwirkende] [Publikationen]
Mathematik-Online-Lexikon:

Trapez-Regel für eine Differentialgleichung

A B C D E F G H I J K L M N O P Q R S T U V W X Y Z Übersicht

Integriert man die Lösung $ (u_1(t),\ldots,u_d(t))^t$ eines Differentialgleichungssystems

$\displaystyle u^\prime = f(t,u) \,, $

so folgt

$\displaystyle u(t+h) = u(t) + \int_t^{t+h}f(s,u(s))\,ds \,. $

Das Integral kann durch die Trapezregel angenähert werden. Ein Schritt $ u(t)\approx v \to w \approx u(t+h)$ des resultierenden Verfahrens hat dann die Form

$\displaystyle w = v + \frac{h}{2}(f(t,v)+f(t+h,w)) \,. $

Die Trapezregel ist ein implizites Einschrittverfahren der Ordnung $ 2$.

Zur Durchführung eines Verfahrensschritts $ v\to w$ muss im allgemeinen ein nichtlineares Gleichungssystem gelöst werden. Näherungsweise kann dies durch einige wenige Schritte einer Fixpunktiteration

$\displaystyle w^{\mathrm{neu}} = v + \frac{h}{2}(f(t,v)+f(t+h,w^{\mathrm{alt}})) $

mit Startwert $ w=v$ geschehen.

Für ein lineares Differentialgleichungssystem

$\displaystyle u^\prime = A(t)u + b(t) $

ist die Trapezregel einfacher realisierbar. Ein Schritt erfordert lediglich die Lösung des linearen Gleichungssystems

$\displaystyle \left(E-\frac{h}{2} A(t+h)\right) w = \left(E+\frac{h}{2} A(t)\right) v + \frac{h}{2}(b(t) + b(t+h)) $

mit $ E$ der $ d\times d$-Einheitsmatrix.

Beispiel:

[Erläuterungen] [Verweise]
  automatisch erstellt am 19.  8. 2013