Integración de Monte Carlo

En matemáticas, se le conoce como integración de Monte Carlo a un método que utiliza números aleatorios para estimar el valor de una integral definida, este método es muy utilizado para evaluar integrales múltiples, es decir, integrales de la forma

siendo ${displaystyle Omega subset mathbb {R} ^{m}}$.

La integración de Monte Carlo forma parte de una familia de algoritmos llamados genéricamente métodos de Monte Carlo, estos algoritmos utilizan números aleatorios para resolver diferentes tipos de problemas matemáticos y reciben su nombre debido al casino de Monte Carlo.

Considere una función ${displaystyle f:[a,b] o mathbb {R} }$ que es integrable en el intervalo abierto ${displaystyle (a,b)}$ y suponga que se desea evaluar la integral

siendo esta una integral complicada de evaluar de forma analítica.

Esta integral puede expresarse en términos de la esperanza de cierta variable aleatoria como sigue

donde ${displaystyle f(X)}$ es la variable aleatoria con ${displaystyle Xsim operatorname {U} (a,b)}$, lo anterior es válido pues

por lo que

siendo

la función de densidad de ${displaystyle Xsim operatorname {U} (a,b)}$.

Entonces para calcular la integral hay que calcular la esperanza de ${displaystyle f(X)}$, lo que hace la Integración de Monte Carlo es estimarla.

Lo que se busca es estimar ${displaystyle operatorname {E} [f(X)]}$ entonces por la Ley fuerte de los Grandes Números

por lo que

Mientras más grande sea ${displaystyle n}$, es decir, si ${displaystyle n o infty }$, más exacta será la aproximación entonces

Si en lugar de generar una muestra de ${displaystyle Xsim operatorname {U} (a,b)}$ se genera de ${displaystyle Xsim operatorname {U} (0,1)}$ entonces uno puede utilizar el Método de la transformada inversa para convertir la muestra al intervalo deseado. Si ${displaystyle {x_{i}}_{i=1}^{n}}$ denota una muestra de tamaño ${displaystyle n}$ de ${displaystyle Xsim operatorname {U} (0,1)}$ entonces para convertirla a una muestra de ${displaystyle Xsim operatorname {U} (a,b)}$ se considera ${displaystyle a+bx_{i}}$ por lo que

Para evaluar una integral de la forma

se sigue el siguiente algoritmo

A través de un ejemplo se ilustrará el método. El proceso consiste en calcular el área encerrada por una línea cerrada cualquiera que está incluida en un cuadrado de lado unitario (y área unitaria).

Al generar puntos al azar (mediante dos números aleatorios) se calcula la fracción que se establece entre la cantidad de puntos que caen dentro del área asociada a la curva y la cantidad total de puntos (o puntos en el cuadrado).

Supongamos que el área a calcular es un cuarto de círculo, de radio unitario, que está dentro de un cuadrado de lado unitario. La fracción será:

(Área del ${displaystyle {}^{1}/_{4}}$ círculo)${displaystyle /}$(Área del cuadrado unitario ${displaystyle [0,1] imes [0,1]}$)=(Puntos en el ${displaystyle {}^{1}/_{4}}$ de círculo)${displaystyle /}$(Puntos en el cuadrado)${displaystyle ={frac {pi r^{2}}{4}}}$

Para generar los puntos utilizamos dos sucesiones de números aleatorios ${displaystyle R_{1}}$ y ${displaystyle R_{2}}$. Si queremos saber si un punto pertenece al cuarto de círculo, establecemos, a partir de la ecuación del círculo, la condición de pertenencia:

Si se verifica la relación anterior, el punto pertenece al cuarto de círculo (y al cuadrado). De lo contrario pertenecerá solo al cuadrado.

Para el caso de una simulación con 25 pares de números aleatorios, es decir, para 25 puntos generados, nos dará una fracción tal como ${displaystyle {frac {21}{25}}=0,84}$, mientras que el área buscada será:

La precisión del método se mejora utilizando una gran cantidad de simulaciones, siendo el error del orden del 0,001 cuando se emplean unos 15.000 puntos simulados.