16.07.2014 Views

martingales discrètes, chaines de Markov, processus de Poisson

martingales discrètes, chaines de Markov, processus de Poisson

martingales discrètes, chaines de Markov, processus de Poisson

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

1.7 Algorithme <strong>de</strong> Robbins-Monro<br />

Exemple : le dosage.<br />

Il s’agit <strong>de</strong> déterminer le dosage optimal d’un produit chimique pour obtenir l’effet voulu ;<br />

on sait, à cet effet, effectuer <strong>de</strong>s tests.<br />

Au dosage z ∈ R, on associe l’effet F (z, η), où η est une variable aléatoire dont la<br />

réalisation dépend du test effectué (pas du dosage). L’effet moyen du dosage z est donné<br />

par f(z) = E[F (z, η)] ; on suppose que f est strictement croissante et continue.<br />

Soit a ∈ R l’effet désiré. On cherche à déterminer l’unique z ∗ tel que f(z ∗ ) = a. Pour cela<br />

on construit une suite <strong>de</strong> variables aléatoires (Z n , n ∈ N) définie par récurrence :<br />

Z n+1 = Z n − γ n [F (Z n , η n+1 ) − a]<br />

où γ est une suite <strong>de</strong> réels positifs convergeant vers 0 et (η n ) est une suite <strong>de</strong> variables<br />

aléatoires indépendantes <strong>de</strong> même loi que η. On va donner <strong>de</strong>s hypothèses pour que cette<br />

suite (Z n , n ∈ N) converge vers z ∗ .<br />

Théorème 1.35 <strong>de</strong> ROBBINS-MONRO : Soit h ∈ C(R, R), (γ n ) une suite <strong>de</strong> réels positifs<br />

et ε n <strong>de</strong>s variables aléatoires <strong>de</strong> carré intégrable telles que si F n = σ(ε 1 , · · · , ε n ),<br />

E(ε n+1 /F n ) = 0. On pose<br />

Z n+1 = Z n + γ n [h(Z n ) + ε n+1 ], Z 0 donnée constante,<br />

et on suppose<br />

(i) ∑ ∑<br />

n γ n = +∞ ; n γn 2 < +∞,<br />

(ii) E(ε 2 n+1/F n ) ≤ c(1 + Zn),<br />

2<br />

(iii) |h(z)| 2 ≤ c(1+z 2 ) et il existe un unique z ∗ ∈ R : h(z ∗ ) = 0, (z−z ∗ )h(z) < 0, ∀z ≠ z ∗ .<br />

Alors, la suite (Z n ) converge presque sûrement vers z ∗ .<br />

Un exemple <strong>de</strong> h : z ↦→ a − f(z) dans l’exemple di dosage ci-<strong>de</strong>ssus.<br />

Preuve en exercice<br />

Première étape : trouver une majoration <strong>de</strong> la forme<br />

E[(Z n+1 − z ∗ ) 2 /F n ] ≤ (Z n − z ∗ ) 2 (1 + αγ 2 n) + γ 2 nβ.<br />

En effet, d’une part les hypothèses montrent par récurrence que Z n est <strong>de</strong> carré intégrable<br />

pour tout n, et d’autre part :<br />

(Z n+1 − z ∗ ) 2 = (Z n − z ∗ ) 2 + 2γ n (Z n − z ∗ )(h(Z n ) + ε n+1 ) + γ 2 n(h(Z n ) + ε n+1 ) 2<br />

que l’on conditionne par F n :<br />

(1) E[(Z n+1 − z ∗ ) 2 /F n ] = (Z n − z ∗ ) 2 + 2γ n (Z n − z ∗ )h(Z n ) + γ 2 n[h(Z n ) 2 + E[ε 2 n+1/F n ]].<br />

19

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!