ATLAS Ouvrir dans l'atlas →

Théorème 1Chaînes de Markov

Construction d'une chaîne à partir de variables uniformes ⋆

⋆ Axiome du choix

Le point (iii) hérite du choix dénombrable dont dépend la mesure de Lebesgue, par la suite de variables uniformes indépendantes du nœud Suites indépendantes ; c'est le seul usage du choix dans ce chapitre, signalé ici une fois. Les points (i) et (ii) n'en font aucun : l'énumération de II est fixée une fois, et les constructions sont explicites. Les fiches suivantes énoncent des propriétés de chaînes données ; plusieurs, comme la récurrence, la mesure γx\gamma_x ou la formule π(x)=1/Ex[Tx]\pi(x) = 1/E_x[T_x], ne portent que sur QQ, mais s'écrivent avec PxP_x, qui suppose une chaîne issue de xx. Une chaîne irréductible donnée en fournit une pour chaque état sans aucun choix (Propriété de Markov, (ii)) ; à défaut, (iii) en donne une, et ces énoncés héritent alors de ce ⋆. Le couplage du nœud Convergence vers l'équilibre, points (i) et (ii), reste sans choix : il se fait sur un ensemble fini, par (ii).

Soient μ\mu une probabilité et QQ une matrice de transition sur II, dont les états sont énumérés y0,y1,…y_0, y_1, \dots Pour x∈Ix \in I et uu réel, on pose f(x,u)=ykf(x, u) = y_k, où kk est le plus petit indice tel que u⩽Q(x,y0)+⋯+Q(x,yk),u \leqslant Q(x, y_0) + \cdots + Q(x, y_k), si 0<u⩽10 < u \leqslant 1 et si un tel indice existe, et f(x,u)=y0f(x, u) = y_0 sinon. On définit de même g(u)g(u) à partir des μ(yk)\mu(y_k).

Énoncé

(i) Soient, sur un espace probabilisé quelconque, U0,U1,…U_0, U_1, \dots des variables indépendantes de loi uniforme sur [0,1][0, 1]. La suite définie par X0=g(U0),Xn+1=f(Xn,Un+1)X_0 = g(U_0), \qquad X_{n+1} = f(X_n, U_{n+1}) est une chaîne (μ,Q)(\mu, Q). Pour tout x∈Ix \in I, la suite X0x=xX^x_0 = x, Xn+1x=f(Xnx,Un+1)X^x_{n+1} = f(X^x_n, U_{n+1}) est une chaîne issue de xx.

(ii) Horizon fini. Pour m∈ℕm \in \mathbb N, la formule p(x0,…,xm)=μ(x0)Q(x0,x1)⋯Q(xm−1,xm)\begin{aligned} &p(x_0, \dots, x_m) \\ &= \mu(x_0)Q(x_0, x_1)\cdots Q(x_{m-1}, x_m) \end{aligned} définit une probabilité sur l'ensemble fini ou dénombrable Im+1I^{m+1}, muni de toutes ses parties, sous laquelle les coordonnées forment une chaîne (μ,Q)(\mu, Q) jusqu'à l'instant m.

(iii) Il existe une chaîne (μ,Q)(\mu, Q) sur [0,1[[0, 1[ muni de la mesure de Lebesgue.

Démonstration

(i) La loi d'un pas. Fixons xx, posons s−1=0s_{-1} = 0 et sk=Q(x,y0)+⋯+Q(x,yk)s_k = Q(x, y_0) + \cdots + Q(x, y_k). Pour k⩾1k \geqslant 1, on a f(x,u)=ykf(x, u) = y_k si et seulement si 0<u⩽10 < u \leqslant 1 et sk−1<u⩽sks_{k-1} < u \leqslant s_k.

Si 0<u<10 < u < 1, un indice kk tel que u⩽sku \leqslant s_k existe, puisque les sks_k croissent vers 11 ; l'ensemble des uu où f(x,u)=y0f(x, u) = y_0 est donc ]0,s0]]0, s_0], augmenté d'une partie de ℝ∖ ]0,1[{\mathbb R} \setminus \,]0, 1[.

Si UU est uniforme sur [0,1][0, 1], P(a<U⩽b)=b−aP(a < U \leqslant b) = b - a pour 0⩽a⩽b⩽10 \leqslant a \leqslant b \leqslant 1 et P(U∉ ]0,1[)=0P(U \notin \,]0, 1[) = 0 ; donc f(x,U)f(x, U) est une variable à valeurs dans I, et P(f(x,U)=yk)=sk−sk−1=Q(x,yk)P(f(x, U) = y_k) = s_k - s_{k-1} = Q(x, y_k). De même, g(U)g(U) a pour loi μ\mu.

La récurrence. Par récurrence sur nn, XnX_n est une variable à valeurs dans II, mesurable pour σ(U0,…,Un)\sigma(U_0, \dots, U_n), puisque {Xn+1=y}=⋃x∈I({Xn=x}∩{f(x,Un+1)=y}).\begin{aligned} &\{X_{n+1} = y\} \\ &= \bigcup_{x \in I}\big(\{X_n = x\} \cap \{f(x, U_{n+1}) = y\}\big). \end{aligned}

L'événement A={X0=x0,…,Xn=xn}A = \{X_0 = x_0, \dots, X_n = x_n\} appartient donc à σ(U0,…,Un)\sigma(U_0, \dots, U_n), qui est indépendante de Un+1U_{n+1} par regroupement (Indépendance, (ii)). Ainsi P(A∩{Xn+1=xn+1})=P(A∩{f(xn,Un+1)=xn+1})=P(A) Q(xn,xn+1),\begin{aligned} &P(A \cap \{X_{n+1} = x_{n+1}\}) \\ &= P\big(A \cap \{f(x_n, U_{n+1}) = x_{n+1}\}\big) \\ &= P(A)\,Q(x_n, x_{n+1}), \end{aligned} et la formule de définition s'obtient par récurrence, à partir de P(X0=x0)=μ(x0)P(X_0 = x_0) = \mu(x_0). Le cas de XxX^x est identique, avec X0x=xX^x_0 = x.

(ii) Les poids sont positifs. En sommant d'abord sur xmx_m, puis sur xm−1x_{m-1}, et ainsi de suite, par paquets (Sommes dans [0, +∞], (iv)), chaque somme sur une ligne de QQ vaut 11, puis ∑μ=1\sum \mu = 1 : la masse totale est 11, et A↦∑w∈Ap(w)A \mapsto \sum_{w \in A}p(w) est une probabilité, σ-additive par sommation par paquets.

Pour n⩽mn \leqslant m, la même sommation sur xn+1,…,xmx_{n+1}, \dots, x_m seulement donne P(X0=x0,…,Xn=xn)=μ(x0)⋯Q(xn−1,xn)P(X_0 = x_0, \dots, X_n = x_n) = \mu(x_0)\cdots Q(x_{n-1}, x_n).

(iii) Le nœud Suites indépendantes, (i), fournit sur [0,1[[0, 1[ une suite de variables indépendantes de loi uniforme ; on applique (i).

□
Remarque

La construction évite tout théorème d'extension de Kolmogorov ou d'Ionescu Tulcea, qui construirait directement une probabilité sur (Iℕ,𝒞)(I^{\mathbb N}, \mathcal C) : la chaîne est une fonction explicite d'une suite de variables uniformes. C'est aussi la façon dont on simule une chaîne, chaque pas tirant un nombre au hasard dans [0,1][0, 1] et cherchant dans quel intervalle de longueurs Q(x,y0),Q(x,y1),…Q(x, y_0), Q(x, y_1), \dots il tombe. Pour des espaces d'états non dénombrables, ces théorèmes deviennent nécessaires ; ils sont admis ici.

Remarque

Les chaînes XxX^x et XyX^y de (i), issues de deux états, sont construites sur les mêmes variables UnU_n. L'erreur serait de les croire indépendantes : dès qu'elles se rencontrent, elles restent égales, puisque chaque pas applique la même fonction ff au même Un+1U_{n+1}. Faire marcher deux chaînes ensemble jusqu'à ce qu'elles se rejoignent est l'idée du couplage, que le nœud Convergence vers l'équilibre mettra en œuvre avec deux chaînes indépendantes jusqu'à leur rencontre.

Ouvrir dans l'atlas interactif →