ATLAS Ouvrir dans l'atlas →

Théorème 12Analyse de Fourier

Théorème central limite sous un moment d'ordre 2

Énoncé

(i) Restes de l'exponentielle. Pour tout réel yy, ∣eiy−1∣⩽∣y∣|e^{iy} - 1| \leqslant |y|, ∣eiy−1−iy∣⩽y22|e^{iy} - 1 - iy| \leqslant \frac{y^2}2, et ∣eiy−1−iy+y22∣⩽min⁡(∣y∣36,y2).\Big|e^{iy} - 1 - iy + \frac{y^2}2\Big| \leqslant \min\Big(\frac{|y|^3}6, y^2\Big).

(ii) Développement à l'ordre 22. Si E[X2]<+∞E[X^2] < +\infty, alors pour tout réel tt, ∣φX(t)−1−itE[X]+t22E[X2]∣⩽t2 η(t),\begin{aligned} &\Big|\varphi_X(t) - 1 - itE[X] + \frac{t^2}2E[X^2]\Big| \\ &\leqslant t^2\,\eta(t), \end{aligned} où η(t)=E[min⁡(∣t∣ ∣X∣36,X2)]\eta(t) = E\big[\min\big(\frac{|t|\,|X|^3}6, X^2\big)\big] tend vers 00 quand tt tend vers 00.

(iii) Théorème central limite. Soient X1,X2,…X_1, X_2, \dots des variables indépendantes, de même loi, avec E[X1]=0E[X_1] = 0 et E[X12]=1E[X_1^2] = 1, et Sn=X1+⋯+XnS_n = X_1 + \cdots + X_n. Si NN suit la loi 𝒩(0,1)\mathcal N(0, 1), alors pour tout réel xx, P(Snn⩽x)⟶P(N⩽x).P\Big(\frac{S_n}{\sqrt n} \leqslant x\Big) \longrightarrow P(N \leqslant x).

Démonstration

(i) Posons r0(y)=eiyr_0(y) = e^{iy}, r1(y)=eiy−1r_1(y) = e^{iy} - 1, r2(y)=eiy−1−iyr_2(y) = e^{iy} - 1 - iy et r3(y)=eiy−1−iy+y22.r_3(y) = e^{iy} - 1 - iy + \frac{y^2}2.

Chaque rk+1r_{k+1} s'annule en 00 et a pour dérivée i rki\,r_k, puisque ddyeiy=ieiy\frac{d}{dy}e^{iy} = ie^{iy} (Exponentielle, (iv)). Le théorème fondamental (Théorème fondamental) et le fait qu'une fonction de dérivée nulle sur un intervalle est constante (Signe de la dérivée, (ii)), appliqués aux parties réelle et imaginaire, donnent rk+1(y)=i∫0yrk(s) ds.r_{k+1}(y) = i\int_0^y r_k(s)\,ds.

Par Transformée de Fourier, (i), si ∣rk(s)∣⩽c ∣s∣p|r_k(s)| \leqslant c\,|s|^p entre 00 et yy, alors ∣rk+1(y)∣⩽c ∣y∣p+1p+1|r_{k+1}(y)| \leqslant \frac{c\,|y|^{p+1}}{p+1}. Partant de ∣r0∣=1|r_0| = 1 (Exponentielle, (iii)), on obtient successivement ∣r1(y)∣⩽∣y∣|r_1(y)| \leqslant |y|, ∣r2(y)∣⩽y22|r_2(y)| \leqslant \frac{y^2}2 et ∣r3(y)∣⩽∣y∣36|r_3(y)| \leqslant \frac{|y|^3}6.

D'autre part, ∣r2(s)∣⩽∣r1(s)∣+∣s∣⩽2∣s∣|r_2(s)| \leqslant |r_1(s)| + |s| \leqslant 2|s|, d'où ∣r3(y)∣⩽y2|r_3(y)| \leqslant y^2.

(ii) La variable XX est intégrable (Espérance, variance, (ii)). Par linéarité de l'espérance, appliquée aux parties réelle et imaginaire, φX(t)−1−itE[X]+t22E[X2]=E[r3(tX)].\begin{aligned} &\varphi_X(t) - 1 - itE[X] + \frac{t^2}2E[X^2] \\ &= E[r_3(tX)]. \end{aligned}

La preuve de Transformée de Fourier, (i), vaut pour toute probabilité, et donne ∣E[r3(tX)]∣⩽E[∣r3(tX)∣]|E[r_3(tX)]| \leqslant E[|r_3(tX)|]. Par (i), ∣r3(tX)∣⩽t2min⁡(∣t∣ ∣X∣36,X2)|r_3(tX)| \leqslant t^2\min\big(\frac{|t|\,|X|^3}6, X^2\big), d'où l'inégalité annoncée.

La fonction η\eta est croissante en ∣t∣|t|, puisque min⁡(∣t∣ ∣X∣36,X2)\min\big(\frac{|t|\,|X|^3}6, X^2\big) l'est en chaque point. Les variables min⁡(∣X∣36m,X2)\min\big(\frac{|X|^3}{6m}, X^2\big) sont dominées par X2X^2, intégrable, et tendent vers 00 en tout point quand m→∞m \to \infty ; par convergence dominée (Convergence dominée, (i)), η(1m)→0\eta(\frac1m) \to 0.

Pour ε>0\varepsilon > 0, soit m⩾1m \geqslant 1 tel que η(1m)⩽ε\eta(\frac1m) \leqslant \varepsilon ; pour ∣t∣⩽1m|t| \leqslant \frac1m, on a 0⩽η(t)⩽η(1m)⩽ε0 \leqslant \eta(t) \leqslant \eta(\frac1m) \leqslant \varepsilon. Donc η(t)→0\eta(t) \to 0 quand t→0t \to 0, sans passer par les suites quelconques, dont l'usage demanderait le choix dénombrable.

(iii) Produit. Notons φ=φX1\varphi = \varphi_{X_1}, qui est aussi la fonction caractéristique de chaque XkX_k, de même loi. La variable Sn−1S_{n-1} est mesurable pour la tribu engendrée par X1,…,Xn−1X_1, \dots, X_{n-1}, indépendante de σ(Xn)\sigma(X_n) par regroupement (Indépendance, (ii)). Par Fonctions caractéristiques, (ii), et récurrence, φSn=φn\varphi_{S_n} = \varphi^n, d'où φSn/n(t)=φ(tn)n.\varphi_{S_n/\sqrt n}(t) = \varphi\Big(\frac t{\sqrt n}\Big)^n.

Comparaison. Fixons tt et n⩾1n \geqslant 1, et posons a=φ(t/n)a = \varphi(t/\sqrt n) et b=e−t2/(2n)b = e^{-t^2/(2n)}. Comme ∣a∣⩽1|a| \leqslant 1 et ∣b∣⩽1|b| \leqslant 1, l'identité an−bn=(a−b)∑k=0n−1akbn−1−ka^n - b^n = (a - b)\sum_{k=0}^{n-1}a^kb^{n-1-k} donne ∣an−bn∣⩽n ∣a−b∣|a^n - b^n| \leqslant n\,|a - b| ; et bn=e−t2/2b^n = e^{-t^2/2} (Exponentielle, (ii)).

Par (ii), avec E[X1]=0E[X_1] = 0 et E[X12]=1E[X_1^2] = 1, ∣a−(1−t22n)∣⩽t2n η(tn).\Big|a - \Big(1 - \frac{t^2}{2n}\Big)\Big| \leqslant \frac{t^2}n\,\eta\Big(\frac t{\sqrt n}\Big).

Pour x⩾0x \geqslant 0, on a 0⩽e−x−1+x⩽x220 \leqslant e^{-x} - 1 + x \leqslant \frac{x^2}2. En effet, e−x−1+x=∫0x(1−e−s) dse^{-x} - 1 + x = \int_0^x (1 - e^{-s})\,ds, et 0⩽1−e−s=∫0se−v dv⩽s0 \leqslant 1 - e^{-s} = \int_0^s e^{-v}\,dv \leqslant s pour s⩾0s \geqslant 0, puisque 0<e−v⩽10 < e^{-v} \leqslant 1 pour v⩾0v \geqslant 0 (Exponentielle, (iv)). Avec x=t22nx = \frac{t^2}{2n}, il vient ∣b−(1−t22n)∣⩽t48n2\big|b - \big(1 - \frac{t^2}{2n}\big)\big| \leqslant \frac{t^4}{8n^2}.

En rassemblant, ∣φSn/n(t)−e−t2/2∣⩽t2 η(tn)+t48n⟶0.\begin{aligned} &\big|\varphi_{S_n/\sqrt n}(t) - e^{-t^2/2}\big| \\ &\leqslant t^2\,\eta\Big(\frac t{\sqrt n}\Big) + \frac{t^4}{8n} \longrightarrow 0. \end{aligned}

Conclusion. La limite e−t2/2e^{-t^2/2} est la fonction caractéristique de NN (Fonctions caractéristiques, (iii)). Le théorème de Lévy (Continuité de Lévy, (ii)) donne la convergence en loi de Sn/nS_n/\sqrt n vers NN.

La fonction de répartition de NN est continue en tout point, puisque P(N=x)=∫{x}e−u2/22π du=0,P(N = x) = \int_{\{x\}} \frac{e^{-u^2/2}}{\sqrt{2\pi}}\,du = 0, la mesure de Lebesgue de {x}\{x\} étant nulle (Mesure de Lebesgue) : la convergence a lieu pour tout réel xx.

□
Remarque

Ce théorème s'applique là où la version démontrée au nœud Théorème central limite, qui suppose E∣X∣3<+∞E|X|^3 < +\infty, ne s'applique pas. La loi de Student à trois degrés de liberté, de densité 63π(3+x2)2\frac{6\sqrt3}{\pi(3 + x^2)^2}, a une variance égale à 33, mais E∣X∣3=+∞E|X|^3 = +\infty, l'intégrande se comportant comme c∣x∣\frac{c}{|x|} à l'infini.

La variable réduite Y=X/3Y = X/\sqrt3 a pour densité h(x)=2π(1+x2)2h(x) = \frac{2}{\pi(1 + x^2)^2}. Soit g(t)=(1+∣t∣)e−∣t∣g(t) = (1 + |t|)e^{-|t|}, paire, continue et intégrable. Les intégrales ∫0∞e−tcos⁡(tξ) dt=11+ξ2\int_0^\infty e^{-t}\cos(t\xi)\,dt = \frac1{1 + \xi^2} et ∫0∞te−tcos⁡(tξ) dt=1−ξ2(1+ξ2)2\int_0^\infty te^{-t}\cos(t\xi)\,dt = \frac{1 - \xi^2}{(1 + \xi^2)^2} donnent g^(ξ)=4(1+ξ2)2=2πh(ξ)\hat g(\xi) = \frac{4}{(1 + \xi^2)^2} = 2\pi h(\xi), qui est intégrable. La formule d'inversion (Inversion), appliquée à gg, continue, donne alors g(t)=∫h(ξ)eitξ dξ=φY(t)g(t) = \int h(\xi)e^{it\xi}\,d\xi = \varphi_Y(t) pour tout tt.

Par la série de l'exponentielle (Exponentielle), g(t)=1−t22+∣t∣33+O(t4).g(t) = 1 - \frac{t^2}2 + \frac{|t|^3}3 + O(t^4).

Supposons qu'il existe un polynôme PP de degré au plus 33 tel que g=P+o(t3)g = P + o(t^3). Alors P−(1−t22+t33)P - (1 - \frac{t^2}2 + \frac{t^3}3), polynôme de degré au plus 33, est un o(t3)o(t^3) quand t→0+t \to 0^+, donc nul, ses coefficients s'annulant l'un après l'autre ; de même P=1−t22−t33P = 1 - \frac{t^2}2 - \frac{t^3}3 en considérant t→0−t \to 0^-.

Le coefficient de t3t^3 vaudrait à la fois 13\frac13 et −13-\frac13, ce qui est absurde.

La fonction φY\varphi_Y n'a donc pas de développement d'ordre 33, ce qui confirme ici l'absence de moment d'ordre 33 (la réciproque est fausse en général) : seul le développement d'ordre 22 de (ii) est disponible.

t 1 n = 1 n = 10 e−t²/2
Pour des copies indépendantes de YY, les fonctions caractéristiques g(t/n)ng(t/\sqrt n)^n de Sn/nS_n/\sqrt n, pour n=1n = 1 en vert et n=10n = 10 en bleu, se rapprochent de e−t2/2e^{-t^2/2}, en rouge. La convergence est lente : en t=1t = 1, l'écart, qu'on ne voit plus à cette échelle, vaut encore environ 0,0190{,}019 pour n=100n = 100, de l'ordre de e−1/23n\frac{e^{-1/2}}{3\sqrt n}, à cause du terme en ∣t∣3|t|^3.
Remarque

Pour des variables valant 11 et −1-1 avec probabilité 12\frac12, φ=cos⁡\varphi = \cos, et cos⁡(1/n)n\cos(1/\sqrt n)^n vaut environ 0,60600{,}6060 pour n=100n = 100, contre e−1/2≈0,6065e^{-1/2} \approx 0{,}6065. L'écart, voisin de e−1/212n\frac{e^{-1/2}}{12n}, est d'ordre 1n\frac1n, bien plus petit que dans l'exemple précédent : la vitesse dépend des moments supérieurs, que le théorème n'utilise pas, et qu'il ne quantifie pas non plus, contrairement à la majoration de Lindeberg.

Repose sur

Sert à

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