Estoy bastante seguro de que entiendo cómo funciona la integración de Monte Carlo, pero no entiendo la formulación de cómo se usa para estimar Pi. Voy por el procedimiento descrito en la 5ta diapositiva de esta presentación http://homepages.inf.ed.ac.uk/imurray2/teaching/09mlss/slides.pdf
Entiendo los pasos preliminares. Pi es igual a 4 veces el área de un cuarto del círculo unitario. Y el área del cuarto superior derecho del círculo unitario centrado en (0,0) es equivalente a la integral de la curva que es el cuarto superior derecho del círculo unitario en y 0 < y < 1 .
Lo que no entiendo es cómo es esta integral
donde se distribuye uniformemente en el cuadrado de la unidad alrededor del cuarto de círculo (es decir, siempre es igual a 1 si 0 < x < 1 y 0 < y < 1 y 0 de lo contrario). Esto significaría que
I ( ( x 2 + y 2 ) < 1 ) P ( x , y )
es la función que es el cuadrante superior derecho del círculo unitario en 0 < x < 1 y
pero no entiendo cómo esto es cierto ya que la función del indicador solo puede ser 1 o 0. Entiendo que probablemente esté escrito de esta manera para facilitar el muestreo de Monte Carlo (es decir, es una expectativa, así que solo muestree de P ( x , y ) y obtenga el promedio de las muestras aplicadas a I ( ( x 2 + y 2 ) < 1 ) ) pero simplemente no tiene sentido intuitivo para mí por qué esa integral representa el área bajo esa curva.
¿Podría alguien proporcionar una explicación intuitiva de esto? ¿Quizás mostrar cómo se derivó esa integral paso a paso?
EDITAR:
Pude obtener una mejor comprensión al relacionar las expectativas con un área. Lo explicaré aquí en caso de que ayude a alguien. Primero comience relacionando Pi con el área del cuadrante superior derecho del círculo unitario
Luego colocamos el cuadrante superior derecho en el cuadrado de la unidad. Y bajo una distribución uniforme sobre el cuadrado de la unidad, el área del cuadrante del círculo es proporcional a la probabilidad de obtener una muestra de él. Se deduce que la siguiente igualdad es válida
Y sustituyendo en la ecuación original
Entonces lo entendí relacionando el área con una probabilidad y luego relacionando esa probabilidad con una expectativa que es equivalente a la integral. Avísame si he cometido algún error.