📚 الرئيسية

مسابقة دكتوراه 2021Université Abou Bekr Belkaïd - Tlemcen — الموضوع 01

مسابقة تخصص · Probabilités & Statistiques · المدة: 2سا

Concours d'accès à la formation de 3ème cycle LMD 2020/2021, Domaine Mathématiques-Informatique, Filière Mathématiques, Épreuve de spécialité : Probabilités et Statistiques, Université Aboubekr Belkaïd – Tlemcen, Faculté des Sciences – Tidjani Haddam, durée 2h00, le samedi 06 mars 2021.

التمرين 1

Problème 1 — Troncature, inégalité de Kolmogorov et loi forte de Marcinkiewicz

#probability-theory#strong-law-of-large-numbers#almost-sure-convergence#kolmogorov-inequality#borel-cantelli

Soient X={Xn,n1}X=\{X_n,\,n\geq 1\} et Y={Yn,n1}Y=\{Y_n,\,n\geq 1\} deux suites de variables aléatoires définies sur un même espace de probabilité complet (Ω,F,P)(\Omega,\mathcal{F},P). On note {an,n1}\{a_n,\,n\geq 1\} une suite croissante de nombres réels strictement positifs avec limnan=+\lim_{n\to\infty}a_n=+\infty. E(X)E(X) désigne l'espérance mathématique de XX et 1(A)\mathbf{1}(A) la fonction indicatrice de l'événement AA. Pour n1n\geq 1, on pose

Sn=k=1nXketXn=Xn1(Xnan).S_n=\sum_{k=1}^n X_k\qquad\text{et}\qquad \overline{X}_n=X_n\,\mathbf{1}(|X_n|\leq a_n).

(p.s. : presque sûrement ; v.a.r. : variable aléatoire réelle.)

Première partie.

  1. Soit ff une fonction réelle borélienne et impaire. Montrer que si UU est une v.a.r. symétrique (UU et U-U ont même loi), alors V:=f(U)V:=f(U) est aussi une v.a.r. symétrique.
  2. Montrer que si k1P(XkXk)<\sum_{k\geq 1}P(X_k\neq \overline{X}_k)\lt\infty et limnSnan=0\lim_{n\to\infty}\dfrac{S_n}{a_n}=0 p.s., alors limn1ank=1nXk=0\lim_{n\to\infty}\dfrac{1}{a_n}\sum_{k=1}^n \overline{X}_k=0 p.s.

Deuxième partie.

  1. Montrer que si X={Xn,n1}X=\{X_n,\,n\geq 1\} est une suite de v.a.r. indépendantes avec XnL2X_n\in L^2 et E(Xn)=0E(X_n)=0 pour n1n\geq 1, alors

n1E(Xn2)<  {Sn,n1} converge dans L2 et p.s.\sum_{n\geq 1}E(X_n^2)\lt\infty\ \Longrightarrow\ \{S_n,\,n\geq 1\}\ \text{converge dans }L^2\ \text{et p.s.}

Troisième partie.

On suppose à présent que X={Xn,n1}X=\{X_n,\,n\geq 1\} est une suite de v.a.r. indépendantes, symétriques et identiquement distribuées, et qu'il existe C>0C\gt 0 et n01n_0\geq 1 tels que

nn0,k=n1ak2Cnan2.\forall n\geq n_0,\qquad \sum_{k=n}^{\infty}\frac{1}{a_k^2}\leq \frac{C\,n}{a_n^2}.

  1. Montrer que si k1P(X1>ak)<\sum_{k\geq 1}P(|X_1|\gt a_k)\lt\infty, alors {k=1nXkak,n1}\Big\{\sum_{k=1}^n \dfrac{\overline{X}_k}{a_k},\,n\geq 1\Big\} converge p.s.
  2. En déduire que si 12<α<1\tfrac12\lt\alpha\lt 1, alors Snnα0\dfrac{S_n}{n^\alpha}\to 0 p.s.     k1P(X1>kα)<\iff \sum_{k\geq 1}P(|X_1|\gt k^\alpha)\lt\infty.
  3. Montrer que cette dernière équivalence ne peut pas être vraie pour α=12\alpha=\tfrac12.

Rappels.

Inégalité de Kolmogorov : soit X={Xn}X=\{X_n\} une suite de v.a.r. indépendantes avec XnL2X_n\in L^2 et E(Xn)=0E(X_n)=0. Alors pour tout ε>0\varepsilon\gt 0 et tout n1n\geq 1,

P(max1knj=1kXj>ε)1ε2j=1nE(Xj2).P\Big(\max_{1\leq k\leq n}\big|\textstyle\sum_{j=1}^k X_j\big|\gt\varepsilon\Big)\leq \frac{1}{\varepsilon^2}\sum_{j=1}^n E(X_j^2).

Lemme de Kronecker : soient {xn}\{x_n\} et {bn}\{b_n\} deux suites réelles avec 0<bn0\lt b_n\uparrow\infty. Alors

k1xkbk<  limn1bnk=1nxk=0.\sum_{k\geq 1}\frac{x_k}{b_k}\lt\infty\ \Longrightarrow\ \lim_{n\to\infty}\frac{1}{b_n}\sum_{k=1}^n x_k=0.

الحل

1.

Comme ff est impaire, f(u)=f(u)f(-u)=-f(u) pour tout uu. Pour tout borélien BB, en utilisant que UU et U-U ont même loi :

P(VB)=P(f(U)B)=P(Uf1(B))=P(Uf1(B))=P(f(U)B)=P(f(U)B)=P(VB).P(V\in B)=P\big(f(U)\in B\big)=P\big(U\in f^{-1}(B)\big)=P\big(-U\in f^{-1}(B)\big)=P\big(f(-U)\in B\big)=P\big(-f(U)\in B\big)=P(-V\in B).

Donc VV et V-V ont même loi :

V=f(U) est une v.a.r. symeˊtrique.\boxed{V=f(U)\ \text{est une v.a.r. symétrique.}}

2.

Comme k1P(XkXk)<\sum_{k\geq 1}P(X_k\neq \overline{X}_k)\lt\infty, le premier lemme de Borel-Cantelli donne P(XkXk i.s.)=0P(X_k\neq\overline{X}_k\ \text{i.s.})=0. Donc, presque sûrement, il existe un rang (aléatoire) N=N(ω)N=N(\omega) tel que Xk=XkX_k=\overline{X}_k pour tout k>Nk\gt N. Pour n>Nn\gt N,

k=1nXkSn=k=1N(XkXk)=:R(ω),\sum_{k=1}^n \overline{X}_k-S_n=\sum_{k=1}^{N}\big(\overline{X}_k-X_k\big)=:R(\omega),

quantité finie et indépendante de nn. Par conséquent

1ank=1nXk=R(ω)an+Snann0+0=0p.s.,\frac{1}{a_n}\sum_{k=1}^n \overline{X}_k=\frac{R(\omega)}{a_n}+\frac{S_n}{a_n}\xrightarrow[n\to\infty]{}0+0=0\quad\text{p.s.},

car R(ω)R(\omega) est fixé, an+a_n\to+\infty et Sn/an0S_n/a_n\to 0 p.s. par hypothèse.

limn1ank=1nXk=0 p.s.\boxed{\lim_{n\to\infty}\frac{1}{a_n}\sum_{k=1}^n \overline{X}_k=0\ \text{p.s.}}

3.

Indépendance et espérance nulle entraînent l'orthogonalité dans L2L^2 : pour m<nm\lt n,

E[(SnSm)2]=E[(k=m+1nXk)2]=k=m+1nE(Xk2),E\big[(S_n-S_m)^2\big]=E\Big[\Big(\sum_{k=m+1}^n X_k\Big)^2\Big]=\sum_{k=m+1}^n E(X_k^2),

les termes croisés étant nuls (E[XiXj]=E[Xi]E[Xj]=0E[X_iX_j]=E[X_i]E[X_j]=0 pour iji\neq j). Comme kE(Xk2)<\sum_k E(X_k^2)\lt\infty, le reste tend vers 00, donc {Sn}\{S_n\} est de Cauchy dans L2L^2 ; cet espace étant complet, SnS_n converge dans L2L^2.

Pour la convergence presque sûre, on applique l'inégalité de Kolmogorov aux variables {Xj}j>m\{X_j\}_{j\gt m} : pour tout ε>0\varepsilon\gt 0,

P(maxm<knSkSm>ε)1ε2j=m+1nE(Xj2).P\Big(\max_{m\lt k\leq n}|S_k-S_m|\gt\varepsilon\Big)\leq \frac{1}{\varepsilon^2}\sum_{j=m+1}^n E(X_j^2).

En faisant nn\to\infty,

P(supk>mSkSmε)1ε2j>mE(Xj2)=:rmε2,rm0.P\Big(\sup_{k\gt m}|S_k-S_m|\geq\varepsilon\Big)\leq \frac{1}{\varepsilon^2}\sum_{j\gt m} E(X_j^2)=:\frac{r_m}{\varepsilon^2},\qquad r_m\to 0.

Posons Mm=supk,lmSkSlM_m=\sup_{k,l\geq m}|S_k-S_l|, suite décroissante en mm, de limite MM. Comme Mm2supkmSkSmM_m\leq 2\sup_{k\geq m}|S_k-S_m|, on a P(Mm>2ε)rm/ε20P(M_m\gt 2\varepsilon)\leq r_m/\varepsilon^2\to 0, d'où M=0M=0 p.s. : la suite {Sn}\{S_n\} est p.s. de Cauchy, donc converge p.s.

n1E(Xn2)<  Sn converge dans L2 et p.s.\boxed{\sum_{n\geq 1}E(X_n^2)\lt\infty\ \Longrightarrow\ S_n\ \text{converge dans }L^2\ \text{et p.s.}}

4.

On a Xk=fk(Xk)\overline{X}_k=f_k(X_k) avec fk(x)=x1(xak)f_k(x)=x\,\mathbf{1}(|x|\leq a_k), fonction impaire. Les XkX_k étant symétriques (question 1), Xk\overline{X}_k est symétrique ; bornée (Xkak|\overline{X}_k|\leq a_k) donc intégrable, une variable symétrique intégrable vérifie E(Xk)=0E(\overline{X}_k)=0. De plus les Xk\overline{X}_k sont indépendantes et dans L2L^2. Posons Zk=XkakZ_k=\dfrac{\overline{X}_k}{a_k} ; alors E(Zk)=0E(Z_k)=0 et, par la question 3, kZk\sum_k Z_k converge p.s. dès que kE(Zk2)<\sum_k E(Z_k^2)\lt\infty.

Or E(Zk2)=1ak2E[X121(X1ak)]E(Z_k^2)=\dfrac{1}{a_k^2}E\big[X_1^2\mathbf{1}(|X_1|\leq a_k)\big] (les XkX_k sont identiquement distribuées). Par le théorème de Tonelli,

k1E(Zk2)=E[X12k:akX11ak2].\sum_{k\geq 1}E(Z_k^2)=E\Big[X_1^2\sum_{k:\,a_k\geq |X_1|}\frac{1}{a_k^2}\Big].

Pour x=X1x=|X_1|, notons ν(x)=min{k1: akx}\nu(x)=\min\{k\geq 1:\ a_k\geq x\}. Alors k:akxak2=kν(x)ak2\sum_{k:\,a_k\geq x}a_k^{-2}=\sum_{k\geq \nu(x)}a_k^{-2}. En séparant xan01x\leq a_{n_0-1} (où X12X_1^2 est bornée et kak2<\sum_k a_k^{-2}\lt\infty) et xx tel que ν(x)n0\nu(x)\geq n_0, la condition donnée fournit

kν(x)1ak2Cν(x)aν(x)2Cν(x)x2(car aν(x)x).\sum_{k\geq \nu(x)}\frac{1}{a_k^2}\leq \frac{C\,\nu(x)}{a_{\nu(x)}^2}\leq \frac{C\,\nu(x)}{x^2}\qquad(\text{car }a_{\nu(x)}\geq x).

Donc, à une constante additive finie près,

k1E(Zk2)cste+CE[ν(X1)].\sum_{k\geq 1}E(Z_k^2)\leq \text{cste}+C\,E\big[\nu(|X_1|)\big].

Comme ν(x)k+1    ak<x\nu(x)\geq k+1\iff a_k\lt x,

E[ν(X1)]=k1P(ν(X1)k)=1+k1P(X1>ak)<.E\big[\nu(|X_1|)\big]=\sum_{k\geq 1}P\big(\nu(|X_1|)\geq k\big)=1+\sum_{k\geq 1}P(|X_1|\gt a_k)\lt\infty.

Ainsi kE(Zk2)<\sum_k E(Z_k^2)\lt\infty et, par la question 3,

k=1nXkak converge p.s.\boxed{\sum_{k=1}^n \frac{\overline{X}_k}{a_k}\ \text{converge p.s.}}

5.

On prend an=nαa_n=n^\alpha avec 12<α<1\tfrac12\lt\alpha\lt 1. Comme 2α>12\alpha\gt 1, knk2αn12α2α1\sum_{k\geq n}k^{-2\alpha}\sim \dfrac{n^{1-2\alpha}}{2\alpha-1}, donc il existe C>0C\gt 0 avec knak2Cn12α=Cn/an2\sum_{k\geq n}a_k^{-2}\leq C\,n^{1-2\alpha}=C\,n/a_n^2 : la condition de la troisième partie est satisfaite.

()(\Longleftarrow) Supposons kP(X1>kα)<\sum_k P(|X_1|\gt k^\alpha)\lt\infty. Les XkX_k étant i.i.d., kP(Xk>ak)=kP(X1>kα)<\sum_k P(|X_k|\gt a_k)=\sum_k P(|X_1|\gt k^\alpha)\lt\infty, donc kP(XkXk)<\sum_k P(X_k\neq \overline{X}_k)\lt\infty et, par Borel-Cantelli, Xk=XkX_k=\overline{X}_k à partir d'un certain rang p.s. Par la question 4, kXk/ak\sum_k \overline{X}_k/a_k converge p.s. ; le lemme de Kronecker (avec bn=an=nαb_n=a_n=n^\alpha\uparrow\infty) donne

1nαk=1nXk0 p.s.\frac{1}{n^\alpha}\sum_{k=1}^n \overline{X}_k\to 0\ \text{p.s.}

Puisque Xk=XkX_k=\overline{X}_k à partir d'un certain rang, Snk=1nXkS_n-\sum_{k=1}^n\overline{X}_k est finalement constant, donc (SnkXk)/nα0\big(S_n-\sum_k\overline{X}_k\big)/n^\alpha\to 0, et finalement Sn/nα0S_n/n^\alpha\to 0 p.s.

()(\Longrightarrow) Supposons Sn/nα0S_n/n^\alpha\to 0 p.s. Alors

Xnnα=Snnα(n1n)αSn1(n1)α0 p.s.,\frac{X_n}{n^\alpha}=\frac{S_n}{n^\alpha}-\Big(\frac{n-1}{n}\Big)^{\alpha}\frac{S_{n-1}}{(n-1)^\alpha}\to 0\ \text{p.s.},

donc P(Xn>nα i.s.)=0P(|X_n|\gt n^\alpha\ \text{i.s.})=0. Les événements {Xn>nα}\{|X_n|\gt n^\alpha\} étant indépendants, le second lemme de Borel-Cantelli impose nP(Xn>nα)<\sum_n P(|X_n|\gt n^\alpha)\lt\infty, soit, par équidistribution, nP(X1>nα)<\sum_n P(|X_1|\gt n^\alpha)\lt\infty.

12<α<1:Snnα0 p.s.    k1P(X1>kα)<.\boxed{\tfrac12\lt\alpha\lt 1:\quad \frac{S_n}{n^\alpha}\to 0\ \text{p.s.}\iff \sum_{k\geq 1}P(|X_1|\gt k^\alpha)\lt\infty.}

6.

Pour α=12\alpha=\tfrac12, la condition de la troisième partie tombe en défaut : knak2=kn1k=+\sum_{k\geq n}a_k^{-2}=\sum_{k\geq n}\tfrac1k=+\infty (série harmonique), donc la méthode précédente ne s'applique plus. On construit un contre-exemple : soit (Xn)(X_n) i.i.d. de Rademacher, P(X1=1)=P(X1=1)=12P(X_1=1)=P(X_1=-1)=\tfrac12 (symétriques, E(X1)=0E(X_1)=0, Var(X1)=1\operatorname{Var}(X_1)=1). Alors

k1P(X1>k)=0<(car X1=1k k1),\sum_{k\geq 1}P(|X_1|\gt \sqrt{k})=0\lt\infty\quad(\text{car }|X_1|=1\leq\sqrt{k}\ \forall k\geq 1),

donc le membre de droite est satisfait. Pourtant, par le théorème central limite Snn L N(0,1)\dfrac{S_n}{\sqrt{n}}\xrightarrow{\ \mathcal{L}\ }\mathcal{N}(0,1) (loi non dégénérée), et par la loi du logarithme itéré lim supnSnn=+\limsup_n \dfrac{|S_n|}{\sqrt{n}}=+\infty p.s., de sorte que Sn/n↛0S_n/\sqrt{n}\not\to 0 p.s. Le membre de gauche est donc faux.

Pour α=12, le membre de droite peut eˆtre vrai sans que Sn/n0: l’eˊquivalence est fausse.\boxed{\text{Pour }\alpha=\tfrac12,\ \text{le membre de droite peut être vrai sans que }S_n/\sqrt{n}\to 0:\ \text{l'équivalence est fausse.}}

التمرين 2

Problème 2 — Chaîne de Markov à 4 états : irréductibilité, loi stationnaire, temps de retour et période

#markov-chain#stationary-distribution#recurrence#hitting-time#periodicity

Soit X=(Xn,nN)X=(X_n,\,n\in\mathbb{N}) une chaîne de Markov à valeurs dans F={1,2,3,4}F=\{1,2,3,4\} de matrice de transition

P=(150252510003515150045150).P=\begin{pmatrix}\tfrac15 & 0 & \tfrac25 & \tfrac25\\ 1 & 0 & 0 & 0\\ \tfrac35 & \tfrac15 & \tfrac15 & 0\\ 0 & \tfrac45 & \tfrac15 & 0\end{pmatrix}.

  1. La chaîne est-elle irréductible ?
  2. Déterminer les états récurrents et les états transients.
  3. Déterminer les lois stationnaires.
  4. Pour xFx\in F on note Tx=inf{n1: Xn=x}T_x=\inf\{n\geq 1:\ X_n=x\} le temps de retour à xx. Calculer Ex[Tx]\mathbb{E}_x[T_x].
  5. Calculer la période de chaque état.
الحل

1.

Les transitions possibles sont : 1{1,3,4}1\to\{1,3,4\}, 2{1}2\to\{1\}, 3{1,2,3}3\to\{1,2,3\}, 4{2,3}4\to\{2,3\}. On vérifie que tous les états communiquent :

1321,142,21,31.1\to 3\to 2\to 1,\qquad 1\to 4\to 2,\qquad 2\to 1,\qquad 3\to 1.

Depuis n'importe quel état on atteint 11, et depuis 11 on atteint 33, 44 (puis 22). Tous les états forment une seule classe de communication :

la chaıˆne est irreˊductible.\boxed{\text{la chaîne est irréductible.}}

2.

L'espace d'états FF est fini et la chaîne est irréductible : tous les états appartiennent à une unique classe finie et fermée, donc sont récurrents (récurrents positifs).

Les eˊtats 1,2,3,4 sont tous reˊcurrents positifs ; aucun eˊtat n’est transient.\boxed{\text{Les états }1,2,3,4\ \text{sont tous récurrents positifs ; aucun état n'est transient.}}

3.

On cherche π=(π1,π2,π3,π4)\pi=(\pi_1,\pi_2,\pi_3,\pi_4) vérifiant πP=π\pi P=\pi et iπi=1\sum_i\pi_i=1. En colonnes :

π1=15π1+π2+35π3,π2=15π3+45π4,π3=25π1+15π3+15π4,π4=25π1.\pi_1=\tfrac15\pi_1+\pi_2+\tfrac35\pi_3,\quad \pi_2=\tfrac15\pi_3+\tfrac45\pi_4,\quad \pi_3=\tfrac25\pi_1+\tfrac15\pi_3+\tfrac15\pi_4,\quad \pi_4=\tfrac25\pi_1.

De la dernière équation π4=25π1\pi_4=\tfrac25\pi_1. La troisième donne 45π3=25π1+15π4=1225π1\tfrac45\pi_3=\tfrac25\pi_1+\tfrac15\pi_4=\tfrac{12}{25}\pi_1, soit π3=35π1\pi_3=\tfrac35\pi_1. La deuxième donne π2=15π3+45π4=325π1+825π1=1125π1\pi_2=\tfrac15\pi_3+\tfrac45\pi_4=\tfrac{3}{25}\pi_1+\tfrac{8}{25}\pi_1=\tfrac{11}{25}\pi_1. La normalisation

π1(1+1125+35+25)=π16125=1  π1=2561.\pi_1\Big(1+\tfrac{11}{25}+\tfrac35+\tfrac25\Big)=\pi_1\cdot\frac{61}{25}=1\ \Longrightarrow\ \pi_1=\frac{25}{61}.

D'où l'unique loi stationnaire (la chaîne est irréductible) :

π=(2561, 1161, 1561, 1061).\boxed{\pi=\Big(\tfrac{25}{61},\ \tfrac{11}{61},\ \tfrac{15}{61},\ \tfrac{10}{61}\Big).}

4.

La chaîne étant irréductible récurrente positive, le temps moyen de retour à un état est l'inverse de sa probabilité stationnaire : Ex[Tx]=1πx\mathbb{E}_x[T_x]=\dfrac{1}{\pi_x}. Donc

E1[T1]=6125,E2[T2]=6111,E3[T3]=6115,E4[T4]=6110.\boxed{\mathbb{E}_1[T_1]=\frac{61}{25},\quad \mathbb{E}_2[T_2]=\frac{61}{11},\quad \mathbb{E}_3[T_3]=\frac{61}{15},\quad \mathbb{E}_4[T_4]=\frac{61}{10}.}

5.

Come P11=15>0P_{11}=\tfrac15\gt 0, l'état 11 admet un retour en une étape, donc sa période vaut gcd{n1: P11n>0}=1\gcd\{n\geq 1:\ P^n_{11}\gt 0\}=1. La chaîne étant irréductible, tous les états ont la même période :

chaque eˊtat est de peˊriode 1 : la chaıˆne est apeˊriodique (donc ergodique).\boxed{\text{chaque état est de période }1\ :\ \text{la chaîne est apériodique (donc ergodique).}}