Trapez-Methode

Trapez-Methode

Das implizite Trapez-Verfahren ist ein Verfahren zur numerischen Lösung eines Anfangswert-Problems

y'(t)  = f\left(t, y(t)\right), \quad
       y(t=0) = y_0

Es lässt sich sowohl den Runge-Kutta-Verfahren als auch den Adams-Moulton-Verfahren zuordnen. Das Trapezverfahren ist A-stabil mit der Besonderheit, dass für die Schwingungslösung y'=iαiy kein Amplitudenfehler auftritt[1]. Das Verfahren lässt sich aus der Trapezregel herleiten:

y_{n+1}=y_n + \frac{h}{2}(f_{n+1}+f_n)

mit

f_n \ := f(t_n,y_n).

Inhaltsverzeichnis

Lösungsmethode

Zur Lösung dieses, in der Regel nichtlinearen Gleichungssystems können verschiedene numerische Verfahren zur Lösung nichtlinearer Gleichungssysteme genutzt werden. Für das quadratisch konvergente Newton-Verfahren ergibt sich konkret:

y^{(k + 1)}_{n+1} =
y^{(k)}_{n+1} - \left( y^{(k)}_{n+1} - y_n - \frac{h}{2}(f^{(k)}_{n+1} + f_n) \right)
\left(I - \frac{h}{2}\frac{\partial f^{(k)}_{n+1}}{\partial y^{(k)}_{n+1}}\right)^{-1}.

Man erhält also ein lineares Gleichungssystem

(I-\frac{h}{2} J^{(k)})y^{(k + 1)}_{n+1} =
-\frac{h}{2} J^{(k)} y^{(k)}_{n+1} +y_n + \frac{h}{2}(f^{(k)}_{n+1} + f_n),

wobei J die Jacobi-Matrix

J^{(k)} := \left(\frac{\partial f}{\partial y}\right)^{(k)}_{n+1},

I die Einheitsmatrix und k der Iterationsschritt ist.

Schrittweite h

Die (variable) Schrittweite kann aus folgender Beziehung berechnet werden:

\vert\frac{F(h\lambda)}{e^{h\lambda}} - 1\vert = \delta;

δ bezeichnet den zugelassenen lokalen Diskretisierungsfehler. Der Ansatz y_{n+1}=y_n+\frac{h}{2}(f_{n+1} + f_n)=:F(h\lambda)y_n liefert für die implizite Trapez-Methode

F(h\lambda)=\frac{2+h\lambda}{2-h\lambda}.

Dabei ist \lambda \, := \max_j{|\lambda_j|} der Betrag des betragsmäßig größten Eigenwerts der Jacobi-Matrix (Spektralradius). Die numerische Bestimmung der Eigenwerte ist sehr zeitaufwendig; für den Zweck der Schrittweitenberechnung ist es im Allgemeinen ausreichend die Gesamtnorm \lambda = N \cdot \max_{i,j} |a_{i j}| heranzuziehen, die immer größer oder gleich der Spektralnorm ist. N ist der Rang der Jacobi-Matrix und aij deren Elemente.

Literatur

  • H. R. Schwarz: Numerische Mathematik, B.G.Teubner Stuttgart; 1986

Einzelnachweise

  1. M. Kloker: Numerische Löser für die Gewöhnliche Modelldifferentialgleichung y'=αy, Universität Stuttgart, 1996

Wikimedia Foundation.

Игры ⚽ Поможем написать реферат

Schlagen Sie auch in anderen Wörterbüchern nach:

  • Implizite Trapez-Methode — Das implizite Trapez Verfahren ist ein Verfahren zur numerischen Lösung eines Anfangswert Problems Es lässt sich sowohl den Runge Kutta Verfahren als auch den Adams Moulton Verfahren zuordnen. Das Trapezverfahren ist A stabil mit der Besonderheit …   Deutsch Wikipedia

  • Runge-Kutta-Methode — Einige Runge Kutta Verfahren im Vergleich. Die s stufigen Runge Kutta Verfahren (nach Carl Runge und Martin Wilhelm Kutta) sind Einschrittverfahren zur näherungsweisen Lösung von Anfangswertproblemen in der numerischen Mathematik. Wenn vom Runge… …   Deutsch Wikipedia

  • Implizites Trapez-Verfahren — Das implizite Trapez Verfahren ist ein Verfahren zur numerischen Lösung eines Anfangswert Problems Es lässt sich sowohl den Runge Kutta Verfahren als auch den Adams Moulton Verfahren zuordnen. Das Trapezverfahren ist A stabil mit der Besonderheit …   Deutsch Wikipedia

  • Ratengleichungen — beschreiben die zeitlichen Konzentrationsverläufe verschiedener Spezies, zum Beispiel bei gekoppelten chemischen Reaktionen, oder bei An und Abregungsprozessen verschiedener Niveaus bei Atomen oder Molekülen: (Gl.1)         die an den Reaktionen… …   Deutsch Wikipedia

  • Trapezregel — Die Trapezregel beschreibt ein mathematisches Verfahren, wie man das Integral einer Funktion f(x) im Intervall [a,b] numerisch annähert (Numerische Quadratur). Dazu ersetzt man die Fläche unter der Kurve y = f(x) im gegebenen Intervall durch ein… …   Deutsch Wikipedia

  • Butcher-Tableau — Einige Runge Kutta Verfahren im Vergleich. Die s stufigen Runge Kutta Verfahren (nach Carl Runge und Martin Wilhelm Kutta) sind Einschrittverfahren zur näherungsweisen Lösung von Anfangswertproblemen in der numerischen Mathematik. Wenn vom Runge… …   Deutsch Wikipedia

  • Runge-Kutta — Einige Runge Kutta Verfahren im Vergleich. Die s stufigen Runge Kutta Verfahren (nach Carl Runge und Martin Wilhelm Kutta) sind Einschrittverfahren zur näherungsweisen Lösung von Anfangswertproblemen in der numerischen Mathematik. Wenn vom Runge… …   Deutsch Wikipedia

  • Trapezmethode — Das implizite Trapez Verfahren ist ein Verfahren zur numerischen Lösung eines Anfangswert Problems Es lässt sich sowohl den Runge Kutta Verfahren als auch den Adams Moulton Verfahren zuordnen. Das Trapezverfahren ist A stabil mit der Besonderheit …   Deutsch Wikipedia

  • Ratengleichung — Ratengleichungen beschreiben die zeitlichen Konzentrationsverläufe verschiedener Spezies, zum Beispiel bei gekoppelten chemischen Reaktionen, oder bei An und Abregungsprozessen verschiedener Niveaus bei Atomen oder Molekülen:       …   Deutsch Wikipedia

  • Rastern von Kreisen — Unter der Rasterung von Kreisen versteht man in der Computergrafik das Zeichnen (Rastern) eines Kreises auf dem Punktraster einer Rastergrafik oder eines Raster Grafikgeräts durch Einfärben entsprechender Pixel. Es gibt hierfür sowohl Algorithmen …   Deutsch Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”