Análise de Monte Carlo
Introdução
Em geral
O método Monte Carlo[1] é um método estatístico não determinístico ou numérico, usado para aproximar expressões matemáticas complexas cuja avaliação precisa é cara. O método recebeu esse nome em referência ao Cassino de Monte Carlo (Mônaco) por ser “a capital do jogo de azar”, já que a roleta é um simples gerador de números aleatórios. O nome e o desenvolvimento sistemático dos métodos de Monte Carlo remontam aproximadamente a 1944 e foram bastante aprimorados com o desenvolvimento do computador.
O uso dos métodos de Monte Carlo em pesquisas científicas teve origem durante o desenvolvimento da bomba atômica na Segunda Guerra Mundial, no Laboratório Nacional de Los Alamos, nos Estados Unidos. Nesse contexto, os pesquisadores utilizaram simulações probabilísticas para estudar fenômenos hidrodinâmicos relacionados à difusão de nêutrons em materiais de fissão, um processo comportamental inerentemente aleatório.
Atualmente, os métodos de Monte Carlo são aplicados em uma ampla variedade de campos, desde física nuclear e estatística computacional até computação gráfica, onde formam a base de numerosos algoritmos de ray tracing usados para a geração de imagens tridimensionais realistas.
Na primeira fase destas investigações, John von Neumann e Stanislaw Ulam refinaram esta roleta e os métodos de divisão de tarefas. No entanto, o desenvolvimento sistemático dessas ideias teve que esperar pelo trabalho de Harris e Herman Kahn em 1948. Aproximadamente no mesmo ano, Enrico Fermi, Nicholas Metropolis e Ulam obtiveram estimadores para os valores característicos da equação de Schrödinger para captura de nêutrons no nível nuclear usando este método.
O método de Monte Carlo fornece soluções aproximadas para uma ampla variedade de problemas matemáticos, possibilitando a realização de experimentos com amostragem de números pseudoaleatórios em um computador. O método é aplicável a qualquer tipo de problema, seja estocástico ou determinístico. Ao contrário dos métodos numéricos que dependem de avaliações em N pontos em um espaço M-dimensional para produzir uma solução aproximada, o método de Monte Carlo tem um erro absoluto de estimativa que diminui conforme o teorema do limite central.