Mathématiques pour l’ingénieur

R. Harrimann

Février 2025
2

0.1 Introduction
Un système d’équations différentielles ordinaires (EDO) est un ensemble d’équations de la
forme :
dy
= f (t, y) (1)
dt
ou y = (y1 , y2 , . . . , yn ) est un vecteur de fonctions inconnues.

0.2 Méthodes de Résolution

0.2.1 Méthode de la Diagonalisation
Si le système peut s’écrire sous la forme matricielle :

dy
= Ay (2)
dt
où A est une matrice constante, on cherche une solution sous la forme y = eAt c.

0.2.2 Méthode de Variation des Paramètres

dy
Si le système est de la forme dt = Ay + b(t), on peut chercher une solution particulière via
la variation des paramètres.

0.2.3 Méthode de Runge-Kutta

Pour les systèmes non linéaires, des méthodes numériques comme Runge-Kutta permettent
d’obtenir des approximations des solutions.

0.2.4 Méthode d’Euler

Une méthode simple pour obtenir une approximation numérique est la méthode d’Euler :

yn+1 = yn + hf (tn , yn ) (3)

où h est le pas de discrétisation.

0.2.5 Méthode de Laplace

La transformée de Laplace permet de convertir un système d’EDO en un système algébrique,
plus simple à résoudre.

0.3 Exemples et Applications

0.3.1 Système Linéaire Homogène
Considérons le système :     
    
y1  2 −1 y1 
d   
 =
 
  (4)
dt   
    
 
    
y2 1 3 y2

On trouve les valeurs propres et vecteurs propres pour résoudre ce système.

0.4. CONCLUSION 3

0.3.2 Système Non Linéaire

Un système comme :
dy1 dy2
= y1 y2 , = y1 − y2 (5)
dt dt
peut être résolu numériquement.

0.3.3 Application en Physique

Les EDO apparaissent en physique dans la dynamique des systèmes, comme l’étude du
mouvement des planètes ou des circuits èlectriques.

0.3.4 Application en Biologie

Les modèles de croissance des populations, comme l’équation de Lotka-Volterra, sont formulés
avec des systèmes d’EDO.

0.3.5 Application en Électricité et Électronique

Les circuits èlectriques contenant des composants RLC obéissent à des EDO. Par exemple,
un circuit RLC en série suit l’équation :

d2 i di i
L 2
+R + = V (t) (6)
dt dt C

où L, R, et C sont respectivement l’inductance, la résistance et la capacitance.

0.3.6 Application en Mécanique

En mécanique, les oscillateurs harmoniques sont modélisés par des EDO. L’équation du
mouvement d’un ressort amorti est donnée par :

d2 x dx
m 2
+b + kx = F (t) (7)
dt dt

où m est la masse, b le coefficient d’amortissement, et k la constante de rappel du ressort.

0.4 Conclusion
La résolution des systèmes d’EDO repose sur des méthodes analytiques et numériques.
Chaque cas demande une approche adaptée en fonction de sa structure.

Résolution d’une équation différentielle avec la transformée de La-

place
Considérons l’équation différentielle linéaire du premier ordre avec une condition initiale :

dy(t)
+ 2y(t) = 3e−t , y(0) = 1
dt
4

Étape 1 : Appliquer la transformée de Laplace

On applique la transformée de Laplace à chaque terme de l’équation. Rappelons que la
transformée de Laplace de y(t) est notée Y (s), et celle de dy(t)
dt est sY (s) − y(0).
 
dy(t)
L + 2L{y(t)} = L{3e−t }
dt
Cela donne :
3
sY (s) − y(0) + 2Y (s) =
s+1

Étape 2 : Substituer la condition initiale

On remplace y(0) = 1 dans l’équation :
3
sY (s) − 1 + 2Y (s) =
s+1

Étape 3 : Résoudre pour Y (s)

On regroupe les termes en Y (s) :
3
Y (s)(s + 2) = +1
s+1
3 + (s + 1) s+4
Y (s)(s + 2) = =
s+1 s+1
On isole Y (s) :
s+4
Y (s) =
(s + 1)(s + 2)

Étape 4 : Décomposer en fractions partielles

On décompose Y (s) en fractions partielles pour faciliter la transformation inverse :

s+4 A B
= +
(s + 1)(s + 2) s+1 s+2
En résolvant pour A et B, on trouve :

A = 3, B = −2
Donc :
3 2
Y (s) = −
s+1 s+2

Étape 5 : Appliquer la transformée de Laplace inverse

On applique la transformée de Laplace inverse à chaque terme pour retrouver y(t) :
   
−1 3 −1 2
y(t) = L −L
s+1 s+2

y(t) = 3e−t − 2e−2t

0.4. CONCLUSION 5

Solution finale
La solution de l’équation différentielle est :

y(t) = 3e−t − 2e−2t

Vérification
On peut vérifier que cette solution satisfait à la fois l’équation différentielle et la condition
initiale y(0) = 1. La transformée de Laplace d’une fonction f (t) est définie par :
Z ∞
L{f (t)} = F (s) = e−st f (t)dt
0
où s est un paramètre complexe tel que l’intégrale converge.
6

0.5 Transformées de Laplace des fonctions usuelles

Fonction f (t) Transformée de Laplace F (s) = L{f (t)} Domaine de convergence

1
1 R(s) > 0
s

n!
tn (n ∈ N) R(s) > 0
sn+1

1
eat R(s) > a
s−a

ω
sin(ωt) R(s) > 0
s2 + ω2

s
cos(ωt) R(s) > 0
s2 + ω2

n!
tn eat (n ∈ N) R(s) > a
(s − a)n+1

ω
sinh(ωt) R(s) > |ω|
s2 − ω2

s
cosh(ωt) R(s) > |ω|
s2 − ω2

2ωs
t sin(ωt) R(s) > 0
(s2 + ω 2 )2

s2 − ω 2
t cos(ωt) R(s) > 0
(s2 + ω 2 )2

ω
eat sin(ωt) R(s) > a
(s − a)2 + ω 2

s−a
eat cos(ωt) R(s) > a
(s − a)2 + ω 2

δ(t) (Dirac) 1 ∀s

n!
tn u(t) (n ∈ N) R(s) > 0
sn+1

1 1
√ √ R(s) > 0
πt s
 
1 at s−a
(e − ebt ) ln R(s) > max(a, b)
t s−b
0.5. TRANSFORMÉES DE LAPLACE DES FONCTIONS USUELLES 7

Démonstration
1. Transformée de Laplace de f (t) = 1

Z ∞
L{f (t)} = f (t)e−st dt
Z0 ∞ Z ∞
−st
= 1×e dt = e−st dt
0 0

En intégrant e−st par rapport à t,

∞
e−st

0 − (−1) 1
L{1} = = = , s>0
−s 0 s s

2. Transformée de Laplace de f (t) = tn

Avant de chercher la transformée de f (t) = tn , faisons un petit rappel sur la fonction Γ

Définition
Soit z un élément de C, la fonction gamma est une fonction spéciale définie par :
Z +∞
∀ R(z) > 0, Γ(z) = tz−1 e−t dt.
0

Quelques propriétés de la fonction Γ

Γ(z + 1) = z · Γ(z)
Z ∞ Z ∞
En effet,par définition, Γ(z + 1) = tz+1−1 e−t dt = tz e−t dt.
0 0

Calculons ce Γ(z + 1)
Z Z
Utilisons l’intégration par parties de la forme : u dv = uv − v du.

Choisissons :
u = tz ⇒ du = ztz−1 dt,
dv = e−t dt ⇒ v = −e−t .
Appliquons l’intégration par parties :
Z ∞ Z ∞
z −t
 z −t ∞
Γ(z + 1) = t e dt = −t e 0 + ztz−1 e−t dt.
0 0
∞
Il est facile de constater que notre premier terme s’annule aux bornes, i-e : −tz e−t 0 = 0.


L’intégrale restante est donc :

Z ∞ Z ∞
z−1 −t
zt e dt = z tz−1 e−t dt = z · Γ(z).
0 0
C’est ainsi que :
Γ(z + 1) = z · Γ(z) .
La relation Γ(z + 1) = z · Γ(z) est fondamentale pour la fonction Gamma. Elle permet de
généraliser la factorielle aux nombres complexes. En effet, pour n ∈ N, on a :
Γ(n + 1) = n!.
8

Démonstration par récurrence

En appliquant la relation Γ(z + 1) = z · Γ(z) de manière répétée, on obtient :

Γ(n + 1) = n · Γ(n) = n · (n − 1) · Γ(n − 1) = · · · = n! · Γ(1).

Or, Γ(1) se calcule directement à partir de la définition de la fonction Gamma :

Z ∞ ∞
e−t dt = −e−t 0 = 1.

Γ(1) =
0

Par conséquent, on a :
Γ(n + 1) = n! · Γ(1) = n! · 1 = n!.

Exemple
Par exemple, pour calculer Γ(5), on prend n = 4 :

Γ(5) = Γ(4 + 1) = 4 · Γ(4) = 4 · 3 · Γ(3) = 4 · 3 · 2 · Γ(2) = 4 · 3 · 2 · 1 · Γ(1) = 4! · 1 = 24.

Ce qui correspond bien à :

4! = 24

.
Revenons maintenant à la recherche de L{tn }
Z ∞
L{f (t)} = L{tn } = tn e−st dt
Z0 ∞
= tn+1−1 e−st dt
0

dT
Effectuons un changement de variable en posant T = st, soit = s ou encore 1s dT = dt
dt
Z ∞
L{f (t)} = L{tn } = tn+1−1 e−st dt
0
Z ∞  n+1−1
T 1
= e−T dT
0 s s
Z ∞
1 1
= · n T n+1−1 e−T dT
s s 0

En utilisant l’intégrale gamma :

Z ∞
1
tn e−st dt = Γ(n + 1), or Γ(n + 1) = n!
0 sn+1

On obtient donc :

n!
L{tn } = , s>0
sn+1
0.5. TRANSFORMÉES DE LAPLACE DES FONCTIONS USUELLES 9

3. Transformée de Laplace de eat

Z ∞
L{eat } = e(a−s)t dt
0
L’intégrale donne :
" #∞
e(a−s)t
L{eat } =
a−s
0
1
= , s>a
s−a

4. Transformée de Laplace de cos(ωt)

Z ∞
L{cos(ωt)} = e−st cos(ωt)dt
0
eiωt +e−iωt
En utilisant la formule d’Euler cos(ωt) = 2 , on obtient :
s
L{cos(ωt)} = , s>0
s2 + ω2

5. Transformée de Laplace de sin(ωt)

De manière similaire,
Z ∞
L{sin(ωt)} = e−st sin(ωt)dt
0
ω
= , s>0
s2 + ω 2

6. Transformée de Laplace de l’impulsion de Dirac δ(t)

Z ∞
L{δ(t)} = e−st δ(t)dt
0
Par la propriété de l’impulsion de Dirac :

L{δ(t)} = 1

Correction des exercices sur la transformée de Laplace

Exercice 1 : Transformée de Laplace de f (t) = e3t
La transformée de Laplace de f (t) = eat est donnée par :
1
L{eat } = , pour s > a.
s−a
Pour f (t) = e3t , on a a = 3, donc :
1
L{e3t } = , pour s > 3.
s−3
10

Exercice 2 : Transformée de Laplace de f (t) = t2 + 3t + 2

La transformée de Laplace est linéaire, donc :

L{t2 + 3t + 2} = L{t2 } + 3L{t} + 2L{1}.

On utilise les formules suivantes :

n! 1
L{tn } = , L{1} = .
sn+1 s
Ainsi :
2! 2 1 1
L{t2 } = = 3, L{t} = , L{1} = .
s3 s s2 s
En combinant ces résultats :
2 3 2
L{t2 + 3t + 2} = + + .
s3 s2 s

Exercice 3 : Transformée de Laplace de f (t) = sin(2t)

La transformée de Laplace de sin(ωt) est donnée par :
ω
L{sin(ωt)} = .
s2 + ω 2
Pour f (t) = sin(2t), on a ω = 2, donc :
2
L{sin(2t)} = .
s2 + 4

Exercice 4 : Transformée de Laplace d’une fonction définie par morceaux

La fonction est définie par :
(
1 si 0 ≤ t < 2,
f (t) =
t si t ≥ 2.

On peut écrire f (t) comme :

f (t) = 1 · 1[0,2) (t) + t · 1[2,∞) (t).

La transformée de Laplace est alors :

Z 2 Z ∞
−st
L{f (t)} = e · 1 dt + e−st · t dt.
0 2

Calculons chaque intégrale séparément :

2
2
−e−st 1 − e−2s
Z 
−st
e dt = = ,
0 s 0 s
et ∞
∞ ∞
−te−st e−st 2e−2s e−2s
Z  Z
−st
te dt = + dt = + 2 .
2 s 2 2 s s s
En combinant ces résultats :
1 − e−2s 2e−2s e−2s
L{f (t)} = + + 2 .
s s s
0.5. TRANSFORMÉES DE LAPLACE DES FONCTIONS USUELLES 11

2s+3
Exercice 5 : Transformée de Laplace inverse de F (s) = s2 +4s+13

On commence par compléter le carré au dénominateur :

s2 + 4s + 13 = (s + 2)2 + 9.

Ainsi :
2s + 3
F (s) = .
(s + 2)2 + 9
On décompose F (s) en deux termes :
s+2 1
F (s) = 2 · − .
(s + 2) + 9 (s + 2)2 + 9
2

En utilisant les propriétés de la transformée de Laplace inverse :

 
−1 s+a
L = e−at cos(ωt),
(s + a)2 + ω 2
et  
−1 ω
L = e−at sin(ωt).
(s + a)2 + ω 2
On obtient :
1
f (t) = 2e−2t cos(3t) − e−2t sin(3t).
3
12

y
Résolution de l’équation différentielle y 0 = ey/x + x

Étape 1 : Substitution
Posons u = xy . Alors, y = ux, et en dérivant par rapport à x, on obtient :

y 0 = u + xu0 .
y
Substituons y 0 et x = u dans l’équation différentielle :

u + xu0 = eu + u.

Simplifions cette équation :

xu0 = eu .

Étape 2 : Séparation des variables

L’équation xu0 = eu peut être réécrite comme :
du dx
u
= .
e x
Cela donne :
1
e−u du = dx.
x

Étape 3 : Intégration
Intégrons les deux côtés : Z Z
−u 1
e du = dx.
x
L’intégration donne :
−e−u = ln |x| + C,
où C est la constante d’intégration.

Étape 4 : Résolution pour u

Isolons e−u :
e−u = − ln |x| − C.
Prenons le logarithme naturel des deux côtés :

−u = ln(− ln |x| − C).

Ainsi,
u = − ln(− ln |x| − C).

Étape 5 : Revenir à y
Rappelons que u = xy . Donc,
y
= − ln(− ln |x| − C).
x
Multiplions par x pour obtenir y :

y = −x ln(− ln |x| − C).

0.5. TRANSFORMÉES DE LAPLACE DES FONCTIONS USUELLES 13

Étape 6 : Solution finale

La solution analytique de l’équation différentielle est :

y(x) = −x ln(− ln |x| − C) ,

où C est une constante d’intégration. Cette solution est valide pour ln |x| + C < 0, c’est-à-dire
pour |x| < e−C .
14

Exercices d’équations différentielles du second ordre

Exercice 1
Équation différentielle :
y 00 − 4y 0 + 4y = e2x .

Indications :
1. Résolvez l’équation homogène associée.
2. Utilisez la méthode de variation des constantes pour trouver une solution particulière.
3. Donnez la solution générale.
—

Exercice 2
Équation différentielle :
y 00 + y = tan(x).

Indications :
1. Résolvez l’équation homogène associée.
2. Utilisez la méthode de variation des constantes pour trouver une solution particulière.
3. Donnez la solution générale.
—

Exercice 3
Équation différentielle :
ex
y 00 − 2y 0 + y = .
x
Indications :
1. Résolvez l’équation homogène associée.
2. Utilisez la méthode de variation des constantes pour trouver une solution particulière.
3. Donnez la solution générale.
—

Exercice 4
Équation différentielle :
y 00 + 4y = sec(2x).

Indications :
1. Résolvez l’équation homogène associée.
2. Utilisez la méthode de variation des constantes pour trouver une solution particulière.
3. Donnez la solution générale.
—
0.5. TRANSFORMÉES DE LAPLACE DES FONCTIONS USUELLES 15

Exercice 5
Équation différentielle :
ex
y 00 − 3y 0 + 2y = .
1 + ex
Indications :
1. Résolvez l’équation homogène associée.
2. Utilisez la méthode de variation des constantes pour trouver une solution particulière.
3. Donnez la solution générale.
—

Méthode de variation des constantes (rappel)

Pour une équation différentielle du second ordre de la forme :

y 00 + p(x)y 0 + q(x)y = f (x),

la solution particulière yp (x) est donnée par :

yp (x) = u1 (x)y1 (x) + u2 (x)y2 (x),

où y1 (x) et y2 (x) sont les solutions de l’équation homogène, et u1 (x), u2 (x) sont déterminées en
résolvant le système : (
u01 (x)y1 (x) + u02 (x)y2 (x) = 0,
u01 (x)y10 (x) + u02 (x)y20 (x) = f (x).