¿Cómo puedo aproximar una integral impropia?

13

Tengo una función f(x,y,z) tal que
R3f(x,y,z)dV
es finita, y quiero aproximar esta integral.

Estoy familiarizado con las reglas de cuadratura y las aproximaciones de las integrales de Monte Carlo, pero veo algunas dificultades para implementarlas en un dominio infinito. En el caso de Monte Carlo, ¿cómo se hace para muestrear una región infinita (especialmente si las regiones que contribuyen más significativamente a la integral son desconocidas)? En el caso de la cuadratura, ¿cómo encuentro los puntos óptimos? ¿Debo simplemente arreglar una región arbitrariamente grande centrada alrededor del origen y aplicar reglas de cuadratura dispersas? ¿Cómo puedo aproximarme a esta integral?

Paul
fuente

Respuestas:

20

En una dimensión, puede asignar su intervalo infinito a un intervalo finito utilizando la integración por sustitución, p. Ej.

abf(x)dx=u1(a)u1(b)f(u(t))u(t)dt

Donde es alguna función que se dispara al infinito en un rango finito, por ejemplo, tan ( x ) :u(x)tan(x)

f(x)dx=2π/2π/2f(tan(t))1cos(2t)+1dt

Luego puede usar cualquier rutina de cuadratura numérica regular para la integral finita modificada.

La sustitución de múltiples variables es un poco más complicada, pero aquí se describe bastante bien .

Pedro
fuente
Eso es muy interesante ... ¡Nunca consideré la posibilidad de sustitución! Pero, ¿la elección de la función tiene algún efecto en la precisión de la aproximación? u(t)
Paul
@Paul: ¡Sí, definitivamente! La función debe ser lo más suave posible para mantener f ( u ( t ) ) lo más suave posible, permitiendo así una integración más precisa. u(t)f(u(t))
Pedro
Eso es cierto, pero lo que tenía en mente era la velocidad a la que u (t) converge al infinito. ¿Eso también afecta la precisión?
Paul
1
@Paul: No sé si entiendo tu pregunta correctamente, pero la función tiene que terminar en el infinito en un momento u otro. Si se toma su tiempo y luego crece bruscamente, esto introducirá algunos gradientes grandes en , lo que dificulta la integración y, por lo tanto, podría afectar la precisión. f(u(t))
Pedro
1
Su derivada para la tangente estaba equivocada; Lo arreglé.
JM
11

La forma estándar de hacerlo es extraer de la expresión para un prefactor exponencial, transformarlo en e - x 2 , y luego usar las reglas de la cuadratura gaussiana (o Gauss Kronrod) con esto como un peso. Si f es suave, esto generalmente da excelentes resultados.f(x)ex2f

En , lo mismo funciona con el peso e - | x | 2 , y se pueden encontrar fórmulas de cubicación apropiadas, por ejemplo, en el libro de Engels, cuadratura numérica y cubicación.R3e|x|2

Las fórmulas en línea están en http://nines.cs.kuleuven.be/ecf/

Arnold Neumaier
fuente
2
Eso funciona bien si su integrando es aproximadamente exp (-x ^ 2). Si su integrando es aproximadamente normal, pero está centrado lejos del origen, este enfoque puede funcionar mal.
John D. Cook
1
@ JohnD.Cook: Por eso escribí '' extrae un prefactor exponencial, transfórmalo a '', que generalmente implica una transformación lineal, combinando una traslación que mueve el centro al origen, y rotaciones y escalas para hacer el el nivel establece aproximadamente esférico. La función en sí puede estar bastante lejos de lo normal. ex2
Arnold Neumaier
7

Para la cuadratura unidimensional, puede consultar el libro en Quadpack (un viejo dorado pero aún muy relevante en la cuadratura unidimensional) y las técnicas utilizadas en el algoritmo QAGI, un integrador automático para un rango infinito.

Otra técnica es la fórmula de cuadratura doble exponencial, bien implementada por Ooura para un intervalo infinito .

Para la cubicación, puede consultar la Enciclopedia de fórmulas de cubicación de Ronald Cools.

GertVdE
fuente
2
Tenga en cuenta que la cuadratura doble exponencial es esencialmente un método de sustitución; realiza una sustitución que transforma su integral de rango infinito en otra integral de rango infinito cuya tasa de descomposición es, bueno, doble exponencial ...
JM
1
@JM Correcto. Y lo hace para obtener lo mejor de la fórmula de suma de Euler-Mclaurin para la regla del trapecio, al igual que la transformación IMT y la transformación TANH. Aquí
GertVdE
6

f(x)f~(x)f~f

f(x)f~(x)=ex2p(x)p(x)f(x)ex2f~(x)dx

Wolfgang Bangerth
fuente
4

Si desea utilizar la integración de Monte Carlo, puede comenzar utilizando el muestreo de importancia con un muestreador que se aproxima aproximadamente a su integrando. Cuanto mejor coincida su muestra con su integrando, menor será la variación en sus estimaciones integrales. No importa que su dominio sea infinito siempre que su muestra tenga el mismo dominio.

John D. Cook
fuente