Résolution numérique d’équations différentielles ordinaires

EDO

Abahchouch Youness

CPGE Agadir
REDA SLAOUI
Département d’informatique
[email protected]

2025
Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Introduction

De nombreux problèmes physiques se modélisent à l’aide des équations

différentielles pour lesquelles on ne dispose pas de solutions analytiques, il
existe alors de nombreuses méthodes numériques pour résoudre d’une façon
approchée ces équations. Ces méthodes numériques sont classées en
différentes catégories : méthodes à un pas, méthodes à multi pas, méthodes
implicites, méthodes explicites ....

1/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Définition 1
Une équation différentielle ordinaire EDO est une équation faisant intervenir une
fonction inconnue d’une seule variable (temps ou espace) avec ses dérivées.

f (t, y(t), y′ (t), ..., y(n) (t)) = 0

L’ordre d’une équation différentielle correspond au degré maximal de dérivation
auquel la fonction inconnue a été soumise.

2/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Définition 2
Un problème de Cauchy est une EDO du premier ordre munie d’une condition
initiale imposée y0 , définie sur un intervalle [t0 , t0 + T ] par :
 ′
y (t) = f (t, y(t))
y(t0 ) = y0

On admettra le théorème de Cauchy-Lipschitz donnant l’existence et

l’unicité de la solution de ce problème sous certaines conditions sur la
fonction f, qu’on supposera vérifiées.

3/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Tout système d’équations différentielles d’ordre plus élevé peut être ramené
à un problème de Cauchy :
Une équation différentielle scalaire d’ordre supérieur se ramène à des
systèmes différentiels couplés d’ordre 1
Un système d’équations différentielles couplées, se ramène à une seule
équation différentielle vectorielle d’ordre

Remarque
Pour une EDO d’ordre n ayant comme inconnue la fonction y(t), le problème
de Cauchy requiert la connaissance au temps t0 de y(t0 ) et de ses n − 1
premières dérivées : y′ (t0 ), ..., yn−1 (t0 )

4/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Exemples : Problème scalaire du premier ordre

Considérons le circuit électrique suivant :

Il est régi par les relations s = e − Ri et i = C ds

dt qui conduisent à
l’équation différentielle : RC ds
dt + s = e.En posant τ = RC on est
amené à résoudre le problème de Cauchy :
 ′
 s (t) = τ1 (e(t) − s(t))
q0
s(0) =

C

5/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Exemples : Problème scalaire du second ordre

Considérons une masse m liée à un point fixe par un fil de
longueur l, et oscillant dans le champ de pesanteur terrestre. En
notant θ l’angle que fait le fil avec la verticale, on obtient
2 g
l’équation du mouvement : ddt2θ + ω 2 sin(θ ) = 0 avec ω 2 = l

6/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Exemples : Problème scalaire du second ordre

Considérons une masse m liée à un point fixe par un fil de
longueur l, et oscillant dans le champ de pesanteur terrestre. En
notant θ l’angle que fait le fil avec la verticale, on obtient
2 g
l’équation du mouvement : ddt2θ + ω 2 sin(θ ) = 0 avec ω 2 = l
En posant y = (y1 , y2 ) = (θ, θ ′ ), notre problème devient
équivalent au problème de Cauchy suivant :

y′ = (θ ′ , θ”) = (θ ′ , −ω 2 sin(θ )) = (y2 , −ω 2 sin(y1 ))

6/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Exemples : Problème vectoriel 2D du premier ordre

Considérons le circuit électrique suivant :

di
Il est régi par les relations s = e − Ri − L dt et i = C ds
dt qui
conduisent au système différentiel
 ′ 1
 i (t) = L (e(t) − s(t) − Ri(t))

s′ (t ) = 1
C i(t)

7/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Exemples : Problème vectoriel 2D du premier ordre

Considérons le circuit électrique suivant :

di
Il est régi par les relations s = e − Ri − L dt et i = C ds
dt qui
conduisent au système différentiel
 ′ 1
 i (t) = L (e(t) − s(t) − Ri(t))
⇒ y′ = (i′ (t), s′ (t)) = ( L1 (e(t) − s(t) − Ri(t)), C1 i(t)) = ( L1 (e(t) − y2 − Ry1 ), C1 y1 )
s′ (t ) = 1
C i(t)


En posant y = (y1 , y2 ) = (i(t), s(t)), notre problème devient

équivalent au problème de Cauchy suivant
8/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Exemples : Problème vectoriel 2D du second ordre

Système masse-ressort en mouvement plan sans frottement

 x”(t) = −k(1 − ll0 )x


p
avec l = x2 + y2
y”(t) = −k(1 − ll0 )y


9/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Exemples : Problème vectoriel 2D du second ordre

Système masse-ressort en mouvement plan sans frottement

 x”(t) = −k(1 − ll0 )x


p
avec l = x2 + y2
y”(t) = −k(1 − ll0 )y


En posant y = (y1 , y2 , y3 , y4 ) = (x(t), y(t), x′ (t), y′ (t)), notre

problème devient équivalent au problème de Cauchy suivant

l0 l0
y′ = (y3 , y4 , −k(1 − q ) y1 , − k ( 1 − q ) y2
y21 + y22 y21 + y22

9/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Principe : L’objectif est d’obtenir une solution approchée au

problème de Cauchy en calculant des valeurs approchées yn
aux valeurs exactes y(tn ) en des instants tn , obtenues par une
discrétisation d’intervalle d’étude.
L’idée consiste à progresser par petits pas, et à calculer à
chaque étape une valeur approchée au pas suivant à partir des
valeurs déjà obtenues.
Ainsi, on appel schéma numérique d’une méthode de
résolution des équations différentielles, les formules de
récurrence associés, qui permettent le calcul de ces valeurs
approchées

10/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Méthodes à un pas

Ces méthodes sont définies d’une façon générale sous la forme récurrente
suivante appelée schéma numérique :

 ∀n ∈ J0, N − 1K yn+1 = yn + hn Φ(tn , yn , hn )



tn+1 = tn + hn


11/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Méthode d’Euler explicite

La méthode d’Euler explicite est une méthode à un pas qui correspond au

schéma suivant :

 ∀n ∈ J0, N − 1K yn+1 = yn + hn f (tn , yn )
tn+1 = tn + hn
y 0 = y ( t0 )


Justification du schéma :
Cette méthode consiste à discrétiser l’intervalle puis approcher
la relation :
Z tn+1 Z tn+1
y ( tn + 1 ) − y ( tn ) = y′ (t)dt = f (t, y(t))dt
tn tn
R tn+1
On approche donc l’intégrale tn
f (t, y(t))dt par la méthode
du rectangle gauche
12/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Méthode d’Euler implicite

La méthode d’Euler implicite est une méthode à un pas qui correspond au

schéma suivant :

 ∀n ∈ J0, N − 1K yn+1 = yn + hn f (tn+1 , yn+1 )
tn+1 = tn + hn
y 0 = y ( t0 )


Justification du schéma :
Cette méthode consiste à discrétiser l’intervalle puis approcher
la relation :
Z tn+1 Z tn+1
y(tn+1 ) − y(tn ) = y′ (t)dt = f (t, y(t))dt
tn tn
R tn+1
On approche donc l’intégrale tn
f (t, y(t))dt par la méthode
du rectangle droit
13/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Ce code utilise la méthode d’Euler explicite pour approximer

dy
numériquement la solution de l’équation différentielle dt = 4y + 3 avec la
condition initiale y(0) = 1, sur l’intervalle [0, 2] divisé en 300 pas. Il calcule
les valeurs approchées de y(t) étape par étape et affiche la courbe
correspondante à l’aide de la bibliothèque matplotlib.

14/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Méthode d’Euler implicite

La méthode d’Euler implicite est une méthode à un pas qui correspond au

schéma suivant :

 ∀n ∈ J0, N − 1K yn+1 = yn + hn f (tn+1 , yn+1 )
tn+1 = tn + hn
y 0 = y ( t0 )


On observe que cette relation ne procure pas une relation explicite de yn+1
puisque ce terme est aussi présent dans le second membre. Pour calculer ce
terme, on peut utiliser le schéma implicite pour avoir une première
approximation de yn+1 du membre droit, ou bien con peut utiliser une
méthode de résolution numérique des équations telle la méthode de
Newton-Raphson(. Pour cette raison, elle se révèle plus coûteuse à mettre en
œuvre.

15/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Méthode d’Euler modifié (point milieu)

La méthode d’Euler modifié est une méthode à un pas qui correspond au

schéma suivant :

h
 ∀n ∈ J0, N − 1K yn+1 = yn + hn f (yn+ 21 , tn + 2n )


 tn + 1 = tn + h n
 y0 = y(t0 )
Justification du schéma :
Rt
La méthode consiste à approcher l’intégrale tnn+1 f (t, y(t))dt par
la méthode du point milieu, c-à-dire :

Z tn+1
y(tn+1 ) + y(tn ) tn+1 + tn
f (t, y(t))dt = (tn+1 − tn )f ( , )
tn 2 2

16/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Méthode d’Euler modifié (point milieu)

La méthode d’Euler modifié est une méthode à un pas qui correspond au

schéma suivant :

h
 ∀n ∈ J0, N − 1K yn+1 = yn + hn f (yn+ 21 , tn + 2n )


 tn + 1 = tn + h n
 y0 = y(t0 )
Cependant, comme on ne connaît pas yn+ 1 , on peut donc utiliser la méthode
2
d’Euler explicite pour l’approximer :

hn
yn+ 1 = yn + f (tn , yn )
2 2

17/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Méthode de Heun

Rt
La méthode consiste à approcher l’intégrale tnn+1 f (t, y(t))dt par la méthode
des trapèzes, ce qui conduit au schéma numérique implicite suivant :

 ∀n ∈ J0, N − 1K yn+1 = yn + h2n (f (yn , tn ) + f (yn+1 , tn+1 ))



tn+1 = tn + hn
y 0 = y ( t0 )


On utilise donc le schéma d’Euler explicite pour approcher le terme yn+1 du

membre droit.

18/19

Résolution numérique d’équations différentielles ordinaires EDO

Introduction Définitions et Notations Systèmes différentiels et équations d’ordre supérieur Résolution numérique des EDO

Méthode de Runge-Kutta 4

Rt
La méthode consiste à approcher l’intégrale tnn+1 f (t, y(t))dt par la méthode
de
 simpson, ce qui conduit au schéma numérique implicite suivant :
h h h
 yn+1 = yn + 6n (f (yn , tn ) + 2f (y1n+ 12 , tn + 2n ) + 2f (y2n+ 21 , tn + 2n ) + f (yn+1 , tn+1 ))


 tn+1 = tn + hn
 y = y(t )
0 0

On utilise donc le schéma d’Euler explicite pour approcher le terme y1n+ 1 ,

2
puis le schéma d’Euler implicite pour approcher le terme y2n+ 1 , enfin le
2
schéma du point milieu pour approcher le terme yn+1 du membre droit.

19/19

Résolution numérique d’équations différentielles ordinaires EDO