EQUATIONS DIFFERENTIELLES :

RESOLUTION THEORIQUE ET NUMERIQUE

Soyibou SY
Enseignant en Mathématiques

EPF École d’Ingénieur-e-s de Cahan (France)

Yaoundé, du 12 au 16 février 2023

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France)) 1 / 28THEORIQUE ET NUMERIQUE
 Plan de la présentation

1 Partie I: Résultats classiques sur les EDO

Motivations
Equations différentielles d’ordre 1
À variables séparables
Linéaires à coefficients constants
Linéaires à coefficients non constants
Equations différentielles d’ordre 2
Linéaires homogènes à coefficients constants
Linéaires non homogènes à coefficients constants
2 Partie II: Approximations numériques
Méthodes d’Euler
Méthodes de Runge Kutta d’ordre 2 et 4
3 Partie III: Applications à la physique
Utilisation du logiciel Scilab ou Matlab
4 Conclusion
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France)) 2 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 1

Motivations
1 Pour un circuit RL, avec la loi des mailles :
di R u
+ i= : i ′ (t) + a i(t) = f (t).
dt L L
2 Chute d’un corps avec une accélération g sans résistance de l’air :
d 2x
= −g .
dt 2
Cela permet d’avoir
1
v = −gt + v0 , et x = − gt 2 + v0 t + x0 .
2
3 Pendule simple :
d 2θ
= −g sin(θ)
dt 2
d 2θ x
Si θ est suffisamment petit : = −g θ = −g .
dt 2 l
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France)) 3 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 1

Équations à variables séparées

On considère l’équation différentielle d’ordre 1 :

f (y (x))y ′ (x) = g (x),

On note F : primitive de f et G : primitive de g . Alors :

y (x) = F −1 (G (x) + C ), C est une constante.

Astuce :
dy
y′ = .
dx

Example (Exercice)
 ′
T (t) = K (T (t) − 25)
Loi de Newton :
T (0) = 75
1 Solution : T (t) = 25 + 50e Kt .
ln(2)
2 Déterminer K si T (5) = 50. T (t) = 25 + 50e − 5 .
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France)) 4 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 1

Exercices séparation de variables

1 Résoudre l’équation différentielle sur R :
2
e −x +y ′
y = x.
2 Résoudre l’équation différentielle R :
y
y′ = .
1 + t2
3 Montrer que
(x 2 + 1)y ′ (x) − xy (x) = 0
a pour solution p
y (x) = C x 2 + 1, C ∈ R.
4 Montrer que *
2
e −x +y ′
y (x) = x
a pour solution
1 2

y (x) = ln e x + C , C ∈ R.
2
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France)) 5 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 1

Équations de la forme : ay ′ (t) + by (t) = f (t).

Les équations différentielles de la forme :

ay ′ (t) + by (t) = f (t), (a, b) ∈ R∗ × R.

Équations différentielles linéaires d’ordre 1 non homogènes à coefficients constants.

Théorème (Méthode de résolution)

1 Résoudre l’équation homogène : ay ′ (t) + by (t) = 0 :

b
yH (t) = Ke rt , r =− , K ∈ R∗ .
a

2 On trouve une solution particulière yP de ay ′ (t) + by (t) = f (t).

3 On en déduit la solution générale sous la forme :

yG = yH + yP .
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France)) 6 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 1

Méthode de calcul de yP
1 Suivant la forme de f (t) :
Proposition
1 f (t) = P(t) −→ yP = P(t) : en particulier P = Cste.
deg (P) = deg (Q) si 0 n’est pas racine de ar + b = 0.
deg (P) = deg (Q) + 1 si 0 est racine simple de ar + b = 0.

2 f (t) = ke λt on cherche yP sous la forme :

yP (t) = ae λt si λ n’est pas racine de ar + b = 0.
yP (t) = ate λt si λ est racine simple de ar + b = 0.

3 f (t) = k cos(wt) ou f (t) = k sin(wt)

yP = a sin(wt) + b cos(wt). avec a(iw ) + b ̸= 0.

2 Méthode de variation de constantes.
2 1 2
y ′ (t) − 2y (t) = (t − 1)e t −→ yG = e t + Ke 2t
2 Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France)) 7 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 1

Exercices d’applications

Résoudre les équations différentielles suivantes :

5 2
1 3y ′ (t) + 2y (t) = 5 −→ yG = + Ce − 3 t
2
5 13 2
2 3y ′ (t) + 2y (t) = 5t + 1. −→ yG = t − + Ce − 3 t
2 4
5 2
3 3y ′ (t) + 2y (t) = 5e 2t . −→ yG = e 2t + Ce − 3 t
8
1 2
4 y ′ (t) + y (t) = cos(2t). −→ yG = cos(2t) + sin(2t) + Ce −t
5 5
5 Principe de superposition : ay ′ + by (t) = g1 (t) + g2 (t).

Example
1 2
y ′ (t) + y (t) = 5 + cos(2t). −→ yG = 5 + cos(2t) + sin(2t) + Ce −t
5 5

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France)) 8 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 1

Équations différentielles linéaires d’ordre 1 à coefficients

non constants
On considère l’équation différentielle :

a(t)y ′ (t) + b(t)y (t) = f (t), a(t) ̸= 0. ←→ y ′ (t) + c(t)y (t) = g (t).

1 On résout y ′ (t) + c(t)y (t) = 0 −→ yH (t) = Ke −C (t) , avec C ′ = c.

2 On applique la méthode de variation de constante −→ yG .

Example
1 Résoudre l’équation différentielle :

cos(x)y ′ (x) + sin(x)y (x) = 0.

2 En déduire que

cos(x)y ′ (x) + sin(x)y (x) = 1, y (0) = 1 : −→ yG = sin(x) + cos(x).

Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France)) 9 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 2

Équations différentielles : ay ′′ (t) + by ′ (t) + cy (t) = 0

On considère l’équation différentielle :
ay ′′ (t) + by ′ (t) + cy (t) = 0, a, b ∈ R.
Équation caractéristique : ar 2 + br + c = 0. On pose : ∆ = b 2 − ac.
Théorème
√ √
−b − ∆
r1 t r2 t −b + ∆
1 Si ∆ > 0 −→ yH (t) = K1 e + K2 e , r1 = et r2 = .
2a 2a
−b
2 Si ∆ = 0 −→ yH (t) = (K1 t + K2 )e r0 t , r0 = .
2a
3 Si ∆ < 0 −→ p
αt

b |∆|
yH (t) = e K1 cos(βt) + K2 sin(βt) , α = − et β = .
2a 2a

Example
Résoudre y ′′ (t) − 4y ′ (t) + 3y (t) = 0 −→: yH (t) = Ke t + K2 e 3t .
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))10 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 2

Équations différentielles : ay ′′ (t) + by ′ (t) + cy (t) = 0

Example
1 4y ′′ (t) + 4y ′ (t) + y (t) = 0 −→

yH (t) = (K1 t + K2 )e −0·5t

2 y ′′ (t) + 6y ′ (t) + 13y (t) = 0 −→

yH (t) = e −3t A cos(2t) + B sin(2t) .

3 2y ′′ (t) − 3y ′ (t) + y (t) = 0 −→

yH (t) = K1 e 0·5t + K2 e t , K1 , K2 ∈ R.

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))11 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 2

Équations différentielles : ay ′′ (t) + by ′ (t) + cy (t) = f (t).

Théorème
Si f (t) = P(t)e αt et EC : ar 2 + br + c = 0.
1 Si α n’est pas racine de EC : yp (t) = Q(t)e αt , deg (Q) = deg (P).

2 Si α racine de EC : yp (t) = Q(t)e αt , deg (Q) = deg (P) + 1.

3 Si α racine double de EC : yp (t) = Q(t)e αt , deg (Q) = deg (P) + 2.

Remarque
On retrouve les résultats des points 1&2 du théorème sur les équations
différentielles d’ordre 1 en posant α = 0.

Example
1 −x
1 y ′′ + 5y ′ − 6y = e −x (P = 1, α = −1) −→ yG = − e + Ae −6x + Be x .
10
2 y ′′ +y ′ −2y = 2x +1 (P = 2x + 1, α = 0) −→ yG = −x − 1 + Ae −2x + Be x .
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))12 / 28THEORIQUE ET NUMERIQUE
 Équations différentielles d’ordre 2

Équations différentielles avec conditions initiales

On considère l’équation différentielle ordinaire avec conditions initiales :

 ′′
 y + 4y ′ + 13y = 0

y (0) = 1, y ′ (0) = 1


1 r 2 + 4r + 13 = 0 −→ ∆ = (6i)2 et α = −2, β = 3.
2

yH (t) = e −2t (A cos(3t) + B sin(3t)).

3 y (0) = 1, y ′ (0) = 1 donc A = B = 1, ainsi :
 
yG (t) = e −2t cos(3t) + sin(3t) .

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))13 / 28THEORIQUE ET NUMERIQUE
 Schéma d’Euler

Présentation
On considère le problème de Cauchy (EDO avec cdts initiales) suivant :
 ′
 y (t) = f (t, y (t)), t ∈ [0, T ]
(1)
y (0) = y0


On subdivise l’intervalle [0, T ] en : 0 = t0 < t1 < · · · < tn = T , avec

T
, ti = t0 + ih, i = 1, · · · , N.
h=
N
En intégrant 1) entre tn et tn+1 :
Z tn+1
y (tn+1 ) − y (tn ) = f (t, y (t)) dt. (2)
tn

On approche l’intégrale dans (2) en utilisant les méthodes :

Rectangle à gauche et rectangle à droite : méthodes d’Euler.
Trapèze : méthode de Heun ou Runge Kutta 2.
Simpson : Runge Kutta 4.
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))14 / 28THEORIQUE ET NUMERIQUE
 Méthodes d’Euler explicite et implicite

Rectangles gauche et droite

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))15 / 28THEORIQUE ET NUMERIQUE
 Méthodes d’Euler explicite et implicite

Méthodes d’Euler explicite

Euler explicite : rectangle à gauche :

Z tn+1
f (t, y (t)) dt ≃ h × f (tn , y (tn )), où h = tn+1 − tn .
tn

Schéma d’Euler explicite :


 yn+1 = yn + hf (tn , yn )

y0 = y (0)


1 Sous-estime l’aire quand f croit au cours du temps.

2 Méthode d’ordre 1.
3 Facile à implémenter, problème de stabilité (choisir le pas h très petit).

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))16 / 28THEORIQUE ET NUMERIQUE
 Méthodes d’Euler explicite et implicite

Méthodes d’Euler implicite

Euler implicite : rectangle à droite :
Z tn+1
f (t, y (t)) dt ≃ h × f (tn+1 , yn+1 ), où h = tn+1 − tn .
tn

Schéma d’Euler implicite :


 yn+1 = yn + hf (tn+1 , yn+1 )

y0 = y (0)


1 Sur-estime l’aire quand f croit au cours du temps.

2 Méthode d’ordre 1.
3 yn+1 est définie de manière implicite −→ résolution d’une équation implicite
(Newton, dichotomie).
4 Schéma très stable (inconvénient compensé).
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))17 / 28THEORIQUE ET NUMERIQUE
 Schémas de Heun et Runge Kutta d’ordre 2

Méthode des trapezes

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))18 / 28THEORIQUE ET NUMERIQUE
 Schémas de Heun et Runge Kutta d’ordre 2

Méthodes de Runge Kutta d’ordre 2 : Heun

Méthode du trapèze :
b
b − ah
Z i
f (x)dx = f (a) + f (b) .
a 2
Z tn+1
hh i
f (t, y (t)) dt ≃ f (tn , yn ) + f (tn+1 , yn+1 ) , où h = tn+1 − tn .
tn 2
Schéma de Heun :

h   k1 = f (tn , yn )
yn+1 = yn + k1 + k2 , avec k2 = f (tn + h, yn + hk1 ) (3)
2
y0 = y (0)


1 Neutralisation de l’erreur d’estimation de l’aire.

2 Méthode d’ordre 2. Assez simple à implémenter.
3 Schéma stable.
4 Algorithme numérique : appliquer les étapes de (3).
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))19 / 28THEORIQUE ET NUMERIQUE
 Schémas de Runge Kutta d’ordre 4

Méthode de Simpson

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))20 / 28THEORIQUE ET NUMERIQUE
 Schémas de Runge Kutta d’ordre 4

Méthodes de Runge Kutta d’ordre 4

b
b − ah
Z a + b i
Méthode du Simpson : f (x)dx = f (a) + 4f + f (b) .
a 6 2
Z tn+1
hh i
f (t, y (t)) dt ≃ f (tn , yn ) + 4f (tn+ 21 , y (tn+ 12 )) + f (tn+1 , yn+1 ) .
tn 6

Schéma de Runge Kutta d’ordre 4 :



 k1 = f (t
 n , yn ) 
hk1
 h
k = f t + , y +

2 n n

 2 2
h    
yn+1 = yn + k1 + 2k2 + 2k3 + k4 , avec k 3 = f t n + h
, y n + hk2 (4)
6 
 2 2
k = f (t + h, y + hk 3)

4 n n



y0 = y (0)

1 Étape supplémentaire dans le raffinement du calcul de l’intégrale.

2 Méthode d’ordre 4. Méthode très précise.

3 Deux fois plus de calcul que la méthode de Heun. Schéma stable.
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))21 / 28THEORIQUE ET NUMERIQUE
 Applications avec Scilab

Exemple d’équation différentielle

On cherche à comparer par la méthode d’Euler explicite la solution exacte et la

solution numérique de l’équation différentielle :
 ′
y (t) = −2ty (t), ∀ t ∈ [0; 2]
y (0) = 1

dont la solution exacte est donnée par la fonction :

2
t 7→ e −t .

Fonctions préliminaires : codes scilab

function g = fonction(t, y ) function y = myfonction(t)

g = −2 ∗ t. ∗ y ; y = exp(−t.2 );
end end

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))22 / 28THEORIQUE ET NUMERIQUE
 Applications avec Scilab

Code Scilab
function y=MethodeEuler(y0,t0,T,h)
Intervalles et vecteurs
t=[t0 : h : T];
N=length(t);
y=zeros(N,1);
sol=zeros(N,1);
y(1)=y0;
Boucle for pour résolution
for k=1:N-1
y(k+1)=y(k)+h*fonction(t(k),y(k));
end
Affichage des solutions
plot(t, myfonction(t), ’b-’, t,y,’r-’)
legend(’solution exacte’,’solution numérique’);
title ’Comparaison de solutions pour h=0.01’;
hold on
end
Yaoundé, du 12 au 16 février 2023
Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))23 / 28THEORIQUE ET NUMERIQUE
 Applications avec Scilab

Illustration graphique Euler h=0.1

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))24 / 28THEORIQUE ET NUMERIQUE
 Applications avec Scilab

Illustration graphique Euler h=0.01

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))25 / 28THEORIQUE ET NUMERIQUE
 Applications avec Scilab

Exercices d’application
1 Écrire à l’aide Scilab les fonctions : ”fonction” (pour le second membre),
”myfonction” (pour la solution exacte) et ”fonctionEuler” (pour le schéma
d’Euler explicite) pour résoudre numérique l’équation différentielle :
 ′
y (t) = 3y (t) − 4e −t ∀ t ∈ [0; 1]
y (0) = 2

On montrera au passage que la solution exacte est donnée par :

yexacte (t) = e 3t + e −t .

2 Écrire un programme Scilab utilisant le schéma d’Euler explicite pour

résoudre numérique l’équation différentielle :

 y ′ (t) = y (t) − 2t , ∀ t ∈ [0; 1]
y (t)
 y (0) = 1

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))26 / 28THEORIQUE ET NUMERIQUE
 Conclusion

Conclusion

1 Résultats classiques sur les équations différentielles d’ordre 1 et 2.

2 Présentation des méthodes numériques pour résoudre les équations
différentielles.
3 Connaissance de la solution analytique. Comparaison des méthodes
numériques.
4 Sans la solution analytique. Approximation de la solution de l’équation.
5 Etude de la stabilité. Meilleur schéma.

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))27 / 28THEORIQUE ET NUMERIQUE
 Merci pour votre attention !

Yaoundé, du 12 au 16 février 2023

Soyibou SY Enseignant en Mathématiques (EPF
EQUATIONS
École d’Ingénieur-e-s
DIFFERENTIELLES
de Cahan
: RESOLUTION
(France))28 / 28THEORIQUE ET NUMERIQUE