La complexité de problèmes réels rencontrés en économie, en biologie, en physique, en finance, voire en mathématiques, fait que les méthodes déterministes s’avèrent parfois insuffisantes. Le recours à des méthodes aléatoires numériques utilisant l’informatique se révèle alors d’un secours incroyablement utile. Parmi celles-ci, on retrouve la méthode dite de Monte-Carlo, qui fut mise au point peu avant 1950 par le physicien Nicholas Constantine Metropolis. Elle se montre efficace en mathématiques notamment pour le calcul d’intégrales à dimensions multiples, portant sur de très larges domaines, et fournit en pratique des résultats numériques plus précis que les méthodes classiques déterministes. Voyons en quoi consiste son principe dans un cas élémentaire.
Reproduire artificiellement le hasard
Pour calculer l’intégrale , on va considérer le réel x (qui varie entre 0 et 1) comme une variable aléatoire pour laquelle il convient de construire une suite de valeurs x1, x2… xn indépendantes, distribuées uniformément sur l’intervalle [0, 1[. On estime alors l’intégrale définie I par la relation suivante (dans laquelle on reconnaît la notion de moyenne) :
Les lois faibles des grands nombres (ou le théorème central limite) donnent une idée assez précise de la vitesse de convergence de la méthode, qui ...
Lire la suite