Introduction à la commande stochastique v.0.9 - Jean-Pierre ...
Introduction à la commande stochastique v.0.9 - Jean-Pierre ...
Introduction à la commande stochastique v.0.9 - Jean-Pierre ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
124 8. PROBLÈMES<br />
En utilisant <strong>la</strong> propriété P0x k = 1six = k modulo E, Pk 0x = 0 sinon,<br />
on obtient<br />
[(n−x)/E]<br />
p n x = ∑<br />
C (kE+x)<br />
n<br />
(1 − λ) n−kE−x λ kE+x .<br />
k=0<br />
3. Les valeurs propres de <strong>la</strong> matrice M sont<br />
k2 jπ/E<br />
(1 − λ) + λe<br />
en utilisant les propriétés spectrales des matrices de permutation. La<br />
matrice M est diagonalisable et admet <strong>la</strong> valeur propre 1 avec <strong>la</strong> multiplicité<br />
1. Les autres valeurs propres ont des modules strictement<br />
inferieurs <strong>à</strong> 1. On voit, en se p<strong>la</strong>çant dans <strong>la</strong> base propre, que M n<br />
converge vers le projecteur spectral sur le vecteur propre associé <strong>à</strong><strong>la</strong><br />
valeur propre 1. Or ce vecteur propre est p = ( 1/E ··· 1/E ) .<br />
p n converge donc vers l’unique mesure de probabilité vérifiant q =<br />
qM c.a.d. p.<br />
4. En faisant x = l dans <strong>la</strong> formule donnant explicitement p n x<br />
le résultat<br />
de convergence de <strong>la</strong> question précédente donne précisément <strong>la</strong> formule<br />
cherchée.<br />
2.2.2. OPTIMISATION DE LA POLITIQUE DE MAINTENANCE.<br />
1. La <strong>commande</strong> u peut prendre deux valeurs u = 0 signifie pas d’entretien,<br />
u = 1 signifie un entretien. La matrice de transition de <strong>la</strong><br />
chaîne de Markov controlée X n est alors<br />
⎛<br />
M u =<br />
⎜<br />
⎝<br />
1 − λ + uλ (1 − u)λ 0 · ·<br />
uλ 1 − λ (1 − u)λ 0 ·<br />
· · · · ·<br />
uλ 0 · 1 − λ (1 − u)λ<br />
λ 0 · 0 1− λ<br />
⎞<br />
⎟ .<br />
⎠<br />
Notant dx u = k x si u = 1 et 0 sinon, on veut résoudre<br />
∞∑<br />
min E 1<br />
U (1 + µ) (c i=0<br />
i+1 X n<br />
+ d U n<br />
X n<br />
).<br />
2. Pour résoudre ce problème on introduit <strong>la</strong> fonction de <strong>la</strong> programmation<br />
dynamique<br />
{ }<br />
∞<br />
v x = min E ∑ 1<br />
[ ]<br />
c<br />
U (1 + µ)<br />
i=0<br />
i+1 Xn + d U n<br />
X n<br />
| X 0 = x .<br />
v vérifie l’équation de <strong>la</strong> programmation dynamique<br />
0 = min<br />
u∈{0,1} [(Mu − (1 + µ)I )v + d u ] + c.<br />
Explicitant cette équation de <strong>la</strong> programmation dynamique, on obtient<br />
:<br />
(λ + µ)v x = min {λv x+1 + c x ,λv 0 + k x } , ∀x = 0, ··· , E − 2.