¿Puedo confiar en esta integral triple numérica de Matlab?

Gente de Ciencias Computacionales:

Me Iniciado esta cuestión en Matemáticas Pila de cambio y alguien comento que yo podría conseguir "mucho mejor" respuestas aquí:

Soy un novato en métodos numéricos y Matlab. Estoy tratando de evaluar la siguiente suma de dos integrales triples (obviamente, se puede escribir de manera más simple, pero aún así no se puede evaluar simbólicamente (?)). Tengo problemas para hacer que el funcione aquí, así que de mala gana lo partí en pedazos aquí: quiero encontrar la suma de $\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),$

\frac{2}{((1 / / 0,3) - 1)^{2}} (\int_{1}^{1 / / 0,3} \int_{1}^{r_{1}} \int_{r_{1} - r_{0 0}}^{r_{1} + r_{0 0}} F_{2} (r_{0 0}, r_{1}, t) Exp (- \frac{(0,3)^{2} t^{2}}{4 4}) re t re r_{0 0} re 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),$

dónde

F_{1} (r_{0 0}, r_{1}, t) = \frac{t^{2} r_{0 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}}$

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)}.$

EDITAR (2 de marzo de 2013): Alguien respondió que consiguió que Mathematica hiciera las integrales simbólicamente. Solo intenté hacer esto (con versiones simplificadas de las integrales) y Mathematica solo pudo hacer los dos exteriores del primero, y me detuve en el segundo. Agradecería algo de ayuda. Aquí esta lo que hice.:

Intenté evaluar

\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$ vía

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

y Mathematica regresa (Tuve problemas con $\LaTeX$ aquí porque el resultado es largo. Lo dividí en dos ecuaciones. Si alguien sabe una buena manera de mostrar esto, dígame):

\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.$

Entonces traté de evaluar

\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$

utilizando

Integrar [(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}]

justo ahora, y Mathematica no ha devuelto una respuesta después de aproximadamente media hora (pero estoy teniendo problemas de red informática en este momento, y pueden ser los culpables).

[FIN DE LA EDICIÓN DEL 2 DE MARZO]

Utilicé el comando "triplequad" de Matlab, sin opciones adicionales. Manejé los límites variables de integración por medio de funciones de cabecera, porque no conocía otra forma de hacerlo. Matlab me dio . $0.007164820144202$

Sé que Matlab es un buen software, pero he oído que las integrales triples numéricas son difíciles de hacer con precisión, y se supone que los matemáticos son escépticos, por lo que quiero alguna forma de verificar la precisión de esta respuesta. Las integrales dan el valor esperado de un determinado experimento (si alguien quiere, puedo editar esta pregunta para describir el experimento): implementé el experimento en Matlab usando números generados al azar de manera apropiada, un millón de veces, y promedié los resultados. Repetí este proceso cuatro veces. Aquí están los resultados (pido disculpas si he usado la palabra "prueba" de forma incorrecta):

Prueba 1: $0.007133292603256$

Prueba 2: $0.007120455071989$

Prueba 3: $0.007062595022049$

Prueba 4: $0.007154940168452$

Prueba 5: $0.007215000289130$

Aunque cada ensayo utilizó un millón de muestras, los valores de simulación solo coinciden en el primer dígito significativo. No están lo suficientemente cerca el uno del otro para determinar si la triple integral numérica es precisa.

Entonces, ¿alguien puede decirme si puedo confiar en el resultado de "triplequad" aquí, y bajo qué circunstancias uno puede confiar en él en general?

Una sugerencia que recibí en Math Stack Exchange fue probar otro software como Mathematica, Octave, Maple y SciPy. ¿Es este un buen consejo? ¿La gente realmente hace trabajo numérico en Mathematica y Maple? Octave es una especie de clon de Matlab, ¿puedo asumir que usa los mismos algoritmos de integración? Ni siquiera he oído hablar de SciPy antes y agradecería cualquier opinión al respecto.

ACTUALIZACIÓN: Alguien de Math Stack Exchange lo hizo en Maple y obtuvo . Eso está de acuerdo con tres cifras significativas. Eso es un buen signo. $0.007163085468$

Además, agradecería sugerencias sobre cómo ingresar expresiones largas de varias líneas en en Stack Exchange. ¿Puedes usar el entorno "alineado" aquí? Lo intenté y no pude hacerlo funcionar. $\LaTeX$

matlab Stefan Smith
fuente

Los resultados de su simulación son perfectamente consistentes con el valor numérico devuelto por Matlab: su media de es solo errores estándar menos de lo que Matlab devolvió. FWIW, Mathematica devuelve . También puede evaluar estas integrales simbólicamente en términos de polinomios y funciones de error.

0.00713726

$0.00713726$

- 1.11

$-1.11$

0.00716308537 \dots

$0.00716308537\ldots$

whuber

@whuber Gracias. Podría jurar que lo intenté simbólicamente en Maple y Maple no pudo hacerlo. Lo intentaré nuevamente en Maple, y si no funciona, lo intentaré en Mathematica. Por cierto, hice una integral similar en Maple, y obtuve una gran respuesta simbólica. Parecía ser una suma y diferencia de números muy grandes cuyo gran total era bastante pequeño. Sospecho que es probable que haya un error de redondeo en la respuesta final. En un problema como este, ¿debería usar la respuesta simbólica o simplemente hacer la integral numéricamente?

Stefan Smith

Las respuestas simbólicas tienen la ventaja de ser combinaciones de funciones que (a menudo) pueden calcularse eficientemente con precisión arbitraria. Por lo general, también, la solución simbólica también se presta a un cálculo rápido cuando los parámetros varían. Por estas razones, a menudo vale la pena buscar una solución simbólica.

whuber

@whuber: Intenté hacer algunas integrales esencialmente equivalentes (cambiar algunas de las constantes y eliminar algunas constantes multiplicativas) en Mathematica, y Mathematica solo pudo hacer las dos integraciones externas de la primera integral, y parece haberse estancado en la segunda. Publiqué mi código y resultados arriba.

Stefan Smith

Re Edición del 2 de marzo: al reducir la triple integral simbólicamente a una sola integral (en la primera mitad de sus integrales), ha logrado mucho. El integrando se comporta muy bien y puede integrarse numéricamente con una precisión extremadamente alta en una fracción de segundo.

whuber

¿Puedo confiar en esta integral triple numérica de Matlab?

Respuestas: