FIABILIT

E
Corriges des exercices du polycopie
49 La duree de bon fonctionnement dun compteur electrique suit la loi exponentielle, de
moyenne 10 mois. Lorsquil est en panne, louvrier met un temps exponentiel de moyenne
15 jours avant de sen apercevoir et alors, la reparation, elle aussi exponentielle, dure en
moyenne 2 mois.
Determiner la disponibilite `a linstant t et en regime stationnaire.
On denit 3 etats : 1 le compteur fonctionne , 2 le compteur est en panne et louvrier ne la pas vu,
3 le compteur est en repartion.
Equations de Kolmogorov :
_
_
_
p

1
(t) = p
1
(t) + p
3
(t)
p

2
(t) = p
2
(t) + p
1
(t)
p

3
(t) = p
3
(t) + p
2
(t)
.
Equations des transformees de Laplace :
_
_
_
sp
1
(s) 1 = p
1
(s) + p
3
(s)
sp
2
(s) = p
2
(s) + p
1
(s)
sp
3
(s) = p
3
(s) + p
2
(s)
, soit
_
_
_
(s + )p
1
(s) = p
3
(s) + 1
(s + )p
2
(s) = p
1
(s)
(s + )p
3
(s) = p
2
(s)
.
On en deduit p
3
(s) =

s+
p
2
(s) =

(s+)(s+)
p
1
(s).
En reportant dans la premi`ere equation:
_
s +

(s+)(s+)
_
p
1
(s) = 1, soit
p
1
(s) =
(s + )(s + )
(s + )(s + )(s + )
=
(s + )(s + )
s(s
2
+ ( + + )s + + + )
.
On note s
1
et s
2
les racines de s
2
+ ( + + )s + + + .
On a alors p
1
(s) =
A
s
+
B
ss1
+
C
ss2
avec :
A =

+ +
; B =
(s
1
+ )(s
1
+ )
s
1
(s
1
s
2
)
; C =
(s
2
+ )(s
2
+ )
s
2
(s
2
s
1
)
.
Enn A(t) = p
1
(t) =

++
+
1
s1s2
_
(s1+)(s1+)
s1
e
s1t

(s2+)(s2+)
s2
e
s2t
_
.
A.N. : =
1
10
= 0, 1,
1

=
1
2
donc = 2 et
1

= 2 donc = 0, 5, = 1, s
1
et s
2
sont les racines de
s
2
+ 2, 6s + 1, 25.
On trouve s
1
0, 64, s
2
1, 96 et

++
=
1
1,25
= 0, 8, puis
A(t) = 0, 8 + 0, 22e
0,64t
0, 02e
1,96t
et lim
t+
A(t) = A = 0, 8
(A est la disponibilite en regime stationnaire).
2 DESS ISMAG
50 On dispose de 2 composants : C
1
est le composant principal, utilise en priorite et
C
2
le moteur de secours, mis en reserve. Les durees de fonctionnement et de reparation
sont exponentielle, les taux de defaillance
i
et de reparation
i
. C
2
ne tombent pas en
panne au repos, mais refuse de demarrer avec la probabilite quand on le lance (son
temps de reparation est suppose le meme sil est tombe en panne durant son service ou
au demarrage). Il y a un reparateur par moteur.
Faire le graphe. Determiner la disponibilite en regime stationnaire et la abilite
en regime transitoire.
A.N. :
1
=
2
= 1,
2
=
1
= 2 et =
1
2
.
On fait le graphe des transitions et on ecrit les equations de balance en considerant les 4 etats suivants
:
(1, 1) : C
1
marche, C
2
en reserve ;
(0, 1) : C
1
en reparation, C
2
marche ;
(0, 0) : C
1
et C
2
en reparation ;
(1, 0) : C
1
marche, C
2
en reparation.
On a les equations :

1
p
1,1
=
1
p
0,1
+
2
p
1,0
(
1
+
2
)p
1,0
=
1
p
0,0
(
2
+
1
)p
0,1
=
1
(1 )p
1,1
+
2
p
0,0
.
En remplacant dans la troisi`eme equation, on obtient :
(
2
+
1
)p
0,1
= (1 )
1
p
0,1
+
(1 )
1

1
+
2
p
0,0
+
2
p
0,0
do` u p
0,1
=
2
2+1
_
1 +
(1)1
1+2
_
p
0,0
et p
1,0
=
1
1+2
p
0,0
,
puis p
1,1
=
_
12
1(2+1)
_
1 +
(1)1
1+2
_
+
12
1(1+2)
_
p
0,0
.
On obtient alors p
0,0
en ecrivant p
0,0
+ p
0,1
+ p
1,0
+ p
1,1
= 1, soit
p
0,0
_
1 +

1

1
+

2
(
1
+
1
)

1
(
2
+
1
)
_
1 +
(1 )
1

1
+
2
__
1
= 1.
Avec
1
=
2
= 1,
2
=
1
= 2 et =
1
2
, on obtient p
0,0
= p
1,0
=
2
9
, p
0,1
=
1
9
et p
1,1
=
4
9
.
Pour determiner la abilite, on rend letat (0, 0) absorbant (plus de `eches sortantes et donc letat
(1, 0) nest plus `a considerer). Cette fois-ci, on ecrit les equations de Kolmororo en regime transitoire,
soit
_
_
_
p

1,1
(t) =
1
p
1,1
(t) +
1
p
0,1
(t)
p

0,1
(t) = (
2
+
1
)p
0,1
(t) +
1
(1 )p
1,1
(t)
p

0,0
(t) =
2
p
0,1
(t) +
1
p
1,1
(t)
et par transformee de Laplace,
_
_
_
sp
1,1
1 =
1
p
1,1
+
1
p
0,1
sp
0,1
= (
2
+
1
)p
0,1
+
1
(1 )p
1,1
sp
0,0
=
2
p
0,1
+
1
p
1,1
.
De la deuxi`eme equation, on tire p
1,1
=
s+2+1
1(1)
p
0,1
, puis en remplacant dans la premi`ere, on obtient
p
0,1
_
(s +
1
)(s +
2
+
1
)

1
(1 )

1
_
= 1
Fiabilite 3
et avec les valeurs numeriques, p
0,1
[2(s + 1)(s + 4) 2] = 1, soit p
0,1
=
1
2(s
2
+5s+3)
et
p
1,1
+ p
0,1
=
2s + 9
2(s
2
+ 5s + 3)
=
1
2(s
1
s
2
)
_
2s
1
+ 9
s s
1

2s
2
+ 9
s s
2
_
o` u s
1
et s
2
sont les racines de s
2
+5s+3 = 0, soit s
i
=
1
2
(5

13) et la abilite est R(t) = p

1,1
(t)+p
0,1
(t).
On a 2(s
1
s
2
) = 2

13 et 2s
i
+ 9 = 4

13, donc
R(t) =
1
2

13
_
(4 +

13)e

t
2
(5

13)
(4

13)e

t
2
(5+

13)
_
.
51 On reprend le podium de lexercice 46 avec N = 2 et = 1 et on sinteresse
maintenant `a sa abilite. Pour cela, on note T linstant de premier eondrement, et F sa
fonction de repartition.
En rendant letat C absorbant, verier que F(t) = p
C
(t).
a) Montrer que F verie lequation dierentielle F

+(2+)F

+F = 1 et resoudre
cette equation lorsque = 1/2. Determiner la loi de T et la duree moyenne avant le
premier eondrement.
b) Retrouver le resultat precedent `a laide des transformees de Laplace.
On rend letat C absorbant en supprimant la reparation du podium car on sinteresse au premier
eondrement. On ecrit les equations en regime transitoire, avec = 1 :
_
_
_
p

0
(t) = p
1
(t) p
0
(t)
p

1
(t) = p
0
(t) (1 + )p
1
(t)
p

C
(t) = p
1
(t)
Si T est linstant du premier eondrement, P([T t]) = P([X
t
= C]), soit F(t) = p
C
(t) .
a) On a alors F

(t) = p
1
(t) et F

(t) = p

1
(t) = p
0
(t) (1 +)p
1
(t) avec p
0
(t) = 1 p
1
(t) p
C
(t),
do` u F

(t) = 1 (2 + )p
1
(t) p
C
(t) et on a bien F

+ (2 + )F

+ F = 1 . Cette equation a pour

solution particuli`ere 1 et pour equation caracteristique associee r
2
+ (2 + )r + 1 = (r + 2)(r + 1/2) = 0
si = 1/2. On a alors F(t) = 1 +
1
e
2t
+
2
e
t/2
. Letat initial etant 0, on a p
C
(0) = F(0) = 0 et
p
1
(0) = F

(0) = 0, ce qui donne

1
+
2
= 1 et 2
1
+
1
2

2
= 0, soit
2
= 4
1
et 3
1
= 1, do` u

1
=
1
3
et
2
=
4
3
, et nalement F(t) = 1 +
1
3
e
2t

4
3
e
t/2
, et f
T
(t) =
2
3
e
t/2

2
3
e
2t
(pour t > 0).
On a alors
E(T) =
_
t f
T
(t) dt =
2
3
_
+
0
te
t/2
dt
2
3
_
+
0
te
2t
dt
=
2
3
_
+
0
4ue
u
du
2
3
_
+
0
v
4
e
v
dv
=
8
3

1
6
=
15
6
=
5
2
= 2, 5
en posant u = t/2 dans la premi`ere integrale et v = 2t dans la deuxi`eme.
b) Avec les transformees de Laplace, on a :
4 DESS ISMAG
_
_
_
s p
0
1 = p
1
p
0
s p
1
= p
0
(1 + ) p
1
s p
C
= p
1
, soit
_
_
_
(s + 1 + ) p
1
= p
0
(s + 1) p
0
= p
1
+ 1 = (s + 1)(s + 1 + ) p
1
s p
C
= p
1
, do` u
p
C
=
1
s
p
1
=
1
s
1
((s + 1)(s + 1 + ) )
=
1
s(s
2
+ (2 + )s + 1)
.
Avec =
1
2
, p
C
=
1
s(s+2)(s+1/2)
=
a
s
+
b
s+2
+
c
s+1/2
o` u a =
1
1/22
= 1, b =
1
2(2+1/2)
=
1
3
et
c =
1
1/2(1/2+2)
=
4
3
, donc p
C
= 1 +
1
3
1
s+2

4
3
1
s+1/2
et, par unicite de la transformee de Laplace, on
obtient p
C
(t) = 1 +
1
3
e
2t

4
3
e
t/2
: on retrouve le meme resultat quau a).
52 Un appareil fonctionne avec une batterie, de duree de charge exponentielle, de
moyenne
1

= 30 mn. La batterie dechargee doit etre rechargee sur secteur pendant un

temps exponentiel de moyenne
1

= 1 h. Pour minimiser lindisponibilite de lappareil, on

dispose de 3 batteries identiques, pouvant eventuellement etre rechargees simultanement.
a) Determiner la disponibilite en regime stationnaire.
b) On suppose maintenant que, si lappareil se retrouve sans batterie, il perd toutes
ses donnees et ne vaut plus rien. En utilisant les transformees de Laplace, determiner la
probabilite quil soit encore valable `a linstant t.
A.N. : t = 1 h.
a) Ici, un etat designe le nombre de batteries `a plat. On fait le graphe pour la disponibilite et on
determine les probabilites en regime stationnaire.
_

_
p
0
= p
1
p
1
= 2p
2
p
2
= 3p
3
p
0
+ p
1
+ p
2
+ p
3
= 1
soit, avec =

,
_

_
p
1
= p
0
p
2
=

2
2
p
0
p
3
=

3
6
3p
0
p
0
+ p
1
+ p
2
+ p
3
= 1
, do` u p
0
=
_
1 + +

2
2
+

3
6
_
1
et, comme A = 1 p
3
, on a
donc A =
1++

2
2
1++

2
2
+

3
6
.
Avec = 2 et = 1, soit = 2, on obtient A =
1+2+2
1+2+2+
4
3
=
15
19
, soit A 0, 79 .
b) Pour la abilite, on rend maintenant letat 3 absorbant et on ecrit les equations de Kolmogorov.
_

_
p

0
= p
0
+ p
1
p

1
= ( + )p
1
+ p
0
+ 2p
1
p

2
= ( + 2)p2 + p
1
p

3
= p
2
;
et pour resoudre ces equations, on passe par les transformees de Laplace :
_

_
(s + )p
0
1 = p
1
(s + + )p
1
= p
0
+ 2p
2
(s + + 2)p
2
= p
1
sp
3
= p
2
qui donne p
3
=

s
p
2
; p
1
=
s++2

p
2
puis, en reportant dans la deuxi`eme
equation, p
2
=

2
(s++)(s++2)2
p
0
puis p
1
=
(s++2)
s
2
+(2+3)s+
2
++2
2
et par la premi`ere
_
s +
(s + + 2)
s
2
+ (2 + 3)s +
2
+ + 2
2
_
p
0
= 1,
Fiabilite 5
soit p
0
=
s
2
+(2+3)s+
2
++2
2
s
3
+3(+)s
2
+(3
2
+3+2
2
)s+
3
et nalement p
3
=

3
s(s
3
+3(+)s
2
+(3
2
+3+2
2
)s+
3
)
.
Avec = 2 et = 1, ceci donne p
3
=
8
s(s
3
+9s
2
+20s+8)

1
s

1,35
s+0,51
+
0,44
s+2,71

0,09
s+5,78
(avec la
calculatrice).
On en en deduit R(t) = 1 p
3
(t) 1, 35e
0,51t
0, 44e
2,71t
+ 0, 09e
5,78t
.
En particulier, R(1) 0, 78.
53 Comparer la abilite `a un instant t quelconque des 2 dispositifs suivants, constitues
`a laide de 2n composants identiques :
1. en serie, n syst`emes formes de 2 composants en parall`ele ;
2. en parall`ele, 2 syst`emes formes de n composants en serie.
(On pourra poser u = 1 R(t), o` u R(t) designe la abilite dun composant `a linstant t).
Pour le montage en serie, on a R
1
= [1 (1 R(t))
2
]
n
= (1 u
2
)
n
.
Pour le montage en parall`ele, on a R
2
= 1 [1 R(t)n]
2
= 2(1 u)
n
(1 u)
2n
.
R
1
R
2
= (1 u)
n
[(1 + u)
n
2 + (1 u)
n
]
= (1 u)
n
_
n

k=0
C
k
n
u
k
2 +
n

k=0
C
k
n
(1)
k
u
k
_
= (1 u)
n
_
n

k=1
C
k
n
(1 + (1)
k
)u
k
_
0
car u ]0, 1[ et 1 + (1)
k
{0, 2}.
Ainsi, le syst`eme 1 est plus able que le syst`eme 2.