ATLAS Ouvrir dans l'atlas →

Proposition 1Martingales

Série harmonique aléatoire et urne de Pólya

Énoncé

(i) Si (εn)(\varepsilon_n) sont indépendantes et valent ±1\pm 1 avec probabilité 12\frac12, la série ∑n⩾1εn/n\sum_{n \geqslant 1} \varepsilon_n/n converge presque sûrement.

(ii) Urne de Pólya. Une urne contient au départ une boule rouge et une bleue ; à chaque étape, on tire une boule uniformément, et on la remet avec une boule supplémentaire de même couleur. La proportion XnX_n de boules rouges après nn tirages est une martingale, qui converge presque sûrement.

Démonstration

(i) Les sommes partielles Mn=∑k⩽nεk/kM_n = \sum_{k \leqslant n}\varepsilon_k/k forment une martingale (Martingales, (i), première affirmation), et E[Mn2]=∑k⩽n1/k2E[M_n^2] = \sum_{k \leqslant n}1/k^2 est borné. On applique Convergence des martingales. Des variables εn\varepsilon_n existent sur [0,1[[0, 1[ : εn=2dn−1\varepsilon_n = 2d_n - 1, où dnd_n est le nn-ième chiffre binaire (Nombres normaux).

(ii) On réalise l'urne à partir d'une suite de variables uniformes indépendantes UnU_n, qui existe sur [0,1[[0, 1[ (Suites indépendantes, (i)) : au tirage n+1n + 1, on tire une rouge si Un+1⩽XnU_{n+1} \leqslant X_n.

Après nn tirages, l'urne contient n+2n + 2 boules, dont RnR_n rouges, et Xn=Rn/(n+2)X_n = R_n/(n + 2) ; on prend ℱn=σ(U1,…,Un)\mathcal F_n = \sigma(U_1, \dots, U_n), et RnR_n est ℱn\mathcal F_n-mesurable, par récurrence.

Comme XnX_n ne prend que les valeurs r/(n+2)r/(n + 2), 1⩽r⩽n+11 \leqslant r \leqslant n + 1, 1{Un+1⩽Xn}=∑r1{Xn=r/(n+2)} 1{Un+1⩽r/(n+2)},\begin{aligned} &\mathbf 1_{\{U_{n+1} \leqslant X_n\}} \\ &= \sum_r \mathbf 1_{\{X_n = r/(n+2)\}}\,\mathbf 1_{\{U_{n+1} \leqslant r/(n+2)\}}, \end{aligned} où chaque 1{Xn=r/(n+2)}\mathbf 1_{\{X_n = r/(n+2)\}} est bornée et ℱn\mathcal F_n-mesurable, et Un+1U_{n+1} est indépendante de ℱn\mathcal F_n par regroupement (Indépendance, (ii)).

Par Espérance conditionnelle, (iv), l'espérance conditionnelle de ce tirage sachant ℱn\mathcal F_n vaut ∑r1{Xn=r/(n+2)} rn+2=Xn,\sum_r \mathbf 1_{\{X_n = r/(n+2)\}}\,\frac r{n+2} = X_n, d'où E[Rn+1∣ℱn]=Rn+XnE[R_{n+1} \mid \mathcal F_n] = R_n + X_n et E[Xn+1∣ℱn]=Rn(n+3)/(n+2)n+3=Xn.\begin{aligned} E[X_{n+1} \mid \mathcal F_n] &= \frac{R_n(n + 3)/(n + 2)}{n + 3} \\ &= X_n. \end{aligned}

Comme 0⩽Xn⩽10 \leqslant X_n \leqslant 1, E[Xn2]⩽1E[X_n^2] \leqslant 1, et la martingale converge presque sûrement (Convergence des martingales).

□
Remarque

La limite X∞X_\infty de l'urne de Pólya suit la loi uniforme sur [0,1][0, 1]. D'abord, R0=1R_0 = 1 ; si RnR_n est uniforme sur {1,…,n+1}\{1, \dots, n + 1\}, le tirage n+1n + 1 est rouge avec la probabilité rn+2\frac r{n+2} sur {Rn=r}\{R_n = r\}, par l'indépendance de Un+1U_{n+1} et de ℱn\mathcal F_n, d'où P(Rn+1=r)=1n+1(r−1n+2+n+2−rn+2)=1n+2\begin{aligned} &P(R_{n+1} = r) \\ &= \frac1{n+1}\Big(\frac{r - 1}{n + 2} + \frac{n + 2 - r}{n + 2}\Big) = \frac1{n+2} \end{aligned} pour 1⩽r⩽n+21 \leqslant r \leqslant n + 2. Par récurrence, RnR_n est uniforme sur {1,…,n+1}\{1, \dots, n + 1\} pour tout nn, et P(Xn⩽x)=min⁡(⌊(n+2)x⌋,n+1)/(n+1)P(X_n \leqslant x) = \min\big(\lfloor (n + 2)x \rfloor, n + 1\big)/(n + 1) tend vers xx pour x∈[0,1]x \in [0, 1].

Soit FF la fonction de répartition de X∞X_\infty, qui est à valeurs dans [0,1][0, 1]. Pour ε>0\varepsilon > 0, l'événement {X∞⩽x−ε}\{X_\infty \leqslant x - \varepsilon\} est contenu dans {Xn⩽x}∪{∣Xn−X∞∣⩾ε}\{X_n \leqslant x\} \cup \{|X_n - X_\infty| \geqslant \varepsilon\}, et {Xn⩽x}\{X_n \leqslant x\} dans {X∞⩽x+ε}∪{∣Xn−X∞∣⩾ε}\{X_\infty \leqslant x + \varepsilon\} \cup \{|X_n - X_\infty| \geqslant \varepsilon\}. La convergence presque sûre entraînant la convergence en probabilité (Modes de convergence, (i)), il vient F(x−ε)⩽x⩽F(x+ε)F(x - \varepsilon) \leqslant x \leqslant F(x + \varepsilon) pour x∈[0,1]x \in [0, 1] ; donc F(y)=yF(y) = y sur [0,1][0, 1], FF est nulle sur ]−∞,0[]{-\infty}, 0[ et vaut 11 sur ]1,+∞[]1, +\infty[, et c'est la fonction de répartition de la loi uniforme, qui est donc la loi de X∞X_\infty (Variables aléatoires, (i)).

Chaque urne converge, mais vers une proportion aléatoire, fixée par les tout premiers tirages ; c'est un modèle de renforcement, où le hasard initial devient une propriété permanente.

0 1/2 1 1 10 100 500 0 n Xn
Six urnes simulées selon la règle de la preuve, sur 500500 tirages, l'axe horizontal étant gradué en échelle logarithmique : les proportions se séparent dès les premiers tirages, puis chacune se stabilise vers sa propre limite.
Remarque

La série harmonique alternée ∑(−1)n+1/n\sum (-1)^{n+1}/n converge, vers une valeur précise, alors que ∑1/n\sum 1/n diverge. Le théorème (i) dit que des signes tirés au hasard suffisent presque sûrement à la convergence ; la somme, elle, est une variable aléatoire non constante, de variance ∑n⩾11/n2\sum_{n \geqslant 1} 1/n^2, puisque la convergence a aussi lieu dans L2{L^2} (Convergence des martingales).

n −1 0 1 1 50 100 150 200
Les sommes partielles ∑k⩽nεk/k\sum_{k \leqslant n}\varepsilon_k/k pour quatre tirages des signes, en rouge, et celles de la série alternée, en pointillé : chaque trajectoire se stabilise, mais vers une limite qui dépend du tirage.

Repose sur

Sert à

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