Die Simpsonregel oder simpsonsche Formel (nach Thomas Simpson) ist ein Verfahren der numerischen Integration, bei dem eine Näherung zum Integral einer in einem Intervall schwer zu integrierenden Funktion berechnet wird, indem man die Funktion durch eine exakt integrierbare Parabel annähert.
Inhaltsverzeichnis
1Regel
2Fehlerabschätzung
3Veranschaulichung durch Rechteckflächen
4Beispiel
5Summierte simpsonsche Formel
5.1Variante 1
5.1.1Beispiel
5.1.2Fehlerabschätzung
5.1.3Beispiel
5.1.4Fehlerschätzung
5.2Variante 2
5.2.1Beispiel
5.2.2Fehlerabschätzung
6Zusammenhang mit anderen Formeln
7Geschichte
7.1Keplersche Fassregel
7.1.1Parabolische Krümmung
8Verwendung als Runge-Kutta-Verfahren
9Literatur
10Weblinks
11Einzelnachweise
Regel
Simpsonsche Formel
Für eine Funktion im Intervall wird eine Parabel als Interpolationspolynom durch die Funktionswerte an den Stellen , und gelegt. Das Integral von nähert man dann durch das Integral der Parabel an. Die Simpsonsregel für das Integral
lautet dann
.
Der Wert ist dann eine Näherung von . Somit ist die Simpsonregel eine abgeschlossene Newton-Cotes-Formel.[1]
Fehlerabschätzung
Das Restglied (alternativ der Quadraturfehler) beschreibt die Differenz des tatsächlichen Integrals und der Näherung durch die Simpsonregel:
Ist viermal stetig differenzierbar in , dann gilt für das Restglied die Abschätzung
Ist zusätzlich noch reellwertig, dann gilt mit einer geeigneten Zwischenstelle aus für das Restglied
Diese Restglieddarstellung wurde 1887 von Giuseppe Peano gefunden. Sie besagt insbesondere, dass die Simpsonregel Polynome vom Grad drei exakt integriert, also einen Grad höher, als man nach Konstruktion erwarten würde. Diese Eigenschaft haben alle abgeschlossenen und offenen Newton-Cotes-Formeln von geradem Grad.[2]
Veranschaulichung durch Rechteckflächen
Simpsonsche Formel Veranschaulichung
Das Integral der Näherungs-Parabel ist gleich der schraffierten Fläche von sechs Rechtecken, deren Breite jeweils des Intervalls ist. Ein Rechteck hat dabei die Höhe , ein Rechteck die Höhe und vier Rechtecke die Höhe .
Hier sieht man auch den Zusammenhang mit der Sehnentrapezformel
Während die zwei äußeren Rechtecke der mit skalierten Sehnentrapezformel entsprechen, entsprechen die übrigen Rechtecke der mit skalierten Tangententrapezformel. Es ergibt sich die Ausgangsformel der simpsonschen Formel:
aus der sich nach Einsetzen und Umformen die bereits bekannte simpsonsche Formel ergibt:
.
Beispiel
Wir betrachten im Folgenden das Beispiel
Bemerkung: Diese Funktion dient der Veranschaulichung ist auf einfache Weise elementar integrierbar. In der Praxis werden die simpsonsche Formel und effizientere Verfahren zur numerischen Integration üblicherweise auf Funktionen angewendet, die nicht oder nicht effektiv elementar integrierbar sind.
Gute Ergebnisse erzielt man, wenn das Intervall in mehrere Teilintervalle zerlegt und die Simpsonregel in jedem dieser Teilintervalle angewendet wird. Diese Summationsvarianten werden nachfolgend vorgestellt und jeweils auf dieses Beispiel angewendet.
Summierte simpsonsche Formel
Um das Integral noch besser annähern zu können, unterteilt man das Intervall in nebeneinanderliegende, gleich große Teilintervalle. In jedem Teilintervall wendet man die simpsonsche Formel für die einzelnen Teilflächen an und addiert danach die entstandenen Näherungen. Damit erhält man die summierte oder zusammengesetzte Simpsonregel. Es gibt unterschiedliche Notationen für die Unterteilung in Teilintervalle, die zu verschiedenen Formulierungen der summierten simpsonschen Formel führen.
Variante 1
Summierte simpsonsche Formel für
Hier unterteilt man das Intervall in nebeneinanderliegende, gleich große Teilintervalle der Länge . In jedem Teilintervall wendet man die simpsonsche Formel
an und addiert danach die entstandenen Näherungen. Mit
und
erhält man:
bzw.
.
Beispiel
Angewandt auf obiges Beispiel:
Sei und somit die Schrittweite . Dann ist
.
Sei und somit die Schrittweite . Dann ist
Dabei ist der Wert der Sehnentrapezregel und der Wert der Tangententrapezregel.
Fehlerabschätzung
Die Fehlerabschätzung für das Restglied lautet:
beziehungsweise für reellwertige Funktionen mit einer geeigneten Zwischenstelle aus dem Intervall
Der Faktor in obiger Formel bedeutet, dass bei einer Halbierung der Schrittweite (Verdoppelung der Intervalle), wie es beim Romberg-Verfahren mit der Romberg-Folge der Fall ist, der Fehler in etwa um den Faktor 16 kleiner wird, wie auch nachfolgendes Beispiel zeigt:
Beispiel
Angewandt auf obiges Beispiel:
Mit folgt
und somit die Fehlerabschätzung
,
die erwartungsgemäß einen größeren Wert ergibt als den exakten Wert
Analog erhält man die Fehlerabschätzung
,
die erwartungsgemäß einen größeren Wert ergibt als den exakten Wert
.
Es gilt
.
Fehlerschätzung
Rechnet man die Simpsonregel zweimal mit 2 verschiedenen Anzahlen von Intervallen , so erhält man folgende Fehlerschätzung:
Speziell bei der Verdoppelung der Intervalle (Halbierung der Schrittweite) erhält man die Fehlerschätzung:
Angewandt auf das obige Beispiel erhält man
Variante 2
Summierte simpsonsche Formel für N=4
Hier unterteilt man das Intervall in nebeneinanderliegende, gleich große Teilintervalle mit Mittelpunkt und Länge mit . Da jetzt gegenüber Variante 1 doppelt so groß ist, ist gegenüber Variante 1 nur halb so groß. Somit muss in allen Formeln von Variante 1 das durch ersetzt werden.
Für jedes gerade wendet man auf das Intervall die simpsonsche Formel
an und addiert danach die entstandenen Näherungen.
Für gerades gilt nun , und und man erhält:
bzw.
Beispiel
Angewandt auf obiges Beispiel:
Sei , und die Schrittweite . Dann ist
.
Das ist das gleiche Resultat wie in Variante 1.
Fehlerabschätzung
Die Fehlerabschätzung für das Restglied lautet nun
beziehungsweise für reellwertige Funktionen mit einer geeigneten Zwischenstelle aus dem Intervall
Zusammenhang mit anderen Formeln
Addiert man zum Näherungswert die Fehlerschätzung für , so erhält man die i. A. bessere Formel:
Das ist die Formel für die 3. Spalte des Rechenschemas der Romberg-Integration bei Verwendung der Romberg-Folge und gleichzeitig das Resultat der Milne-Regel (Abgeschlossene Newton-Cotes-Formel mit Genauigkeitsgrad 5) bei Anwendung auf Teilintervalle von
Angewandt auf obiges Beispiel erhält man mit
eine bessere Näherung für das exakte Integral
als mit oder
bei gleicher Anzahl auszuwertender Funktionswerte wie bei , nämlich 13 Stück.
Geschichte
Die Formel wurde erstmals von dem 1608 geborenen Evangelista Torricelli benutzt, ist aber nach dem 1710 geborenen englischen Mathematiker Thomas Simpson benannt.
Keplersche Fassregel
Die Anwendung der Simpsonregel auf Rotationskörper entspricht der Keplerschen Fassregel, die Johannes Kepler bereits 1615 aufstellte. Über die Entstehungsgeschichte berichtet Kepler in der Widmung der späteren Veröffentlichung. Nachdem 1611 Keplers erste Frau in Prag gestorben war, heiratete er – nun in Linz arbeitend – 1613 wieder. Er kaufte für die Hochzeit einige Fässer Wein. Als der Wein eingekellert war, kam der Verkäufer mit einer Messrute und bestimmte den Inhalt für alle Fässer ohne Überlegung oder Rechnung nach der gleichen Methode. Die Messrute wurde mit ihrer metallenen Spitze durch das Spundloch quer bis zu den Rändern der beiden Böden eingeführt und die Marke am Spundloch ergab den Rauminhalt. Kepler wunderte sich, dass eine Diagonale durch die Fasshälfte ein Maß für den Rauminhalt abgeben sollte, und bezweifelte die Richtigkeit dieser Methode, da ein sehr niedriges Fass mit etwas breiteren Böden und daher sehr viel kleinerem Rauminhalt die gleiche Visierlänge besitzen könnte.
Kepler verfasste daraufhin die Schrift Nova Stereometria doliorum vinariorum 1615 (Neue Inhaltsberechnung von Weinfässern), in der er nach überprüfbaren Methoden zur Inhaltsberechnung von Weinfässern suchte. Eine dieser Methoden bestand darin, die Krümmung des Fasses durch eine Parabel anzunähern, da Inhaltsberechnungen mit Hilfe von Parabeln seit Archimedes exakt durchgeführt werden konnten.
Der Name Fassregel lässt sich durch die folgende Anwendung begründen: Zur Berechnung des Volumens eines Weinfasses sei die Querschnittsfläche quer zur Längsachse in der Entfernung vom Boden des Fasses. Sie lässt sich durch Bestimmung des Umfanges leicht ausrechnen. Ist die Höhe des Fasses, so ist das Volumen gleich
Die Keplersche Fassregel gibt nun
als Näherungswert für das Volumen eines Körpers, dessen Querschnitt an drei Stellen bekannt ist. Ist der Körper ein Rotationskörper, so gilt bei Rotation der (Radius-)Funktion um die x-Achse:
Ist der Umfang von Boden und Deckel und der Umfang in der Mitte des Fasses, so ergibt sich daraus der Näherungswert
Also (selbstverständlich Innenmaße):
oder
.
Parabolische Krümmung
Fass mit parabolischer Krümmung
Hat das Fass eine parabolische Krümmung, so erhält man das Fass durch Rotation der Funktion des Radius um die x-Achse. Legt man zur Vereinfachung das Achsenkreuz in die Mitte des Fasses, so gilt
An diesem Beispiel kann man die Gültigkeit der oben angegebenen Formel
gut verifizieren. Hier ist und
ein Polynom vom Grad 4 mit konstanter 4. Ableitung:
Für den Fehler erhält man
somit den gleichen Wert wie oben.
Der Fehler ist kleiner oder gleich Null. Damit ist die Näherung größer oder gleich dem exakten Volumen. Der Fehler ist umso größer, je mehr sich und unterscheiden, je gewölbter das Fass ist. Der Fehler ist genau dann 0, wenn , das Fass also ein Zylinder ist, in Übereinstimmung mit obiger Aussage, dass die Formel für Zylinder exakt ist.
Hans R. Schwarz, Norbert Köckler: Numerische Mathematik. 6. Auflage, Teubner, Stuttgart 2006, ISBN 3-519-42960-8, S. 311–316.
Johannes Kepler: Neue Stereometrie der Fässer. Aus dem Lateinischen übersetzt und herausgegeben von R. Klug. W. Engelmann. Leipzig, 1908.
Weblinks
Commons: Simpson's rule – Sammlung von Bildern, Videos und Audiodateien
Website des Museumsraum zu Johannes Kepler des Bundesrealgymnasium Kepler, Graz
Einzelnachweise
↑Josef Stoer: Einführung in die Numerische Mathematik I. Unter Berücksichtigung von Vorlesungen von F. L. Bauer (= Heidelberger Taschenbücher. Band105). 9. Auflage. Springer-Verlag, Berlin / Heidelberg / New York / Tokyo 2005, ISBN 3-540-21395-3, S.157.
↑Josef Stoer: Einführung in die Numerische Mathematik I. Unter Berücksichtigung von Vorlesungen von F. L. Bauer (= Heidelberger Taschenbücher. Band105). 9. Auflage. Springer-Verlag, Berlin / Heidelberg / New York / Tokyo 2005, ISBN 3-540-21395-3, S.161ff.