ATLAS Ouvrir dans l'atlas →

Théorème 3Probabilités

Loi faible des grands nombres et théorème de Weierstrass

Énoncé

(i) Loi faible. Soient X1,X2,…X_1, X_2, \dots des variables deux à deux indépendantes, de même loi, de carré intégrable, d'espérance mm et de variance σ2\sigma^2, et Sn=X1+⋯+XnS_n = X_1 + \cdots + X_n. Pour tout ε>0\varepsilon > 0, P(∣Snn−m∣⩾ε)⩽σ2nε2⟶0.P\Big(\Big|\frac{S_n}{n} - m\Big| \geqslant \varepsilon\Big) \leqslant \frac{\sigma^2}{n\varepsilon^2} \longrightarrow 0.

(ii) Weierstrass, par Bernstein. Pour f:[0,1]→ℝf : [0, 1] \to {\mathbb R} continue, les polynômes Bnf(x)=∑k=0nNn,k f(kn) xk(1−x)n−kB_nf(x) = \sum_{k=0}^n N_{n,k}\, f\Big(\frac kn\Big)\, x^k (1 - x)^{n-k} convergent uniformément vers ff sur [0,1][0, 1]. En particulier, toute fonction continue sur un segment est limite uniforme de polynômes.

Démonstration

(i) E[Sn/n]=mE[S_n/n] = m et Var⁡(Sn/n)=σ2/n\operatorname{Var}(S_n/n) = \sigma^2/n (Espérance, variance, (iv)) ; on applique l'inégalité de Bienaymé-Tchebychev.

(ii) Fixons x∈[0,1]x \in [0, 1]. Sur l'ensemble fini Ω={0,1}n\Omega = \{0, 1\}^n, muni de toutes ses parties, posons P({ω})=x∣ω∣(1−x)n−∣ω∣P(\{\omega\}) = x^{|\omega|}(1 - x)^{n - |\omega|}, où ∣ω∣|\omega| est le nombre de 11 ; en développant ∏i(x+(1−x))=1\prod_i \big(x + (1 - x)\big) = 1, ces poids ont pour somme 11.

Les coordonnées ωi\omega_i sont des variables de Bernoulli de paramètre xx, indépendantes, car P(ω1=ε1,…,ωn=εn)P(\omega_1 = \varepsilon_1, \dots, \omega_n = \varepsilon_n) est le produit des xεi(1−x)1−εix^{\varepsilon_i}(1 - x)^{1 - \varepsilon_i}. Leur somme SnS_n vérifie P(Sn=k)=Nn,kxk(1−x)n−kP(S_n = k) = N_{n,k}x^k(1 - x)^{n-k}, d'où Bnf(x)=E[f(Sn/n)].B_nf(x) = E[f(S_n/n)].

0 1 x x − δ x + δ n = 10 hors fenêtre : 0,300 borne x(1 − x)/(nδ2) = 0,933 0 1 x x − δ x + δ n = 100 hors fenêtre : 0,0015 borne x(1 − x)/(nδ2) = 0,0933
Pour x=0,3x = 0{,}3, les poids P(Sn=k)P(S_n = k) sont placés aux points k/nk/n ; ceux de la fenêtre ∣k/n−x∣<δ|k/n - x| < \delta, avec δ=0,15\delta = 0{,}15, sont en bleu, les autres en rouge. La masse rouge, 0,3000{,}300 pour n=10n = 10 et 0,00150{,}0015 pour n=100n = 100, est celle que la preuve majore par Bienaymé-Tchebychev, par x(1−x)/(nδ2)x(1 - x)/(n\delta^2), puis par 1/(4nδ2)1/(4n\delta^2), qui ne dépend plus de xx.

Soit ε>0\varepsilon > 0. Par le théorème de Heine (Théorème de Heine), il existe δ>0\delta > 0 tel que ∣u−v∣<δ|u - v| < \delta entraîne ∣f(u)−f(v)∣<ε|f(u) - f(v)| < \varepsilon. Soit ∥f∥=sup⁡∣f∣\|f\| = \sup |f|, fini par le théorème des bornes atteintes (Bornes atteintes).

En séparant les cas ∣Sn/n−x∣<δ|S_n/n - x| < \delta et ⩾δ\geqslant \delta, et par l'inégalité de Bienaymé-Tchebychev, comme en (i), avec Var⁡(Sn/n)=x(1−x)/n⩽14n\operatorname{Var}(S_n/n) = x(1 - x)/n \leqslant \frac1{4n}, ∣Bnf(x)−f(x)∣⩽E[∣f(Sn/n)−f(x)∣]⩽ε+2∥f∥4nδ2.\begin{aligned} |B_nf(x) - f(x)| &\leqslant E\big[|f(S_n/n) - f(x)|\big] \\ &\leqslant \varepsilon + \frac{2\|f\|}{4n\delta^2}. \end{aligned}

La majoration ne dépend pas de xx : pour nn assez grand, ∣Bnf−f∣<2ε|B_nf - f| < 2\varepsilon sur [0,1][0, 1]. Un segment [a,b][a, b], avec a<ba < b, se ramène à [0,1][0, 1] par la bijection s↦a+(b−a)ss \mapsto a + (b - a)s, qui transforme les polynômes en polynômes.

□
Remarque

La convergence est lente. Pour f(x)=∣x−12∣f(x) = |x - \frac12|, on calcule Bnf(12)=E[∣Sn/n−12∣]≈0,123B_nf(\frac12) = E\big[|S_n/n - \frac12|\big] \approx 0{,}123 pour n=10n = 10 et ≈0,040\approx 0{,}040 pour n=100n = 100 : l'erreur au point de non-dérivabilité décroît comme 1/n1/\sqrt n, alors que f(12)=0f(\frac12) = 0.

1/2 1 1/2 0,123 0,040 f B10f B100f
Au point 12\frac12, où ff n'est pas dérivable, B10fB_{10}f (en bleu) et B100fB_{100}f (en vert) dépassent ff de 0,1230{,}123 et 0,0400{,}040 environ, longueurs des segments rouges : multiplier nn par 1010 ne divise l'erreur que par un peu plus de 33, à peu près 10\sqrt{10}.
Remarque

La démonstration de Bernstein (1912) est l'un des premiers exemples d'une preuve probabiliste d'un théorème d'analyse qui ne parle pas de hasard. L'aléa n'y est qu'un moyen de calcul : Bnf(x)B_nf(x) est une moyenne de valeurs de ff pondérée par des poids qui se concentrent près de xx. L'analyse fonctionnelle refait le calcul sans probabilités, par les seules identités du binôme, et en tire l'estimation L2n\frac{L}{2\sqrt n} pour une fonction LL-lipschitzienne (Polynômes de Bernstein).

Repose sur

Sert à

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