Integración de Monte Carlo

(Redirigido desde «Integración de Montecarlo»)

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 editar

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 editar

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 editar

Para evaluar una integral de la forma

 

se sigue el siguiente algoritmo

  1. Generar una muestra de tamaño   de  .
  2. Evaluar cada elemento de la muestra en la función  .
  3. Calcular
 

Ejemplo editar

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 editar