Integración de Monte Carlo
![](http://upload.wikimedia.org/wikipedia/commons/thumb/b/b0/MonteCarloIntegrationCircle.png/220px-MonteCarloIntegrationCircle.png)
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.
Generalidades
Considere una función 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 es la variable aleatoria con , lo anterior es válido pues
por lo que
siendo
la función de densidad de .
Entonces para calcular la integral hay que calcular la esperanza de , lo que hace la Integración de Monte Carlo es estimarla.
Estimador
Lo que se busca es estimar entonces por la Ley fuerte de los Grandes Números
por lo que
Mientras más grande sea , es decir, si , más exacta será la aproximación entonces
Si en lugar de generar una muestra de se genera de entonces uno puede utilizar el Método de la transformada inversa para convertir la muestra al intervalo deseado. Si denota una muestra de tamaño de entonces para convertirla a una muestra de se considera por lo que
Algoritmo
Para evaluar una integral de la forma
se sigue el siguiente algoritmo
- Generar una muestra de tamaño de .
- Evaluar cada elemento de la muestra en la función .
- Calcular
Ejemplo
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 círculo)(Área del cuadrado unitario )=(Puntos en el de círculo)(Puntos en el cuadrado)
Para generar los puntos utilizamos dos sucesiones de números aleatorios y . 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.
Véase también
- Método de Monte Carlo
- Teorema del límite central
- Ley de los grandes números
- Distribución uniforme continua
- Números pseudoaleatorios
Datos: Q39879