Metodos Numericos

INTEGRACIÓN X MÉTODOS NUMÉRICOS

Cálculo de integrales por Monte-Carlo Los métodos de integración de Montecarlo son algoritmos para encontrar una evaluación aproximada de una integral definida, normalmente de integrales múltiples.

1.- Sea la integral definida en un intervalo de , siendo una función de mucha dificultad para integrar; se hace el siguiente cambio de variable , y se sustituye en la función a integrar definida quedando una función a integral definida en un intervalo de la forma , con .

2.- Luego se generan las variables aleatorias , de veces se quiera por medio de la siguiente fórmula (generador congruencial de números pseudo-aleatorios): donde es la semilla, es el incremento, es el multiplicativo y es el módulo; con . Estos valores se toman al azar y se van obteniendo los , las veces que se requiera.

3.- Cuando tengamos los se comienza a generar los por esta fórmula , donde inicialmente , es decir , , así sucesivamente. Los números obtenidos son números aleatorios en un intervalo [0,1], que se van a sustituir en la fórmula siguiente:

Mientras más grande sea , más exacta es la aproximación. Es decir:

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 relación pitagórica, 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á sólo 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

...