Correction de TD:

Résolution d’équations différentielles ordinaires

Unité pédagogique: Mathématiques

logo.eps

UP Mathématique 1 26 Avril 2020 1 / 26

Plan

1 Exercice 1: Schéma d’Euler implicite.

2 Correction exercice 1: Schéma d’Euler implicite.
3 Exercice 2: Schéma d’Euler explicite et implicite.
4 Correction exercice 2: Schéma d’Euler explicite et implicite.
5 Exercice 3: Schéma d’Euler explicite .
6 Correction exercice 3: Schéma d’Euler explicite.

logo.eps

UP Mathématique 2 26 Avril 2020 2 / 26

Plan

1 Exercice 1: Schéma d’Euler implicite

2 Correction Exercice 1: Schéma d’Euler implicite

3 Exercice 2: Schéma d’Euler explicite et implicite

4 Correction exercice 2: Schéma d’Euler explicite et implicite

5 Exercice 3: Schéma d’Euler explicite

6 Correction exercice 3: Schéma d’Euler explicite

logo.eps

UP Mathématique 3 26 Avril 2020 3 / 26

Exercice 1

On considère le problème de Cauchy suivante

1
 0
 x (t) = − x(t)
(P C) : 1 + t2

x(0) = 5

1 Vérifier que la solution analytique du problème de Cauchy (P C) est donnée par:

x(t) = 5 exp(− arctan(t)) pour t ≥ 0.
2 Donner le schéma d’Euler implicite (régressif) avec un pas de temps h constant.
3 En déduire que pour h = 12 , la solution numérique xn+1 (approchant la solution exacte x
au point de discrétisation tn+1 , n ≥ 0) du problème de Cauchy (P C) trouvée par la
méthode d’Euler implicite vérifie la relation suivante :
4 + (n + 1)2
xn+1 = xn , ∀n ≥ 0.
6 + (n + 1)2
4 Appliquer le schéma itératif de la question (3) pour résoudre numériquement (P C) sur
l’intervalle [0, 2].
logo.eps
5 Calculer l’erreur commise par la méthode d’Euler implicite au point t = 2.

UP Mathématique 4 26 Avril 2020 4 / 26

Plan

1 Exercice 1: Schéma d’Euler implicite

2 Correction Exercice 1: Schéma d’Euler implicite

3 Exercice 2: Schéma d’Euler explicite et implicite

4 Correction exercice 2: Schéma d’Euler explicite et implicite

5 Exercice 3: Schéma d’Euler explicite

6 Correction exercice 3: Schéma d’Euler explicite

logo.eps

UP Mathématique 5 26 Avril 2020 5 / 26

Correction

Question 1
Pour que x(t) = 5 exp(− arctan(t)) soit la solution analytique de (P C) pour t ≥ 0 if faut que:
t 7→ x(t) de classe C 1 sur [0, +∞[.
t 7→ x(t) vérifie le système (P C).
On a
Il est évident que t 7→ 5 exp(− arctan(t)) est de classe C 1 sur [0, +∞[.
On a
0 0 1 1
x (t) = −5(arctan(t)) exp(− arctan(t)) = −5 exp(− arctan(t)) = − x(t).
1 + t2 1 + t2
Donc x(t) = 5 exp(− arctan(t)) vérifie la première équation du système (P C).
De plus pour t = 0 on a x(0) = 5 exp(− arctan(0)) = 5 exp(0) = 5. Donc
x(t) = 5 exp(− arctan(t)) vérifie la deuxième équation du système (P C).
Ainsi x(t) = 5 exp(− arctan(t)) est la solution analytique de (P C) pour t ≥ 0.
Question 2
On commence d’abord par la discrétisation de l’intervalle de temps avec un pas de temps h on
pose: t0 = 0, t1 = t0 + h = h,...,tn = t0 + nh = nh,....
On note x(tn ) = xn .
logo.eps

UP Mathématique 6 26 Avril 2020 6 / 26

Correction
Le schéma d’Euler implicite (régressif) avec un pas de temps h pour le système (P C) consiste à
écrire:
0 1
x (tn+1 ) = f (tn+1 , xn+1 ) avec f (t, x) = − x.
1 + t2
0
On utilisant la méthode de dérivation numérique pour x (t) en t = tn+1 par:
0 x(tn+1 ) − x(tn ) x(tn+1 ) − x(tn )
x (tn+1 ) = = .
tn+1 − tn h
Donc le schéma d’Euler implicite:
x(tn+1 ) − x(tn )
= f (tn+1 , xn+1 ), ∀n ≥ 0 ⇔
h
x(tn+1 ) − x(tn ) 1
=− x(tn+1 ), ∀n ≥ 0 ⇔
h 1 + t2n+1
1
x(tn+1 ) = x(tn ) − h x(tn+1 ), ∀n ≥ 0 ⇔
1 + t2n+1
1
xn+1 = xn − h xn+1 , ∀n ≥ 0.
1 + t2n+1
Donc
schéma d’Euler implicite pour (P C)
1 logo.eps
∀n ≥ 0 on a xn+1 = xn − h xn+1 , x0 = 5.
1 + t2n+1
UP Mathématique 7 26 Avril 2020 7 / 26
Correction
Question 3
1 1 n n+1
Pour h = alors on a t0 = 0, t1 = t0 + h = ,...,tn = t0 + nh = ainsi tn+1 = . si on
2 2 2 2
remplace dans le schéma d’Euler implicite pour (P C) alors on a:

1
∀n ≥ 0 on a xn+1 = xn − h xn+1 ⇔
1 + t2n+1
1
∀n ≥ 0 on a xn+1 = xn − xn+1 ⇔
2(1 + t2n+1 )
1
∀n ≥ 0 on a xn+1 = xn − xn+1 ⇔
2(1 + ( n+1
2
)2 )
1
∀n ≥ 0 on a xn+1 = xn − xn+1 ⇔
(n+1)2
2(1 + 4
)
1
∀n ≥ 0 on a xn+1 (1 + ) = xn ⇔
(n+1)2
2(1 + 4
)
2
∀n ≥ 0 on a xn+1 (1 + ) = xn ⇔
4 + (n + 1)2
6 + (n + 1) 2
∀n ≥ 0 on a xn+1 ( ) = xn ⇔
4 + (n + 1)2
4 + (n + 1)2 logo.eps
∀n ≥ 0 on a xn+1 = xn ( ).
6 + (n + 1)2

UP Mathématique 8 26 Avril 2020 8 / 26

Correction
schéma d’Euler implicite pour (P C)
4 + (n + 1)2
∀n ≥ 0 on a xn+1 = xn ( ), x0 = 5.
6 + (n + 1)2

Question 4
1 3
On cherche la solution du problème (P C) au point t0 = 0, t1 = , t2 = 1, t3 = et t4 = 2
2 2
Pour t0 = 0 ,x0 = 5.
1
Pour Pour n = 0 tn+1 = t0+1 = t1 = on a
2
4 + (0 + 1)2 5 25
x1 = x0 ( ) = 5( ) = = 3.571428571
6 + (0 + 1)2 7 7
Pour Pour n = 1 tn+1 = t1+1 = t2 = 1 on a
4 + (1 + 1)2 25 8 20
x2 = x1 ( )= ( )= = 2.857142857
6 + (1 + 1)2 7 10 7
3
Pour Pour n = 2 tn+1 = t2+1 = t3 = on a
2
4 + (2 + 1)2 20 13 52
x3 = x2 ( )= ( )= = 2.476190476
6 + (2 + 1)2 7 15 21
Pour Pour n = 3 tn+1 = t3+1 = t4 = 2 on a logo.eps
4 + (3 + 1)2 52 10 520
x4 = x3 ( )= ( )= = 2.251082251.
6 + (3 + 1)2 21 11 231
UP Mathématique 9 26 Avril 2020 9 / 26
Correction

1 3
La solution du problème (P C) au point t0 = 0, t1 = , t2 = 1, t3 = et t4 = 2 est donnée
2 2
par le tableau suivant:

ti t0 t1 t2 t3 t4
xi 5 3.571428571 2.857142857 2.476190476 2.251082251

Table : Valeur de xi à l’instant ti .

Question 5
L’erreur commise par la méthode d’Euler implicite au point t = 2 est donnée par:

E(t = 2) = |x(t = 2) − x4 | or
x(2) = 5 exp(− arctan(2)) = 1.652499838 donc
E(t = 2) = |1.652499838 − 2.251082251| = 0.598582413.

logo.eps

UP Mathématique 10 26 Avril 2020 10 / 26

Plan

1 Exercice 1: Schéma d’Euler implicite

2 Correction Exercice 1: Schéma d’Euler implicite

3 Exercice 2: Schéma d’Euler explicite et implicite

4 Correction exercice 2: Schéma d’Euler explicite et implicite

5 Exercice 3: Schéma d’Euler explicite

6 Correction exercice 3: Schéma d’Euler explicite

logo.eps

UP Mathématique 11 26 Avril 2020 11 / 26

Exercice

On considère le problème de Cauchy défini par :

x0 = −λx,
(
λ, t ≥ 0
(P C)
x(0) = x0 .

1 Donner la solution analytique x(t) de (P C).

2 Calculer lim x(t).
t→+∞
3 Montrer que les solutions numériques du problème de Cauchy (P C) trouvées par la
méthode d’Euler explicite xE I
n et par la méthode d’Euler implicite xn , sont données
respectivement par :
xE n
n = x0 (1 − λh) ,

xIn = x0 (1 + λh)−n ,

où h désigne le pas de discrétisation et n > 0.

4 Calculer les limites lim xE
n et lim xIn .
n→+∞ n→+∞
logo.eps

UP Mathématique 12 26 Avril 2020 12 / 26

Plan

1 Exercice 1: Schéma d’Euler implicite

2 Correction Exercice 1: Schéma d’Euler implicite

3 Exercice 2: Schéma d’Euler explicite et implicite

4 Correction exercice 2: Schéma d’Euler explicite et implicite

5 Exercice 3: Schéma d’Euler explicite

6 Correction exercice 3: Schéma d’Euler explicite

logo.eps

UP Mathématique 13 26 Avril 2020 13 / 26

Correction

Question 1 Z
La solution analytique de (P C) est x(t) = K exp(A(t)) avec A(t) = −λdt = −λt et K
constante.
Donc x(t) = K exp(−λt) or
D’une part pour t = 0 on a x(0) = x0 .
D’autre part x(0) = K exp(−λ ∗ 0) = K exp(0) = K.
ainsi K = x0 d’ou la solution analytique de (P C) est x(t) = x0 exp(−λt).
Question 2
On a lim x(t) = lim x0 exp(−λt) = 0.
t→+∞ t→+∞
Question 3
On note les solutions numériques du problème de Cauchy (P C) par:
xE
n la solution par la méthode d’Euler explicite.
xIn la solution par la méthode d’Euler implicite.

logo.eps

UP Mathématique 14 26 Avril 2020 14 / 26

Correction

on commence d’abord par la discrétisation de l’intervalle de temps avec un pas de temps h on

pose: t0 = 0, t1 = t0 + h = h,...,tn = t0 + nh = nh,....
La méthode d’Euler explicite:
On note x(tn ) = xE
n.
0
En intégrant l’équation différentielle x (t) = f (t, x(t)) avec f (t, x(t)) = −λx(t) entre tn
et tn+1 on obtient:
Z tn+1 Z tn+1
0
x (t)dt = f (t, x(t))dt.
tn tn
Z tn+1
On approcher f (t, x(t))dt par la méthode des rectangles à gauche:
tn
x(tn+1 ) − x(tn ) = (tn+1 − tn )f (tn , x(tn )), ∀n ≥ 0 ⇔
x(tn+1 ) − x(tn ) = −hλx(tn ), ∀n ≥ 0 ⇔
xE E E
n+1 − xn = −hλxn , ∀n ≥ 0 ⇔
E E
xn+1 = (1 − hλ)xn , ∀n ≥ 0.

logo.eps

UP Mathématique 15 26 Avril 2020 15 / 26

Correction

On a xE E
n+1 = (1 − hλ)xn , ∀n ≥ 0.
Montrer par récurrence que:

∀n ≥ 0 on a xE n E
n = (1 − hλ) x0 .

Pour n = 0 on a (1 − hλ)0 xE E
0 = x0 donc vrai pour n = 0.
Supposons pour n fixée on a xE n E
n = (1 − hλ) x0 .
Montrons que xEn+1 = (1 − hλ)
n+1 xE :
0
On a xEn+1 = (1 − hλ)x E or d’après la supposition on xE = (1 − hλ)n xE .
n n 0
Ainsi xE n E
n+1 = (1 − hλ)(1 − hλ) x0 = (1 − hλ)
n+1 xE . D’où on a montrer par récurrence
0
que:
∀n ≥ 0 on a xE n E
n = (1 − hλ) x0 .

logo.eps

UP Mathématique 16 26 Avril 2020 16 / 26

Correction

La méthode d’Euler implicite:

On note x(tn ) = xIn .
0
En intégrant l’équation différentielle x (t) = f (t, x(t)) avec f (t, x(t)) = −λx(t) entre tn
et tn+1 on obtient:
Z tn+1 Z tn+1
0
x (t)dt = f (t, x(t))dt.
tn tn
Z tn+1
On approcher f (t, x(t))dt par la méthode des rectangles à droite:
tn
x(tn+1 ) − x(tn ) = (tn+1 − tn )f (tn+1 , x(tn+1 )), ∀n ≥ 0 ⇔
x(tn+1 ) − x(tn ) = −hλx(tn+1 ), ∀n ≥ 0 ⇔
xIn+1 − xIn = −hλxIn+1 , ∀n ≥ 0 ⇔
xIn+1 (1 + hλ) = xIn , ∀n ≥ 0 ⇔
1
xIn+1 = xI , ∀n ≥ 0.
(1 + hλ) n

logo.eps

UP Mathématique 17 26 Avril 2020 17 / 26

Correction
1
On a xIn+1 = xI , ∀n ≥ 0.
(1 + hλ) n
Montrer par récurrence que:

1
∀n ≥ 0 on a xIn = xI .
(1 + hλ)n 0

1
Pour n = 0 on a xI = xI0 donc vrai pour n = 0.
(1 + hλ)0 0
1
Supposons pour n fixée on a xIn = xI .
(1 + hλ)n 0
1
Montrons que xIn+1 = xI :
(1 + hλ)n+1 0
1 1
On a xIn+1 = xI or d’après la supposition on xIn = xI .
(1 + hλ) n (1 + hλ)n 0
1 1 1
Ainsi xIn+1 = xI = xI . D’où on a montrer par
(1 + hλ) (1 + hλ)n 0 (1 + hλ)n+1 0
récurrence que:
1
∀n ≥ 0 on a xIn = xI . logo.eps
(1 + hλ)n 0

UP Mathématique 18 26 Avril 2020 18 / 26

Correction

Question 4
1
On a lim xIn = lim xI =?.
n→+∞ + hλ)n 0
n→+∞ (1
On a h > 0 et λ > 0 donc 1 + hλ > 1 donc lim (1 + hλ)n = +∞ alors
n→+∞
1
lim = 0.
n→+∞ (1 + hλ)n
I
Donc lim xn = 0.
n→+∞

lim xE
n = lim (1 − hλ)n xE
0 =?
n→+∞ n→+∞

Si | 1 − hλ |< 1 ⇔ −1 < 1 − hλ < 1 ⇔ −2 < −hλ < 0 ⇔ 0 < hλ <

2
2 ⇔ 0 < h < alors: lim (1 − hλ)n = 0 donc lim xE n = 0.
λ n→+∞ n→+∞
Si | 1 − hλ |> 1 alors lim (1 − hλ) = ∞ donc lim xE
n
n = ∞.
n→+∞ n→+∞
Si | 1 − hλ |= 1 ⇔ hλ = 0 ⇔ h = 0ouλ = 0 imposible car h > 0 et
λ > 0.
logo.eps

UP Mathématique 19 26 Avril 2020 19 / 26

Plan

1 Exercice 1: Schéma d’Euler implicite

2 Correction Exercice 1: Schéma d’Euler implicite

3 Exercice 2: Schéma d’Euler explicite et implicite

4 Correction exercice 2: Schéma d’Euler explicite et implicite

5 Exercice 3: Schéma d’Euler explicite

6 Correction exercice 3: Schéma d’Euler explicite

logo.eps

UP Mathématique 20 26 Avril 2020 20 / 26

Exercice

Le problème de Cauchy défini par :

2t
(
x0 = x + e , t≥0
(P C) (1)
x(0) = 2.

t 2t
possède la solution analytique suivante x(t) = e + e .
1 Donner le schéma d’Euler explicite (progressif) avec un pas de temps h constant.
2 En prenant h = 0.1:
a) Faire 3 itérations de la méthode d’Euler explicite.
b) Calculer l’erreur commise sur x3 , la valeur approchée de x(0.3).
3 En prenant h = 0.05:
1 Faire 6 itérations de la méthode d’Euler explicite et
2 Calculer l’erreur commise sur x6 , la valeur approchée de x(0.3).
4 Interpréter les résultats obtenus.
logo.eps

UP Mathématique 21 26 Avril 2020 21 / 26

Plan

1 Exercice 1: Schéma d’Euler implicite

2 Correction Exercice 1: Schéma d’Euler implicite

3 Exercice 2: Schéma d’Euler explicite et implicite

4 Correction exercice 2: Schéma d’Euler explicite et implicite

5 Exercice 3: Schéma d’Euler explicite

6 Correction exercice 3: Schéma d’Euler explicite

logo.eps

UP Mathématique 22 26 Avril 2020 22 / 26

Correction
Question 1
On commence d’abord par la discrétisation de l’intervalle de temps avec un pas de temps h on
pose: t0 = 0, t1 = t0 + h = h,...,tn = t0 + nh = nh,....
On note x(tn ) = xn .
Le schéma d’Euler explicite (progressif) avec un pas de temps h pour le système (P C) consiste
à écrire:
0
x (tn ) = f (tn , xn ) avec f (t, x(t)) = x(t) + exp(2t).
0
On utilisant la méthode de dérivation numérique pour x (t) en t = tn par:
0 x(tn+1 ) − x(tn ) x(tn+1 ) − x(tn )
x (tn ) = = .
tn+1 − tn h
Donc le schéma d’Euler explicite:
x(tn+1 ) − x(tn )
= f (tn , xn ), ∀n ≥ 0 ⇔
h
x(tn+1 ) − x(tn )
= x(tn ) + exp(2tn ), ∀n ≥ 0 ⇔
h
x(tn+1 ) = x(tn ) + h(xn + exp(2tn )), ∀n ≥ 0 ⇔
xn+1 = (1 + h)xn + h exp(2tn ), ∀n ≥ 0.
Donc
schéma d’Euler explicite pour (P C) logo.eps

∀n ≥ 0 on a xn+1 = (1 + h)xn + h exp(2tn ), ∀n ≥ 0, x0 = 2.

UP Mathématique 23 26 Avril 2020 23 / 26
Correction

Question 2
En prenant h = 0.1: donc le schéma d’euler explicite s’écrit:

∀n ≥ 0 on a xn+1 = 1.1 ∗ xn + 0.1 exp(2tn ), ∀n ≥ 0, x0 = 2.

Question a):
1 Itération 1 pour n = 0 on x1 = 1.1 ∗ x0 + 0.1 exp(2t0 ) or x0 = 2 et t0 = 0. Donc
x1 = 2.3.
2 Itération 2 pour n = 1 on x2 = 1.1 ∗ x1 + 0.1 exp(2t1 ) or x1 = 2.3 et t1 = 0.1. Donc
x2 = 2.6521.
3 Itération 3 pour n = 2 on x3 = 1.1 ∗ x2 + 0.1 exp(2t2 ) or x2 = 2.6521 et t2 = 0.2. Donc
x3 = 3.0665.
Question b): L’erreur commise sur x3 , la valeur approchée de x(0.3) est donnée par:
E3 =| x3 − x(0.3) |=| 3.0665 − exp(0.3) − exp(0.6) |= 0.1055.

logo.eps

UP Mathématique 24 26 Avril 2020 24 / 26

Correction
Question 3
En prenant h = 0.05: donc le schéma d’euler explicite s’écrit:

∀n ≥ 0 on a xn+1 = 1.05 ∗ xn + 0.05 exp(2tn ), ∀n ≥ 0, x0 = 2.

Question a):
1 Itération 1 pour n = 0 on x1 = 1.05 ∗ x0 + 0.05 exp(2t0 ) or x0 = 2 et t0 = 0. Donc
x1 = 2.15.
2 Itération 2 pour n = 1 on x2 = 1.05 ∗ x1 + 0.05 exp(2t1 ) or x1 = 2.15 et t1 = 0.05. Donc
x2 = 2.3128.
3 Itération 3 pour n = 2 on x3 = 1.05 ∗ x2 + 0.05 exp(2t2 ) or x2 = 2.3128 et t2 = 0.1.
Donc x3 = 2.4895.
4 Itération 4 pour n = 3 on x4 = 1.05 ∗ x3 + 0.05 exp(2t3 ) or x3 = 2.4895 et t3 = 0.15.
Donc x4 = 2.6815.
5 Itération 5 pour n = 4 on x5 = 1.05 ∗ x4 + 0.05 exp(2t4 ) or x4 = 2.6815 et t4 = 0.2.
Donc x5 = 2.8902.
6 Itération 6 pour n = 5 on x6 = 1.05 ∗ x5 + 0.05 exp(2t5 ) or x5 = 2.8902 et t5 = 0.25.
Donc x6 = 3.1171.
Question b): L’erreur commise sur x6 , la valeur approchée de x(0.3) est donnée par:
E6 =| x6 − x(0.3) |=| 3.1171 − exp(0.3) − exp(0.6) |= 0.0549.
Question 4
logo.eps
On a E6 < E3 donc lorsque le nombre de points est augmenté c’est-à-dire la valeur de h
diminue la solution de schéma explicite converge vers la solution exacte.
UP Mathématique 25 26 Avril 2020 25 / 26
 MERCI POUR VOTRE ATTENTION

logo.eps

UP Mathématique 26 26 Avril 2020 26 / 26