Numerische Quadratur mit Derivaten


19

Die meisten numerischen Quadraturmethoden behandeln den Integranden als Black-Box-Funktion. Was ist, wenn wir mehr Informationen haben? Welchen Nutzen können wir insbesondere daraus ziehen, die ersten Ableitungen des Integranden zu kennen? Welche anderen Informationen könnten wertvoll sein?

Insbesondere für Ableitungen: Fehlerschätzungen für die Grundquadratur (Rechteck- / Trapez- / Simpsons-Regeln) sind eng miteinander verbunden. Vielleicht gibt es eine Möglichkeit, die Abtastauflösung vorzuwählen, anstatt sich auf dynamische Adaptivität zu verlassen?

Ich interessiere mich sowohl für den univariaten als auch für den mehrdimensionalen Fall.


3
Nur eine kleine Korrektur: Das Rechteck, das Trapez und die Simpsonsche Regel sind Newton-Cotes-Regeln, keine Gaußschen Quadraturen.
Pedro

Antworten:


20

Ich denke, das ist nicht ganz das, was Sie sich vorgestellt haben, aber der Vollständigkeit halber wollen wir mit einigen Grundlagen beginnen. Die meisten Quadraturformeln wie Newton-Cotes und Gauß auf der Idee , dass, um das Integral einer Funktion annähernd zu bewerten, können Sie die Funktion durch, zum Beispiel annähernd können ein Polynom , dass Sie dann genau integrieren:

abf(x)dxabjcjpj(x)dx=jcjabpj(x)dx.

Newton-Cotes und Gauss basieren auf der Lagrange-Interpolation , dh, Sie interpolieren die angegebene Funktion mit ihren Werten auf einer Menge von Knoten (die für Newton-Cotes einheitlich beabstandet und für Gauss in gewissem Sinne optimal gewählt sind). In diesem Fall c j = f ( x j ) , und die Integrale über die Polynom - Grundfunktionen Knoten p j sind genau die Quadraturgewichte.xjcj=f(xj)pj

Der gleiche Ansatz funktioniert mit Hermite-Interpolation , dh Interpolation unter Verwendung der Werte einer Funktion und ihrer Ableitungen bis zu einer bestimmten Reihenfolge auf einer Menge von Knoten. Nur bei Funktions- und Erstableitungswerten haben Sie (Es gibt eineMatlab-Implementierung

einbf(x)dxeinbjf(xj)pj(x)+f(xj)qj(x)dx=jf(xj)wj+f(xj)w¯j.
davon, wenn Sie sehen möchten, wie es funktioniert.)

Dies hängt mit einer Variante der Gauß-Quadratur zusammen, die als Gauß-Legendre-Quadratur bezeichnet wird und bei der die Knoten genau so ausgewählt werden, dass die Gewichte verschwinden (was eine weitere Erklärung für die Tatsache ist, dass die Gauß-Quadratur mit N Knoten genau 2 N ist - 1 ). Ich denke, dies beantwortet Ihre Frage im zweiten Absatz zumindest teilweise. Aus diesem Grund wird in der Regel anstelle der Hermite-Interpolation die Gauß-Quadratur verwendet, da Sie die gleiche Reihenfolge mit der gleichen Anzahl von Punkten erhalten, jedoch keine abgeleiteten Informationen benötigen.w¯jN2N-1

Bei mehrdimensionaler Quadratur besteht das Problem, dass die Anzahl der zu bewertenden Derivate (einschließlich gemischter Derivate) mit zunehmender Reihenfolge sehr schnell zunimmt.

Zurück zu Ihrer Frage: Eine einfache Möglichkeit, abgeleitete Informationen zu nutzen, besteht darin, eine Unterteilung Ihrer Integrationsdomäne zu verwenden und für jede Division eine separate Quadratur zu verwenden. Wenn Sie wissen, dass die Ableitungen Ihrer Funktion in einem Teil der Domäne groß sind, verwenden Sie entweder kleinere Domänen (im Endeffekt eine summierte Quadraturformel) oder eine höhere Quadraturordnung. Dies hängt mit der h- bzw. p-Adaptivität in Finite-Elemente-Methoden zusammen.


6

Es gibt eine Reihe von "korrigierten" Integrationsregeln, die Ableitungen der Endpunkte aufrufen. Ein einfaches Beispiel ist die korrigierte Trapezregel. Angenommen, wir möchten das Integral approximieren

einbf(x)dx.

nh=(b-ein)/n

T=h2(f(ein)+2f(ein+h)+2f(ein+2h)++2f(ein+(n-1)h)+f(b))

h2

T=T-h212(f(b)-f(ein))

erhöht die Genauigkeit erheblich. Zum Beispiel betrachten

ich=01e-x2dx

n=8ich

0,74682413281243

TT

0,7458656148457,0,74682363422375

beziehungsweise. Die Fehler sind

|ich-T|=9.5851796673207534×10-4

und

|ich-T|=4.9858868145236102×10-7

zeigt eine bemerkenswerte Steigerung der Genauigkeit. Es gibt weitere Korrekturen mit höheren Ableitungen oder ausgehend von anderen Newton-Cotes-Regeln oder Gaußschen Regeln.


5

Polynom×Gewichtsfunktion


4

Obwohl dieser Thread ziemlich alt ist, dachte ich, dass es nützlich sein könnte, einen Verweis auf ein von Fachleuten begutachtetes Papier zu haben, um einige allgemeine Quadraturregeln zu verallgemeinern.

Nenad Ujevic, "Eine Verallgemeinerung der modifizierten Simpsons Regel und Fehlergrenzen", ANZIAM Journal, Vol. 47, 2005.

http://journal.austms.org.au/ojs/index.php/ANZIAMJ/article/view/2/1268

Ich dachte, es wäre nützlich, eine gute Referenz zu geben, die frei zugänglich ist und Verweise auf andere Artikel enthält.

Wie Alasdair bereits erwähnt hat, kann das Einbeziehen von Ableitungen der Endpunkte die Genauigkeit erheblich verbessern. Zum Beispiel zeigten Ujevic und Roberts, dass das Hinzufügen der ersten Ableitungen zu Simpsons Regel den Fehler auf die sechste Ordnung im Gitterabstand reduziert, wohingegen es sich um die vierte Ordnung ohne die Ableitungen handelt. Das Ujevic-Papier zeigt, dass noch engere Fehlergrenzen gefunden werden können.

N. Ujevic und AJ Roberts, Eine korrigierte Quadraturformel und Anwendungen, ANZIAM J., 45 (E), (2004), E41 – E56. http://anziamj.austms.org.au/V45/E051

(Christian Clason schlug vor, dass ich einen Kommentar in die Antwort einfüge, weil er meinte, dass die Referenzen, die ich gebe, gut sind und sie möglicherweise verloren gehen, wenn Kommentare irgendwann gelöscht werden.)


Können Sie die im Artikel präsentierten Ergebnisse kommentieren?
Nicoguaro

Ich kann jetzt, dass ich genug Wiederholungspunkte habe! Ich dachte, es wäre nützlich, eine gute Referenz zu geben, die frei zugänglich ist und Verweise auf andere Artikel enthält. Wie Alasdair bereits erwähnt hat, kann das Einbeziehen von Ableitungen der Endpunkte die Genauigkeit erheblich verbessern. Zum Beispiel zeigten Roberts und Ujevic in Referenz 6 der Veröffentlichung, dass das Hinzufügen erster Ableitungen zu Simpsons Regel den Fehler auf die 6. Ordnung im Gitterabstand reduziert, wohingegen es sich um die 4. Ordnung ohne die Ableitungen handelt. Das Ujevic-Papier zeigt, dass noch engere Fehlergrenzen gefunden werden können.
Lysistrata

1
@ Lysistrata Das ist eine nette Referenz. Können Sie Ihre Kommentare in der Antwort selbst bearbeiten? Kommentare können verschwinden, und es wäre schade, sie zu verlieren.
Christian Clason
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.