COURS Chapitre4

Dynamique de populations en intéractions
2. Modèles discrets
Pr: Noha EL KHATTABI
Université Mohammed-V Agdal - Rabat

Faculté des sciences
17 mars 2020
Pr: Noha EL KHATTABI (COURS) FSR 17 mars 2020 1 / 13

Etude d’un système de deux équations en temps
discret
Point fixe et linéarisation.
Un système de deux équations en temps discret se définit de fao̧n générale par le
système suivante : 
 x = f (xn , yn )
n+1
 y = g (xn , yn )
n+1
où les variables xn et yn sont calculées à des intervalles de temps discrets :

0, 1, 2, ..., n, ...
Un point fixe du système est une solution constante du système précédent. Elle
vérifie donc pour tout n, les équations suivantes :

 x
n+1 = xn
 y = y
n+1 n

Ce qui nous ramènent à résoudre le système

 x = f (x, y )
 y = g (x, y )
Soit(x̂, ŷ ) un tel point fixe. Posons un = xn − x̂ et vn = yn − ŷ . La linéarisation du

système au voisinage de ce point fixe conduit au système suivant :

 u δf δf
n+1 = δx (x̂, ŷ )un + δy (x̂, ŷ )vn
 v δg δg
n+1 = δx (x̂, ŷ )un + δy (x̂, ŷ )vn
que l’on peut écrire    

un+1 un
  = A 
vn+1 vn
où A est la matrice jacobienne au point (x̂, ŷ ).

Cas de deux valeurs propres réelles distinctes.
Dans la nouvelle base formé e des vecteurs propres le système s’écrit :
   
un+1 un
 =J 
vn+1 vn
avec  
λ1 0
J = P −1 AP =  .
0 λ2
La solution évidente de ce système est :
un = (λ1 )n u0 , vn = (λ2 )n v0

On distingue trois cas :
|λ1 | > 1 et |λ2 | > 1. Les deux valeurs propres sont de module supérieur à 1.
Lorsque n → ∞, les solutions un et vn tendent vers ±∞. Toutes les
trajectoires s’éloignent du point fixe (0, 0) qui est par conséquent instable. Il
s’agit d’un noeud instable ou encore source (voir figure 1) .
|λ1 | < 1 et |λ2 | < 1. Les deux valeurs propres sont de module supéieur à 1.
Lorsque n → ∞, les solutions un et vn tendent vers 0. Toutes les trajectoires
tendent vers le point fixe (0, 0) qui est par conséquent stable. Il s’agit d’un
noeud stable ou encore puits. (voir figure 2)
|λ1 | < 1 et |λ2 | > 1 ou inversement, il s’agit d’un point selle ou point col.
Lorsque n → ∞, une des solutions tend vers 0 et l’autre tend vers ±∞. Le
point fixe (0, 0) est instable (voir figure 3)

Cas d’une valeur propre double.
Dans la nouvelle base formée des vecteurs propres le système s’ecrit :
   
un+1 un
 =J 
vn+1 vn
avec  
λ0 1
J = P −1 AP =  .
0 λ0
Et la solution s’écrit :
      
un u 0 λn nλn−1 u0
  = Jn  = 0 0  .
vn v0 0 λn0 v0
|λ0 | < 1. Lorsque n → ∞, les solutions un et vn tendent vers 0. Le point fixe

est un noeud stable.
|λ0 | > 1. Le point fixe est un noeud instable.

Cas de deux valeurs propres complexes conjuguées. Les valeurs propres sont
de la forme λ1 = α + iβ et λ2 = α − iβ. Dans la nouvelle base formée des
vecteurs propres le système s’écrit :
      
un+1 un α −β un
 =J =  .
vn+1 vn β α vn
En posant α = ρ cos ω et β = ρ sin ω, la solution s’écrit :

      
u u cos nω − sin nω u
 n  = J n  0  = ρn   0 
vn v0 sin nω cos nω v0
ρ < 1. Lorsque n → ∞, les solutions un et vn tendent vers 0. Le point fixe est

un noeud stable (voir figure 4).
ρ > 1. Le point fixe est noeud instable(voir figure 5).
ρ = 1, il s’agit d’un centre (voir figure 6).

On peut résumer ainsi ce qui précède :
(x ∗ , y ∗ ) est localement asymptotiquement stable ⇐⇒ |λi | < 1, i = 1, 2
On montre que cette condition peut s’exprimer en fonction de la trace et du
déterminant de l’équation caractéristique, et donc : (x ∗ , y ∗ ) est localement
asymptotiquement stable ⇐⇒ 2 > 1 + detA > |trA|. Cette dernière condition est
très utilisée en pratique.

Modèle de Nicholson-Bailey
C’est un modèle qui décrit la dynamique d’un système hôte parasitȯide d’une
génération à la suivante. Il s’écrit :

 x = λxn exp(−ayn )
n+1
 y
n+1 = cxn (1 − exp(−ayn )),
xn étant le nombre d’hôtes et yn le nombre de parasitȯides à la génération n. λ est

le taux de croissance intrinsèque des hôtes. La fonction exp(−ayn ) < 1 représente
la fraction des hôtes qui sont non parasités et vont donc contribuer à la génération
suivante des hôtes. a est un paramètre strictement positif qui quantifie l’impact du
parasitoide sur son hôte. (1 − exp(−ayn )) représente la fraction des hôtes qui sont
parasités et vont donc contribuer à la prochaine génération des parasitoides. c est
le nombre moyen de parasitoides viables issus d’un seul hôte parasité.
Les points fixes de ce système sont solutions du système suivant :

 x ∗ = λx ∗ exp(−ay ∗ )
 y ∗ = cx ∗ (1 − exp(−ay ∗ )),
L’origine est point fixe ainsi que le point non trivial (x ∗ , y ∗ ) tel que :
λln(λ)
x∗ =
ac(λ − 1)
ln(λ)
y∗ = ,
a
Le taux de croissance des hôtes doit être supérieur à 1 pour que ce point ait un
sens biologiquement parlant. Il s’agit d’une hypothèse réaliste puisqu’elle signifie
qu’en l’absence de parasitoide la population hôte se développerait. Pour étudier la
stabilité locale des points fixes, calculons la matrice Jacobienne :
   
∗ ∗ ∗ ∗
λ exp(−ay ) −aλx exp(−ay ) 1 −ax
J(x ∗ , y ∗ ) =  =
∗
,
c(1 − exp(−ay ∗ )) acλx ∗ exp(−ay ∗ ) c(1 − λ1 ) acx
λ

dont le déterminant est :
λlnλ
detJ(x ∗ , y ∗ ) = .
λ−1
On vérifie aisément que ce déterminant est strictement supérieur à un et le point
fixe (x ∗ , y ∗ ) est donc instable. Voir figure avec a = 0.15,λ = 2 et c = 1.

COURS Chapitre4

Transféré par

Droits d'auteur :

Formats disponibles

COURS Chapitre4

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

COURS Chapitre4

Transféré par

Droits d'auteur :

Formats disponibles

Dynamique de populations en intéractions

Pr: Noha EL KHATTABI

Université Mohammed-V Agdal - Rabat

Pr: Noha EL KHATTABI (COURS) FSR 17 mars 2020 1 / 13

où les variables xn et yn sont calculées à des intervalles de temps discrets :

Pr: Noha EL KHATTABI (COURS) FSR 17 mars 2020 2 / 13

Soit(x̂, ŷ ) un tel point fixe. Posons un = xn − x̂ et vn = yn − ŷ . La linéarisation du

que l’on peut écrire    

Pr: Noha EL KHATTABI (COURS) FSR 17 mars 2020 3 / 13

Pr: Noha EL KHATTABI (COURS) FSR 17 mars 2020 4 / 13

Pr: Noha EL KHATTABI (COURS) FSR 17 mars 2020 5 / 13

|λ0 | < 1. Lorsque n → ∞, les solutions un et vn tendent vers 0. Le point fixe

|λ0 | > 1. Le point fixe est un noeud instable.

En posant α = ρ cos ω et β = ρ sin ω, la solution s’écrit :

ρ < 1. Lorsque n → ∞, les solutions un et vn tendent vers 0. Le point fixe est

ρ > 1. Le point fixe est noeud instable(voir figure 5).

ρ = 1, il s’agit d’un centre (voir figure 6).

Pr: Noha EL KHATTABI (COURS) FSR 17 mars 2020 10 / 13

xn étant le nombre d’hôtes et yn le nombre de parasitȯides à la génération n. λ est

Pr: Noha EL KHATTABI (COURS) FSR 17 mars 2020 12 / 13

Pr: Noha EL KHATTABI (COURS) FSR 17 mars 2020 13 / 13

Vous aimerez peut-être aussi