Exercice Avancr un Mod idologique SEIRV

onc l’exercice
Considns le mod SEIRV suivant pour la propagation d’une maladie infectieuse :
dS SI
= µN − β − µS − ρS
dt N
dE SI
=β − (σ + µ)E
dt N
dI
= σE − (γ + µ)I
dt
dR
= γI − µR
dt
dV
= ρS − µV
dt
avec :
— S(t), E(t), I(t), R(t), V (t) : les populations sensibles, expos, infectieuses, rblies et
vaccin instant t ;
— N = S + E + I + R + V : la population totale supposonstante ;
— µ : taux de mortalitaissance ;
— β : taux de transmission ;
— σ : taux de progression de l’exposition infection ;
— γ : taux de guson ;
— ρ : taux de vaccination.

Travail demand
1. Montrez que le mod est biologiquement cohnt : les solutions restent positives et
born si les conditions initiales sont positives.
2. Drminez le nombre de reproduction de base R0 .
3. Identifiez les points d’ilibre, notamment le point d’ilibre endque si R0 > 1.
4. udiez la stabilitcale du point d’ilibre endque aide de la matrice Jacobienne.
5. Proposez une fonction de Lyapunov pour dntrer la stabilitobale de l’ilibre endque.
6. Drminez le taux de vaccination minimal ρc nssaire pour garantir l’extinction de la
maladie (i.e. R0eff < 1).

1
Corrigtaill
1. Cohnce biologique :
— Le syst est un syst diffntiel autonome efficients continus. Par le th de Cauchy-
Lipschitz, il admet une solution unique.
— Si les conditions initiales sont positives, un raisonnement par contradiction et
le th d’invariance de Nagumo permettent de montrer que S, E, I, R, V ≥ 0 pour
tout t ≥ 0.
— Puisque dN dt
= µN − µN = 0, on a N (t) = N (0) constant, donc toutes les
variables sont born.
2. Nombre de reproduction de base R0 :
— En prnce de vaccination, le taux effectif de transmission est diminur la propor-
tion vaccin
— L’expression de R0 devient :
βσ µ
R0 = ·
(µ + σ)(µ + γ) µ + ρ
3. Points d’ilibre :
— uilibre sans maladie (ESM) : E = I = 0, donc

µN ρ ∗
S∗ = , V∗ = S , R∗ = N − S ∗ − V ∗ .
µ+ρ µ

— uilibre endque (si R0 > 1) : on rut le syst stationnaire complet avec I ∗ > 0,
ce qui donne aprubstitutions :

(µ + γ)(µ + σ)
S∗ = N,
βσ
(γ + µ) ∗
E∗ = I ,
σ
µ
I ∗ = (R0 − 1), avec R0 > 1.
β
4. Stabilitcale de lilibre endque :
— Calcul de la matrice Jacobienne J du syst au point d’ilibre E ∗ = (S ∗ , E ∗ , I ∗ , R∗ , V ∗ ).
— On restreint lde aux ations S, E, I, les autres nt dupl.
— La matrice Jacobienne restreinte est :
∗ ∗ 
−µ − ρ − β IN −β SN

0
∗ ∗
J = β IN −σ − µ β SN 
0 σ −γ − µ

— On vfie que les coefficients de Hurwitz (ou on montre par calcul direct des valeurs
propres) sont strictement positifs si R0 > 1 ⇒ stabilitcaleasymptotique.
5. Stabilitobale (fonction de Lyapunov) :
— Fonction candidate :
   
∗ ∗ E ∗ ∗ I
V (E, I) = a E − E − E ln ∗ + b I − I − I ln ∗
E I

2
 dV
— Calcul de dt
le long des trajectoires du syst :

E ∗ dE I ∗ dI
   
dV
=a 1− +b 1−
dt E dt I dt
= (expression ntive sous condition que a, b bien choisis)

— En utilisant le th de LaSalle, on conclut que lilibre endque est globalement

asymptotiquement stable.
6. Taux de vaccination minimal ρc pour que R0eff < 1 :
— On impose :
βσ µ
R0eff = · <1
(µ + σ)(µ + γ) µ + ρ
— Ce qui donne :  
βσ
ρ > ρc = µ −1
(µ + σ)(µ + γ)
— Interprtion : vacciner suffisamment de la population rit efficacement le rrvoir
de susceptibles et permet de contenir lde.

Dil de la stabilitcale (Crit de Routh-Hurwitz)

La matrice Jacobienne restreinte aux variables (S, E, I) au point d’ilibre endque est :
∗ ∗
−µ − ρ − β IN −β SN
 
0
∗ ∗
J = β IN −(σ + µ) β SN 
0 σ −(γ + µ)

Polynme caractstique
Le polynme caractstique P (λ) = det(J − λI) s’it :

P (λ) = λ3 + a1 λ2 + a2 λ + a3 = 0
avec les coefficients :
I∗
a1 = tr(J) = 3µ + ρ + σ + γ + β > 0
N
X I∗ I∗
a2 = minors principaux = (µ + ρ + β )(σ + µ) + (µ + ρ + β )(γ + µ)
N N
S∗
+ (σ + µ)(γ + µ) − σβ
N
a3 = − det(J) = Expression positive lorsque R0 > 1

Conditions de stabilit
D’apre crit de Routh-Hurwitz, les conditions pour la stabilitnt :
1. a1 > 0 (toujours vfi
2. a3 > 0 (vfiorsque R0 > 1)
3. a1 a2 − a3 > 0 (vfious la m condition)

3
Stabilitobale par fonction de Lyapunov
Construction de la fonction
On dnit la fonction de Lyapunov candidate :
   
∗ ∗ E σ+µ ∗ ∗ I
V (E, I) = E − E − E ln ∗ + I − I − I ln ∗
E σ I
| {z } | {z }
V1 V2

Dve long des trajectoires

En dvant et en utilisant les ations du syst :

E ∗ dE σ + µ I ∗ dI
   
dV
= 1− + 1−
dt E dt σ I dt
∗
  
E SI
= 1− β − (σ + µ)E
E N
I∗
 
σ+µ
+ 1− (σE − (γ + µ)I)
σ I

Simplification ilibre
∗ ∗
En utilisant les relations d’ilibre β SNI = (σ + µ)E ∗ et σE ∗ = (γ + µ)I ∗ , on obtient
apralculs :

dV (E − E ∗ )2 (σ + µ)(γ + µ) (I − I ∗ )2
= −(σ + µ) − ≤0
dt E σ I

Conclusion
La dvst ntive semi-dnie et s’annule uniquement en (E ∗ , I ∗ ). Par le th de LaSalle, l’ilibre
endque est globalement asymptotiquement stable dans R+ ∗.