Método Montecarlo y los físico

Mi hermano es físico y actualmente está dando un curso de Física Computacional con sus amigos. Me plático sus inquietudes del método Montecarlo y se las respondí aunque al final no se las expusieron a sus alumnos.

El Método Montecarlo es un método numérico para simular sistemas físicos o matemáticos la virtud de esté método radica en un eficacia mayor a los usuales debido que hace uso de métodos aleatorios.

Uno de los principales usos de este método es la integración consiste en tratar de obtener el valor

para lo cual se generan variables aleatorias independientes con densidad , de donde se sigue que y por lo tanto por la ley de los Grandes Números se sigue que

,

por lo cual si quiero obtener el valor basta con que promedie un número suficientemente grande de variables; de donde se desprende la primera pregunta y la más importante para entregar resultados: ¿qué tan grande debe ser para obtener un resultado bueno?. Es en este sentido que me planteo su duda mi hermano y es aquí dónde los físicos que conozco no saben la respuesta correcta.

El primer acercamiento a esta pregunta, según los físicos, es calcular la varianza de de donde se tiene que

,

de donde implican que lo único importante es hacer sufiecientemente grande para que , pero esto no siempre es suficiente para tener una buena aproximación ya que la varianza de una variable no siempre nos dice cosas sobre como es la variable aleatoria (vease casos en los que es inútil la desigualdad de Chebyshev), es aquí donde el otro gran teorema de la probabilidad entra a resolver esta duda.

El Teorema del Límite Central nos dice que

,

esto es

donde es la distribución normal. Pero qué quiere decir esto al problema de decir cual es la que debemos poner para obtener una buena estimación.

Notemos que es nuestro error en la aproximación por lo cual usando el Teorema del Límite Central se tiene que

.

Aplicando esto se tiene por ejemplo que si tengo quiero que con un 95% de eficacia el error sea menor que 1/1000 tenemos

de donde, usando un tabla para , se tiene que y ahora se tiene que se tiene que elegir tal que , y es así como se resuelve la pregunta tan molesta de que tan grande debe ser $n$. Como se puede ver también en este método además del error de la aproximación se tiene un error de confiabilidad.

Queda un problema en el aire el cual ya no discutiré pero debe de tomarse en cuenta este es: ¿cómo calcular ?, a veces será más fácil que calcular la esperanza pero la gran mayoría de las veces será igual o más complicado que calcular esta.

Me queda una duda que esperaría que me resolviera algún físico que me lea: ¿La computación cuántica, lo que eso quiera decir, resolvería el problema de generación de aleatoriedad en las computadoras?. Si la respuesta fuera afirmativa creo este tipo de algoritmos serían de una aplicación perfecta y con una eficacia óptima.

Anuncio publicitario