*** Welcome to piglix ***

Monte Carlo integration


In mathematics, Monte Carlo integration is a technique for numerical integration using random numbers. It is a particular Monte Carlo method that numerically computes a definite integral. While other algorithms usually evaluate the integrand at a regular grid, Monte Carlo randomly choose points at which the integrand is evaluated. This method is particularly useful for higher-dimensional integrals.

There are different methods to perform a Monte Carlo integration, such as uniform sampling, stratified sampling, importance sampling, Sequential Monte Carlo (a.k.a. particle filter), and mean field particle methods.

In numerical integration, methods such as the Trapezoidal rule use a deterministic approach. Monte Carlo integration, on the other hand, employs a approach: each realization provides a different outcome. In Monte Carlo, the final outcome is an approximation of the correct value with respective error bars, and the correct value is within those error bars.

The problem Monte Carlo integration addresses is the computation of a multidimensional definite integral

where Ω, a subset of Rm, has volume

The naive Monte Carlo approach is to sample points uniformly on Ω: given N uniform samples,

I can be approximated by

This is because the law of large numbers ensures that

Given the estimation of I from QN, the error bars of QN can be estimated by the sample variance using the unbiased estimate of the variance.


...
Wikipedia

...