Voir la source de Méthode de Monte-Carlo

{{Voir homonymes|Monte-Carlo (homonymie)}}
{{confusion| algorithme de Monte-Carlo}}
{{sources à lier|date=février 2015}}

Une '''méthode de Monte-Carlo''', ou '''méthode Monte-Carlo''', est une méthode [[algorithmique]] visant à calculer une valeur numérique approchée en utilisant des [[Processus stochastique|procédés aléatoires]], c'est-à-dire des techniques probabilistes.

Les méthodes de Monte-Carlo sont particulièrement utilisées pour calculer des [[Intégration (mathématiques)|intégrales]] en dimensions plus grandes que 1 (en particulier, pour calculer des [[Surface (géométrie analytique)|surface]]s et des [[volume]]s). Elles sont également couramment utilisées en [[physique des particules]], où des [[simulation de phénomènes|simulations]] probabilistes permettent d'estimer la forme d'un signal ou la sensibilité d'un détecteur. La comparaison des données mesurées à ces simulations peut permettre de mettre en évidence des caractéristiques inattendues, par exemple de nouvelles particules.

La méthode de simulation de Monte-Carlo permet aussi d'introduire une approche statistique du [[Risque financier|risque dans une décision financière]]. Elle consiste à isoler des variables clés du projet, telles que le [[chiffre d'affaires]] ou la marge, et à leur affecter une [[Loi de probabilité|loi de probabilités]]. Pour chacun de ces facteurs, un grand nombre de tirages aléatoires, suivant les lois de probabilité déterminées précédemment, est effectué, afin de trouver la probabilité d'occurrence de chacun des résultats. À titre d'exemple, le choix de mode de gestion d'une collectivité territoriale dans le cadre d'un [[partenariat public-privé]] (PPP) peut s'analyser par la méthode de Monte-Carlo, afin de prendre en compte la répartition des risques entre acteurs publics et privés. On parle alors de « risques valorisés » ou « valeurs à risque ».

Le véritable développement des méthodes de Monte-Carlo s'est effectué sous l'impulsion de [[John von Neumann]] et [[Stanislaw Marcin Ulam|Stanislaw Ulam]] notamment, lors de la [[Seconde Guerre mondiale]], et des recherches sur la fabrication de la [[Bombe A|bombe atomique]]. Ils ont en particulier utilisé ces méthodes probabilistes pour résoudre des [[Équation aux dérivées partielles|équations aux dérivées partielles]] dans le cadre de la ''[[Monte-Carlo N-Particle transport]]'' (MCNP).

Le nom de ces méthodes, qui fait allusion aux [[jeu de hasard|jeux de hasard]] pratiqués au [[casino de Monte-Carlo]], a été inventé en 1947 par [[Nicholas Metropolis]]<ref>{{Article|auteur=Nicholas Metropolis|texte=http://library.lanl.gov/cgi-bin/getfile?00326866.pdf|titre=The Beginning of the Monte Carlo Method|journal=Los Alamos Science|date=1987|no=15|pages=125-130}}.</ref>, et apparaît pour la première fois en 1949 dans un article coécrit avec [[Stanislaw Marcin Ulam|Stanislaw Ulam]]<ref>{{Article|doi=10.2307/2280232|url=http://www.chemie.unibas.ch/~meuwly/download/ulam.mc.pdf|date=septembre 1949|vol=44|auteur=Nicholas Metropolis et Stanislaw Ulam|titre=The Monte Carlo Method|journal=Journal of the [[American Statistical Association]]|no=247|pages=335-341}}.</ref>.

== Théorie ==
De manière générale, le problème que l'on cherche à résoudre par les méthodes de Monte-Carlo est celui de l'estimation de l'[[Espérance mathématique|espérance]] d'une [[variable aléatoire]] <math>X</math>, que l'on note généralement <math>\mathbb E[X]</math>.

La méthode de Monte-Carlo la plus simple consiste à générer un échantillon de variables aléatoires [[Variables indépendantes et identiquement distribuées|indépendantes et identiquement distribuées]] (iid) <math>X_1,\dots,X_n</math> suivant la même loi que <math>X</math>. Ensuite on estime l'espérance <math>\mathbb E[X]</math> avec l'[[Estimateur (statistique)|estimateur]], dit de Monte-Carlo ou encore de la moyenne empirique, 

<math display="block">\mu_n := \frac 1n \sum_{k=1}^n X_k.</math>

Ainsi <math>\mu_n</math> est censé être une bonne approximation de la valeur que l'on recherche, à savoir <math>\mathbb E[X]</math>. Cet estimateur est non [[Biais (statistique)|biaisé]] dans le sens où <math>\mathbb E[\mu_n] =\mathbb E[X]</math>. Ce qui justifie la pertinence de cet estimateur est la [[loi forte des grands nombres]] qui nous dit que, si <math>X</math> admet une espérance finie (ce qui est en pratique très souvent le cas puisque l'on cherche justement à estimer cette espérance qui a priori est finie), alors on a la [[convergence presque sûre]] <math>\mu_n \rightarrow \mathbb E[X]</math> lorsque la taille de l'échantillon <math>n</math> tend vers l'infini. 

Il est même possible de connaître l'ordre de grandeur de l'erreur commise par cet estimateur grâce au [[théorème central limite]]. En effet, si <math>X</math> admet une variance finie <math>\sigma^2<\infty</math>, alors le théorème central limite nous dit que la convergence

<math display="block">\frac{\mu_n-\mathbb E[X]}{\sigma/\sqrt n} \rightarrow \mathcal N(0,1) </math>

a lieu en loi quand <math>n </math> tend vers l'infini, où <math>\mathcal N(0,1) </math> désigne une [[loi normale]] centrée réduite. On peut donc en déduire que, si <math>X</math> admet une variance finie, alors l'erreur de l'estimateur de Monte-Carlo est de l'ordre de <math>1/\sqrt n </math> et l'on voit bien que plus la taille de l'échantillon est grande, plus l'erreur se rapproche de zéro. Par exemple, multiplier par 100 la taille de l'échantillon devrait, à peu près, diviser par 10 l'erreur commise.

En utilisant le théorème central limite, on peut déduire des [[Intervalle de confiance|intervalles de confiance]] asymptotiques pour l'estimateur de Monte-Carlo. Soit un réel <math>\alpha\in[0,1] </math> qui correspond au niveau de confiance désiré et notons <math>q_{\beta} </math> le quantile d'ordre <math>\beta </math> de la loi normale centrée réduite. C'est-à-dire que, si <math>Z \sim \mathcal N(0,1) </math> alors <math>\mathbb P(Z \leq q_\beta)=\beta </math> (par exemple pour <math>\beta=97.5% </math> on a <math>q_{0,975} \approx 1,96 </math>). Alors on a la convergence

<math display="block">\mathbb P\left(\mathbb E[X] \in \left[\mu_n-\frac{\sigma q_{1-\alpha/2}}{\sqrt n},\mu_n+\frac{\sigma q_{1-\alpha/2}}{\sqrt n}\right]\right) \rightarrow1-\alpha </math>

lorsque <math>n </math> tend vers l'infini. Il faut garder en tête que cet intervalle de confiance est bien asymptotique et pas exact car il repose sur le théorème central limite qui est un résultat de convergence et non une égalité. En pratique, on connaît rarement la variance <math>\sigma^2 </math> de <math>X</math> ce qui est gênant dans le calcul de l'intervalle de confiance asymptotique ci-dessus. On peut résoudre ce problème en approximant cette variance via l'estimateur

<math display="block">\sigma_n^2 := \frac 1n \sum_{k=1}^n (X_k-\mu_n)^2.</math>

Encore une fois, par la loi forte des grands nombres, cet estimateur converge vers la variance <math>\sigma^2</math> lorsque cette dernière est finie.

Il est à noter que pour appliquer la méthode de Monte-Carlo décrite ci-dessus, il est supposé qu'il est possible de générer des variables aléatoires suivant la loi de <math>X</math>. En pratique, cela peut constituer une difficulté à part entière, de plus il n'est pas toujours possible de générer facilement de telles variables aléatoires avec exactitude. Dans ce dernier cas, on cherchera alors à générer des variables aléatoires suivant une loi approximant au mieux la loi théorique de <math>X</math>.

La méthode de Monte-Carlo présentée ci-dessus est l'une des versions les plus simples et épurées qui soient. Il existe de nombreuses variantes, plus précises ou plus rapides ou adaptées à d'autres circonstances, dont certaines sont mentionnées dans la section dédiées plus bas.

== Applications ==

=== Estimer l'aire d'une surface ou un volume ===
Une grande catégorie de problèmes pouvant se résoudre par l'approche de Monte-Carlo est celle de l'estimation de l'aire d'une surface ou plus généralement d'un volume en dimension quelconque. Le principe général consiste à générer des points uniformément au hasard dans une zone bornée mais assez large pour contenir la surface ou le volume à estimer. Ensuite, il suffit de compter le nombre de points tombant dans la surface ou le volume et de diviser par le nombre total de points générés pour estimer l'aire ou le volume. Quelques exemples sont présentés ci-dessous.

==== Détermination de la superficie d'un lac ====
Cet exemple est un classique en vulgarisation de la méthode de Monte-Carlo. Soit une zone rectangulaire ou carrée dont les côtés sont de longueur connue. Au sein de cette aire se trouve un lac dont la superficie est inconnue. Grâce aux mesures des côtés de la zone, on connaît l'aire du rectangle. Pour trouver l'aire du lac, on demande à une armée de tirer {{Formule|''X''}} coups de canon de manière aléatoire sur cette zone. On compte ensuite le nombre {{Formule|''N''}} de boulets qui sont restés sur le terrain ; on peut ainsi déterminer le nombre de boulets qui sont tombés dans le lac : {{Formule|''X''−''N''}}. Il suffit ensuite d'établir un rapport entre les valeurs :

:<math>\frac{\mathrm{superficie}_{~\mathrm{terrain}}}{\mathrm{superficie}_{~\mathrm{lac}}} = \frac{X}{X-N}</math>

:<math>\Longrightarrow \qquad \mathrm{superficie}_{~\mathrm{lac}} = \frac{(X-N)}{X} \ \times \ \mathrm{superficie}_{~\mathrm{terrain}}</math>

Par exemple, si le terrain fait {{unité|1000|m|2}}, que l'armée tire 500 boulets et que 100 projectiles sont tombés dans le lac, alors une estimation de la superficie du plan d'eau est de : 1000×100÷500 = {{unité|200|m|2}}.

[[Image:Montecarlo.svg|center|vignette|redresse=1.2|Estimation de la surface du lac grâce à des tirs d'artillerie aléatoires.]]

La qualité de l'estimation s'améliore (lentement) en augmentant le nombre de tirs et en s'assurant que les artilleurs ne visent pas toujours le même endroit mais couvrent bien la zone, de manière uniforme. Cette dernière remarque est à mettre en parallèle avec la qualité du [[Générateur de nombres aléatoires|générateur aléatoire]] qui est primordiale pour avoir de bons résultats dans la méthode de Monte-Carlo. Un générateur biaisé est comme un canon qui tire toujours au même endroit : les informations qu'il apporte sont réduites.

==== Détermination de la valeur de {{math|π}} ====
Cette méthode est proche de l'expérience de l'[[aiguille de Buffon]].

Soit un point ''M'' de coordonnées {{math|(''x'', ''y'')}}, où {{math|0 < ''x'' < 1}} et {{math|0 < ''y'' < 1}}. On tire aléatoirement les valeurs de {{math|''x''}} et {{math|''y''}} entre 0 et 1 suivant une loi uniforme. Le point ''M'' appartient au disque de centre (0,0) de rayon {{math|''R'' {{=}} 1}} si et seulement si {{math|''x''{{2}} + ''y''{{2}} ≤ 1}}. La probabilité que le point ''M'' appartienne au disque est {{sfrac|π|4}}, puisque le quart de disque est de surface {{math|σ{{=}}{{sfrac|π ''R''{{2}}|4}} {{=}} {{sfrac|π|4}}}}, et le carré qui le contient est de surface {{math|''S'' {{=}} ''R''{{2}}{{=}}1}} : si la [[loi de probabilité]] du tirage de point est [[Loi uniforme continue|uniforme]], la probabilité de tomber dans le quart de disque vaut {{math|{{sfrac|σ|''S''}} {{=}} {{sfrac|π|4}}}}.

En faisant le rapport du nombre de points dans le disque au nombre de tirages, on obtient une approximation du nombre {{sfrac|π|4}} si le nombre de tirages est grand.

[[Image:Montecarlo-valeur-pi.svg|center|vignette|redresse=1.5|Représentation du calcul de la valeur de {{math|π}} par rapport du nombre de points aléatoires étant contenus dans un quart de cercle, l'ensemble des possibles étant un carré de côté ''R''.]]

==== Calcul d'une intégrale ====
L'intégrale d'une fonction correspond à l'aire algébrique sous sa courbe. Estimer l'intégrale d'une fonction revient donc à estimer l'aire d'une surface et on peut donc appliquer une méthode de Monte-Carlo. 

Par exemple supposons que l'on veuille estimer l'intégrale suivante <math display="block">\int_0^1 f(u)du</math>

où <math>f:[0,1]\rightarrow\R</math> est une [[Continuité (mathématiques)|fonction continue]]. Cette intégrale peut être vue comme une espérance. En effet, soit <math>U</math> une variable aléatoire de loi uniforme sur <math>[0,1]</math>, alors par la [[Théorème de transfert|formule du transfert]] on a<math display="block">\mathbb E[f(U)] =\int_0^1 f(u)du.</math>

On peut alors appliquer une méthode simple de Monte-Carlo et générer <math>n</math> variables aléatoires iid <math>U_1,\dots,U_n</math> uniformément au hasard dans <math>[0,1]</math> et approximer l'intégrale ci-dessus par <math display="block">\frac 1n \sum_{k=1}^n f(U_k).</math>Ce procédé peut s'appliquer à une intégrale en dimension quelconque et pas nécessairement en dimension une comme ici. D'ailleurs en grande dimension, les méthodes de Monte-Carlo sont bien plus efficaces que d'autres méthodes classiques comme celle des sommes de Riemann par exemple.

==== Recouvrement de courbes et méthode contrainte-résistance ====
[[Fichier:Contrainte resistance 2d eloigne.svg|vignette|redresse=1.2|Méthode contrainte-résistance : la fréquence de défaillance correspond à l'aire sous l'intersection des courbes rouge et bleue (haut).]]

La méthode de Monte-Carlo peut être utilisée pour déterminer l'aire sous l'intersection de deux courbes, qui n'est qu'une surface particulière.

Les courbes peuvent être les courbes représentatrices des [[densité de probabilité|densités de probabilité]] de deux lois. C'est par exemple utilisé dans la [[méthode contrainte-résistance]] :
* un système est soumis à une contrainte — une sollicitation quelle qu'elle soit (effort mécanique, variation de température, passage de [[courant électrique]], …) — dont l'intensité est une variable aléatoire S ;
* le système est conçu pour résister à cette contrainte, sa résistance est exprimée par une valeur, une variable aléatoire R ;
* le système est validé si la résistance est supérieure à la contrainte — R &gt; S — dans un certain nombre de cas (typiquement 99 % ou 99,9 %).
La probabilité complémentaire — les cas de défaillance, celle pour laquelle R ≤ S — est l'aire sous l'intersection des deux courbes représentant les lois.

On peut déterminer la probabilité P(R &gt; S) en faisant des tirages aléatoires sur R et S et en dénombrant les cas pour lesquels « R &gt; S » est vrai.

=== Application au modèle d'Ising ===

{{Article détaillé|Modèle d'Ising}}
{{...}}

=== Estimation de la valeur d'un coup au go ===
{{article détaillé|Jeu de go (informatique)#Méthode de Monte-Carlo}}
Aux échecs, comme dans beaucoup de jeux de plateau,  il est possible de mesurer la valeur d'une position, et donc des coups y conduisant, en évaluant quantitativement la position obtenue : nombre de pièces sur l'échiquier, valeurs des pièces (1 point par pion, 5 par tour...), position relative des pièces entre elles, et en pondérant la valeur trouvée par les libertés, les protections des pièces, etc. Cette évaluation basée sur l'analyse et l'expertise est d'autant plus rapide à mesurer qu'on avance dans la partie, car le nombre de pièces diminue.

Dans le jeu de go, l'évaluation d'une position globale reste très difficile avec des méthodes d'analyses classiques du fait de l’enchevêtrement et de la complexité des positions locales et du nombre virtuellement infini de suites de coups possibles. En 2006, le mathématicien Rémi Coulom a fait progresser de manière très sensible cette fonction d'évaluation et l'efficience des logiciels de jeu de go en utilisant la méthode de Monte-Carlo : on joue ''"au hasard"'' un grand nombre de fins de parties réalistes à partir de la position "en cours d'évaluation" et on comptabilise la proportion de parties gagnantes/perdantes. Cette estimation statistique s'affine en biaisant le hasard par élimination de coups ''a priori'' stupides. Cette méthode s'avère très efficace<ref>{{Lien web|langue = |titre = le jeu de go, le seul jeu où l'ordinateur ne bat pas l'homme|url = http://www.slate.fr/life/86963/go-seul-jeu-ordinateur-pas-battre-homme|site =[[Slate (magazine)|Slate]] |date = 5 mai 2014|consulté le = 15 septembre 2015}}.</ref>{{,}}<ref>{{Lien web|langue = |titre = le jeu de go et la révolution de  monte-carlo|url = https://interstices.info/jcms/c_43860/le-jeu-de-go-et-la-revolution-de-monte-carlo|site = interstices.info|date = 28 mai 2009|consulté le = 15 septembre 2015}}.</ref>. Elle est utilisée en particulier par les programmes [[AlphaGo]] et [[AlphaZero]]<ref>{{lien web|langue=en|auteur1=David Silver|auteur2=[[Demis Hassabis]]|url=http://googleresearch.blogspot.com/2016/01/alphago-mastering-ancient-game-of-go.html|titre=AlphaGo: Mastering the ancient game of Go with Machine Learning |date=27 janvier 2016 |site=Google Research Blog}}.</ref>.

=== Estimation de la valeur d'un coup aux échecs ===
Parfois, pour savoir si un coup ambigu doit être fait (échange de pièces par exemple), lorsqu'on manque d'information {{Lesquels}}, ou pour choisir entre plusieurs coups menant tous à des pertes de matériel, il est possible de lancer plusieurs parties rapides au hasard{{C'est-à-dire}}, pour savoir quelle est la moins mauvaise ou la meilleure des solutions{{référence nécessaire}}.

=== Probabilité de la performance en bourse ===
Selon l'[[Hypothèse des marchés financiers efficients|hypothèse des marchés efficients]], les performances boursières sont aléatoires et suivent une loi normale. Dans ces conditions, il est possible de réaliser des milliers de tirages aléatoires pour déterminer les probabilités d'atteindre certaines performances boursières dans le futur<ref>{{Lien web|langue=fr-FR|nom1=Petit|prénom1=Edouard|titre=Prévoir la bourse et votre patrimoine grâce aux simulations de Monte-Carlo |url=https://www.epargnant30.fr/prevoir-bourse-simulation-monte-carlo/|site=Epargnant 3.0|date=2018-04-08|consulté le=2020-01-16}}.</ref>. Cette caractéristique des marchés financiers permet d'utiliser la méthode de Monte-Carlo pour valoriser des [[Call|options]], ou encore des [[Special Purpose Acquisition Company|SPAC]]. 

=== Espérance du gain d'une option dans le modèle de Black et Scholes ===
{{article détaillé|Modèle Black-Scholes}}
On note <math>S_t</math> le prix d'un actif fixé au temps <math>t</math>. Le modèle de Black et Scholes consiste à dire que le prix de cet actif répond à l'équation différentielle stochastique suivante : <math>dS_t = S_t(\mu dt + \sigma dW_t)</math> où <math>\mu</math> est un réel (appelé parfois dérive), <math>\sigma</math> est un réel positif (appelé volatilité) et <math>W</math> désigne un [[mouvement brownien]] standard. 

On suppose que l'on veuille estimer l'espérance du gain perçu par la détention d'une option européenne d'achat (call) de maturité <math>T</math> et de prix d'exercice (strike) <math>K</math>. Si le prix au temps <math>T</math>, à savoir <math>S_T</math>, est plus grand que <math>K</math>, l'option est exercée et le gain est de <math>S_T-K</math>, sinon le gain est nul. On cherche donc à estimer la quantité <math>\mathbb E[(S_T-K)_+]</math> où <math>x_+ := \max(0,x)</math>.

Pour estimer cette espérance, on peut avoir recours à une méthode de Monte-Carlo. Pour cela, il faut alors pouvoir générer des variables aléatoires suivant la loi de <math>S_T</math>. On peut utiliser directement la forme de la solution de l'équation différentielle stochastique, à savoir 

<math display="block">S_t = S_0\exp\left(\mu t - \frac{\sigma^2}{2}t +\sigma W_t\right)\qquad\forall t\geq0.</math>

Il suffit alors de générer des variables aléatoires iid de même loi que <math>W_T</math> qui est simplement une gaussienne centrée de variance <math>T</math>. On approxime ensuite l'espérance du gain de l'option par<math display="block">\frac 1n \sum_{k=1}^n (S_T^k-K)_+</math>

où <math>S_T^1,\dots,S_T^n</math> sont <math>n</math> variables iid de même loi que <math>S_T = S_0\exp(\mu t -\sigma^2/2 t +\sigma W_T)</math>.

Maintenant on suppose que l'on veuille étudier cette fois-ci une option américaine et non plus européenne. Dans ce cas il nous faudra générer tout le processus <math>(S_t)_{t\leq T}</math> et non plus simplement la variable <math>S_T</math>. Une solution consiste à utiliser la [[méthode d'Euler]] à partir de l'équation différentielle stochastique. On commence par découper l'intervalle de temps <math>[0,T]</math> en <math>M</math> parties égales et on pose <math>\Delta t := T/M</math>. Ensuite on génère une suite de variables aléatoires <math display="block">\left\{\begin{array}{lll}
S_0^1 & = & S_0 \\
S_{(i+1)\Delta t}^1 & = &S_{i\Delta t}^1(1+\mu \Delta t +\sigma Z_i^1 \sqrt{\Delta t}) \qquad\forall i < M
\end{array}\right.</math>

où <math>Z_0^ 1,\dots,Z_{M-1}^1</math> sont des gaussiennes centrées réduites indépendantes. On complète ensuite par [[interpolation linéaire]]. Pour <math>M</math> choisi très grand, le processus aléatoire ainsi construit <math>(S_t^1)_{t\leq T}</math> devrait avoir une loi assez proche de celle de <math>(S_t)_{t\leq T}</math>. On peut répéter la construction précédente pour obtenir <math>n</math> processus aléatoires <math>(S_t^1)_{t\leq T} ,\dots,(S_t^n)_{t\leq T}</math> iid de loi proche de celle de <math>(S_t)_{t\leq T}</math>.

== Méthodes accélérées ==

Le nombre de simulations <math>n</math> nécessaires pour atteindre une marge d'erreur voulue est parfois trop important avec la méthode Monte-Carlo. C'est le cas notamment pour estimer la probabilité d'un évènement rare. Dans ce cas la méthode Monte-Carlo est déconseillée, car le nombre de simulations nécessaires pour obtenir une bonne précision est énorme, ce qui donne des temps calcul beaucoup trop long...

Diverses méthodes, dites techniques de [[réduction de la variance]], permettent d'améliorer la précision de l'estimation — ou de diminuer le nombre de simulations nécessaires — en remplaçant <math>X</math> par une autre [[variable aléatoire]]. Ces techniques entrent en général dans l'une des classes suivantes : l'[[échantillonnage préférentiel]], les [[variable de contrôle|variables de contrôle]], le [[conditionnement (Monte-Carlo)|conditionnement]], la [[variable antithétique]], la [[stratification (Monte-Carlo)|stratification]]. L'[[échantillonnage préférentiel]] est théoriquement la méthode qui peut réduire le plus la variance, mais en pratique, cette méthode est parfois difficile a mettre en œuvre. Le cas échéant, les autres méthodes peuvent donner de meilleur résultat. Dans le cas où la variable aléatoire <math>X</math> est une trajectoire aléatoire  (suite de variable aléatoire, ou trajectoire d'un processus), on peut aussi utiliser des [[Filtre particulaire|méthodes particulaires]], ou si l'on dispose d'une densité pour la trajectoire <math>X</math> , l'[[échantillonnage préférentiel]] est aussi une possibilité. 

== Notes et références ==
{{Références}}

== Voir aussi ==
=== Bibliographie ===
* Emmanuel Gobet, ''[http://www.editions.polytechnique.fr/?afficherfiche=188 Méthodes de Monte-Carlo et processus stochastiques - du linéaire au non linéaire]'', Éditions de l'École polytechnique, 2013
* Carl Graham, Denis Talay, ''[http://www.editions.polytechnique.fr/?afficherfiche=165 Simulation stochastique et méthodes de Monte-Carlo]'', Éditions de l'École polytechnique, 2011
* {{en}} [http://cel.archives-ouvertes.fr/cel-00092936/ Introduction to Monte Carlo Algorithm] Werner Krauth - CNRS, Laboratoire de Physique Statistique, École Normale Supérieure
* {{en}} Michael Mascagni, ''Advanced Monte Carlo Methods I & II'', Cours du ETH de Zurich (2005/2006). {{pdf}} [http://www.cs.fsu.edu/~mascagni/Lecture_Notes.html Notes de cours]. On pourra également consulter la [http://www.cs.fsu.edu/~mascagni/Advanced_Monte_Carlo_Methods.html page de présentation du cours], qui contient de nombreuses références disponibles au format en pdf.
* Simon Léger, ''Monte Carlo pour les nuls'', 2006 {{pdf}} {{lire en ligne|lien=http://homepages.nyu.edu/~sl1544/MonteCarloNuls.pdf}}
* [http://www.ina.fr/archivespourtous/index.php?vue=notice&from=fulltext&full=physicien&num_notice=7&total_notices=21 Une explication de la méthode de Monte-Carlo par le physicien Pierre Auger]
*{{Ouvrage|prénom1=J.|nom1=Morio|prénom2=M.|nom2=Balesdent|titre=Estimation of Rare Event  Probabilities in Complex Aerospace and Other Systems|sous-titre=A Practical Approach|lieu=Cambridge|éditeur=Elsevier Science|année=2015|pages totales=216|isbn=978-0-08-100091-5}}
* {{en}} [http://www.ipgp.fr/~tarantola/Files/Professional/Teaching/Paris/index.html Cours de problèmes inverses par méthode de Monte Carlo] - A. Tarantola, Institut de Physique du Globe de Paris
* {{Ouvrage|langue=en|prénom1=Christian|nom1=Robert|lien auteur1=Christian Robert (statisticien)|prénom2=George|nom2=Casella|lien auteur2=George Casella|titre=Monte Carlo Statistical Methods|éditeur=[[Springer Science+Business Media|Springer-Verlag]]|collection=Springer Texts in Statistics|année=2010|isbn=}}
* {{Ouvrage|langue=en|auteur1=Christian Robert (statisticien)|auteur2=George Casella|titre=Introducing Monte Carlo Methods with R|éditeur=[[Springer Science+Business Media|Springer-Verlag]]|collection=Use R! Series|année=2007|pages totales=283|isbn=978-1-4419-1575-7|lire en ligne=https://books.google.com/books?id=WIjMyiEiHCsC&printsec=frontcover}}

=== Articles connexes ===
* [[Simulation informatique]],
* Simuler des processus se produisant à des taux connus : [[Méthode de Monte-Carlo cinétique]].
* les techniques les plus efficaces de [[planification de mouvement]] utilisent des algorithmes probabilistes.
* Méthodes de réduction de la variance de l'estimateur de Monte Carlo : [[Importance sampling]]
* [[Physique statistique]]
* [[Filtre particulaire (statistique)|Filtre particulaire]]
* [[Algorithme probabiliste]]
* [[Algorithme de Las Vegas]]
* [[Algorithme de Monte-Carlo]]
* [[Méthode de Monte-Carlo par chaînes de Markov]]
* [[Prévision d'ensembles]] en météorologie
* [[Recherche arborescente Monte-Carlo]]

'''Codes de simulation utilisant des méthodes de Monte-Carlo'''
* [[Geant4]]
* [[Monte-Carlo N-Particle transport|MCNP]] et MCNPX
* [[Tripoli-4]]

=== Liens externes ===
* {{en}} [http://www.vertex42.com/ExcelArticles/mc/MonteCarloSimulation.html Un exemple de simulation Monte-Carlo clair et didactique] L'exemple se base sur la fabrication d'une pièce mécanique à partir de ses composants.
* {{en}} [http://fr.mathworks.com/matlabcentral/fileexchange/56509-monte-carlo-simulator MATLAB Monte-Carlo simulator with several examples]

{{Portail|probabilités|physique|finance}}

{{DEFAULTSORT:Methode de Monte-Carlo}}

[[Catégorie:Analyse numérique]]
[[Catégorie:Algorithme numérique]]
[[Catégorie:Algorithme probabiliste|Monte Carlo]]
[[Catégorie:Physique statistique]]
[[Catégorie:Science des données]]
[[Catégorie:Évaluation d'actif]]