Transformée de Box-Muller

Transformée de Box-Muller

Méthode de Box-Muller

Représentation graphique de la transformation : les cercles de départ, répartis uniformément autour de l'origine, deviennent un nouvel ensemble de cercles centrés, dont la répartition est proche de l'origine puis s'étiole rapidement. Les plus grands cercles de départ correspondent aux plus petits cercles d'arrivée, et vice-versa.

La méthode de Box-Muller (George Edward Pelham Box et Mervin Edgar Muller, 1958) consiste à générer des paires de nombres aléatoires à distribution normale centrée réduite, à partir d'une source de nombres aléatoires de loi uniforme.

La transformation prend communément deux formes.

  • La forme simple transforme des coordonnées cartésiennes uniformément distribuées dans le cercle unité en des coordonnées normalement distribuées.
  • La forme polaire transforme des coordonnées polaires uniformément distribuées en des coordonnées cartésiennes normalement distribuées.

On peut également utiliser la méthode de la transformée inverse pour générer des nombres normalement distribués; la méthode de Box-Muller a été mise au point pour être algorithmiquement plus efficiente[1]. On peut également envisager la Méthode Ziggourat qui est aussi très efficiente.

Sommaire

Forme cartésienne

Soient x et y choisis indépendamment et uniformément dans [−1,1], et s = x2 + y2. Si s > 1, rejetons-le et choisissons à nouveau un couple (x, y), jusqu'à ce que s appartienne à ]0,1]. Pour ces points "filtrés", calculons ensuite:

z_0 = x \cdot \sqrt{\frac{-2 \ln s}{s}}

et

z_1 = y \cdot \sqrt{\frac{-2 \ln s}{s}}.

Forme polaire

Soient U1 et U2 deux variables aléatoires indépendantes uniformément distribuées dans ]0,1].

Soient

Z_0 = R \cos(\Theta) =\sqrt{-2 \ln U_1} \cos(2 \pi U_2)\,

et

Z_1 = R \sin(\Theta) = \sqrt{-2 \ln U_1} \sin(2 \pi U_2).\,

Alors Z0 et Z1 sont des variables aléatoires indépendantes suivant une loi normale de variance 1.

Cette transformation [2] vient du fait que, dans un système cartésien à deux dimensions où les coordonnées X et Y suivent deux variables aléatoires indépendantes normales, les variables aléatoires R2 et Θ (ci-dessus) sont également indépendantes et peuvent s'écrire

R^2 = -2\cdot\ln U_1\,

et

\Theta = 2\pi U_2.\,

Comparaison entre les deux formes

La forme cartésienne est une méthode d'échantillonnage à rejet, qui n'utilise qu'une partie des nombres générés par la source aléatoire, mais elle est en pratique plus rapide que la forme polaire car elle est plus simple à calculer:

  • la forme cartésienne n'utilise pas de fonctions trigonométriques, coûteuses en temps de calcul
  • la génération de nombres aléatoires est plutôt rapide, il n'est donc pas gênant d'en gaspiller une partie. En moyenne la part de points rejetés est (1-π/4) ≈ 21.46%. On génère donc 4/π ≈ 1.2732 nombres aléatoires uniformes pour obtenir chaque nombre aléatoire normal.

Notes

  1. ↑ Kloeden and Platen, Numerical Solutions of Stochastic Differential Equations, p. 11-12
  2. ↑ Sheldon Ross, A First Course in Probability, (2002), p.279-81
Ce document provient de « M%C3%A9thode de Box-Muller ».

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Transformée de Box-Muller de Wikipédia en français (auteurs)

Игры ⚽ Нужно решить контрольную?

Regardez d'autres dictionnaires:

  • Box-Muller — Méthode de Box Muller Représentation graphique de la transformation : les cercles de départ, répartis uniformément autour de l origine, deviennent un nouvel ensemble de cercles centrés, dont la répartition est proche de l origine puis s… …   Wikipédia en Français

  • Méthode de Box-Muller — Représentation graphique de la transformation : les cercles de départ, répartis uniformément autour de l origine, deviennent un nouvel ensemble de cercles centrés, dont la répartition est proche de l origine puis s étiole rapidement. Les… …   Wikipédia en Français

  • Methode de Box-Muller — Méthode de Box Muller Représentation graphique de la transformation : les cercles de départ, répartis uniformément autour de l origine, deviennent un nouvel ensemble de cercles centrés, dont la répartition est proche de l origine puis s… …   Wikipédia en Français

  • Méthode De Box-Muller — Représentation graphique de la transformation : les cercles de départ, répartis uniformément autour de l origine, deviennent un nouvel ensemble de cercles centrés, dont la répartition est proche de l origine puis s étiole rapidement. Les… …   Wikipédia en Français

  • Méthode de box-muller — Représentation graphique de la transformation : les cercles de départ, répartis uniformément autour de l origine, deviennent un nouvel ensemble de cercles centrés, dont la répartition est proche de l origine puis s étiole rapidement. Les… …   Wikipédia en Français

  • Transformation de Box-Muller — Méthode de Box Muller Représentation graphique de la transformation : les cercles de départ, répartis uniformément autour de l origine, deviennent un nouvel ensemble de cercles centrés, dont la répartition est proche de l origine puis s… …   Wikipédia en Français

  • Projet:Mathématiques/Liste des articles de mathématiques — Cette page n est plus mise à jour depuis l arrêt de DumZiBoT. Pour demander sa remise en service, faire une requête sur WP:RBOT Cette page recense les articles relatifs aux mathématiques, qui sont liés aux portails de mathématiques, géométrie ou… …   Wikipédia en Français

  • Loi Uniforme Continue — Uniforme Densité de probabilité / Fonction de masse Fonction de répartition …   Wikipédia en Français

  • Loi uniforme continue — Uniforme Densité de probabilité / Fonction de masse Fonction de répartition …   Wikipédia en Français

  • Liste des articles de mathematiques — Projet:Mathématiques/Liste des articles de mathématiques Cette page recense les articles relatifs aux mathématiques, qui sont liés aux portails de mathématiques, géométrie ou probabilités et statistiques via l un des trois bandeaux suivants  …   Wikipédia en Français

Share the article and excerpts

Direct link
… Do a right-click on the link above
and select “Copy Link”