Numerische Integration
aus Wikipedia, der freien Enzyklopädie
In der numerischen Mathematik bezeichnet numerische Quadratur bzw. numerische Integration die näherungsweise Berechnung von Integralen. Oft kann man Integrale nicht geschlossen lösen, d.h. man kann keine Stammfunktion zu f(x) angeben. Deshalb versucht man, Näherungswerte zu ermitteln.
Wir bezeichnen mit
das Integral der Funktion f(x) im Intervall [a,b]. Dies wird hier dargestellt als der Wert einer Quadraturformel Q(f) plus dem Fehler E(f).
Dazu unterteilt man die gesuchte Fläche in senkrechte Streifen und nähert jede dieser so erhaltenen Teilflächen durch einfache geometrische Figuren (z.B. Trapez) oder einfache Funktionen (z.B. Polynome) an.
Für die Flächenberechnung dieser einfachen Figuren benötigt man den Wert der Funktion f(x) an den so genannten Stützstellen x0,...xm. Die Summe über diese Teilflächen ergibt eine Näherung Q(f) des Integrals. Je schmaler man die einzelnen Teilflächen wählt desto genauer wird die Näherung. Von Interesse ist dann noch die Frage, wie groß der Fehler ist, der sich durch die Näherung ergibt. Dieser Fehler wird durch das Restglied E(f) beschrieben. Um die Anzahl der Funktionsauswertungen zu minimieren, bei gleichzeitiger Möglichkeit den Fehler zu kontrollieren, verwendet man oft das Rombergsche Extrapolationsverfahren. Hierbei werden die Integralwerte von immer kleiner werdenden 'Streifen' zu einer verschwindenden Breite hin extrapoliert.
Inhaltsverzeichnis |
[Bearbeiten] Spezielle Quadraturformeln
Man hat nun verschiedene Möglichkeiten, die einzelnen Teilflächen durch spezielle einfachere Flächen anzunähern. Die Anwendung der allgemeinen Quadraturformeln auf diese speziellen Flächen liefert einige bekannte und wichtige spezielle Quadraturformeln.
[Bearbeiten] Sehnentrapezformel
Man ersetzt die Kurve f(x) durch die Verbindungsgerade zwischen den Punkten (a,f(a)) und (b,f(b)) - also durch die Sehne - und erhält somit ein Trapez.
Mit
- m = 1
- x0 = a
- x1 = b
erhält man die Koeffizienten βj
- z0 = 0
- z1 = 1
- β0 = 1
- β1 = 0,5
- β2 = -1/6
und daraus schließlich die Sehnentrapezformel
Oder:
Fläche A von einem Trapez kann folgendermaßen berechnet werden: (aus der Geometrie)
In unserem Beispiel schneiden wir ein Trapez im Intervall [a,b] aus. Daraus folgt:
- Seite1 = f(a)
- Seite2 = f(b)
- Hoehe = b − a
Ergebnis:
[Bearbeiten] Tangententrapezformel
Man legt an die Kurve f(x) im Punkt c in der Mitte des Intervalls [a,b] die Tangente und erhält so wieder ein Trapez.
Mit
- m = 1
- x0 = c
- x1 = c
erhält man die Koeffizienten βj
- z0 = 0,5
- z1 = 0,5
- β0 = 1
- β1 = 0
- β2 = 1/12
und daraus schließlich die Tangententrapezformel
Oder:
Die Fläche A von dem Trapez, welches durch die Tangente begrenzt wird, lässt sich folgendermaßen berechnen:
Flaeche = (Funktionswert an der Stelle c) Hoehe
(dadurch wird dieses Trapez zu einem Rechteck eingeebnet)
Die Fläche vom Trapez (mit den Seiten f(a) und f(b)) entspricht der Fläche von einem Rechteck mit der Seite f(c). c befindet sich in der Mitte zwischen a und b.
- Hoehe = b − a
Ergebnis:
[Bearbeiten] Simpsonsche Formel oder Keplersche Fassregel
Interpoliert man die Funktion f(x) mittels eines quadratischen Polynoms in den Punkten a,b,c (c liegt in der Mitte von [a,b]), dann erhält man die simpsonsche Formel.
Mit
- m = 3
- x0 = a
- x1 = b
- x2 = c
- x3 = c
erhält man die Koeffizienten βj
- z0 = 0
- z1 = 1
- z2 = 0,5
- z3 = 0,5
- β0 = 1
- β1 = 1/2
- β2 = -1/6
- β3 = 0
- β4 = -1/120
und damit die Simpsonsche Formel
[Bearbeiten] Hermitsche Formel
Mit
- m = 3
- x0 = a
- x1 = b
- x2 = a
- x3 = b
ergeben sich
- z0 = 0
- z1 = 1
- z2 = 0
- z3 = 1
- β0 = 1
- β1 = 1/2
- β2 = -1/6
- β3 = -1/12
- β4 = -1/30
[Bearbeiten] Allgemeine Quadraturformel
Mit Hilfe von Interpolationspolynomen und deren Lagrange-Darstellung kann man die folgende allgemeine Quadraturformel und das zugehörige Restglied herleiten.
Die allgemeine Quadraturformel für eine Teilfläche lautet
mit den Koeffizienten
Das Restglied beträgt
Ist die Funktion f im Intervall [a,b] (m + 1)-mal stetig differenzierbar ("reellwertig" wird nicht gefordert), dann lässt sich das Restglied nach oben abschätzen durch
Wenn noch zusätzlich für alle Stützstellen im Intervall [a,b] gilt oder alternativ , dann hat der Integrand keinen Vorzeichenwechsel in [a,b] und man kann zeigen:
Daraus folgt dann die Restgliedabschätzung
Ist die Funktion f zusätzlich noch reellwertig in [a,b], dann kann man mit Hilfe des Mittelwertsatzes der Integralrechnung folgende Darstellung für das Restglied herleiten:
mit einer Zwischenstelle ζ im Intervall [a,b].
[Bearbeiten] Summierte Quadraturformeln
Um das Integral noch besser annähern zu können unterteilt man das Intervall [a,b] in N nebeneinanderliegende Teilintervalle .
Die Teilintervalle müssen zunächst nicht die gleiche Länge haben.
In jedem Teilintervall wendet man im Folgenden die gleiche Näherung für die einzelnen Flächen an und addiert danach die entstandenen Näherungen.
Es gilt für jede Teilfläche
Daraus folgt für das gesamte Integral
mit
Sei f nun (m + 1)-mal stetig differenzierbar im Gesamtintervall [a,b]. Ferner sollen ab jetzt alle Teilintervalle die gleiche Länge h haben, also
Dann gilt für die einzelnen Restglieder (siehe oben)
Summierung über die einzelnen Restglieder ergibt die Abschätzung für das gesamte Restglied
mit Nh = b − a.
Ist die Funktion f zudem auf [a,b] reellwertig, dann kann man für das Restglied analog herleiten: