📚 الرئيسية

مسابقة دكتوراه 2016Université des Sciences et de la Technologie Houari Boumediène (USTHB) — الموضوع 04

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

Concours Doctorat LMD en Probabilités, Statistiques et Applications, épreuve de Probabilités (Sujet 2), Faculté des Mathématiques, Département de Probabilités et Statistique, Université des Sciences et de la Technologie Houari Boumediène (USTHB), année universitaire 2016/2017, 3 novembre 2016, durée 2 heures.

التمرين 1

Exercice 1 — Vrai/Faux et mélange Poisson-exponentiel

#markov-chain#brownian-motion#ar1-process#conditional-expectation#poisson-mixture

A/ Les affirmations suivantes sont-elles vraies ou fausses ? Justifiez vos réponses et donnez l'affirmation correcte correspondante dans le cas où l'affirmation proposée est fausse.

  1. Une chaîne de Markov homogène irréductible possède une distribution stationnaire unique.
  2. Si {Wt}t0\{W_t\}_{t\ge0} est un mouvement brownien standard, alors sa covariance est Cov(Wt,Ws)=(ts)2\mathrm{Cov}(W_t,W_s) = (t-s)^2.
  3. L'espérance d'un processus autorégressif AR(1)AR(1) : Xt=ϕ0+ϕ1Xt1+εtX_t = \phi_0 + \phi_1 X_{t-1} + \varepsilon_t, avec ϕ1<1|\phi_1|<1 et εt\varepsilon_t un bruit blanc centré, est égale à ϕ0+ϕ1\phi_0 + \phi_1.
  4. Le temps de séjour d'une chaîne de Markov à temps continu dans un état suit une distribution exponentielle.

B/ Soit (Ω,F,P)(\Omega,\mathcal{F},P) un espace probabilisé, A\mathcal{A} une sous-σ\sigma-algèbre de F\mathcal{F} et XX une variable aléatoire intégrable (E(X)<E(|X|)<\infty) sur (Ω,F,P)(\Omega,\mathcal{F},P).

  1. Définir l'espérance conditionnelle sachant une σ\sigma-algèbre E(XA)E(X\mid\mathcal{A}).
  2. Supposons que XX est distribuée selon une loi de Poisson de paramètre λ\lambda. Le paramètre λ\lambda est lui-même une v.a. de distribution exponentielle de moyenne 1/μ1/\mu. Trouver la distribution de XX.
الحل

A.1.

Fausse. L'irréductibilité seule ne suffit pas : une chaîne irréductible mais récurrente nulle ou transiente (par exemple la marche aléatoire simple symétrique sur Z\mathbb{Z}) ne possède aucune distribution stationnaire.

Correction : une chaîne de Markov homogène irréductible et récurrente positive possède une unique distribution stationnaire (en particulier, toute chaîne finie irréductible en possède une unique).

A.2.

Fausse. Pour sts\le t, on écrit Wt=Ws+(WtWs)W_t = W_s + (W_t-W_s) avec WtWsW_t-W_s indépendant de WsW_s, donc Cov(Wt,Ws)=Var(Ws)=s\mathrm{Cov}(W_t,W_s)=\mathrm{Var}(W_s)=s. En général :

Cov(Wt,Ws)=min(t,s)\boxed{\mathrm{Cov}(W_t,W_s)=\min(t,s)}

A.3.

Fausse. En prenant l'espérance de Xt=ϕ0+ϕ1Xt1+εtX_t=\phi_0+\phi_1 X_{t-1}+\varepsilon_t, avec la stationnarité (E[Xt]=E[Xt1]=mE[X_t]=E[X_{t-1}]=m) et E[εt]=0E[\varepsilon_t]=0, on obtient m=ϕ0+ϕ1mm=\phi_0+\phi_1 m, d'où

E[Xt]=ϕ01ϕ1\boxed{E[X_t]=\frac{\phi_0}{1-\phi_1}}

A.4.

Vraie. Par la propriété de Markov (absence de mémoire), le temps de séjour dans un état ii d'une chaîne de Markov à temps continu est sans mémoire, donc de loi exponentielle, de paramètre qiq_i (le taux de sortie de l'état ii).

B.1.

E(XA)E(X\mid\mathcal A) est l'unique (à une égalité presque sûre près) variable aléatoire YY telle que :

  • YY est A\mathcal A-mesurable et intégrable ;
  • pour tout AAA\in\mathcal A, AYdP=AXdP\displaystyle\int_A Y\,dP=\int_A X\,dP, c'est-à-dire E[Y1A]=E[X1A]E[Y\,\mathbf{1}_A]=E[X\,\mathbf{1}_A].

Son existence et son unicité découlent du théorème de Radon-Nikodym.

B.2.

On a XλP(λ)X\mid\lambda\sim\mathcal P(\lambda) et λE(μ)\lambda\sim\mathcal E(\mu), de densité f(λ)=μeμλf(\lambda)=\mu e^{-\mu\lambda} (λ>0\lambda>0), de moyenne 1/μ1/\mu. Par la formule des probabilités totales, pour kNk\in\mathbb N :

P(X=k)=0eλλkk!μeμλdλ=μk!0λke(1+μ)λdλ=μk!k!(1+μ)k+1P(X=k)=\int_0^{\infty} e^{-\lambda}\frac{\lambda^k}{k!}\,\mu e^{-\mu\lambda}\,d\lambda=\frac{\mu}{k!}\int_0^{\infty}\lambda^k e^{-(1+\mu)\lambda}\,d\lambda=\frac{\mu}{k!}\cdot\frac{k!}{(1+\mu)^{k+1}}

P(X=k)=μ(1+μ)k+1=μ1+μ(11+μ)kP(X=k)=\frac{\mu}{(1+\mu)^{k+1}}=\frac{\mu}{1+\mu}\left(\frac{1}{1+\mu}\right)^{k}

C'est une loi géométrique sur {0,1,2,}\{0,1,2,\dots\} de paramètre p=μ1+μp=\dfrac{\mu}{1+\mu} :

P(X=k)=p(1p)k,p=μ1+μ,k=0,1,2,\boxed{P(X=k)=p(1-p)^k,\quad p=\frac{\mu}{1+\mu},\quad k=0,1,2,\dots}

التمرين 2

Exercice 2 — Chaîne de Markov : classes, absorption et lois limites

#markov-chain#absorption-probability#stationary-distribution#recurrent-classes#limiting-distribution

Soit une chaîne de Markov homogène d'espace d'états S={1,2,3,4,5}S=\{1,2,3,4,5\} dont la matrice de transition en une étape PP est :

P=(13α000β34000000100013230θ00λ0)P=\begin{pmatrix} \tfrac13 & \alpha & 0 & 0 & 0 \\ \beta & \tfrac34 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 & 0 \\ 0 & 0 & \tfrac13 & \tfrac23 & 0 \\ \theta & 0 & 0 & \lambda & 0 \end{pmatrix}

  1. Pour quelles valeurs de α\alpha, β\beta, θ\theta et λ\lambda la matrice ci-dessus est-elle cohérente avec celle d'une chaîne de Markov, sachant que θ=2β\theta=2\beta ?
  2. Étudier la nature de la chaîne.
  3. Calculer les probabilités d'absorption dans les classes récurrentes pour les différents états transitoires.
  4. Déterminer les distributions stationnaires concentrées dans les classes récurrentes.
  5. Calculer les probabilités limites limnpij(n)\lim_{n\to\infty} p_{ij}^{(n)}, i,jSi,j\in S.
الحل

1.

Chaque ligne de PP doit être formée de coefficients positifs sommant à 11 :

  • Ligne 1 : 13+α=1α=23\tfrac13+\alpha=1\Rightarrow \alpha=\tfrac23.
  • Ligne 2 : β+34=1β=14\beta+\tfrac34=1\Rightarrow \beta=\tfrac14.
  • Lignes 3 et 4 : déjà sommées à 11.
  • Ligne 5 : θ+λ=1\theta+\lambda=1 avec θ=2β=12\theta=2\beta=\tfrac12, d'où λ=12\lambda=\tfrac12.

α=23,β=14,θ=12,λ=12\boxed{\alpha=\tfrac23,\quad \beta=\tfrac14,\quad \theta=\tfrac12,\quad \lambda=\tfrac12}

La matrice devient

P=(1/32/30001/43/400000010001/32/301/2001/20)P=\begin{pmatrix} 1/3 & 2/3 & 0 & 0 & 0 \\ 1/4 & 3/4 & 0 & 0 & 0 \\ 0 & 0 & 0 & 1 & 0 \\ 0 & 0 & 1/3 & 2/3 & 0 \\ 1/2 & 0 & 0 & 1/2 & 0 \end{pmatrix}

2.

  • {1,2}\{1,2\} est une classe fermée (aucune sortie), irréductible et apériodique (car P11=13>0P_{11}=\tfrac13>0) : classe récurrente.
  • {3,4}\{3,4\} est une classe fermée, irréductible et apériodique (car P44=23>0P_{44}=\tfrac23>0) : classe récurrente.
  • {5}\{5\} conduit vers 11 et vers 44, mais aucun état n'y revient : état transitoire.

La chaîne est donc réductible, avec deux classes récurrentes C1={1,2}C_1=\{1,2\} et C2={3,4}C_2=\{3,4\}, et un état transitoire {5}\{5\}.

3.

Depuis l'état transitoire 55, un pas mène à 1C11\in C_1 avec probabilité 12\tfrac12 et à 4C24\in C_2 avec probabilité 12\tfrac12 ; une fois dans une classe récurrente fermée, la chaîne y reste. D'où les probabilités d'absorption :

P(absorption dans C1X0=5)=12,P(absorption dans C2X0=5)=12\boxed{P(\text{absorption dans }C_1\mid X_0=5)=\tfrac12,\qquad P(\text{absorption dans }C_2\mid X_0=5)=\tfrac12}

(Pour iC1i\in C_1 : absorption dans C1C_1 avec probabilité 11 ; pour iC2i\in C_2 : absorption dans C2C_2 avec probabilité 11.)

4.

Classe C1={1,2}C_1=\{1,2\} : on résout πP=π\pi P=\pi avec π1+π2=1\pi_1+\pi_2=1 :

π1=13π1+14π2  23π1=14π2  π2=83π1\pi_1=\tfrac13\pi_1+\tfrac14\pi_2\ \Rightarrow\ \tfrac23\pi_1=\tfrac14\pi_2\ \Rightarrow\ \pi_2=\tfrac83\pi_1

d'où π(1)=(311,811,0,0,0)\pi^{(1)}=\left(\tfrac{3}{11},\tfrac{8}{11},0,0,0\right).

Classe C2={3,4}C_2=\{3,4\} :

π3=13π4,π4=π3+23π4  π3=14, π4=34\pi_3=\tfrac13\pi_4,\qquad \pi_4=\pi_3+\tfrac23\pi_4\ \Rightarrow\ \pi_3=\tfrac14,\ \pi_4=\tfrac34

d'où π(2)=(0,0,14,34,0)\pi^{(2)}=\left(0,0,\tfrac14,\tfrac34,0\right).

π(1)=(311,811,0,0,0),π(2)=(0,0,14,34,0)\boxed{\pi^{(1)}=\left(\tfrac{3}{11},\tfrac{8}{11},0,0,0\right),\qquad \pi^{(2)}=\left(0,0,\tfrac14,\tfrac34,0\right)}

Toute distribution stationnaire est une combinaison convexe cπ(1)+(1c)π(2)c\,\pi^{(1)}+(1-c)\,\pi^{(2)}, c[0,1]c\in[0,1] ; l'état 55 a une probabilité stationnaire nulle.

5.

Les deux classes récurrentes étant apériodiques, les limites existent.

  • Pour iC1i\in C_1 : limnpi1(n)=311\lim_n p_{i1}^{(n)}=\tfrac{3}{11}, limnpi2(n)=811\lim_n p_{i2}^{(n)}=\tfrac{8}{11}, les autres nulles.
  • Pour iC2i\in C_2 : limnpi3(n)=14\lim_n p_{i3}^{(n)}=\tfrac14, limnpi4(n)=34\lim_n p_{i4}^{(n)}=\tfrac34, les autres nulles.
  • Pour i=5i=5 (pondération par les probabilités d'absorption 12,12\tfrac12,\tfrac12) :

limnp51(n)=12311=322,limnp52(n)=12811=411\lim_n p_{51}^{(n)}=\tfrac12\cdot\tfrac{3}{11}=\tfrac{3}{22},\qquad \lim_n p_{52}^{(n)}=\tfrac12\cdot\tfrac{8}{11}=\tfrac{4}{11}

limnp53(n)=1214=18,limnp54(n)=1234=38\lim_n p_{53}^{(n)}=\tfrac12\cdot\tfrac14=\tfrac18,\qquad \lim_n p_{54}^{(n)}=\tfrac12\cdot\tfrac34=\tfrac38

  • Pour tout ii : limnpi5(n)=0\lim_n p_{i5}^{(n)}=0 (l'état 55 est transitoire).

limnp5(n)=(322, 411, 18, 38, 0)\boxed{\lim_n p_{5\cdot}^{(n)}=\left(\tfrac{3}{22},\ \tfrac{4}{11},\ \tfrac18,\ \tfrac38,\ 0\right)}