MODÉLISATION ET SIMULATION

DESSIN QUALITATIF DES PORTRAITS DE PHASE

(IV/IV)

1. linéarisation de systèmes non-linéaires

– Soit le système


ẋ1 = f1 (x1 , x2 )



ẋ2 = f2 (x1 , x2 )


et la matrice  
∂f1 ∂f1
∂x
 1 ∂x2 
J=
∂f2 ∂f2
∂x1 ∂x2

(matrice jacobienne).
– La matrice jacobienne, évaluée en un point stationnaire
[x̄1 , x̄2 ] donne une approximation linéaire du comporte-
ment du système au voisinage de ce point.

1
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
2

2. Exercice 1 : dessin qualitatif d’un système

non-linéaire

Soit le système non-linéaire :



 ẋ = x

1 1

 ẋ2 = x21 + x22 − 1




Les étapes à suivre pour tracer l’espace des phases sont les
suivantes:

(1) Calcul des isoclines et partition du plan en régions.

(2) Calcul des points stationnaires (par annulation des

dérivées).

(3) Calcul de la matrice jacobienne.

(4) Pour chaque point stationnaire :

(a) Appliquer la matrice stationnaire au point.

(b) Etudier le système linéaire qui en découle, et qui

donnera une approximation du comportement du
système au voisinage de ce point stationnaire.

(5) Tracer des trajectoires qualitativement.

MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
3

Solution :

(1) Calcul des isoclines

La première isocline est x1 = 0 (et sur cette droite
ẋ2 = x22 − 1).
La deuxième isocline est le cercle ayant pour équation
x21 + x22 = 1 (cercle de rayon 1 autour de l’origine).

(2) Calcul des points stationnaires

On annule les dérivées:


ẋ1 = x1 = 0


ẋ2 = x21 + x22 − 1 = 0




Ce système a deux solutions possibles: [0, 1], [0 − 1].

(3) Calcul de la matrice jacobienne

   
∂f1 ∂f1
∂x ∂x2  1 0
J=  1 = 
∂f2 ∂f2
∂x1 ∂x2 2x1 2x2
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
4

(4) Linéarisation autour du point stationnaire [0, 1]

 
1 0
J = 
0 2

qui a deux valeurs propres positives et les vecteurs pro-

pres v1 = [0, 1]T et v2 = [1, 0]T .
Le point est donc un noeud instable.

(5) Linéarisation autour du point stationnaire [0, −1]

 
1 0
J = 
0 −2

qui a une valeur propre réelles distinctes: l’une positive

et l’autre negative, et les vecteurs propres v1 = [0, 1]T
et v2 = [1, 0]T . Le point est donc une selle.
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
5

(6) Dessin qualitatif des trajectoires

dx1 /dt=0
−+ ++

−− +−

dx2 /dt=0

−+
++
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
6

(7) Solution calculée avec Matlab

1
2

0
x

−1

−2

−3
−3 −2 −1 0 1 2 3
x1
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
7

3. Exercice 2: dessin qualitatif d’un système

non-linéaire

Soit le système non-linéaire :



 ẋ = x

1 2

 ẋ2 = x1 (1 − x21 ) + x2



Les étapes à suivre pour tracer l’espace des phases sont les
suivantes:

(1) Calcul des isoclines et partition du plan en régions.

(2) Calcul des points stationnaires (par annulation des

dérivées).

(3) Calcul de la matrice jacobienne.

(4) Pour chaque point stationnaire :

(a) Appliquer la matrice stationnaire au point.

(b) Etudier le système linéaire qui en découle, et qui

donnera une approximation du comportement du
système au voisinage de ce point stationnaire.

(5) Tracer des trajectoires qualitativement.

MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
8

Solution :

(1) Calcul des isoclines

– Première isocline (ẋ1 = 0 ) :

– x2 = 0.

1
x2

−1

−2

−3
−3 −2 −1 0 1 2 3
x1
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
9

– Sens des vecteurs vitesses verticaux sur cette iso-

cline (signe de ẋ2 ):
Sur cette droite


> 0 pour x1 < −1 et 0 < x1 < 1

ẋ2 =

< 0 pour − 1 < x1 < 0 et x1 > 1


1
x2

−1

−2

−3
−3 −2 −1 0 1 2 3
x1

– Partage du plan en régions ẋ1 < 0 et ẋ1 > 0:

ẋ1 est positive au-dessus de la droite et négative
en-dessous.
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
10

– Deuxième isocline (ẋ2 = 0) :

– x2 = x31 − x1 .

1
x2

−1

−2

−3
−3 −2 −1 0 1 2 3
x1
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
11

– Sens des vecteurs vitesses horizontaux sur cette

isocline (signe de ẋ1 ):
Sur cette courbe ẋ1 = x2 , donc le signe de ẋ1 est
égal au signe de x2

1
x2

−1

−2

−3
−3 −2 −1 0 1 2 3
x1

– Partage du plan en régions ẋ2 < 0 et ẋ2 > 0:

ẋ2 est positive au-dessus de la courbe et négative
en-dessous.
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
12

(2) Calcul des points stationnaires

On annule les dérivées:


ẋ = x = 0

1 2

ẋ2 = x1 (1 − x21 ) + x2 = 0



Ce système a trois solutions possibles: [0, 0], [1, 0] et

[−1, 0].

(3) Calcul de la matrice jacobienne

   
∂f1 ∂f1
∂x ∂x2  0 1
J=  1 = 
∂f2 ∂f2
∂x1 ∂x2 1 − 3x2 1

(4) Linéarisation autour du point stationnaire [0, 0]

 
0 1
J = 
1 1

qui a pour valeurs propres

√
1± 5
λ1,2 =
2

et les vecteurs propres v1 = [1, λ1 ]T et v2 = [1, λ2 ]T .

Le point est donc une selle.
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
13

(5) Linéarisation autour du point stationnaire [1, 0]

 
0 1
J = 
−2 1

qui a pour valeurs propres

√
1±i 7
λ1,2 =
2

Le point est donc un feu instable.

(6) Linéarisation autour du point stationnaire [−1, 0] :

idem que pour [1, 0]

(7) Dessin qualitatif des trajectoires

MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
14

−−
+−

dx2 /dt=0
dx1 /dt=0

dx2 /dt=0
++

−−
−+
++
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
15

(8) Solution calculée avec Matlab

1
x2

−1

−2

−3
−3 −2 −1 0 1 2 3
x1
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
16

4. Exercice 3: dessin qualitatif d’un système

non-linéaire (Lotka-Voltera)

Soit le système non-linéaire :



 ẋ1 = 0.2x1 − 0.08x1 x2



 ẋ2 = 0.1x1 x2 − 0.2x2


Les étapes à suivre pour tracer l’espace des phases sont les
suivantes:

(1) Calcul des isoclines et partition du plan en régions.

(2) Calcul des points stationnaires (par annulation des

dérivées).

(3) Calcul de la matrice jacobienne.

(4) Pour chaque point stationnaire :

(a) Appliquer la matrice stationnaire au point.

(b) Etudier le système linéaire qui en découle, et qui

donnera une approximation du comportement du
système au voisinage de ce point stationnaire.

(5) Tracer des trajectoires qualitativement.

MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
17

(1) Calcul des isoclines

– Première isocline (ẋ1 = 0) :

– x2 = 0.2/0.08 = 2.5.
– Sens des vecteurs vitesses verticaux sur cette iso-
cline (signe de ẋ2 ) :
Sur cette droite


> 0 pour x1 > 2

ẋ2 =

< 0 pour x1 < 2


– Partage du plan en régions ẋ1 < 0 et ẋ1 > 0 :

ẋ1 est négative au-dessus de la droite et positive
en-dessous.
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
18

– Deuxième isocline (ẋ2 = 0 ) :

– x1 = 0.2/0.1 = 2.
– Sens des vecteurs vitesses horizontaux sur cette
isocline (signe de ẋ1 ) :
Sur cette droite


> 0 pour x2 < 2.5

ẋ1 =

< 0 pour x2 > 2.5


– Partage du plan en régions ẋ2 < 0 et ẋ2 > 0 :

ẋ2 est négative (descente) à gauche de la droite
et positive (montée) à sa droite.

(2) Calcul des points stationnaires

On annule les dérivées :


ẋ1 = 0.2x1 − 0.08x1 x2 = 0



ẋ2 = 0.1x1 x2 − 0.2x2 = 0


Ce système a deux solutions possibles : [0, 0], [2, 2.5].

Le point stationnaire [2, 2.5] est à l’intersection des
isoclines.
MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
19

(3) Calcul de la matrice jacobienne

   
∂f1 ∂f1
∂x ∂x2  0.2 − 0.08x2 −0.08x1
J=  1 = 
∂f2 ∂f2
∂x1 ∂x2 0.1x2 0.1x1 − 0.2

(4) Linéarisation autour du point stationnaire [0, 0]

 
0.2 0
J = 
0 −0.2

qui a pour valeurs propres

λ1,2 = ±0.2

Le point est donc (localement) une selle.

(5) Linéarisation autour du point stationnaire [2, 2.5]

 
0 −0.16
J = 
0.25 0

qui a pour valeurs propres

√
λ1,2 = ±i 0.04

Le point est donc un centre.

MODÉLISATION ET SIMULATIONDESSIN QUALITATIF DES PORTRAITS DE PHASE (IV/IV)
20

(6) Solution calculée avec Matlab

4
x2

0
0 1 2 3 4 5 6
x1