28.01.2013 Aufrufe

Runge-Kutta Verfahren

Runge-Kutta Verfahren

Runge-Kutta Verfahren

MEHR ANZEIGEN
WENIGER ANZEIGEN

Erfolgreiche ePaper selbst erstellen

Machen Sie aus Ihren PDF Publikationen ein blätterbares Flipbook mit unserer einzigartigen Google optimierten e-Paper Software.

<strong>Runge</strong>-<strong>Kutta</strong><br />

<strong>Runge</strong> <strong>Kutta</strong> <strong>Verfahren</strong><br />

Ordnungstheorie der<br />

<strong>Runge</strong>-<strong>Kutta</strong><br />

<strong>Runge</strong> <strong>Kutta</strong> <strong>Verfahren</strong>


Übersicht bersicht<br />

1. Wiederholung<br />

2. Motivation und Vorgehensweise<br />

3. Definition der RK-Gewichte<br />

RK Gewichte<br />

4. Die Butcher-Reihe<br />

Butcher Reihe eines<br />

RK-<strong>Verfahren</strong>s<br />

RK <strong>Verfahren</strong>s<br />

5. Die Ordungsgleichungen<br />

6. Bemerkungen<br />

7. Zusammenfassung


1. Wiederholung<br />

�� Ein s-stufiges s stufiges <strong>Verfahren</strong> zur Lösung L sung von<br />

ist eine Abb. der Form:<br />

�� Die Taylor-Entwicklung<br />

Taylor Entwicklung<br />

des Flusses mittels elementarer<br />

Differentiale heißt hei t Butcher-Reihe<br />

Butcher Reihe


2. Motivation & Vorgehensweise<br />

�� Ziel:<br />

Identifizieren der Butcher Reihe von<br />

�� Vorgehensweise:<br />

- Definition der <strong>Runge</strong> <strong>Kutta</strong> Gewichte<br />

- allg. Butcher Reihe der RK-<strong>Verfahren</strong><br />

RK <strong>Verfahren</strong><br />

�� Vergleich von und<br />

�� Definition der Ordungsgleichungen


2.1 Bsp.: Die Butcher-Reihe<br />

Butcher Reihe des<br />

impliziten Euler-<strong>Verfahren</strong>s<br />

Euler <strong>Verfahren</strong>s<br />

�� Sei<br />

die Lösung L sung von<br />

Dann gilt:


3. Definition der RK-Gewichte<br />

RK Gewichte<br />

�� Das Butcher-Schema<br />

Butcher Schema<br />

sei gegeben.<br />

Zum Baum ρτ wähle hle λρτ = (V,E,r ( V,E,r) ) є ρτ mit<br />

V={1,…,n}.<br />

V={1, ,n}.<br />

Definiere die RK-Gewichte:<br />

RK Gewichte:


3.1 Bemerkung<br />

�� Für r eine Kante (α,β) ( ) mit einem Endknoten<br />

β kann eine Summation ausgeführt<br />

ausgef hrt<br />

werden und liefert<br />

�� Es folgt die anschauliche Konstruktion des<br />

Gewichtes Φ: :<br />

• hefte an die Wurzel eine Kopie von b<br />

• an jede Kante eine Kopie von A<br />

• für r "Endkanten" kann diese zu einer Kopie von<br />

c vereinfacht werden<br />

• dann multipliziere alles und addiere.


3.2 Beispiele


3. Butcher Reihe eines<br />

RK-<strong>Verfahren</strong>s<br />

RK <strong>Verfahren</strong>s<br />

�� Für r die RK-Abbildung<br />

RK Abbildung<br />

zur Lösung L sung von gilt :


5. Ordnungsgleichungen<br />

�� Der Vergleich mit dieser Abb. ergibt:


��RK RK-<strong>Verfahren</strong> <strong>Verfahren</strong> hat Ordnung p genau dann<br />

wenn für f r alle Bäume B ume ρτ mit |ρτ | ρτ|


Die Ordnungsgleichungen<br />

bis zu<br />

Ordnung 4


6. Bemerkungen<br />

I. Sei = Anzahl aller Bäume B ume ρτ mit<br />

genau p Knoten.<br />

Die Anzahl der Gleichungen<br />

steigt schnell mit der gewünschten<br />

gew nschten<br />

Ordnung p:


II. Die Konsistenzbedingung<br />

garantiert Konsistenz/Konvergenz.<br />

Die "Büschel" "B schel"<br />

liefern die Quadraturbedingungen<br />

k=1,…,p k=1, ,p


III. Bei der Anwendung von RK-<strong>Verfahren</strong><br />

RK <strong>Verfahren</strong><br />

auf skalare Gleichungen ,<br />

fallen einige elementare Differentiale<br />

zusammen, so daß da im Vergleich der<br />

Butcher-Reihen<br />

Butcher Reihen von und nicht für f r<br />

jeden Baum getrennt<br />

gefordert zu werden braucht, z.B.


7. Zusammenfassung<br />

�� Die allgemeine Definition der Butcher<br />

Reihe für f r ist :<br />

�� Der führende f hrende Fehlerterm hat eine<br />

Abweichung von:<br />

für r ein RK-<strong>Verfahren</strong> RK <strong>Verfahren</strong> der Ordnung p

Hurra! Ihre Datei wurde hochgeladen und ist bereit für die Veröffentlichung.

Erfolgreich gespeichert!

Leider ist etwas schief gelaufen!