ATLAS Ouvrir dans l'atlas →

Théorème 7Probabilités

Suites de variables indépendantes de lois données

On reprend l'espace Ω=[0,1[\Omega = [0, 1[, muni de la mesure de Lebesgue, et les chiffres binaires dkd_k du nœud Nombres normaux. Pour une fonction de répartition FF et u∈ ]0,1[u \in \,]0, 1[, on pose G(u)=inf⁡{x∈ℝ:F(x)⩾u}G(u) = \inf\{x \in {\mathbb R} : F(x) \geqslant u\} : c'est la fonction quantile de FF.

Énoncé

(i) Pour m∈ℕm \in {\mathbb N}, posons Um=∑j⩾0d2m(2j+1) 2−(j+1)U_m = \sum_{j \geqslant 0} d_{2^m(2j+1)}\,2^{-(j+1)}. Les UmU_m sont des variables indépendantes, de loi uniforme sur [0,1][0, 1].

(ii) Si UU est une variable de loi uniforme sur [0,1][0, 1], la variable XX égale à G(U)G(U) sur {0<U<1}\{0 < U < 1\} et à 00 ailleurs a pour fonction de répartition FF.

(iii) Pour toute suite (μm)(\mu_m) de probabilités sur ℝ\mathbb R, il existe sur Ω\Omega une suite de variables indépendantes XmX_m de lois μm\mu_m. En particulier, il existe des suites de variables indépendantes et de même loi, de loi quelconque.

Démonstration

(i) Les sous-suites de chiffres. Tout entier k⩾1k \geqslant 1 s'écrit de manière unique 2m(2j+1)2^m(2j + 1) avec m,j∈ℕm, j \in \mathbb N (ℚ est dénombrable, (i)) : les ensembles Im={2m(2j+1):j∈ℕ}I_m = \{2^m(2j + 1) : j \in \mathbb N\} sont disjoints et recouvrent ℕ∖{0}\mathbb N \setminus \{0\}.

Chaque UmU_m est mesurable, comme limite simple de ses sommes partielles, qui sont des sommes finies de fonctions mesurables (Fonctions mesurables, (ii) et (iii)).

De plus, chaque dkd_k, k∈Imk \in I_m, est mesurable pour la tribu ℬm=σ(⋃k∈Imσ(dk)),\mathcal B_m = \sigma\Big(\bigcup_{k \in I_m}\sigma(d_k)\Big), et le même argument, appliqué sur (Ω,ℬm)(\Omega, \mathcal B_m), montre que UmU_m est ℬm\mathcal B_m-mesurable : σ(Um)⊂ℬm\sigma(U_m) \subset \mathcal B_m.

Les dkd_k étant indépendants, les ℬm\mathcal B_m le sont par regroupement (Indépendance, (ii)), donc aussi les UmU_m.

La loi. Fixons mm et posons ej=d2m(2j−1)e_j = d_{2^m(2j - 1)} pour j⩾1j \geqslant 1 : ce sont des variables de Bernoulli de paramètre 12\frac12 indépendantes, et V=∑jej2−j=UmV = \sum_j e_j2^{-j} = U_m, à valeurs dans [0,1][0, 1].

Pour n⩾1n \geqslant 1 et 0⩽M<2n0 \leqslant M < 2^n, soit EME_M l'événement « les nn premiers eje_j sont les chiffres binaires de MM », de probabilité 2−n2^{-n} par indépendance ; sur EME_M, V∈[M2−n,(M+1)2−n]V \in [M2^{-n}, (M + 1)2^{-n}], la queue de la série étant comprise entre 00 et 2−n2^{-n}.

Soit x=N2−nx = N2^{-n} un dyadique de [0,1][0, 1]. On a ⋃M<NEM⊂{V⩽x}\bigcup_{M < N} E_M \subset \{V \leqslant x\}, d'où P(V⩽x)⩾xP(V \leqslant x) \geqslant x.

Si V⩽xV \leqslant x sur EME_M avec M⩾NM \geqslant N, alors V=xV = x ; et, pour tout n′⩾nn' \geqslant n, l'égalité V=xV = x impose aux n′n' premiers chiffres d'être ceux de N2n′−nN2^{n'-n} ou de N2n′−n−1N2^{n'-n} - 1, d'où P(V=x)⩽2⋅2−n′P(V = x) \leqslant 2 \cdot 2^{-n'}, puis P(V=x)=0P(V = x) = 0. Ainsi P(V⩽x)=xP(V \leqslant x) = x pour tout dyadique xx de [0,1][0, 1].

Pour x∈[0,1]x \in [0, 1] quelconque, des dyadiques x′⩽x⩽x′′x' \leqslant x \leqslant x'' arbitrairement proches de xx et la croissance de FVF_V donnent FV(x)=xF_V(x) = x ; enfin FV=0F_V = 0 sur ]−∞,0[]{-\infty}, 0[ et FV=1F_V = 1 sur [1,+∞[[1, +\infty[. C'est la fonction de répartition de la loi uniforme, qui est donc la loi de VV (Variables aléatoires, (i)).

(ii) Pour u∈ ]0,1[u \in \,]0, 1[, l'ensemble Au={x:F(x)⩾u}A_u = \{x : F(x) \geqslant u\} est non vide et minoré, puisque FF tend vers 1 en +∞+\infty et vers 00 en −∞-\infty ; il a donc une borne inférieure G(u)G(u), par la propriété de la borne supérieure appliquée aux opposés.

Comme FF est croissante, AuA_u contient toute demi-droite [x,+∞[[x, +\infty[ issue de l'un de ses points ; il contient donc G(u)+1kG(u) + \frac1k pour tout k⩾1k \geqslant 1, et, FF étant continue à droite, F(G(u))=lim⁡kF(G(u)+1k)⩾uF(G(u)) = \lim_k F(G(u) + \frac1k) \geqslant u : G(u)∈AuG(u) \in A_u. Donc, pour tout réel xx, G(u)⩽xG(u) \leqslant x si et seulement si u⩽F(x)u \leqslant F(x).

Par suite X=h(U)X = h(U), où h(u)=G(u)h(u) = G(u) pour 0<u<10 < u < 1 et h(u)=0h(u) = 0 sinon, et {h⩽x}=Bx\{h \leqslant x\} = B_x, avec Bx= ]0,1[ ∩ ]−∞,F(x)]B_x = \,]0, 1[ \,\cap\, ]{-\infty}, F(x)] si x<0x < 0 et Bx=(]0,1[ ∩ ]−∞,F(x)])∪(ℝ∖ ]0,1[)\begin{aligned} B_x = {}&(]0, 1[ \,\cap\, ]{-\infty}, F(x)]) \\ &\cup (\mathbb R \setminus \,]0, 1[) \end{aligned} si x⩾0x \geqslant 0 : ces ensembles sont boréliens, donc hh est borélienne (Fonctions mesurables, (i)), XX est mesurable et σ(X)⊂σ(U)\sigma(X) \subset \sigma(U).

Comme P(U∈ ]0,1[)=1P(U \in \,]0, 1[) = 1, P(X⩽x)=P(0<U⩽F(x))=F(x),P(X \leqslant x) = P(0 < U \leqslant F(x)) = F(x), F(x)F(x) étant dans [0,1][0, 1] : la fonction de répartition de XX est FF.

(iii) Soit FmF_m la fonction de répartition de μm\mu_m, GmG_m sa fonction quantile, et XmX_m construite à partir de UmU_m comme en (ii). La loi de XmX_m a pour fonction de répartition FmF_m, donc c'est μm\mu_m (Variables aléatoires, (i)). Chaque Xm=hm(Um)X_m = h_m(U_m) est une fonction borélienne de UmU_m seul, et σ(Xm)⊂σ(Um)\sigma(X_m) \subset \sigma(U_m) par (ii) ; des sous-tribus de tribus indépendantes étant indépendantes, les XmX_m le sont.

□
Remarque

La répartition des chiffres est fixée par l'écriture k=2m(2j+1)k = 2^m(2j + 1) : U0U_0 utilise les chiffres d'indices impairs 1,3,5,…1, 3, 5, \dots, U1U_1 ceux d'indices 2,6,10,…2, 6, 10, \dots, U2U_2 ceux d'indices 4,12,20,…4, 12, 20, \dots

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 … indice k U0 U1 U2 U3 U4
Les vingt-quatre premiers chiffres binaires de ω\omega, colorés selon la variable qui les utilise : le chiffre dkd_k va à UmU_m, où 2m2^m est la plus grande puissance de 22 qui divise kk. Chaque variable reçoit une infinité de chiffres, et deux variables n'en partagent aucun.

Ainsi, pour ω=0,101101…\omega = 0{,}101101\ldots en base 22, U0U_0 commence par les chiffres d1,d3,d5=1,1,0d_1, d_3, d_5 = 1, 1, 0, et U1U_1 par d2,d6=0,1d_2, d_6 = 0, 1, soit U0=0,110…U_0 = 0{,}110\ldots et U1=0,01…U_1 = 0{,}01\ldots en base 22.

Remarque

L'erreur la plus répandue est de croire que des variables uniformes construites sur le même ω\omega sont automatiquement indépendantes. Les variables ω↦2mω−⌊2mω⌋\omega \mapsto 2^m\omega - \lfloor 2^m\omega \rfloor, qui décalent les chiffres de ω\omega, sont toutes uniformes, mais très dépendantes : la deuxième se lit sur la première, dont elle est le double modulo 11. C'est le partage des chiffres en paquets disjoints qui produit l'indépendance.

Remarque

Le théorème donne un sens aux hypothèses des lois des grands nombres et du théorème central limite, qui portent sur des suites indépendantes de même loi : de telles suites existent. Pour une famille non dénombrable de variables, ou des lois sur des espaces plus généraux, il faut construire une probabilité produit infinie (théorème des produits infinis d'Andersen et Jessen, ou d'Ionescu Tulcea, ou théorème d'extension de Kolmogorov pour des espaces polonais), admis ici ; sur [0,1[[0, 1[ lui-même, une famille non dénombrable de variables indépendantes dont aucune n'est presque sûrement constante est impossible.

En effet, si XiX_i n'est pas presque sûrement constante, soit qiq_i le premier rationnel, dans une énumération fixée de ℚ{\mathbb Q}, tel que pi=P(Xi⩽qi)p_i = P(X_i \leqslant q_i) vérifie 0<pi<10 < p_i < 1 : il existe, car si P(Xi⩽q)P(X_i \leqslant q) valait 00 ou 11 pour tout rationnel qq, XiX_i serait presque sûrement égale à la borne supérieure des rationnels qq où cette probabilité est nulle. Les variables Yi=(1{Xi⩽qi}−pi)/pi(1−pi)Y_i = (\mathbf 1_{\{X_i \leqslant q_i\}} - p_i)/\sqrt{p_i(1 - p_i)} sont centrées, de variance 11 et indépendantes, donc E[(Yi−Yj)2]=2E[(Y_i - Y_j)^2] = 2 pour i≠ji \neq j (Espérance, variance, (iv)).

Or les fonctions en escalier à valeurs rationnelles sur les intervalles dyadiques [k2−n,(k+1)2−n[[k2^{-n}, (k + 1)2^{-n}[ forment une suite (sl)(s_l) dense dans L2([0,1[)L^2([0, 1[) : les fonctions continues sont denses (Densité, (iii), en analyse fonctionnelle), et chacune est uniformément approchée par ces fonctions en escalier, par le théorème de Heine (Théorème de Heine). Associer à ii le premier ll tel que E[(Yi−sl)2]<12E[(Y_i - s_l)^2] < \frac12 est alors injectif, par l'inégalité triangulaire : la famille est au plus dénombrable, et aucun choix n'intervient.

Les étapes propres à cette construction, la répartition des chiffres et la fonction quantile, n'utilisent aucun choix ; seul intervient, en amont, le choix dénombrable sur lequel repose la mesure de Lebesgue.

Repose sur

Sert à

Aucune fiche de ce graphe ne l'utilise encore.
Ouvrir dans l'atlas interactif →