Finite-Differenzen-Methode


Finite-Differenzen-Methode

Die Finite-Differenzen-Methode ist das einfachste numerische Verfahren zur Lösung gewöhnlicher und partieller Differentialgleichungen.

Zunächst wird das Gebiet, für das die Gleichung gelten soll, in eine endliche (finite) Zahl von Gitterpunkten eines Gitters von senkrecht aufeinander stehenden Linien eingeteilt. Den Gitterpunkten entsprechen dann die Kreuzungspunkte. Die Ableitungen an den Gitterpunkten werden dann durch Differenzen approximiert (siehe dazu Numerische Differentiation). Die partiellen Differentialgleichungen werden so in ein System von Differenzengleichungen umformuliert und mittels verschiedener Algorithmen entweder implizit oder explizit gelöst.

Verfahren dieser Art finden verbreitete Anwendung unter anderem bei fluiddynamischen Simulationen, zum Beispiel in der Meteorologie und der Astrophysik.

Differenzenquotient

Ein naheliegender Ansatz ist die Verwendung des Vorwärtsdifferenzenquotienten

D^+(x_0)= \frac{f(x_0+h)-f(x_0)}{h},

bzw. des Rückwärtsdifferenzenquotienten

D^-(x_0)= \frac{f(x_0)-f(x_0-h)}{h}.

Dabei ist jedoch die Näherung im Vergleich zur Auslöschung relativ schlecht. Eine bessere Näherung erhält man durch Verwendung des zentralen Differenzenquotienten

D(x_0)= \frac{f(x_0+h)-f(x_0-h)}{2h}.

Beispiel

Wir diskretisieren die lineare Differenzialgleichung f\,''(x) = 2 in Ω = [0;1] mit der Randbedingung f(0) = f(1) = 3 auf einem Gitter mit der Maschenweite  h = \frac{1}{n+1} . Zur Diskretisierung der zweiten Ableitung nehmen wir den zentralen Differenzenquotienten der zweiten Ableitung

D^+(D^-(x_0)) = \frac{f(x_0 + h) - 2f(x_0) + f(x_0 - h)}{h^2}

Unser Gebiet Ω wird also im Inneren durch (n-1) Knoten diskretisiert. Wir nummerieren diese Knoten von x1 bis xn − 1 durch. Für jeden dieser Knoten müssen wir die zweite Ableitung berechnen. Wir erhalten also (n-1) Gleichungen mit (n-1) Unbekannten.

 \frac{1}{h^2}( f(x_{2}) - 2f(x_1) ) = 2 - 3\frac{1}{h^2} und
 \frac{1}{h^2}(f(x_{i+1}) - 2f(x_i) + f(x_{i-1})) = 2 für i = {2,..n − 2} und
 \frac{1}{h^2}( - 2f(x_{n-1}) + f(x_{n-2})) = 2 - 3\frac{1}{h^2}

Daraus erhalten wir ein Gleichungssystem Ax = b mit einer im Allgemeinen dünnbesetzten Matrix A. Solche Gleichungssysteme lassen sich effizient mit iterativen Lösern berechnen.

Weblinks

 Commons: Finite-Differenzen-Methode – Sammlung von Bildern, Videos und Audiodateien

Wikimedia Foundation.

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

  • Finite-Differenzen — Die Finite Differenzen Methode ist das einfachste numerische Verfahren zur Lösung gewöhnlicher und partieller Differentialgleichungen. Zunächst wird das Gebiet, für das die Gleichung gelten soll, in eine endliche (finite) Zahl von Gitterpunkten… …   Deutsch Wikipedia

  • Finite-Differenzen-Verfahren — Die Finite Differenzen Methode ist das einfachste numerische Verfahren zur Lösung gewöhnlicher und partieller Differentialgleichungen. Zunächst wird das Gebiet, für das die Gleichung gelten soll, in eine endliche (finite) Zahl von Gitterpunkten… …   Deutsch Wikipedia

  • Finite Differenzen — Die Finite Differenzen Methode ist das einfachste numerische Verfahren zur Lösung gewöhnlicher und partieller Differentialgleichungen. Zunächst wird das Gebiet, für das die Gleichung gelten soll, in eine endliche (finite) Zahl von Gitterpunkten… …   Deutsch Wikipedia

  • Finite Difference Time Domain — (FDTD, englisch für Finite Differenzen Methode im Zeitbereich) ist ein mathematisches Verfahren zur direkten Integration zeitabhängiger Differentialgleichungen. Vor allem zur Berechnung der Lösungen der Maxwell Gleichungen wird dieses Verfahren… …   Deutsch Wikipedia

  • Finite-Differenz — Die Finite Differenzen Methode ist das einfachste numerische Verfahren zur Lösung gewöhnlicher und partieller Differentialgleichungen. Zunächst wird das Gebiet, für das die Gleichung gelten soll, in eine endliche (finite) Zahl von Gitterpunkten… …   Deutsch Wikipedia

  • Crank-Nicolson-Methode — Das Crank Nicolson Verfahren ist in der numerischen Mathematik eine Finite Differenzen Methode zur Lösung der Wärmeleitungsgleichung und ähnlicher partieller Differentialgleichungen.[1] Es ist ein implizites Verfahren 2. Ordnung und numerisch… …   Deutsch Wikipedia

  • Randelement-Methode — Die Randelementmethode (REM, englisch boundary element method, BEM, gelegentlich auch Momentenmethode, englisch MoM , method of moments) ist ein Diskretisierungsverfahren zur Berechnung von Anfangs und Randwertproblemen mit partiellen… …   Deutsch Wikipedia

  • Transmission Line Matrix Methode — Viele Probleme der modernen Mikrowellentechnik lassen sich nur durch numerische Verfahren zufriedenstellend lösen. Z. B. lassen sich die technischen Eigenschaften einer Handyantenne wie Abstrahlverhalten nicht mehr mit geschlossenen… …   Deutsch Wikipedia

  • Level Set Methode — Die Level Set Methode oder Niveaumengenmethode ist ein mathematisches Verfahren um geometrische Objekte und deren Bewegung numerisch zu berechnen. Der Vorteil der Level Set Methode liegt darin, dass man Kurven und Oberflächen auf einem räumlich… …   Deutsch Wikipedia

  • Transmission-Line-Matrix-Methode — Viele Probleme der modernen Mikrowellentechnik lassen sich nur durch numerische Verfahren zufriedenstellend lösen. Zum Beispiel lassen sich die technischen Eigenschaften einer Handyantenne wie Abstrahlverhalten nicht mehr mit geschlossenen… …   Deutsch Wikipedia