RÉPUBLIQUE DÉMOCRATIQUE DU CONGO

MINISTÈRE DE L’ENSEIGNEMENT SUPÉRIEUR ET UNIVERSITAIRE

Université Loyola du Congo

Faculté d’Ingénierie ULC-Icam
Département de Maintenance et Génie Industriel

Mission catholique KIMWENZA

Commune de Mont-Ngafula
BP 3724 KINSHASA – GOMBE

Cours de math

Travail pratique du 7 décembre 2025

Rédigé par : ISANGANINO M’LEMA Moise

Étudiant en A3 Génie industriel

Dirigé par : Dr. Demco MUKONGO

Année académique : 2025 – 2026

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

y ′′′ + 2y ′′ − y ′ − 2y = 0.

1. Transformation en système d’ordre 1

On pose :
y1 = y, y2 = y ′ , y3 = y ′′ .
Alors :
y1′ = y2 , y2′ = y3 , y3′ = −2y3 + y2 + 2y1 .
Le système s’écrit sous la forme matricielle :
 ′   
y1 0 1 0 y1
y2′  = 0 0 1  y2  .
y3′ 2 1 −2 y3

On note :    
y1 0 1 0
Y = y2  ,
 Y ′ = AY, A = 0 0 1  .
y3 2 1 −2

2. Valeurs propres du système

On calcule :
−λ 1 0
det(A − λI) = 0 −λ 1 .
2 1 −2 − λ
Développement suivant la première ligne :

−λ 1 0 1
det(A − λI) = (−λ) −1 .
1 −2 − λ 2 −2 − λ

On calcule les mineurs :

−λ 1 0 1
= λ2 + 2λ − 1, = −2.
1 −2 − λ 2 −2 − λ

D’où :
det(A − λI) = −λ(λ2 + 2λ − 1) + 2 = −(λ3 + 2λ2 − λ − 2).
L’équation caractéristique est donc :

λ3 + 2λ2 − λ − 2 = 0.

On factorise :
(λ + 2)(λ + 1)(λ − 1) = 0.
Ainsi, les valeurs propres sont :

λ1 = 1, λ2 = −1, λ3 = −2.

1
3. Vecteurs propres
Pour λ1 = 1 :  
−1 1 0
A − I =  0 −1 1  .
2 1 −3
Un vecteur propre associé est :  
1
v1 = 1 .

1

Pour λ2 = −1 :  
1 1 0
A + I = 0 1 1  .
2 1 −1
Un vecteur propre associé est :  
1
v2 = −1 .
1

Pour λ3 = −2 :  
2 1 0
A + 2I = 0 2 1 .
2 1 0
Un vecteur propre associé est :  
1
v3 = −2 .
4

4. Solution générale du système

Nous avons montré que la matrice A admet trois valeurs propres réelles et distinctes :
λ1 = 1, λ2 = −1, λ3 = −2,
avec pour vecteurs propres associés :
     
1 1 1
v1 = 1 ,
  v2 = −1 ,
 v3 = −2 .
1 1 4
Une propriété fondamentale des systèmes différentiels linéaires Y ′ = AY est la suivante :
— Si (λ, v) est un couple valeur propre / vecteur propre de A (c’est-à-dire Av = λv), alors
la fonction vectorielle
Y (x) = eλx v
est une solution du système Y ′ = AY .
Ainsi, à chaque valeur propre λi et son vecteur propre associé vi correspond une solution
fondamentale
Y (i) (x) = eλi x vi .
Nous obtenons donc les trois solutions fondamentales suivantes :

2
Première solution (associée à λ1 = 1).
   x
1 e
(1) λ1 x x 
Y (x) = e v1 = e 1 = ex  .

1 ex

Deuxième solution (associée à λ2 = −1).

  −x 

1 e
(2) λ2 x −x 
Y (x) = e v2 = e −1 = − e−x  .
 
1 e−x

Troisième solution (associée à λ3 = −2).

   −2x 
1 e
Y (3) (x) = eλ3 x v3 = e−2x −2 = −2e−2x  .
4 4e−2x
Comme le système Y ′ = AY est d’ordre 3 et que nous avons trouvé trois solutions vec-
torielles Y (1) , Y (2) , Y (3) linéairement indépendantes, la solution générale du système s’écrit
comme une combinaison linéaire de ces trois solutions :
Y (x) = C1 Y (1) (x) + C2 Y (2) (x) + C3 Y (3) (x),
où C1 , C2 , C3 sont des constantes réelles arbitraires.
En remplaçant explicitement, on obtient :
 
x −x −2x
 x
e
 −x 
e
 −2x 
e C 1 e + C2 e + C 3e
x −x  −2x 
 x −x

−2x 
Y (x) = C1  e + C2  − e + C3 −2e = C 1 e − C 2 e − 2C 3 e .
e−x 4e−2x

ex C ex + C e−x + 4C e−2x
1 2 3

5. Retour à la fonction y(x)

Par construction du système, nous avions posé :
   
y1 (x) y(x)
Y (x) = y2 (x) =  y ′ (x)  .
y3 (x) y ′′ (x)
Ainsi, la première composante de Y (x) est précisément la solution cherchée de l’équation
différentielle :
y(x) = y1 (x) = Y1 (x).
D’après l’expression obtenue ci-dessus pour Y1 (x), on a donc :
y(x) = C1 ex + C2 e−x + C3 e−2x .
On en déduit que la solution générale de l’équation différentielle
y ′′′ + 2y ′′ − y ′ − 2y = 0
est donnée par :
y(x) = C1 ex + C2 e−x + C3 e−2x .

3
Exercice 2
On considère l’équation différentielle linéaire homogène à coefficients constants :

y ′′ + 3y ′ + 2y = 0.

1. Réduction à un système d’ordre 1

On introduit les variables :
y1 = y, y2 = y ′ .
Alors les dérivées s’écrivent :
y1′ = y2 ,
et, comme y2 = y ′ ⇒ y2′ = y ′′ , on utilise l’équation différentielle

y ′′ + 3y ′ + 2y = 0 ⇒ y ′′ = −3y ′ − 2y.

Donc, en remplaçant y par y1 et y ′ par y2 :

y2′ = −3y2 − 2y1 .

On obtient ainsi le système :

( ′
y1 (x) = y2 (x),
y2′ (x) = −2y1 (x) − 3y2 (x).

Sous forme matricielle, on pose

   
y1 (x) y(x)
Y (x) = = ,
y2 (x) y ′ (x)

et le système devient
 
′ 0 1
Y (x) = A Y (x), avec A = .
−2 −3

2. Calcul des valeurs propres de A

Les valeurs propres λ de A sont les solutions de :

det(A − λI) = 0,

où  
−λ 1
A − λI = .
−2 −3 − λ
On calcule le déterminant :
−λ 1
det(A − λI) = = (−λ)(−3 − λ) − (1)(−2).
−2 −3 − λ

On développe :
(−λ)(−3 − λ) = λ(3 + λ) = λ2 + 3λ,

4
donc
det(A − λI) = λ2 + 3λ + 2.
L’équation caractéristique est donc :
λ2 + 3λ + 2 = 0.
On factorise ce polynôme :
λ2 + 3λ + 2 = (λ + 1)(λ + 2).
Les valeurs propres sont alors :
λ1 = −1, λ2 = −2.

3. Calcul des vecteurs propres

Vecteur propre associé à λ1 = −1. On résout :
(A − λ1 I)v1 = 0 ⇐⇒ (A + I)v1 = 0.
On a :      
0 1 1 0 1 1
A+I = + = .
−2 −3 0 1 −2 −2
 
x
On cherche un vecteur non nul v1 = ̸= 0 tel que (A + I)v1 = 0.
y
Le système associé est : (
1 · x + 1 · y = 0,
−2 · x − 2 · y = 0.
Les deux équations sont équivalentes à x + y = 0.
On prend par exemple x = 1, alors y = −1. Un vecteur propre associé à λ1 = −1 est donc :
 
1
v1 = .
−1

Vecteur propre associé à λ2 = −2. On résout :

(A − λ2 I)v2 = 0 ⇐⇒ (A + 2I)v2 = 0.
On a :      
0 1 2 0 2 1
A + 2I = + = .
−2 −3 0 2 −2 −1
 
x
On cherche v2 = ̸= 0 tel que (A + 2I)v2 = 0.
y
Le système associé est : (
2x + y = 0,
−2x − y = 0.
Les deux équations sont équivalentes à 2x + y = 0.
On prend par exemple x = 1, alors y = −2. Un vecteur propre associé à λ2 = −2 est donc :
 
1
v2 = .
−2

5
4. Solutions fondamentales du système Y ′ = AY
Pour chaque paire (λi , vi ), une solution fondamentale du système est :

Y (i) (x) = eλi x vi .

Solution associée à λ1 = −1.

   −x 
(1) λ1 x −x 1 e
Y (x) = e v1 = e = .
−1 − e−x

Solution associée à λ2 = −2.

   −2x 
(2) λ2 x −2x 1 e
Y (x) = e v2 = e = .
−2 −2e−2x

5. Solution générale du système et retour à y(x)

La solution générale du système Y ′ = AY est la combinaison linéaire des solutions fon-
damentales :
Y (x) = C1 Y (1) (x) + C2 Y (2) (x),
où C1 et C2 sont des constantes réelles arbitraires.
En remplaçant Y (1) (x) et Y (2) (x), on obtient :
!
e−x C1 e−x + C2 e−2x
   −2x 
e
Y (x) = C1 −x + C2 = .
−e −2e−2x −C1 e−x − 2C2 e−2x

Or, par définition,    

y1 (x) y(x)
Y (x) = = .
y2 (x) y ′ (x)
On identifie donc :
y(x) = y1 (x) = C1 e−x + C2 e−2x .

6. Solution générale de l’équation différentielle

La solution générale de l’équation différentielle

y ′′ + 3y ′ + 2y = 0

est donc :
y(x) = C1 e−x + C2 e−2x ,
où C1 et C2 sont des constantes réelles arbitraires.

6
Exercice 3
On considère l’équation différentielle linéaire homogène :

y ′′ + 2y ′ − 24y = 0.

1. Réduction à un système d’ordre 1

On introduit les variables :
y1 = y, y2 = y ′ .
Alors :
y1′ = y2 .
Et, comme y2 = y ′ , on a y2′ = y ′′ . En utilisant l’équation différentielle :

y ′′ + 2y ′ − 24y = 0 ⇒ y ′′ = −2y ′ + 24y.

En remplaçant y par y1 et y ′ par y2 :

y2′ = −2y2 + 24y1 .

On obtient le système : ( ′
y1 (x) = y2 (x),
y2′ (x) = 24y1 (x) − 2y2 (x).
On pose le vecteur d’état :    
y1 (x) y(x)
Y (x) = = .
y2 (x) y ′ (x)
Le système s’écrit sous forme matricielle :
 
′ 0 1
Y (x) = A Y (x), avec A = .
24 −2

2. Valeurs propres de la matrice A

Les valeurs propres λ de A sont solutions de :

det(A − λI2 ) = 0.

On a :      
0 1 λ 0 −λ 1
A − λI2 = − = .
24 −2 0 λ 24 −2 − λ
Le déterminant est :
−λ 1
det(A − λI2 ) = = (−λ)(−2 − λ) − 1 · 24.
24 −2 − λ

On développe :
(−λ)(−2 − λ) = λ(2 + λ) = λ2 + 2λ.
Donc :
det(A − λI2 ) = λ2 + 2λ − 24.

7
L’équation caractéristique est :
λ2 + 2λ − 24 = 0.
On cherche deux nombres dont le produit vaut −24 et la somme 2 : ce sont 6 et −4. On
factorise :
λ2 + 2λ − 24 = (λ + 6)(λ − 4).
Les valeurs propres sont donc :

λ1 = 4, λ2 = −6.

3. Vecteurs propres
Vecteur propre associé à λ1 = 4. On résout :

(A − 4I2 )v1 = 0.

On a :      
0 1 4 0 −4 1
A − 4I2 = − = .
24 −2 0 4 24 −6
 
x
On cherche v1 = ̸= 0 tel que (A − 4I2 )v1 = 0.
y
Le système associé est : (
−4x + y = 0,
24x − 6y = 0.
La première équation donne y = 4x. La seconde devient alors 24x − 6 · 4x = 24x − 24x = 0,
vérifiée.
On peut choisir x = 1, donc y = 4. Un vecteur propre associé à λ1 = 4 est :
 
1
v1 = .
4

Vecteur propre associé à λ2 = −6. On résout :

(A + 6I2 )v2 = 0.

On a :      
0 1 6 0 6 1
A + 6I2 = + = .
24 −2 0 6 24 4
 
x
On cherche v2 = ̸= 0 tel que (A + 6I2 )v2 = 0.
y
Le système associé est : (
6x + y = 0,
24x + 4y = 0.
La première équation donne y = −6x. La seconde devient 24x + 4(−6x) = 24x − 24x = 0,
vérifiée.
On peut choisir x = 1, donc y = −6. Un vecteur propre associé à λ2 = −6 est :
 
1
v2 = .
−6

8
4. Solution générale du système Y ′ = AY
Pour chaque paire (λi , vi ), une solution fondamentale est :

Y (i) (x) = eλi x vi .

On obtient :
   4x     −6x 
(1) 4x 1 e (2) −6x 1 e
Y (x) = e = , Y (x) = e = .
4 4e4x −6 −6e−6x

La solution générale du système est la combinaison linéaire :

Y (x) = C1 Y (1) (x) + C2 Y (2) (x),

soit !
−6x C1 e4x + C2 e−6x
 4x
  
e e
Y (x) = C1 + C2 = .
4e4x −6e−6x 4C1 e4x − 6C2 e−6x

5. Retour à la solution y(x)

Par définition,    
y1 (x) y(x)
Y (x) = = .
y2 (x) y ′ (x)
Donc la première composante de Y (x) donne la solution recherchée :

y(x) = y1 (x) = C1 e4x + C2 e−6x .

6. Solution générale de l’équation différentielle

La solution générale de l’équation

y ′′ + 2y ′ − 24y = 0

est donc :
y(x) = C1 e4x + C2 e−6x ,
où C1 et C2 sont des constantes réelles arbitraires.

9
Exercice 4
On considère le système :
( ′
y1 (t) = y1 (t) + y2 (t) + 10 cos t,
y2′ (t) = 3y1 (t) − y2 (t) − 10 sin t.

1. Mise sous forme matricielle

On pose le vecteur d’état  
y1 (t)
Y (t) = .
y2 (t)
Alors le système peut s’écrire :

Y ′ (t) = A Y (t) + F (t),

avec    
1 1 10 cos t
A= , F (t) = .
3 −1 −10 sin t
On va d’abord résoudre le système homogène associé Yh′ (t) = A Yh (t), puis chercher une
solution particulière Yp (t) du système complet, et enfin Y (t) = Yh (t) + Yp (t).

2. Système homogène Yh′ = AYh

2.1 Valeurs propres de A

Les valeurs propres λ de A sont solutions de

det(A − λI2 ) = 0.

On a :  
1−λ 1
A − λI2 = ,
3 −1 − λ
donc
1−λ 1
det(A − λI2 ) = = (1 − λ)(−1 − λ) − 3 · 1.
3 −1 − λ
On développe :
(1 − λ)(−1 − λ) = −1 − λ + λ + λ2 = −1 + λ2 .
Ainsi :
det(A − λI2 ) = (−1 + λ2 ) − 3 = λ2 − 4.
L’équation caractéristique est donc :

λ2 − 4 = 0 =⇒ (λ − 2)(λ + 2) = 0.

Les valeurs propres sont :

λ1 = 2, λ2 = −2.

10
2.2 Vecteurs propres

Vecteur propre associé à λ1 = 2. On résout :

(A − 2I2 )v1 = 0.
On a :    
1−2 1 −1 1
A − 2I2 = = .
3 −1 − 2 3 −3
 
x
On cherche v1 = ̸= 0 tel que (A − 2I2 )v1 = 0.
y
Le système associé est : (
−x + y = 0,
3x − 3y = 0.
Les deux équations donnent y = x.
On choisit x = 1, alors y = 1. Un vecteur propre est donc :
 
1
v1 = .
1

Vecteur propre associé à λ2 = −2. On résout :

(A + 2I2 )v2 = 0.
On a :    
1+2 1 3 1
A + 2I2 = = .
3 −1 + 2 3 1
 
x
On cherche v2 = ̸= 0 tel que (A + 2I2 )v2 = 0.
y
Le système associé est : (
3x + y = 0,
3x + y = 0.
On a donc y = −3x.
On prend par exemple x = 1, alors y = −3. Un vecteur propre est donc :
 
1
v2 = .
−3

2.3 Solution générale du système homogène

Pour chaque paire (λi , vi ), une solution fondamentale de Yh′ = AYh est
(i)
Yh (t) = eλi t vi .
On obtient :    
(1) 1
2t (2) −2t 1
Yh (t) =e , Yh (t) =e .
1 −3
La solution générale du système homogène est donc :
   
2t 1 −2t 1
Yh (t) = C1 e + C2 e ,
1 −3
où C1 et C2 sont des constantes réelles arbitraires.

11
3. Recherche d’une solution particulière Yp (t)
Le second membre F (t) est trigonométrique. On essaie donc une solution particulière de
la forme :
Yp (t) = a cos t + b sin t,
où a et b sont deux vecteurs constants :
   
a1 b
a= , b= 1 .
a2 b2

3.1 Calcul de Yp′ (t)

On dérive Yp (t) :
Yp′ (t) = −a sin t + b cos t.

3.2 Substitution dans l’équation Y ′ = AY + F (t)

On impose que Yp satisfasse :

Yp′ (t) = A Yp (t) + F (t).

En remplaçant :
 
10 cos t
−a sin t + b cos t = A(a cos t + b sin t) + .
−10 sin t

On sépare les parties en cos t et en sin t. À droite :

A(a cos t + b sin t) = (Aa) cos t + (Ab) sin t.

Donc :
   
10
0

A(a cos t + b sin t) + F (t) = Aa + cos t + Ab + sin t.
0 −10

À gauche :
−a sin t + b cos t = b cos t + (−a) sin t.
En identifiant les coefficients devant cos t et sin t, on obtient le système vectoriel :
 !
 10
b = Aa + 0 ,



!
 0
−a = Ab + −10 .




3.3 Passage en coordonnées

On écrit Aa et Ab explicitement avec

     
1 1 a1 b
A= , a= , b= 1 .
3 −1 a2 b2

12
On a :    
a1 + a2 b1 + b2
Aa = , Ab = .
3a1 − a2 3b1 − b2
Les deux équations vectorielles deviennent donc, composante par composante :


 b1 = a1 + a2 + 10,


b2 = 3a1 − a2 ,




 −a1 = b1 + b2 ,


−a2 = 3b1 − b2 − 10.


3.4 Résolution du système pour a1 , a2 , b1 , b2

À partir des deux premières équations :

b1 = a1 + a2 + 10, b2 = 3a1 − a2 .
On calcule d’abord b1 + b2 :
b1 + b2 = (a1 + a2 + 10) + (3a1 − a2 ) = 4a1 + 10.
La troisième équation −a1 = b1 + b2 devient :
−a1 = 4a1 + 10 =⇒ −5a1 = 10 =⇒ a1 = −2.
On remplace dans les expressions de b1 et b2 :
b1 = a1 + a2 + 10 = −2 + a2 + 10 = a2 + 8,
b2 = 3a1 − a2 = 3(−2) − a2 = −6 − a2 .
On utilise maintenant la quatrième équation :
−a2 = 3b1 − b2 − 10.
On calcule 3b1 − b2 − 10 :
3b1 − b2 − 10 = 3(a2 + 8) − (−6 − a2 ) − 10 = 3a2 + 24 + 6 + a2 − 10 = 4a2 + 20.
Ainsi :
−a2 = 4a2 + 20 =⇒ −5a2 = 20 =⇒ a2 = −4.
On en déduit :
b1 = a2 + 8 = −4 + 8 = 4, b2 = −6 − a2 = −6 − (−4) = −2.
Donc finalement :    
−2 4
a= , b= .
−4 −2

3.5 Expression de Yp (t)

La solution particulière est donc

!
−2 cos t + 4 sin t
   
−2 4
Yp (t) = a cos t + b sin t = cos t + sin t = .
−4 −2 −4 cos t − 2 sin t

13
4. Solution générale du système complet
La solution générale du système

Y ′ (t) = A Y (t) + F (t)

est
Y (t) = Yh (t) + Yp (t),
c’est-à-dire :
     
2t1 −2t 1 −2 cos t + 4 sin t
Y (t) = C1 e + C2 e + .
1 −3 −4 cos t − 2 sin t

En écrivant composante par composante :


y1 (t) = C1 e2t + C2 e−2t − 2 cos t + 4 sin t,
y (t) = C e2t − 3C e−2t − 4 cos t − 2 sin t.
2 1 2

14
Exercice 5
On considère le système :
( ′
y1 (t) = 4y1 (t) − 8y2 (t) + 2 cosh t,
y2′ (t) = 2y1 (t) − 6y2 (t) + cosh t + 2 sinh t.

1. Mise sous forme matricielle

On pose le vecteur d’état :  
y1 (t)
Y (t) = .
y2 (t)
Le système peut s’écrire :
Y ′ (t) = A Y (t) + F (t),
avec    
4 −8 2 cosh t
A= , F (t) = .
2 −6 cosh t + 2 sinh t
On va résoudre le système homogène associé Yh′ = AYh , chercher une solution particulière Yp
du système complet, et conclure que Y = Yh + Yp .

2. Système homogène Yh′ = AYh

2.1 Valeurs propres de A

Les valeurs propres λ de A satisfont :

det(A − λI2 ) = 0.

On a :  
4−λ −8
A − λI2 = .
2 −6 − λ
Le déterminant vaut :
4−λ −8
det(A − λI2 ) = = (4 − λ)(−6 − λ) − (−8) · 2.
2 −6 − λ
On développe d’abord :

(4 − λ)(−6 − λ) = −24 − 4λ + 6λ + λ2 = λ2 + 2λ − 24.

Puis on ajoute 16 :

det(A − λI2 ) = (λ2 + 2λ − 24) + 16 = λ2 + 2λ − 8.

L’équation caractéristique est donc :

λ2 + 2λ − 8 = 0.

On factorise :
λ2 + 2λ − 8 = (λ − 2)(λ + 4).
Les valeurs propres sont :
λ1 = 2, λ2 = −4.

15
2.2 Vecteurs propres

Vecteur propre associé à λ1 = 2. On résout

(A − 2I2 )v1 = 0.
On a :    
4−2 −8 2 −8
A − 2I2 = = .
2 −6 − 2 2 −8
 
x
On cherche v1 = ̸= 0 tel que (A − 2I2 )v1 = 0.
y
Le système associé est : (
2x − 8y = 0,
2x − 8y = 0.
La relation est 2x − 8y = 0 ⇒ x = 4y.
On choisit y = 1, alors x = 4. Un vecteur propre associé à λ1 = 2 est donc :
 
4
v1 = .
1

Vecteur propre associé à λ2 = −4. On résout

(A + 4I2 )v2 = 0.
On a :    
4+4 −8 8 −8
A + 4I2 = = .
2 −6 + 4 2 −2
 
x
On cherche v2 = ̸= 0 tel que (A + 4I2 )v2 = 0.
y
Le système associé : (
8x − 8y = 0,
2x − 2y = 0,
donne x = y.
On choisit x = 1, donc y = 1. Un vecteur propre associé à λ2 = −4 est :
 
1
v2 = .
1

2.3 Solution générale du système homogène

Pour chaque (λi , vi ), une solution fondamentale est

(i)
Yh (t) = eλi t vi .
On obtient :    
(1) 2t4 (2) −4t 1
Yh (t) =e , Yh (t) =e .
1 1
La solution générale du système homogène est donc :
   
2t 4 −4t 1
Yh (t) = C1 e + C2 e ,
1 1
où C1 et C2 sont des constantes réelles.

16
3. Recherche d’une solution particulière Yp (t)
Le second membre F (t) contient cosh t et sinh t. On essaie donc une solution particulière
de la forme :
Yp (t) = a cosh t + b sinh t,
où a et b sont deux vecteurs constants :
   
a1 b
a= , b= 1 .
a2 b2

3.1 Dérivée de Yp (t)

On utilise :
d d
(cosh t) = sinh t, (sinh t) = cosh t.
dt dt
Ainsi :
Yp′ (t) = a sinh t + b cosh t.

3.2 Substitution dans Y ′ = AY + F (t)

On impose que Yp vérifie :

Yp′ (t) = A Yp (t) + F (t).

En remplaçant :
 
2 cosh t
a sinh t + b cosh t = A(a cosh t + b sinh t) + .
cosh t + 2 sinh t

On développe le terme A(a cosh t + b sinh t) :

A(a cosh t + b sinh t) = (Aa) cosh t + (Ab) sinh t.

Donc :    
2
0

a sinh t + b cosh t = Aa + cosh t + Ab + sinh t.
1 2
En identifiant les coefficients de cosh t et sinh t, on obtient le système vectoriel :
 !
 2
b = Aa + 1 ,



!
 0
a = Ab + 2 .




3.3 Passage en coordonnées et résolution

Avec      
4 −8 a1 b
A= , a= , b= 1 ,
2 −6 a2 b2
on a :    
4a1 − 8a2 4b1 − 8b2
Aa = , Ab = .
2a1 − 6a2 2b1 − 6b2

17
Les équations vectorielles deviennent donc, composante par composante :


 b1 = 4a1 − 8a2 + 2,


b2 = 2a1 − 6a2 + 1,




 a1 = 4b1 − 8b2 ,


a2 = 2b1 − 6b2 + 2.


On résout ce système. On trouve après simplification :

a1 = 0, a2 = 0, b1 = 2, b2 = 1.

Ainsi :    
0 2
a= , b= .
0 1

3.4 Expression de Yp (t)

La solution particulière est donc :

     
0 2 2 sinh t
Yp (t) = a cosh t + b sinh t = cosh t + sinh t = .
0 1 sinh t

4. Solution générale du système complet

La solution générale du système

Y ′ (t) = A Y (t) + F (t)

est :
Y (t) = Yh (t) + Yp (t).
On regroupe :      
2t4 −4t 1 2 sinh t
Y (t) = C1 e + C2 e + .
1 1 sinh t
En écrivant composante par composante :

y1 (t) = 4C1 e2t + C2 e−4t + 2 sinh t,
y (t) = C e2t + C e−4t + sinh t.
2 1 2

18
Exercice 6
On considère le système

y ′ (t) = −3y1 (t) − 4y2 (t) + 5et ,
 1


y2′ (t) = 5y1 (t) + 6y2 (t) − 6et ,


y1 (0) = 19, y2 (0) = −23.


1. Mise sous forme matricielle

On pose le vecteur d’état :  
y1 (t)
Y (t) = .
y2 (t)
Alors le système s’écrit :
Y ′ (t) = A Y (t) + F (t),
avec      
−3 −4 5et t 5
A= , F (t) = =e .
5 6 −6et −6
On écrit la solution générale sous la forme Y (t) = Yh (t) + Yp (t), où Yh est la solution du
système homogène Yh′ = AYh et Yp une solution particulière du système complet.

2. Système homogène Yh′ = AYh

2.1 Valeurs propres de A

Les valeurs propres λ sont solutions de :

det(A − λI2 ) = 0.

On calcule :  
−3 − λ −4
A − λI2 = .
5 6−λ
Le déterminant vaut :
−3 − λ −4
det(A − λI2 ) = = (−3 − λ)(6 − λ) − (−4) · 5.
5 6−λ

On développe :

(−3 − λ)(6 − λ) = −18 + 3λ − 6λ + λ2 = λ2 − 3λ − 18.

Ainsi :
det(A − λI2 ) = λ2 − 3λ − 18 + 20 = λ2 − 3λ + 2.
L’équation caractéristique est donc :

λ2 − 3λ + 2 = 0 =⇒ (λ − 1)(λ − 2) = 0.

Les valeurs propres sont :

λ1 = 1, λ2 = 2.

19
2.2 Vecteurs propres

Vecteur propre associé à λ1 = 1. On résout :

(A − I2 )v1 = 0.
   
−3 − 1 −4 −4 −4
A − I2 = = .
5 6−1 5 5
 
x
On cherche v1 = ̸= 0 tel que (A − I2 )v1 = 0.
y
Le système associé : (
−4x − 4y = 0,
=⇒ x + y = 0.
5x + 5y = 0
On peut choisir x = −1, y = 1. Un vecteur propre associé est donc :
 
−1
v1 = .
1

Vecteur propre associé à λ2 = 2. On résout :

(A − 2I2 )v2 = 0.
   
−3 − 2 −4 −5 −4
A − 2I2 = = .
5 6−2 5 4
 
x
On cherche v2 = ̸= 0 tel que (A − 2I2 )v2 = 0.
y
Le système : (
−5x − 4y = 0,
=⇒ −5x − 4y = 0.
5x + 4y = 0
On prend y = 5, alors x = −4. Un vecteur propre associé est donc :
 
−4
v2 = .
5

2.3 Solution générale du système homogène

Pour chaque (λi , vi ), une solution fondamentale est

(i)
Yh (t) = eλi t vi .

Donc :    
(1) t −1 (2) 2t−4
Yh (t) =e , Yh (t) =e .
1 5
La solution générale homogène est :
   
t −1 2t −4
Yh (t) = C1 e + C2 e ,
1 5

où C1 , C2 ∈ R.

20
3. Recherche d’une solution particulière Yp (t)
Le second membre est  
t 5
F (t) = e b, b= .
−6
Comme et correspond à la valeur propre λ1 = 1 de A (résonnance), on utilise l’ansatz :

Yp (t) = et Ct + D ,

où C et D sont des vecteurs constants

   
c d1
C= 1 , D= .
c2 d2

3.1 Calcul de Yp′ (t)

On a :
Yp (t) = et (Ct + D).
Par dérivation (produit) :

Yp′ (t) = et (Ct + D) + et C = et Ct + D + C .

3.2 Substitution dans l’équation

L’équation Y ′ (t) = AY (t) + F (t) devient pour Yp :

et Ct + D + C = A et (Ct + D) + et b = et ACt + AD + b .

On simplifie par et :
Ct + D + C = ACt + AD + b.

Coefficient en t :
C = AC ⇒ (A − I2 )C = 0.
Donc C appartient au noyau de (A − I2 ), c’est-à-dire à l’espace propre associé à λ = 1. On
choisit C proportionnel à v1 = (−1, 1)T . Pour simplifier :
 
1
C= = −v1 .
−1

Terme constant :
D + C = AD + b ⇒ (I2 − A)D = b − C.
On calcule :      
1 0 −3 −4 4 4
I2 − A = − = .
0 1 5 6 −5 −5
On a aussi      
5 1 4
b−C = − = .
−6 −1 −5
On cherche D = (d1 , d2 )T tel que :
    
4 4 d1 4
= .
−5 −5 d2 −5

21
Le système associé est : (
4d1 + 4d2 = 4,
−5d1 − 5d2 = −5.
Les deux équations donnent la même relation d1 + d2 = 1.
On choisit par simplicité d1 = 1, d2 = 0, d’où
 
1
D= .
0

3.3 Expression de Yp (t)

On a finalement :    
1 1
C= , D= .
−1 0
Donc    
t t+1 (t + 1)et
t
Yp (t) = e (Ct + D) = e = .
−t −tet

4. Solution générale du système

La solution générale du système complet est :
Y (t) = Yh (t) + Yp (t).
Soit      
t −1 2t −4 (t + 1)et
Y (t) = C1 e + C2 e + .
1 5 −tet
En écrivant composante par composante :
(
y1 (t) = −C1 et − 4C2 e2t + (t + 1)et ,
y2 (t) = C1 et + 5C2 e2t − tet .

5. Utilisation des conditions initiales

On impose :
y1 (0) = 19, y2 (0) = −23.
Comme e0 = 1, on a :
y1 (0) = −C1 − 4C2 + 1 = 19,
y2 (0) = C1 + 5C2 = −23.
On obtient le système : (
−C1 − 4C2 = 18,
C1 + 5C2 = −23.
En additionnant les deux équations :
(−C1 − 4C2 ) + (C1 + 5C2 ) = 18 + (−23) ⇒ C2 = −5.
En remplaçant dans C1 + 5C2 = −23 :
C1 + 5(−5) = −23 ⇒ C1 − 25 = −23 ⇒ C1 = 2.

22
6. Solution particulière avec conditions initiales
En substituant C1 = 2 et C2 = −5 dans les expressions de y1 (t) et y2 (t), on obtient :

y1 (t) = −2et − 4(−5)e2t + (t + 1)et = (t − 1)et + 20e2t .

y2 (t) = 2et + 5(−5)e2t − tet = (2 − t)et − 25e2t .

Ainsi, la solution vérifiant les conditions initiales est :

y1 (t) = (t − 1)et + 20e2t ,

y2 (t) = (2 − t)et − 25e2t .

23
Exercice 7
On considère le système :

y ′ (t) = y2 (t) − 5 sin t,
 1


y2′ (t) = −4y1 (t) + 17 cos t,


y1 (0) = 5, y2 (0) = 2.


1. Réduction à une équation d’ordre 2 pour y1

À partir de la première équation :
y1′ (t) = y2 (t) − 5 sin t ⇒ y2 (t) = y1′ (t) + 5 sin t.
On dérive cette relation :
y2′ (t) = y1′′ (t) + 5 cos t.
Mais d’après la deuxième équation du système :
y2′ (t) = −4y1 (t) + 17 cos t.
On égalise les deux expressions de y2′ (t) :
y1′′ (t) + 5 cos t = −4y1 (t) + 17 cos t.
On regroupe les termes :
y1′′ (t) + 4y1 (t) = 12 cos t.
On obtient donc l’équation différentielle d’ordre 2 pour y1 :
y1′′ (t) + 4y1 (t) = 12 cos t.

2. Conditions initiales pour y1

On connaît y1 (0) directement :
y1 (0) = 5.
Pour y1′ (0), on utilise la première équation du système à t = 0 :
y1′ (0) = y2 (0) − 5 sin 0 = 2 − 0 = 2.
Donc :
y1 (0) = 5, y1′ (0) = 2.

3. Résolution de l’équation y1′′ + 4y1 = 12 cos t

3.1 Solution de l’équation homogène

On considère d’abord l’équation homogène associée :

y1′′ + 4y1 = 0.
L’équation caractéristique est :
r2 + 4 = 0 ⇒ r2 = −4 ⇒ r = ±2i.
La solution générale homogène est donc :
y1h (t) = C1 cos(2t) + C2 sin(2t),
où C1 et C2 sont des constantes réelles.

24
3.2 Recherche d’une solution particulière

Le second membre est 12 cos t. On essaie une solution particulière de la forme :

y1p (t) = a cos t + b sin t,

′′
avec a, b ∈ R. On calcule y1p (t) :
′ ′′
y1p (t) = −a sin t + b cos t, y1p (t) = −a cos t − b sin t.

On remplace dans l’équation :

′′
y1p + 4y1p = (−a cos t − b sin t) + 4(a cos t + b sin t) = 3a cos t + 3b sin t.

On impose :
3a cos t + 3b sin t = 12 cos t.
Par identification des coefficients :
(
3a = 12,
⇒ a = 4, b = 0.
3b = 0.

Ainsi, une solution particulière est :

y1p (t) = 4 cos t.

3.3 Solution générale de l’équation de y1

La solution générale de
y1′′ + 4y1 = 12 cos t
est donc :
y1 (t) = y1h (t) + y1p (t) = C1 cos(2t) + C2 sin(2t) + 4 cos t.

4. Utilisation des conditions initiales

On impose y1 (0) = 5 et y1′ (0) = 2.

Condition sur y1 (0). On évalue y1 (t) en t = 0 :

y1 (0) = C1 cos(0) + C2 sin(0) + 4 cos 0 = C1 + 4.

On impose y1 (0) = 5 :
C1 + 4 = 5 ⇒ C1 = 1.

Condition sur y1′ (0). On calcule d’abord y1′ (t) :

y1′ (t) = −2C1 sin(2t) + 2C2 cos(2t) − 4 sin t.

En t = 0 :
y1′ (0) = −2C1 sin 0 + 2C2 cos 0 − 4 sin 0 = 2C2 .
On impose y1′ (0) = 2 :
2C2 = 2 ⇒ C2 = 1.

25
4.1 Expression finale de y1 (t)

On a donc C1 = 1 et C2 = 1, d’où :

y1 (t) = cos(2t) + sin(2t) + 4 cos t.

5. Calcul de y2 (t) à partir de la première équation du système

On rappelle la première équation :

y1′ (t) = y2 (t) − 5 sin t ⇒ y2 (t) = y1′ (t) + 5 sin t.

On avait :
y1′ (t) = −2 sin(2t) + 2 cos(2t) − 4 sin t.
Donc :
 
y2 (t) = −2 sin(2t) + 2 cos(2t) − 4 sin t + 5 sin t = −2 sin(2t) + 2 cos(2t) + sin t.

On obtient :
y2 (t) = −2 sin(2t) + 2 cos(2t) + sin t.

6. Vérification des conditions initiales pour y2

On vérifie que y2 (0) = 2. En t = 0 :

y2 (0) = −2 sin(0) + 2 cos(0) + sin(0) = 0 + 2 + 0 = 2.

Les conditions initiales sont donc bien satisfaites.

7. Résumé de la solution
La solution du système 
y ′ = y2 − 5 sin t,
 1


y2′ = −4y1 + 17 cos t,


y1 (0) = 5, y2 (0) = 2


est :
y1 (t) = cos(2t) + sin(2t) + 4 cos t,
y2 (t) = −2 sin(2t) + 2 cos(2t) + sin t.

26