Espero que este sea el lugar correcto para preguntar, si no es así, siéntase libre de moverlo a un foro más apropiado.
Hace tiempo que me pregunto cómo tratar las funciones integrables no cuadradas con la integración de Monte Carlo. Sé que MC todavía da una estimación adecuada, pero el error es irrealizable (¿divergente?) Para ese tipo de funciones.
Restringámonos a una dimensión. La integración de Monte Carlo significa que aproximamos la integral
usando la estimación
con puntos aleatorios distribuidos uniformemente. La ley de los grandes números asegura que .La varianza muestraE ≈ I
aproxima la varianza de la distribución inducida por . Sin embargo, si no es integrable al cuadrado, es decir, la integral de la función al cuadrado diverge, esto implica f f
lo que significa que también la varianza diverge.
Un ejemplo simple es la función
para el cual y .σ2=∫10dx
Si es finito, se puede aproximar el error de la media por , pero qué si no es integrable al cuadrado? E S f(x)