Kann ich diesem numerischen Tripelintegral von Matlab vertrauen?

Computational Science Personen:

Ich habe diese Frage ursprünglich bei Math Stack Exchange gepostet und jemand sagte, dass ich hier möglicherweise "viel bessere" Antworten bekomme:

Ich bin ein Anfänger in numerischen Methoden und Matlab. Ich versuche, die folgende Summe von zwei Dreifachintegralen zu bewerten (es kann offensichtlich einfacher geschrieben werden, aber Sie können es immer noch nicht symbolisch bewerten (?)). Ich habe Probleme, das zum zu bringen, deshalb habe ich es hier widerwillig in Stücke zerbrochen: Ich möchte die Summe von finden $\LaTeX$

\frac{2}{((1 / 0,3) - 1)^{2}} (\int_{1}^{1 / 0,3} \int_{1}^{r_{1}} \int_{0}^{r_{1} - r_{0}} F_{1} (r_{0}, r_{1}, t) \exp (- \frac{(0.3)^{2} t^{2}}{4}) d t d r_{0} d r_{1}),

$\frac{2}{((1/0.3) - 1)^2}\left(\int_1^{1/0.3}\int_1^{r_1}\int_0^{r_1-r_0}F_1(r_0,r_1,t)\exp(-\frac{(0.3)^2 t^2}{4})\,dt\,dr_0\,dr_1 \right),$

und

\frac{2}{((1 / 0.3) - 1)^{2}} (\int_{1}^{1 / 0.3} \int_{1}^{r_{1}} \int_{r_{1} - r_{0}}^{r_{1} + r_{0}} F_{2} (r_{0}, r_{1}, t) \exp (- \frac{(0.3)^{2} t^{2}}{4}) d t d r_{0} d r_{1}),

$\frac{2}{((1/0.3) - 1)^2}\left(\int_1^{1/0.3}\int_1^{r_1}\int_{r_1-r_0}^{r_1+r_0} F_2(r_0,r_1,t)\exp(-\frac{(0.3)^2 t^2}{4})\,dt\,dr_0\,dr_1 \right),$

F_{1} (r_{0}, r_{1}, t) = \frac{t^{2} r_{0}^{3} * (0.3)^{3}}{2 r_{1}^{3} \sqrt{π}}

$F_1(r_0,r_1,t)=\frac{t^2 r_0^3*(0.3)^3}{2r_1^3\sqrt{\pi}}$

und

F_{2} (r_{0}, r_{1}, t) = \frac{(0,3)^{3} π^{3 / 2} (r_{0} + r_{1} - t)^{4} (t^{2} + 2 t (r_{0} + r_{1}) - 3 (r_{1} - r_{0})^{2})^{2}}{288 (\frac{4}{3} π r_{0}^{3}) (\frac{4}{3} π r_{1}^{3})} .

$F_2(r_0,r_1,t)=\frac{(0.3)^3\pi^{3/2}(r_0+r_1-t)^4 (t^2+2t(r_0+r_1)-3(r_1-r_0)^2)^2}{288(\frac{4}{3}\pi r_0^3)(\frac{4}{3}\pi r_1^3)}.$

BEARBEITEN (2. März 2013): Jemand antwortete, dass Mathematica die Integrale symbolisch ausführen soll. Ich habe gerade versucht, dies zu tun (mit vereinfachten Versionen der Integrale), und Mathematica konnte nur die beiden äußeren der ersten ausführen und kam bei der zweiten ins Stocken. Ich würde mich über Hilfe freuen. Folgendes habe ich getan:

Ich habe versucht zu bewerten

\int_{1}^{2} \int_{1}^{r_{2}} \int_{0}^{r_{2} - r_{1}} \frac{r_{1}^{3} t^{2} \exp (- t^{2})}{r_{2}^{3}} d t d r_{1} d r_{2}

$\int_1^2 \int_1^{r_2} \int_0^{r_2-r_1} \frac{r_1^3 t^2 \exp(-t^2)}{r_2^3}\,dt\,dr_1\,dr_2$ über

Integriere [r1 ^ 3 / r2 ^ 3 * t ^ 2 * Exp (-t ^ 2), {t, 0, r2 - r1}, {r1, 1, r2}, {r2, 1, 2}]

und Mathematica kehrt zurück (Ich hatte Probleme mit $\LaTeX$ hier, weil das Ergebnis lang ist. Ich habe es in zwei Gleichungen aufgeteilt. Wenn jemand eine gute Möglichkeit kennt, dies anzuzeigen, sag es mir bitte):

\int_{1}^{2} \frac{1}{64 r 2^{2}} e^{- 1 - r 2^{2}} (2 e^{2 r 2} (25 + r 2 (19 + 2 r 2 (1 + r 2))) -

$\int_1^2 \frac{1}{64r2^2} e^{-1-r2^2}(2e^{2r2}(25+r2(19+2r2(1+r2)))-$

e^{1 + r 2^{2}} (32 r 2 (2 + r 2^{2})) + \sqrt{π} (11 + 4 r 2^{2} (9 + r 2^{2})) Erf [1 - r 2]) d r 2.

$e^{1+r2^2}(32r2(2+r2^2)) +\sqrt{\pi}(11+4r2^2(9+r2^2))\operatorname{Erf}[1-r2])\,dr2.$

Dann habe ich versucht zu bewerten

\int_{1}^{2} \int_{1}^{r_{2}} \int_{r_{2} - r_{1}}^{r_{2} + r_{1}} \dots

$\int_1^2\int_1^{r_2}\int_{r_2-r_1}^{r_2+r_1} \ldots \qquad \qquad \qquad$

\dots \frac{\exp (- t^{2}) (r_{1} + r_{2} - t)^{4} (t^{2} + 2 t (r_{1} + r_{2}) - 3 (r_{2} - r_{1})^{2})^{2}}{r_{1}^{3} r_{2}^{3}} d t d r_{} d r_{2}

$\ldots\frac{\exp(-t^2)(r_1+r_2-t)^4(t^2+2t(r_1+r_2)-3(r_2-r_1)^2)^2}{r_1^3 r_2^3}\,dt\,dr_\,dr_2$

mit

Integriere [(r1 + r2 - t) ^ 4 * (t ^ 2 + 2 * t * (r1 + r2) - 3 * (r2 - r1) ^ 2) ^ 2 * Exp [-t ^ 2] / r1 ^ 3 / r2 ^ 3, {r2, 1, 2}, {r1, 1, r2}, {t, r2-r1, r2 + r1}]

Gerade jetzt, und Mathematica hat nach ungefähr einer halben Stunde keine Antwort zurückgegeben (aber ich habe im Moment Probleme mit dem Computernetzwerk und sie können schuld sein).

[ENDE DES 2. MÄRZES]

Ich habe den Matlab-Befehl "triplequad" ohne zusätzliche Optionen verwendet. Die variablen Integrationsgrenzen habe ich mit Hilfe von heaviside-Funktionen behandelt, weil ich keine andere Möglichkeit dazu kannte. Matlab hat mir . $0.007164820144202$

Ich weiß, dass Matlab eine gute Software ist, aber ich habe gehört, dass es schwierig ist, numerische Dreifachintegrale genau zu machen, und Mathematiker sollten skeptisch sein. Deshalb möchte ich einen Weg finden, um die Richtigkeit dieser Antwort zu überprüfen. Die Integrale geben den erwarteten Wert eines bestimmten Experiments an (wenn jemand möchte, kann ich diese Frage bearbeiten, um das Experiment zu beschreiben): Ich habe das Experiment in Matlab mit entsprechend zufällig generierten Zahlen millionenfach implementiert und die Ergebnisse gemittelt. Ich habe diesen Vorgang viermal wiederholt. Hier sind die Ergebnisse (ich entschuldige mich, wenn ich das Wort "Studie" falsch verwendet habe):

Versuch 1: $0.007133292603256$

Versuch 2: $0.007120455071989$

Versuch 3: $0.007062595022049$

Versuch 4: $0.007154940168452$

Versuch 5: $0.007215000289130$

Obwohl für jeden Versuch eine Million Proben verwendet wurden, stimmen die Simulationswerte nur in der ersten signifikanten Stelle überein. Sie sind nicht nahe genug beieinander, um zu bestimmen, ob das numerische Dreifachintegral korrekt ist.

Kann mir also jemand sagen, ob ich dem Ergebnis von "Triplequad" hier vertrauen kann und unter welchen Umständen man ihm überhaupt vertrauen kann?

Ein Vorschlag, den ich bei Math Stack Exchange erhalten habe, war, andere Software wie Mathematica, Octave, Maple und SciPy auszuprobieren. Ist das ein guter Rat? Arbeiten Menschen in Mathematica und Maple tatsächlich numerisch? Octave ist eine Art Matlab-Klon. Kann ich also davon ausgehen, dass er dieselben Integrationsalgorithmen verwendet? Ich habe noch nie von SciPy gehört und würde mich über jede Meinung dazu freuen.

UPDATE: Jemand von Math Stack Exchange hat es in Maple gemacht und . Das ist Übereinstimmung mit drei bedeutenden Figuren. Das ist ein gutes Zeichen. $0.007163085468$

Ich würde mich auch über Vorschläge freuen, wie Sie einen langen, mehrzeiligen Ausdruck in in Stack Exchange eingeben können . Können Sie hier die "ausgerichtete" Umgebung verwenden? Ich habe es versucht und konnte es nicht zum Laufen bringen. $\LaTeX$

matlab Stefan Smith
quelle

Ihre Simulationsergebnisse stimmen perfekt mit dem von Matlab zurückgegebenen numerischen Wert : Ihr Mittelwert von ist nur Standardfehler weniger als der von Matlab zurückgegebene Wert . FWIW, Mathematica gibt . Sie kann diese Integrale auch symbolisch nach Polynomen und Fehlerfunktionen auswerten.

0.00713726

$0.00713726$

- 1.11

$-1.11$

0.00716308537 \dots

$0.00716308537\ldots$

Whuber

@whuber Danke. Ich könnte schwören, ich habe es symbolisch in Maple versucht und Maple konnte es nicht. Ich werde es in Maple erneut versuchen, und wenn es nicht funktioniert, werde ich es in Mathematica versuchen. Übrigens habe ich ein ähnliches Integral in Maple erstellt und eine große symbolische Antwort erhalten. Es schien eine Summe und Differenz von sehr großen Zahlen zu sein, deren Gesamtsumme ziemlich klein war. Ich vermute Rundungsfehler war wahrscheinlich in der endgültigen Antwort. Sollten Sie in einem Problem wie diesem die symbolische Antwort verwenden oder das Integral nur numerisch ausführen?

Stefan Smith

Symbolische Antworten haben den Vorteil, dass sie Kombinationen von Funktionen darstellen, die (oft) mit beliebiger Genauigkeit effizient berechnet werden können. Auch in der Regel bietet sich die symbolische Lösung bei Variation von Parametern für eine schnelle Neuberechnung an. Aus diesen Gründen lohnt es sich oft, eine symbolische Lösung zu suchen.

Whuber

@whuber: Ich habe in Mathematica versucht, einige im Wesentlichen äquivalente Integrale zu erstellen (einige der Konstanten zu ändern und einige multiplikative Konstanten zu entfernen), und Mathematica konnte nur die beiden äußeren Integrationen des ersten Integrals ausführen und scheint beim zweiten festgefahren zu sein. Ich habe meinen Code und die Ergebnisse oben gepostet.

Stefan Smith

Re The March 2 edit: Indem Sie das Dreifachintegral symbolisch auf ein einziges Integral (in der ersten Hälfte Ihrer Integrale) reduzieren, haben Sie viel erreicht. Der Integrand verhält sich sehr gut und kann innerhalb von Sekundenbruchteilen mit extrem hoher Genauigkeit numerisch integriert werden.

Whuber

Kann ich diesem numerischen Tripelintegral von Matlab vertrauen?

Antworten: