Numerische Integration - LinkFang.de





Numerische Integration


In der numerischen Mathematik bezeichnet numerische Integration (traditionell auch als numerische Quadratur bezeichnet) die näherungsweise Berechnung von Integralen. Oft kann man Integrale nicht geschlossen lösen, weil für den Integranden keine Stammfunktion angegeben werden kann oder er nur durch diskrete Werte, etwa Messungen, gegeben ist. Dann versucht man, Näherungswerte zu ermitteln.

Dazu wird das Integral einer Funktion [math]f[/math] über dem Intervall [math][a,b][/math] dargestellt als Summe aus dem Wert [math]Q(f)[/math] einer Quadraturformel [math]Q[/math] und dem Fehler [math]E(f)[/math]:

[math]\int_{a}^{b}f(x)\,dx= Q(f) + E(f)\,.[/math]

Eine Quadraturformel besteht dabei im Allgemeinen aus einer gewichteten Summe von Funktionswerten

[math]Q(f) = (b-a)\sum_{i=0}^n w_i f(x_i).[/math]

Die Stellen [math]x_0,\ldots, x_n[/math] heißen Stützstellen und die Zahlen [math]w_0,\ldots,w_n[/math] Gewichte. Es existieren verschiedene Ansätze, wie Stützstellen und Gewichte so gewählt werden können, dass der Quadraturfehler [math]E(f)[/math] möglichst klein wird.

Eine Quadraturformel hat den Genauigkeitsgrad (oder auch Exaktheitsgrad) [math]n[/math], wenn sie alle Polynome bis zum Höchstgrad [math]n[/math] exakt integriert, und [math]n[/math] die größtmögliche natürliche Zahl mit dieser Eigenschaft ist.

Ebenso wie das Integral sind Quadraturformeln lineare Operatoren.

Interpolatorische Quadraturformel

Eine wichtige Klasse von Quadraturformeln ergibt sich durch die Idee, die Funktion [math]f(x)[/math] durch ein Interpolationspolynom [math]p_n(x)[/math] vom Grad [math]n[/math] zu approximieren und dieses dann zu integrieren. Die Gewichte ergeben sich dann als die Integrale der Lagrange-Polynome zu den gegebenen Stützstellen. Nach Konstruktion haben diese Quadraturformeln mindestens den Genauigkeitsgrad [math]n[/math]. Die Quadraturformel lautet also

[math]\int_a^b f(x) dx \approx \int_a^b p_n(x) dx = (b-a) \sum_{i=0}^n w_i f(x_i)[/math]

mit den Gewichten

[math]w_i = \frac{1}{(b-a)} \int_a^b L_{in}(x) dx [/math]

und den Lagrangepolynomen

[math]L_{in}(x) = \frac{(x-x_0)\cdots(x-x_{i-1})(x-x_{i+1})\cdots(x-x_n)}{(x_i-x_0)\cdots(x_i-x_{i-1})(x_i-x_{i+1})\cdots(x_i-x_n)}.[/math]

Falls die Integrationsgrenzen Stützstellen sind, spricht man von abgeschlossenen Quadraturformeln, sonst von offenen. Werden die Stützstellen äquidistant gewählt, so ergeben sich unter anderen die Newton-Cotes-Formeln. Zu den abgeschlossenen Newton-Cotes-Formeln gehören die Sehnentrapezregel und die Simpson-Regel, zu den offenen gehört die Tangententrapezregel. Die Newton-Cotes-Formeln für gerades [math]n[/math] haben sogar den Genauigkeitsgrad [math]n+1[/math]. Zu den offenen Quadraturformeln gehören auch die Gauß-Quadraturformeln.

Fehlerabschätzung

Mit [math][c,d][/math] sei das kleinste Intervall bezeichnet, das die Stützstellen [math]x_i[/math] und das Intervall [math][a,b][/math] enthält. Ferner sei [math]f[/math] (n+1)-mal stetig differenzierbar auf [math][c,d][/math]. Gemäß der Interpolationsgüte des Interpolationspolynoms gibt es ein [math]\xi(x) \in [c,d] [/math], so dass gilt:

[math]f(x) - p_n(x) = \frac{f^{(n+1)}(\xi(x))}{(n+1)!} \prod_{i=0}^n (x-x_i).[/math]

Durch Integration erhält man die Fehlerformel für die numerische Quadratur

[math]E(f) = \int_a^bf(x) dx - \int_a^bp_n(x) dx = \frac{1}{(n+1)!} \int_a^b f^{(n+1)}(\xi(x))\prod_{i=0}^n (x-x_i)dx.[/math]

Falls [math]f^{(n+1)}(x)=0[/math] für alle [math]x \in [c,d][/math] gilt, ist der Quadraturfehler gleich 0. Da das für alle Polynome bis zum Grad [math]n[/math] der Fall ist, ist der Genauigkeitsgrad dieser Quadraturformeln mindestens [math]n[/math].

Aus dieser Fehlerformel folgt die Fehlerabschätzung

[math]\left| E(f) \right| \le {1 \over (n+1)!}\ \max_{c\le x \le d} {\left| f^{(n+1)}(x) \right|} \int_a^b \left|\prod_{i=0}^n (x-x_i)\right|dx.[/math]

Falls die Funktion [math]\prod_{i=0}^n (x-x_i)[/math] im Intervall [math][a,b][/math] ihr Vorzeichen nicht wechselt, d.h. wenn keine Stützstelle im Intervall [math](a,b)[/math] liegt, kann man mit Hilfe des Mittelwertsatzes der Integralrechnung folgende Darstellung für das Restglied herleiten:

[math]E(f) = \frac{f^{(n+1)}(\zeta)}{(n+1)!} \int_a^b \prod_{i=0}^n (x-x_i)dx.[/math]

mit einer Zwischenstelle [math]\zeta\in[c,d].[/math]

Ähnliche Formeln für den Quadraturfehler erhält man auch bei speziellen Verteilungen der Stützstellen im Intervall [math][a,b][/math], etwa für die Newton-Cotes-Formeln oder die Gauß-Quadraturformeln.

Ist die Funktion [math]f[/math] nur stetig, so gelten obige Aussagen nicht, der Fehler kann sehr groß werden.

Weitere Quadraturformeln

Der Versuch, die Fehlerordnung der Quadraturformel zu minimieren, führt auf die Gauß-Quadratur. Diese nutzen die Theorie orthogonaler Polynome, um Formeln zu erhalten, die den Genauigkeitsgrad [math]2n[/math] haben, wobei [math]n[/math] die Anzahl der genutzten Funktionsauswertungen ist.

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 Streifenbreite hin extrapoliert.

Summierte Quadraturformeln

Um das Integral noch besser annähern zu können, unterteilt man das Intervall [math][a,b][/math] in [math]N[/math] nebeneinanderliegende Teilintervalle. Die Teilintervalle müssen 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. Von besonderem Interesse sind adaptive Formeln, die keine weitere Unterteilung eines Intervalls vornehmen, wenn der in dem Intervall geschätzte Fehler unterhalb einer Schranke liegt.

Monte-Carlo-Integration

Ein Verfahren, das nicht versucht, eine Näherungsformel für die zu integrierende Funktion heranzuziehen, ist die Monte-Carlo-Integration. Anschaulich gesagt wird hierbei das Integral dadurch bestimmt, dass [math]n[/math] zufällige Punkte [math]x_1,\dots,x_n[/math] im Integrationsintervall [math][a,b][/math] (horizontal) erzeugt werden. Dann ergibt sich eine Näherung des Integrals als Durchschnitt der Funktionswerte dieser Stellen

[math]S_n(f)=\frac {b-a}n \sum_{i=1}^n f(x_i).[/math]

Der Vorteil ist die vergleichsweise einfache Implementierung sowie die relativ einfache Erweiterbarkeit auf Vielfachintegrale. Der Rechenaufwand ist etwas höher im Vergleich zu Quadrationsformeln.

Literatur

  • Hans R. Schwarz, Norbert Köckler: Numerische Mathematik. 6. Auflage, Teubner, Stuttgart 2006, ISBN 3-519-42960-8

Weblinks

 Commons: Graphen zur numerischen Integration  – Sammlung von Bildern, Videos und Audiodateien

Kategorien: Numerische Mathematik

Quelle: Wikipedia - http://de.wikipedia.org/wiki/Numerische Integration (Vollständige Liste der Autoren des Textes [Versionsgeschichte])    Lizenz: CC-by-sa-3.0

Änderungen: Alle Bilder mit den meisten Bildunterschriften wurden entfernt. Ebenso alle zu nicht-existierenden Artikeln/Kategorien gehenden internen Wikipedia-Links (Bsp. Portal-Links, Redlinks, Bearbeiten-Links). Entfernung von Navigationsframes, Geo & Normdaten, Mediadateien, gesprochene Versionen, z.T. ID&Class-Namen, Style von Div-Containern, Metadaten, Vorlagen, wie lesenwerte Artikel. Ansonsten sind keine Inhaltsänderungen vorgenommen worden. Weiterhin kann es durch die maschinelle Bearbeitung des Inhalts zu Fehlern gerade in der Darstellung kommen. Darum würden wir jeden Besucher unserer Seite darum bitten uns diese Fehler über den Support mittels einer Nachricht mit Link zu melden. Vielen Dank!

Stand der Informationen: August 201& - Wichtiger Hinweis: Da die Inhalte maschinell von Wikipedia übernommen wurden, ist eine manuelle Überprüfung nicht möglich. Somit garantiert LinkFang.de nicht die Richtigkeit und Aktualität der übernommenen Inhalte. Sollten die Informationen mittlerweile fehlerhaft sein, bitten wir Sie darum uns per Support oder E-Mail zu kontaktieren. Wir werden uns dann innerhalb von spätestens 10 Tagen um Ihr Anliegen kümmern. Auch ohne Anliegen erfolgt mindestens alle drei Monate ein Update der gesamten Inhalte.