Équations Diérentielles Ordinaires du Second Ordre

A. MEDELFEF
Faculté de Physique
Université des Sciences et de la Technologie Houari Boumediene

3ème Année Licence de Physique Energétique

 Sommaire

1 Dénitions et concepts de base

2 Équations diérentielles linéaires

Principe de superposition
Le Wronskien

3 Équations homogènes à coecients constants

Problème aux valeurs initiales (IVP problem)
Étude du polynôme caractéristique

4 L'oscillateur mécanique : un exemple récapitulatif

5 Équations non homogènes à coecients constant

Quelques cas particuliers
Les phénomènes de Résonance
La variation de la constante

A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 2 / 17

 Dénitions et concepts de base

Soit la fonction y(t) de la variable indépendante t 1 . Une EDO du second ordre est
une équation liant la fonction y(t) à ses dérivées première y0 (t) et seconde y00 (t)
ayant la forme standard suivante :
y00 (t) = F (t; y(t); y0 (t)

(1)
Dans ce cours, on va s'intéresser uniquement au cas où F (t; y(t); y0 (t) peut être


écrit sous la forme xp(t)y0 (t) x q(t)y(t) + g (t), auquel cas, la forme standard (1)
devient :
y00 (t) + p(t)y0 (t) + q(t)y(t) = g (t) (2)
ici on distingue les propriétés suivantes :
a Si p(t ), q (t ) et g (t ) ne dépendent pas de y (t ), alors l'EDO est dite linéaire
b Si g (t ) = 0, alors l'EDO est dite homogène
c Si les coecients p(t ) et q (t ) sont des constantes, c.-à-d. p(t ) = p et
q(t) = q, alors l'EDO est dite à coecients constants

1. Cette variable peut être x pour y(x), ou toute autre variable indépendante !
A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 3 / 17
 Équations diérentielles linéaires Principe de superposition

Principe de superposition
Si y1 (t) et y2 (t) sont deux solutions indépendantes de l'EDO (2) homogène, alors
y( t) = C1 y1 (t) + C2 y2 (t) est aussi une solution.

Il sut de remplacer dans l'EDO (2) pour trouver

C1 [y00 (t) + p(t)y0 (t) + q(t)y1 (t)] + C2 [y00 (t) + p(t)y0 (t) + q(t)y2 (t)] = 0
1 1 2 2

Le même principe s'applique au cas d'une EDO non-homogène. En eet, on peut

aisément montrer que si yh (t) = C1 y1 (t)+ C2 y2 (t) est solution du problème homo-
gène, et que yp (t) est une solution particulière du cas général (g (t) ä 0), alors la
solution générale de l'EDO est
y(t) = yh (t) + yp (t)
De même, si le terme g (t) est constitué de plusieurs fonctions, par exemple g1 (t) =
cos(! t) et g2 (t) = exp( t) soit g (t) = cos(! t) + exp( t), alors yp (t) est donnée
par
yp (t) = yp1 (t) + yp2 (t)
et de ce fait, la solution générale est tout simplement la somme de toutes les
solutions.
Comment déterminer si y1 (t) et y2 (t) sont des solutions indépendantes ?
A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 4 / 17
 Équations diérentielles linéaires Le Wronskien

Indépendance des solutions : le Wronskien

Si y1 (t) et y2 (t) sont deux solutions de l'EDO (2), leur wronskien est déni de la
manière suivante :
y1 (t) y2 (t)
W fy1 (t); y2 (t)g = = y (t)y0 (t) x y0 (t)y (t) (3)
y10 (t) y20 (t) 1 2 1 2

Si W fy1 (t); y2 (t)g ä 0, alors les solutions y1 (t) et y2 (t) sont indépendantes.
Quelques exemples :
a Montrer que y1 (t ) = cos (!t ) et y2 (t ) = sin(!t ) sont indépendantes.
On calcule le Wronskien
W fcos(!t); sin(!t)g = cos(!t) !cos(!t) x (x!sin(!t)) cos(!t)
= ! cos (!t) + sin (!t) ä 0


2 2

donc y (t) = cos(!t) et y (t) = sin(!t) sont indépendantes et on pourra

écrire yh (t) = C y (t) + C y (t).
1 2

1 1 2 2

b Sous quelle condition y (t) = exp( t) et y (t) = exp( t) peuvent former une
1 2

solution ?
A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 5 / 17
 Équations diérentielles linéaires Le Wronskien

b (suite) On calcule toujours le Wronskien :

W fexp( t); exp( t)g = exp( t) exp( t) x exp( t) exp( t)
= ( x ) exp([ + ]t)
Pour que y1 (t) et y2 (t) soient indépendantes, il faut que
W fexp( t); exp( t)g ä 0, c.-a-d., ä .

Exercice :
Montrer que
1 y1 (t) = exp( t)sin(!t) et y2 (t) = exp( t)cos(!t)
2 y1 (t) = exp( t) et y2 (t) = t exp( t)
constituent des familles de solutions indépendantes.

A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 6 / 17

 Équations homogènes à coecients constants Problème aux valeurs initiales (IVP problem)

Avant de commencer, il est utile de savoir ce qu'un IVP !

Problème aux valeurs initiales (IVP problem)
On appelle problème aux valeurs initiales le problème mathématique suivant
y00 (t) + py0 (t) + qy(t) = g (t)


y(t0 ) = y0 ; y0 (t0 ) = v0

Nous allons traiter le cas simple d'un IVP homogène pour lequel nous avons :
p(t) = p g (t) = 0 q(t) = q et les CIs : y(0) = y0 y0 (0) = v0
c'est à dire :
y00 (t) + py0 (t) + qy(t) = 0


y(0) = y0 ; y0 (0) = v0
(4)

Dans ce cas, nous cherchons une solution de type exponentielle y(t)  exp( t) .
En reportant dans (4), on obtient une équation, dite polynôme caractéristique, pour
les valeurs du paramètre .
2 + p + q = 0 (5)
Selon la valeur du discriminant
= p2 x 4q, trois cas sont à distinguer.
A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 7 / 17
 Équations homogènes à coecients constants Étude du polynôme caractéristique

= p2 x 4q > 0
Dans ce cas, nous avons deux racines réelles distinctes 1 et 2 telles que
xp ~ p2 x 4q
p

1;2 =
2
De ce fait, nous obtenons deux familles de solutions y1 (t) = exp(1 t) et y2 (t) =
exp(2 t).
La solutions générale est alors une combinaison linéaire de y1 (t) et y2 (t) :
y(t) = C1 exp(1 t) + C2 exp(2 t) (6)
où C1 et C2 sont déterminés en appliquant les CIs y(0) = y0 et y0 (0) = v0 .
Ce faisant, on obtient le système d'équations suivant pour C1 et C2 :
C1 + C2 = y0


1 C1 + 2 C2 = v0
dont la solution permet d'obtenir y(t) :
v0 x y0 2 v0 x y0 1
   

y(t) = exp(1 t) x exp(2 t)

1 x 2 1 x 2
A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 8 / 17
 Équations homogènes à coecients constants Étude du polynôme caractéristique

= p2 x 4q < 0
Dans ce cas, nous avons deux racines complexes conjuguées distinctes  + j! et
 x j! telles que

= xp ~ j 24q x p =  ~ j!
p
2

z 1; 2
De ce fait, nous obtenons deux familles de solutions y1 (t) = exp(t)cos(!t) et
y2 (t) = exp(t)sin(!t).
La solutions générale est alors une combinaison linéaire de y1 (t) et y2 (t) :
y(t) = C1 exp(t)cos(!t) + C2 exp(t)sin(!t) (7)
où C1 et C2 sont déterminés en appliquant les CIs y(0) = y0 et y0 (0) = v0 .
Ce faisant, on obtient le système d'équations suivant pour C1 et C2 :
C1 = y0


C1 + !C2 = v0
dont la solution permet d'obtenir y(t) :
v0 x y0 
 

y(t) = y0 exp(t)cos(!t) + exp(t)sin(!t)

!
A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 9 / 17
 Équations homogènes à coecients constants Étude du polynôme caractéristique

Pourquoi y1 (t) = exp(t)cos(!t) et y2 (t) = exp(t)sin(!t) ?

nous avons
z(t) = exp([ + j!]t) = exp(t)exp(j!t)


z (t) = exp([ x j!]t) = exp(t)exp(xj!t)

qui sont des fonctions complexes. Or, nous cherchons des solutions réelles !
En se servant de l'identité d'Euler :
exp(~j ) = cos( ) ~ jsin( )

et du principe de superposition, il vient :

y1 (t) = (z(t) + z (t))=2 = Refz(t)g = exp(t)cos(!t)


y2 (t) = (z(t) x z (t))=2 = I mfz(t)g = exp(t)sin(!t)

qui sont bel et bien des fonctions réelles.

A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 10 / 17

 Équations homogènes à coecients constants Étude du polynôme caractéristique

= p2 x 4q = 0
Dans ce cas, nous avons une racine double !  = xp=2.
Question : Peut-on obtenir y1 (t) et y2 (t) directement ?
Réponse : Non ! Car y1 (t) et y2 (t) doivent être indépendantes ce qui n'est pas le
cas d'une racine double (W fexp(t; exp(t)g = 0)
Deux issues nous sont proposées :
a Changer la forme de la solution en cherchant y (t ) = Pn (t ) exp(t ) où Pn (t )
est un polynôme à déterminer
En eet, en injectant y(t) dans l'EDO, il vient :
P 00 (t) = 0 =) Pn (t) = C1 + C2 t =) y(t) = (C1 + C2 t) exp(t)
n

L'application des CIs, y(0) = y0 et y0 (0) = v0 , permet d'obtenir :

C1 = y0


C1 + C2 = v0
dont la solution y(t) est :
y(t) = y0 exp(t) + (v0 x y0 ) t exp(t) (8)
A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 11 / 17
 Équations homogènes à coecients constants Étude du polynôme caractéristique

b Observer le fait suivant : le cas d'une racine double n'est qu'un cas limite des
racines complexes conjuguées lorsque la partie imaginaire de celles-ci tend
vers zéro. C'est à dire
v0 x y0 
   

y(t) = lim y0 exp(t)cos(!t) + exp(t)sin(!t)

! !0 !
Calculons cette limite :
   
v xy 
y(t) = lim y exp(t)cos(!t) + exp(t)sin(!t) 0 0
!! !
0
0  
= y exp(t)lim
0
!!
cos(!t) + (v x y ) exp(t)lim sin!!t
0 !!
0 0
0
( )

sachant que :
lim
!! 
0
cos(!t) = 1

lim
!!0
sin !t
!
(
= lim
)

!!
t cos(!t) = t
0
(régle de L'Hôpital)

On obtient la solution nale :

y(t) = y0 exp(t) + (v0 x y0 ) t exp(t) (9)
A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 12 / 17
 L'oscillateur mécanique : un exemple récapitulatif

Un oscillateur mécanique constitué d'une masse m

dont le centre d'inertie G, repéré par la position x(t),
oscille verticalement autour de sa position d'équilibre
xeq . L'oscillateur est relié à un support xe par l'inter-
médiaire d'une ressort de raideur k et d'un amortisseur
de viscosité  . En appliquant la RFD et en introdui-
sant la fonction X (t) = x(t) x xeq , on obtient l'équa-
tion diérentielle régissant le mouvement mécanique
de l'oscillateur :
X 00 (t) + 2!0 X 0 (t) + !20 X (t) = 0


X (0) = X0 et X 0 (0) = V0
(10)
p p
où !0 = k=m et  = =2 mk désignent, respectivement, la pulsation propre et le
facteur d'amortissement du système.
Résoudre le problème dans les cas suivants :
1  =0
2 0< <1
3  =1
4  >1
Discuter les résultats physiques dans chaque cas.
A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 13 / 17
 Équations non homogènes à coecients constant

séance du 10/11/2021

A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 14 / 17

 Équations non homogènes à coecients constant Quelques cas particuliers

séance du 10/11/2021

A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 15 / 17

 Équations non homogènes à coecients constant Les phénomènes de Résonance

séance du 10/11/2021

A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 16 / 17

 Équations non homogènes à coecients constant La variation de la constante

séance de 10/11/2021

A. MEDELFEF (FPH-USTHB) Chapitre I.2 - EDO du Second Ordre 2021-2022 17 / 17