ATLAS Ouvrir dans l'atlas →

ExerciceMaths Sup – MPSI, chapitre 21 — Dénombrement et probabilités · Inégalités probabilistes et loi faible des grands nombres

Exercice 101

Calcul

L'inégalité de Hoeffding pour la loi binomiale. La proposition 24 du cours traite le cas d'une pièce équilibrée ; on étend ici le résultat à une pièce quelconque. Soit p∈]0 ,1[p \in \left]0\,,1\right[ et SnS_{n} une variable de loi ℬ(n,p)\mathcal{B}(n,p).

1. Soit a∈]p ,1]a \in \left]p\,,1\right]. Montrer que, pour tout réel s⩾0s \geqslant 0, P(Snn⩾a)⩽e−nsa(1−p+p es)n.\mathrm{P}\Big(\frac{S_{n}}{n} \geqslant a\Big) \leqslant \mathrm{e}^{-nsa}\big(1 - p + p\,\mathrm{e}^{s}\big)^{n}.

2. Pour a∈]p ,1[a \in \left]p\,,1\right[, déterminer la meilleure de ces majorations, et l'exprimer à l'aide de D(a,p)=(1−a)ln⁡(1−a1−p)+aln⁡(ap).D(a, p) = (1 - a)\ln\Big(\frac{1 - a}{1 - p}\Big) + a\ln\Big(\frac{a}{p}\Big).

3. Pour 0<q<r<10 < q < r < 1, exprimer D(r,q)D(r, q) comme l'intégrale sur [q ,r][q\,,r] de la fonction t↦r−tt(1−t)t \mapsto \frac{r - t}{t(1 - t)}, et en déduire que D(r,q)⩾2(r−q)2D(r, q) \geqslant 2(r - q)^{2}.

4. En déduire que P(∣Snn−p∣⩾ε)⩽2 e−2nε2\mathrm{P}\Big(\Big|\dfrac{S_{n}}{n} - p\Big| \geqslant \varepsilon\Big) \leqslant 2\,\mathrm{e}^{-2n\varepsilon^{2}} pour tout ε>0\varepsilon > 0, et comparer avec la proposition 24.

Corrigé

1. Soit s⩾0s \geqslant 0. Si s=0s = 0, le majorant vaut 11, et l'inégalité est évidente. Si s>0s > 0, la fonction x↦esxx \mapsto \mathrm{e}^{sx} est strictement croissante, et les événements (Snn⩾a)\big(\frac{S_{n}}{n} \geqslant a\big) et (esSn⩾ensa)\big(\mathrm{e}^{sS_{n}} \geqslant \mathrm{e}^{nsa}\big) sont égaux. L'inégalité de Markov, appliquée à la variable positive esSn\mathrm{e}^{sS_{n}}, donne P(Snn⩾a)⩽e−nsa E(esSn).\mathrm{P}\Big(\frac{S_{n}}{n} \geqslant a\Big) \leqslant \mathrm{e}^{-nsa}\,\mathrm{E}\big(\mathrm{e}^{sS_{n}}\big). La formule de transfert (théorème 10) et la formule du binôme donnent E(esSn)=∑k=0nesk(nk)pk(1−p)n−k=∑k=0n(nk)(p es)k(1−p)n−k=(1−p+p es)n,\mathrm{E}\big(\mathrm{e}^{sS_{n}}\big) = \sum_{k=0}^{n} \mathrm{e}^{sk}\binom{n}{k}p^{k}(1 - p)^{n-k} = \sum_{k=0}^{n} \binom{n}{k}\big(p\,\mathrm{e}^{s}\big)^{k}(1 - p)^{n-k} = \big(1 - p + p\,\mathrm{e}^{s}\big)^{n},

P(Snn⩾a)⩽e−nsa(1−p+p es)npour tout s⩾0\boxed{\mathrm{P}\Big(\frac{S_{n}}{n} \geqslant a\Big) \leqslant \mathrm{e}^{-nsa}\big(1 - p + p\,\mathrm{e}^{s}\big)^{n} \quad\text{pour tout } s \geqslant 0}

Piège. L'hypothèse s⩾0s \geqslant 0 est indispensable : pour s<0s < 0, l'exponentielle est décroissante et l'inclusion change de sens. Pour n=10n = 10, p=12p = \frac{1}{2} et a=14a = \frac{1}{4}, l'énoncé étendu à s<0s < 0 donnerait un majorant d'environ 0,270{,}27, alors que P(S10⩾2,5)≈0,945\mathrm{P}(S_{10} \geqslant 2{,}5) \approx 0{,}945 ✓.

2. Minimiser le majorant revient à minimiser son logarithme divisé par nn, h(s)=−sa+ln⁡(1−p+p es),h′(s)=−a+p es1−p+p es=−a+1−1−p1−p+p es.h(s) = -sa + \ln\big(1 - p + p\,\mathrm{e}^{s}\big), \qquad h'(s) = -a + \frac{p\,\mathrm{e}^{s}}{1 - p + p\,\mathrm{e}^{s}} = -a + 1 - \frac{1 - p}{1 - p + p\,\mathrm{e}^{s}}. La fonction h′h' est strictement croissante, h′(0)=p−a<0h'(0) = p - a < 0, et h′(s)h'(s) tend vers 1−a>01 - a > 0 en +∞+\infty. Elle s'annule en un unique point s∗s^{*}, défini par p es∗(1−a)=a(1−p)p\,\mathrm{e}^{s^{*}}(1 - a) = a(1 - p) : es∗=a(1−p)p(1−a)>1,\mathrm{e}^{s^{*}} = \frac{a(1 - p)}{p(1 - a)} > 1, puisque a(1−p)−p(1−a)=a−p>0a(1 - p) - p(1 - a) = a - p > 0. La fonction hh décroît sur [0 ,s∗][0\,,s^{*}] et croît ensuite : son minimum sur ℝ+\mathbb{R}_{+} est h(s∗)h(s^{*}). Or 1−p+p es∗=(1−p)(1+a1−a)=1−p1−a1 - p + p\,\mathrm{e}^{s^{*}} = (1 - p)\big(1 + \frac{a}{1 - a}\big) = \frac{1 - p}{1 - a}, d'où h(s∗)=−aln⁡a(1−p)p(1−a)+ln⁡1−p1−a=−aln⁡ap−(1−a)ln⁡1−a1−p=−D(a,p).h(s^{*}) = -a\ln\frac{a(1 - p)}{p(1 - a)} + \ln\frac{1 - p}{1 - a} = -a\ln\frac{a}{p} - (1 - a)\ln\frac{1 - a}{1 - p} = -D(a, p).

la meilleure majoration est P(Snn⩾a)⩽e−nD(a,p)\boxed{\text{la meilleure majoration est } \mathrm{P}\Big(\frac{S_{n}}{n} \geqslant a\Big) \leqslant \mathrm{e}^{-nD(a, p)}}

3. Pour t∈[q ,r]t \in [q\,,r], la décomposition en éléments simples (chapitre 13) donne r−tt(1−t)=rt+r−11−t,\frac{r - t}{t(1 - t)} = \frac{r}{t} + \frac{r - 1}{1 - t}, ce que l'on vérifie en réduisant au même dénominateur : r(1−t)+(r−1)t=r−tr(1 - t) + (r - 1)t = r - t. Comme 0<q⩽t⩽r<10 < q \leqslant t \leqslant r < 1, ∫qrr−tt(1−t) dt=rln⁡rq+(r−1)[−ln⁡(1−t)]qr=rln⁡rq+(1−r)ln⁡1−r1−q,\int_{q}^{r} \frac{r - t}{t(1 - t)}\,\mathrm{d} t = r\ln\frac{r}{q} + (r - 1)\big[-\ln(1 - t)\big]_{q}^{r} = r\ln\frac{r}{q} + (1 - r)\ln\frac{1 - r}{1 - q},

D(r,q)=∫qrr−tt(1−t) dt\boxed{D(r, q) = \int_{q}^{r} \frac{r - t}{t(1 - t)}\,\mathrm{d} t} Pour t∈[q ,r]t \in [q\,,r], on a 0<t(1−t)⩽140 < t(1 - t) \leqslant \frac{1}{4} et r−t⩾0r - t \geqslant 0, donc r−tt(1−t)⩾4(r−t)\frac{r - t}{t(1 - t)} \geqslant 4(r - t). Par croissance de l'intégrale (chapitre 18), les bornes étant dans l'ordre croissant, D(r,q)⩾4∫qr(r−t) dt=2(r−q)2\boxed{D(r, q) \geqslant 4\int_{q}^{r} (r - t)\,\mathrm{d} t = 2(r - q)^{2}}

4. Soit ε>0\varepsilon > 0 et a=p+εa = p + \varepsilon. Montrons d'abord que P(Snn−p⩾ε)⩽e−2nε2\mathrm{P}\big(\frac{S_{n}}{n} - p \geqslant \varepsilon\big) \leqslant \mathrm{e}^{-2n\varepsilon^{2}}.

• Si a>1a > 1, l'événement est impossible, puisque Sn⩽nS_{n} \leqslant n.

• Si a<1a < 1, les questions 2 et 3, avec r=ar = a et q=pq = p, donnent e−nD(a,p)⩽e−2nε2\mathrm{e}^{-nD(a, p)} \leqslant \mathrm{e}^{-2n\varepsilon^{2}}.

• Si a=1a = 1, l'événement est (Sn=n)(S_{n} = n), de probabilité pn=e−nln⁡(1/p)p^{n} = \mathrm{e}^{-n\ln(1/p)}, et le calcul de la question 3 avec r=1r = 1 s'applique encore : ln⁡1p=∫p1dtt⩾∫p14(1−t) dt=2(1−p)2\ln\frac{1}{p} = \int_{p}^{1} \frac{\mathrm{d} t}{t} \geqslant \int_{p}^{1} 4(1 - t)\,\mathrm{d} t = 2(1 - p)^{2}, puisque 1t=1−tt(1−t)⩾4(1−t)\frac{1}{t} = \frac{1 - t}{t(1 - t)} \geqslant 4(1 - t).

Pour l'autre côté, la variable n−Snn - S_{n} suit la loi ℬ(n,1−p)\mathcal{B}(n,1 - p), puisque P(n−Sn=k)=P(Sn=n−k)=(nk)(1−p)kpn−k.\mathrm{P}(n - S_{n} = k) = \mathrm{P}(S_{n} = n - k) = \binom{n}{k}(1 - p)^{k}p^{n-k}. Le résultat précédent, qui ne dépend que de la loi, appliqué à n−Snn - S_{n} et au paramètre 1−p1 - p, donne P(Snn−p⩽−ε)=P(n−Snn−(1−p)⩾ε)⩽e−2nε2.\mathrm{P}\Big(\frac{S_{n}}{n} - p \leqslant -\varepsilon\Big) = \mathrm{P}\Big(\frac{n - S_{n}}{n} - (1 - p) \geqslant \varepsilon\Big) \leqslant \mathrm{e}^{-2n\varepsilon^{2}}. L'événement (∣Snn−p∣⩾ε)\big(\big|\frac{S_{n}}{n} - p\big| \geqslant \varepsilon\big) est la réunion des deux, et la sous-additivité conclut : P(∣Snn−p∣⩾ε)⩽2 e−2nε2\boxed{\mathrm{P}\Big(\Big|\frac{S_{n}}{n} - p\Big| \geqslant \varepsilon\Big) \leqslant 2\,\mathrm{e}^{-2n\varepsilon^{2}}} Pour p=12p = \frac{1}{2}, c'est exactement la proposition 24, obtenue ici sans la majoration du cosinus hyperbolique, et le résultat vaut maintenant pour tout pp. La forme intermédiaire e−nD(a,p)\mathrm{e}^{-nD(a, p)} est même bien meilleure loin de 12\frac{1}{2}, où t(1−t)t(1 - t) est petit devant 14\frac{1}{4}. Pour n=1 000n = 1\,000, p=0,1p = 0{,}1 et ε=0,05\varepsilon = 0{,}05, la probabilité exacte vaut environ 4,5×10−74{,}5 \times 10^{-7} ; la somme des deux majorants e−nD\mathrm{e}^{-nD} vaut environ 4,9×10−64{,}9 \times 10^{-6}, la borne 2 e−2nε22\,\mathrm{e}^{-2n\varepsilon^{2}} environ 0,0130{,}013, et la borne p(1−p)nε2\frac{p(1 - p)}{n\varepsilon^{2}} du théorème 16 vaut 0,0360{,}036 ✓.

À retenir. La méthode est celle de la proposition 24 : l'inégalité de Markov appliquée à esSn\mathrm{e}^{sS_{n}}, puis l'optimisation en ss. Le choix de l'exponentielle transforme la somme de variables indépendantes en produit, dont l'espérance se calcule par la proposition 20 ou, comme ici, par la formule du binôme.

Ouvrir dans l'atlas interactif →