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
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.
{\displaystyle f:[a,b]\to \mathbb {R} }
que es integrable en el intervalo abierto
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
, lo que hace la Integración de Monte Carlo es estimarla.
Lo que se busca es estimar
entonces por la Ley fuerte de los Grandes Números por lo que Mientras más grande sea
, más exacta será la aproximación entonces Si en lugar de generar una muestra de
entonces uno puede utilizar el Método de la transformada inversa para convertir la muestra al intervalo deseado.
denota una muestra de tamaño
entonces para convertirla a una muestra de
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.
(Área del cuadrado unitario
Para generar los puntos utilizamos dos sucesiones de números aleatorios
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
, 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.
