Théorème central limite
Le théorème central limite (aussi appelé théorème limite central, théorème de la limite centrale ou théorème de la limite centrée) établit la convergence en loi de la somme d'une suite de variables aléatoires vers la loi normale. Intuitivement, ce résultat affirme qu'une somme de variables aléatoires indépendantes et identiquement distribuées tend (le plus souvent) vers une variable aléatoire gaussienne.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/3/37/Standard_deviation_diagram_%28decimal_comma%29.svg/440px-Standard_deviation_diagram_%28decimal_comma%29.svg.png)
Ce théorème et ses généralisations offrent une explication de l'omniprésence de la loi normale dans la nature : de nombreux phénomènes sont dus à l'addition d'un grand nombre de petites perturbations aléatoires.
Histoire
modifierHistoire de la démonstration
modifierLa première démonstration de ce théorème, publiée en 1809, est due à Pierre-Simon de Laplace[1],[2], mais le cas particulier où les variables suivent la loi de Bernoulli de paramètre p = 0,5 était connu depuis les travaux de De Moivre[3], en 1733.
À propos de la dénomination de ce théorème
modifierLa dénomination « théorème central limite » fait référence à un document scientifique écrit par George Pólya en 1920, intitulé Über den zentralen Grenzwertsatz der Wahrscheinlichkeitsrechnung und das Momentenproblem[4] (Sur le théorème [ayant rapport à la notion de] limite central du calcul probabiliste et le problème des moments). Historiquement, et conformément à la traduction du titre, c'est donc bien le théorème qui est central, d'où l'appellation « théorème central limite ».
Cependant, dans la littérature mathématique française, on peut trouver d'autres dénominations, comme « théorème limite central »[5], « théorème de la limite centrale » ou « théorème de la limite centrée ». Une justification avancée par certains auteurs est que l'adjectif « central » s'applique au centre de la distribution, par opposition à sa queue[4].
Illustrations
modifierCe théorème est évident si n variables aléatoires suivent une loi normale d'espérance μ : en effet leur somme suit une loi normale de paramètre nμ. Dans le cas de variables ne suivant pas une loi normale, le théorème peut sembler étonnant au premier abord.
Pile ou face
modifier![](http://upload.wikimedia.org/wikipedia/commons/thumb/0/0f/Pile_ou_face.png/110px-Pile_ou_face.png)
Donnons d'abord une illustration avec le jeu de pile ou face, illustration ne nécessitant pas de connaissance particulière en statistiques, mais uniquement en dénombrement.
Considérons le jeu de pile ou face et mettons des valeurs sur les faces de la pièce, par exemple 0 pour pile et 1 pour face ; on s'intéresse à la somme de n tirages. La pièce est équilibrée, chaque face a une chance sur deux d'être tirée. Si l'on fait un seul tirage, nous avons donc le tirage n°1 (et aucun autre), et son résultat peut être 0 ou 1 ; nous faisons la somme d'une seule valeur.
Résultat tirage n°1 |
Somme |
---|---|
0 | 0 |
1 | 1 |
Nous avons donc n = 2 possibilités pour la valeur de la somme, apparaissant avec les fréquences suivantes :
Valeurs de la somme |
Nombre d'apparitions |
Fréquence |
---|---|---|
0 | 1 | 1/2 = 0,5 (50 %) |
1 | 1 | 1/2 = 0,5 (50 %) |
Avec deux tirages, chaque tirage peut donner 0 ou 1, ce qui donne le tableau suivant :
Résultat tirage n°1 |
Résultat tirage n°2 |
Somme |
---|---|---|
0 | 0 | 0 |
0 | 1 | 1 |
1 | 0 | 1 |
1 | 1 | 2 |
nous avons n = 4 possibilités, soit le tableau des fréquences.
Valeurs de la somme |
Nombre d'apparitions |
Fréquence |
---|---|---|
0 | 1 | 1/4 = 0,25 (25 %) |
1 | 2 | 2/4 = 0,5 (50 %) |
2 | 1 | 1/4 = 0,25 (25 %) |
Et ainsi de suite :
Résultat tirage n°1 |
Résultat tirage n°2 |
Résultat tirage n°3 |
Somme | Valeurs de la somme |
Nombre d'apparitions |
Fréquence | |
---|---|---|---|---|---|---|---|
0 | 0 | 0 | 0 | 0 | 1 | 0,125 (12,5 %) | |
0 | 0 | 1 | 1 | 1 | 3 | 0,375 (37,5 %) | |
0 | 1 | 0 | 1 | 2 | 3 | 0,375 (37,5 %) | |
0 | 1 | 1 | 2 | 3 | 1 | 0,125 (12,5 %) | |
1 | 0 | 0 | 1 | ||||
1 | 0 | 1 | 2 | ||||
1 | 1 | 0 | 2 | ||||
1 | 1 | 1 | 3 |
Graphiquement, on constate que plus le nombre de tirages augmente, plus la courbe de fréquence se rapproche d'une courbe en cloche symétrique, caractéristique de la densité de probabilité de la loi normale.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/2/29/Somme_n_tirages_pile_ou_face_1_a_4.svg/220px-Somme_n_tirages_pile_ou_face_1_a_4.svg.png)
![](http://upload.wikimedia.org/wikipedia/commons/thumb/5/50/Somme_n_tirages_pile_ou_face_1_a_12.svg/220px-Somme_n_tirages_pile_ou_face_1_a_12.svg.png)
Dé à six faces
modifierOn obtient un résultat similaire en jetant plusieurs dés à six faces (d6) et en faisant la somme des nombres apparus sur leurs faces supérieures, mais le dénombrement est plus fastidieux (il y a six valeurs par dé)[6].
![](http://upload.wikimedia.org/wikipedia/commons/thumb/c/c3/Somme_n_tirages_d6_1_a_5.svg/220px-Somme_n_tirages_d6_1_a_5.svg.png)
Cas avec des lois différentes
modifierDans toutes les situations ci-dessus, on a des lois uniformes ; et pourtant, la somme d'un grand nombre d'événements tend graphiquement vers une courbe en cloche symétrique. Et cela reste vrai même lorsque les lois sont différentes (cas des dés polyédriques non réguliers).
![](http://upload.wikimedia.org/wikipedia/commons/thumb/1/1e/Somme_tirages_des_polyedriques.svg/220px-Somme_tirages_des_polyedriques.svg.png)
En effet, on ne s'intéresse pas au tirage en lui-même, mais à la somme du tirage. De ce point de vue, plusieurs tirages sont équivalents, donc une valeur de somme peut être obtenue par plusieurs tirages ; par exemple, pour deux dés à six faces (2d6), on peut obtenir 7 par 1+6, 2+5, 3+4, 4+3, 5+2 et 6+1, il y a six tirages équivalents. Or, il y a toujours plus de combinaisons permettant d'obtenir une valeur moyenne qu'une valeur extrême, ce qui donne la courbe en cloche.
![alternative textuelle](http://upload.wikimedia.org/wikipedia/commons/thumb/7/79/IllustrationTheoremeCentral.png/550px-IllustrationTheoremeCentral.png)
Le théorème central limite
modifierÉnoncé
modifierSoit une suite de variables aléatoires réelles définies sur le même espace de probabilité, indépendantes et identiquement distribuées suivant la même loi . Supposons que l'espérance et l'écart-type de existent et soient finis avec .
Considérons la somme
Alors
- l'espérance de est et
- son écart-type vaut
De plus, quand est assez grand, la loi normale est une bonne approximation de la loi de .
Afin de formuler mathématiquement cette approximation, nous allons poser
et
- ,
de sorte que l'espérance et l'écart-type de Zn valent respectivement 0 et 1 : la variable est ainsi dite centrée et réduite.
Le théorème central limite énonce alors que la suite de variables aléatoires converge en loi vers une variable aléatoire , définie sur le même espace probabilisé, et de loi normale centrée réduite lorsque n tend vers l'infini.
Cela signifie que si est la fonction de répartition de , alors pour tout réel :
Démonstration du théorème central limite
modifierPour un théorème d'une telle importance en statistiques et en probabilité appliquée, il existe une démonstration particulièrement simple utilisant les fonctions caractéristiques. Cette démonstration ressemble à celle d'une des lois des grands nombres.
Développement limité des variable centré réduite
modifierPour une variable aléatoire d'espérance 0 et de variance 1. Comme sa fonction caractéristique est de classe et par propriété de la fonction caractéristique [8]
donc on a et
Ainsi, admet un développement limité au voisinage de 0 de la forme
Calcul de la fonction caractéristique de la variable Zn
modifierOr si on prend les comme les variable centrée réduite des c'est-à-dire : , on a simplement :
- .
- Soit ( ), et pour tout ( ), comme () sont mutuellement indépendantes et identiquement distribuées, on a :[9]
Convergence et identification
modifierDonc la fonction caractéristique de est lorsque .
Mais cette limite est la fonction caractéristique de la loi normale centrée réduite , d'où l'on déduit le théorème central limite grâce au théorème de convergence de Lévy, qui affirme que la convergence simple des fonctions caractéristiques implique la convergence en loi.
Convergence vers la limite
modifier![](http://upload.wikimedia.org/wikipedia/commons/thumb/f/ff/Approximation_d%27une_distribution_normale.gif/220px-Approximation_d%27une_distribution_normale.gif)
La convergence de la fonction de répartition de Zn est uniforme, en vertu du deuxième théorème de Dini. Si le moment centré d'ordre 3, existe et est fini, alors la vitesse de convergence est au moins d'ordre (voir le théorème de Berry-Esseen).
Images d'une loi lissées par sommation qui montrent la distribution de la loi originale et trois sommations successives (obtenues par convolution) :
![](http://upload.wikimedia.org/wikipedia/commons/thumb/a/a6/Central_limit_thm_1.png/220px-Central_limit_thm_1.png)
![](http://upload.wikimedia.org/wikipedia/commons/thumb/f/f8/Central_limit_thm_2.png/220px-Central_limit_thm_2.png)
![](http://upload.wikimedia.org/wikipedia/commons/thumb/6/69/Central_limit_thm_3.png/220px-Central_limit_thm_3.png)
![](http://upload.wikimedia.org/wikipedia/commons/thumb/8/85/Central_limit_thm_4.png/220px-Central_limit_thm_4.png)
Dans les applications pratiques, ce théorème permet en particulier de remplacer une somme de variables aléatoires en nombre assez grand mais fini par une approximation normale, généralement plus facile à manipuler. Il est donc intéressant de voir comment la somme s'approche de la limite. Les termes utilisés sont expliqués dans l'article Variable aléatoire.
Une somme de variables continues est une variable continue dont on peut comparer la densité de probabilité à celle de la limite normale.
![](http://upload.wikimedia.org/wikipedia/commons/thumb/e/ec/Theoreme_de_la_limite_centrale.png/220px-Theoreme_de_la_limite_centrale.png)
Avec une somme de variables discrètes, il est parfois commode de définir une pseudo-densité de probabilité mais l'outil le plus efficace est la fonction de probabilité représentée par un diagramme en bâtons. On peut constater graphiquement une certaine cohérence entre les deux diagrammes, difficile à interpréter. Dans ce cas, il est plus efficace de comparer les fonctions de répartition.
D'autre part, l'approximation normale est particulièrement efficace au voisinage des valeurs centrales. Certains disent même qu'en matière de convergence vers la loi normale, l'infini commence souvent à six[réf. nécessaire].
La précision se dégrade à mesure qu'on s'éloigne de ces valeurs centrales. C'est particulièrement vrai pour une somme de variables positives par nature : la loi normale fait toujours apparaître des valeurs négatives avec des probabilités faibles mais non nulles. Même si c'est moins choquant, cela reste vrai en toutes circonstances : alors que toute grandeur physique est nécessairement bornée, la loi normale qui couvre un intervalle infini n'est qu'une approximation utile.
Enfin, pour un nombre donné de termes de la somme, l'approximation normale est d'autant meilleure que la distribution est plus symétrique.
Application à la statistique mathématique
modifierCe théorème de probabilités possède une interprétation en statistique mathématique. Cette dernière associe une loi de probabilité à une population. Chaque élément extrait de la population est donc considéré comme une variable aléatoire et, en réunissant un nombre n de ces variables supposées indépendantes, on obtient un échantillon. La somme de ces variables aléatoires divisée par n donne une nouvelle variable nommée la moyenne empirique. Celle-ci, une fois réduite, tend vers une variable normale réduite lorsque n tend vers l'infini.
Autres formulations du théorème
modifierDensités de probabilité
modifierLa densité de probabilité de la somme de plusieurs variables indépendantes s'obtient par convolution de leurs densités (si celles-ci existent). Ainsi on peut interpréter le théorème central limite comme une formulation des propriétés des densités de probabilité soumises à une convolution : sous les conditions établies précédemment, la convoluée d'un certain nombre de densités de probabilité tend vers la densité normale lorsque leur nombre croît indéfiniment.
Comme la fonction caractéristique d'une convolution est le produit des fonctions caractéristiques des variables en cause, le théorème central limite peut se formuler d'une manière différente : sous les conditions précédentes, le produit des fonctions caractéristiques de plusieurs densités de probabilité tend vers la fonction caractéristique de la loi normale lorsque le nombre de variables croît indéfiniment.
Produits de variables aléatoires
modifierLe théorème central limite nous dit à quoi il faut s'attendre en matière de sommes de variables aléatoires indépendantes ; mais qu'en est-il des produits ? Eh bien, le logarithme d'un produit (à facteurs strictement positifs) est la somme des logarithmes des facteurs, de sorte que le logarithme d'un produit de variables aléatoires (à valeurs strictement positives) tend vers une loi normale, ce qui entraîne une loi log-normale pour le produit lui-même.
Bon nombre de grandeurs physiques (en particulier la masse et la longueur, c'est une question de dimension, ne peuvent être négatives) sont le produit de différents facteurs aléatoires, de sorte qu'elles suivent une loi log-normale.
Généralisations du théorème central limite
modifierLe théorème central limite admet plusieurs généralisations qui donnent la convergence de sommes de variables aléatoires sous des hypothèses beaucoup plus faibles. Ces généralisations ne nécessitent pas des lois identiques mais font appel à des conditions qui assurent qu'aucune des variables n'exerce une influence significativement plus importante que les autres. Telles sont la condition de Lindeberg et la condition de Lyapounov. D'autres généralisations autorisent même une dépendance « faible ». De plus, une généralisation due à Gnedenko et Kolmogorov énonce que la somme d'un certain nombre de variables aléatoires avec une queue de distribution décroissante selon avec (ayant donc une variance infinie) tend vers une loi de Lévy tronquée symétrique et stable quand le nombre de variables augmente.
Condition de Liapounov
modifier![](http://upload.wikimedia.org/wikipedia/commons/thumb/9/9b/Somme_tirage_1_a_6d6_et_loi_normale.svg/220px-Somme_tirage_1_a_6d6_et_loi_normale.svg.png)
On peut, au prix d'une formulation un peu moins simple, supprimer l'hypothèse selon laquelle les variables sont de même loi. Les variables restent toutefois indépendantes : soit donc une suite de variables aléatoires définies sur le même espace de probabilité, indépendantes. Supposons que, pour , ait une espérance finie et un écart-type fini , et posons
et
- .
Supposons que pour un certain la condition de Liapounov
soit satisfaite, alors la somme normalisée des converge vers une loi normale centrée réduite, c'est-à-dire :
- .
Condition de Lindeberg
modifierAvec les mêmes définitions et les mêmes notations que précédemment, nous pouvons remplacer la condition de Liapounov par la suivante qui est plus faible[10].
Théorème (Lindeberg, 1920) — Si, pour tout ε > 0
où 1{…} est la fonction indicatrice, alors la loi de converge vers la loi normale centrée réduite .
Cas des variables dépendantes
modifierIl existe quelques théorèmes qui traitent le cas de sommes de variables aléatoires réelles dépendantes, par exemple le théorème central limite pour les suites m-dépendantes, le théorème central limite pour les martingales et le théorème central limite pour les processus mélangeants.
Cas des vecteurs aléatoires
modifierIl existe une généralisation à des vecteurs aléatoires indépendants et de même loi, dont les composantes sont de carrés intégrables, la limite étant alors un vecteur gaussien. Une première version de ce théorème central limite vectoriel, due à Pierre-Simon de Laplace, parait en 1812[11]. Parmi les nombreuses conséquences de ce théorème, on compte par exemple la convergence vers la loi du χ², cruciale, par exemple, pour ses applications en statistiques, ou encore la convergence des marches aléatoires vers le mouvement Brownien.
Notes et références
modifierRéférences
modifier- Pierre-Simon Laplace, « Mémoire sur les approximations des formules qui sont fonctions de très-grands nombres, et sur leur application aux probabilités », Mémoires de la Classe des sciences mathématiques et physiques de l'Institut de France, , p. 353-415 (lire en ligne [PDF], consulté le )
- Pierre-Simon Laplace, « Supplément au mémoire sur les approximations de formules qui sont fonctions de très-grands nombres », Mémoires de la Classe des sciences mathématiques et physiques de l'Institut de France, , p. 559-565 (lire en ligne [PDF], consulté le )
- (en) Stephen M. Stigler, The History of Statistics : The Measurement of Uncertainty before 1900, Harvard, Belknap Press of Harvard University Press, , 1re éd., 432 p. (ISBN 978-0-674-40341-3 et 067440341X, lire en ligne), chap. 2 (« Probabilists and the measurement of uncertainty »). Le cas particulier des variables de Bernoulli est appelé théorème de Moivre-Laplace. Sa démonstration par De Moivre, dans le cas p = 0,5, n'a été possible qu'à travers la démonstration, toujours par De Moivre, de la formule de Stirling.
- (en) Lucien Le Cam, « The central limit theorem around 1935 », Statistical Science, vol. 1, no 1, , p. 78-91 (DOI 10.1214/ss/1177013818, lire en ligne)
- Jean-Yves Ouvrard, Probabilités 1 licence capes, Paris, Cassini, 244 p. (ISBN 978-2-84225-130-7)
- On obtient également une courbe en cloche avec tout dé ayant un nombre de faces différent de 6.
- Illustration extraite du livre Calcul d'incertitudes
- Mathieu Mansuy, Cours Magistral 5 de l'Université de Reims Champagne Ardenne:Fonctions caractéristiques, Université de Reims Champagne Ardenne, (lire en ligne)
- Parreaux 2018.
- (en) William Feller, An Introduction to Probability Theory and Its Applications, vol. 2, New York/Chichester/Brisbane etc., Wiley, , 2e éd., 669 p. (ISBN 0-471-25709-5 et 978-0471257097), p. 262-263.
- Pierre-Simon de Laplace, Théorie analytique des Probabilités, 2e édition, 1812, Livre II, chapitre IV, section 21. [1]
Bibliographie
modifierCours
modifierJulie Parreaux, Théorème central limite, (lire en ligne [PDF]).
Voir aussi
modifierArticle connexe
modifierLiens externes
modifier- Théorème Central Limit, Java
- Central Limit Theorem Simulation interactive pour faire des expériences utilisant plusieurs paramètres.
- Bibliothèque AtelieR pour le logiciel libre R Permet de découvrir le théorème central limite par simulation.
- Article de Jean-Pierre Kahane