PCSI5 Chapitre6
PCSI5 Chapitre6
PCSI5 Chapitre6
Équations différentielles
Mathieu Mansuy - Professeur de Mathématiques en supérieures PCSI au Lycée Saint Louis (Paris)
[email protected]
PCSI5 Lycée Saint Louis, Paris
On dit que f est solution sur I de (E) si f est dérivable sur I et si pour tout t ∈ I, f 0 (t)+a(t)f (t) = b(t).
Remarques.
Soit f une solution de (E) sur I, f est donc dérivable sur I. Elle est en fait de classe C 1 sur I puisque
f 0 (t) = b(t) − a(t)f (t) qui est par hypothèse continue sur I.
Si on a une équation différentielle de la forme α(t)y 0 (t) + β(t)y(t) = γ(t), on peut la ramener à la forme
précédente (dite normalisée) en divisant par α(t). Il faut en revanche se placer sur un intervalle où la
fonction α ne s’annule pas.
La tension U aux bornes d’un condensateur de capacité C en série avec une résistance R et un générateur de
force électromotrice E(t) vérifie l’équation différentielle :
RCU 0 + U = E(t).
L’intensité I dans un circuit RL constitué d’une résistance R, d’une bobine d’inductance L et d’un générateur
de force électromotrice E(t) vérifie l’équation différentielle :
LI 0 + RI = E(t).
La proportion y de carbone 14 dans le carbone total des êtres vivants est constante. Après la mort, elle
diminue de 1/8000 par an et vérifie donc l’équation différentielle suivante (où le temps est exprimé en années) :
1
y0 + y = 0.
8000
On peut ainsi dater la mort d’un être vivant grâce à cette relation (datation au carbone 14).
On appelle équation différentielle homogène (ou sans second membre) associée à (E), l’équation
Dans ce paragraphe, on étudie et on résout l’équation homogène (E0 ) associée à l’équation différentielle linéaire
du premier ordre (E).
2
PCSI5 Lycée Saint Louis, Paris
Propriété 1
∀λ, µ ∈ K, ∀f, g ∈ S0 , λf + µg ∈ S0 .
Preuve.
Propriété 2
d at
∀t ∈ R, eat (y 0 (t) + ay(t)) == (e y(t)) = 0.
dt
Ceci est équivalent à l’existence de λ ∈ K tel que eat y(t) = λ, ou y(t) = λeat .
Propriété 3
t 7→ y(t) = λe−A(t)
Preuve. L’idée de la démonstration est la même que précédemment : on multiplie l’équation y 0 + a(t)y = 0
par eA(t) (6= 0 sur R) :
d
∀t ∈ R, eA(t) (y 0 (t) + a(t)y(t)) == (eA(t) y(t)) = 0.
dt
Ceci est équivalent à l’existence de λ ∈ K tel que eA(t) y(t) = λ, ou y(t) = λe−A(t) .
Remarques.
Cette équation est une équation différentielle linéaire d’ordre un, homogène (ou sans second membre), non
normalisée. Comme sin ne s’annule pas sur I, on peut la normalisée : cette équation équivaut donc à
cos t
y0 − y = 0.
sin t
Une primitive de t 7→ − cos t
sin t est donnée par t 7→ − ln(|sin(t)|) sur I. On en déduit que les solutions de l’équation
sont de la forme :
t 7→ λeln(sin(t)) = λ sin t avec λ ∈ R.
3
PCSI5 Lycée Saint Louis, Paris
Cette équation est une équation différentielle linéaire d’ordre un, homogène (ou sans second membre), non
normalisée. On la résout sur R∗+ ou sur R∗− . Plaçons nous sur l’un de ces intervalles. On peut alors la
normalisée :
1
y 0 − y = 0.
2t
1
Une primitive de t 7→ − 2t est donnée par t 7→ − 12 ln(|t|). On en déduit que les solutions de l’équation sont de
la forme :
1
t 7→ λe 2 ln(|t|) .
√ √
Ainsi sur R∗+ les solutions sont de la forme t 7→ λ1 t. Sur R∗− les solutions sont de la forme t 7→ λ2 −t.
Soit (E) : y 0 + a(t)y = b(t) une équation différentielle linéaire du premier ordre avec a, b : I → K
continues.
Alors toute solution de (E) s’obtient en ajoutant à une solution particulière yP de (E) une solution
de l’équation homogène associée (E0 ).
Ainsi l’ensemble S des solutions de (E) satisfait :
S = yp + S0 .
Il s’agit d’une équation différentielle linéaire d’ordre 1 à coefficients continus non constants, avec second membre.
On montre que l’équation homogène associée a pour solution générale x 7→ λe− arctan x , λ ∈ R. La fonction
7 1 + λe− arctan x ,
constante égale à 1 est solution particulière évidente, donc la solution générale de (E) est x →
λ ∈ R.
On donne maintenant des méthodes générales pour obtenir une solution particulière d’une équation différentielle
linéaire.
4
PCSI5 Lycée Saint Louis, Paris
(λf1 + µf2 )0 + a(t)(λf1 + µf2 ) = λ(f10 + a(t)f1 ) + µ(f20 + a(t)f2 ) = λb1 + µb2 .
I Cette proposition peut être utile dans la pratique : pour résoudre une équation différentielle linéaire, on peut
éventuellement scinder le second membre afin de se ramener à plusieurs équations différentielles dont les second
membres seront plus simples à traiter.
Alors l’équation (E) a au moins une solution particulière de la forme t → eλt tm Q(t) où Q est une
fonction polynomiales à coefficients dans K telle que deg(Q) = deg(P ) et :
m = 0 si λ + a 6= 0.
m = 1 si λ + a = 0.
Soit y 0 (t) + a(t)y(t) = b(t) une équation différentielle linéaire du premier ordre.
Les solutions de l’équation homogène sont du type t 7→ λe−A(t) où λ est une constante.
Les solutions de l’équation complète sont du type t 7→ λ(t)e−A(t) où λ est de classe C 1 sur I.
I On utilisera donc la méthode suivante, due à Lagrange, pour résoudre pratiquement l’équation différentielle
y 0 (t) + a(t)y(t) = b(t) :
On résout l’équation homogène y 0 + a(t)y = 0 : les solutions sont t 7→ λ exp(−A(t)) où λ ∈ K est une
constante quelconque et A une primitive de a.
On résout l’équation complète y 0 (t) + a(t)y(t) = b(t) en cherchant une solution particulière sous la forme
t 7→ λ(t) exp(−A(t)), où λ est une fonction de classe C 1 sur I.
5
PCSI5 Lycée Saint Louis, Paris
Propriété 8
Pour (t0 , y0 ) ∈ I × K, le problème de Cauchy associé à la condition initiale y(t0 ) = y0 admet une
unique solution sur I.
Remarque. On a montré en particulier dans la preuve de ce résultat que toute solution de (E) est bien de la
forme t 7→ λ(t)e−A(t) , ce qui justifie la méthode de variation de la constante.
√
Exemple. Résoudre y 0 − 2x 1
y = x sur I = R∗+ , avec pour condition initiale y(1) = 1.
Par la méthode
√ de la variation de la constante, on cherche une solution particulière
√ √de (E) sous la forme
x 7→ λ(x) x, avec λ dérivable sur I. En reportant dans l’équation, il vient√λ0 (x) x = x donc λ0 (x) = 1 puis
on peut prendre λ√: x 7→ x et une solution particulière de (E) est x 7→ x x. La solution générale de (E) est
donc x 7→ (λ + x) x, λ ∈ R.
Pour trouver l’unique solution du problème considéré, il√faut résoudre l’équation (λ + 1) = 1 i.e. λ = 0. Ainsi
la solution du problème de Cauchy considéré est x 7→ x x.
Propriété 9
(1) Une seule courbe intégrale définie sur I passe par tout point (t0 , y0 ) ∈ I × R.
6
PCSI5 Lycée Saint Louis, Paris
Preuve. Le premier point est une reformulation du résultat précédent. Montrons le deuxième point : supposons
que deux courbes intégrales se coupent en M0 (t0 , y0 ). Notons y1 et y2 les solutions de (E) sur I correspondantes.
Elles sont donc toutes deux solutions du même problème de Cauchy :
(
(E)
y(t0 ) = y0
Remarque. La seule solution de (E0 ) : y 0 + a(t)y = 0 telle que y(t0 ) = 0 est la fonction nulle. En effet, la
fonction nulle est solution de (E0 ). Et par la proposition précédente, c’est l’unique solution telle que y(t0 ) = 0.
En particulier, la seule solution de (E0 ) qui s’annule est la fonction nulle.
où a, b, c sont continues sur un intervalle I, mais où a est susceptible de s’annuler sur I. Nous allons voir sur
des exemples que les résultats précédents, comme le théorème de Cauchy, peuvent alors tomber en défaut.
Supposons par exemple que a s’annule en un unique point t0 de I. On procèdera alors comme suit.
On commence par résoudre l’équation sur les intervalles Ig =] − ∞, t0 [∩I et Id =]t0 , +∞[∩I où a ne
s’annule pas.
On cherche ensuite une solution f de (E) sur I. On remarque que f est solution de (E) sur Ig . On en
connait donc la forme explicite f = fg . De même sur l’intervalle(Id , on connait la forme fd de f . La
fg (t) si t ∈ Ig
fonction f étant solution de (E) sur I, la fonction t ∈ I \ {t0 } 7→ doit être :
fd (t) si t ∈ Id
On va voir sur des exemples que ces conditions sont satisfaites dans certains cas, pas dans d’autres.
C’est une équation différentielle linéaire d’ordre un à coefficients variables et continus, non normalisée. On
commence par la résoudre sur des intervalles sur lesquels t 6= 0, soit sur R∗+ et sur R∗− . On montre alors alors
que :
Ici on a :
lim λt2 + t3 = 0 et lim− µt2 + t3 = 0.
x→0+ x→0
7
PCSI5 Lycée Saint Louis, Paris
On étudie la dérivabilité :
f (t) − f (0) λt2 + t3
lim = lim =0
t→0+ t t→0 + t
f (t) − f (0) µt2 + t3
lim− = lim− = 0.
t→0 t t→0 t
Là encore, cette condition est bien satisfaite pour tout λ, µ ∈ R.
Finalement, on obtient que les solutions de (E) sur R sont les fonctions de la forme :
(
λt2 + t3 si t ≥ 0
t 7→
µt2 + t3 si t < 0
pour tout λ, µ ∈ R.
(
(E)
Remarque. Le problème de Cauchy a donc une infinité de solutions. À l’opposé, le problème de
y(0) = 0
(
(E)
Cauchy n’a aucune solution. On voit ici que le théorème de Cauchy tombe en défaut en t = 0, là
y(0) = 1
où le coefficient t de y 0 s’annule.
C’est une équation différentielle linéaire d’ordre un à coefficients variables et continus, non normalisée. On
commence par la résoudre sur des intervalles sur lesquels 1 − t 6= 0, soit sur ]1, +∞[ et sur ] − ∞, 1[. On montre
alors alors que :
λ t2 µ t2
S]1,+∞[ = {t → + , λ ∈ R} et S]−∞,1[ = {t → + , λ ∈ R}.
1 − t 2(1 − t) 1 − t 2(1 − t)
On cherche à présent les solutions de (E) sur R. Soit f une solution l’une d’elles. On a alors :
λ t2
∃λ ∈ R, f (t) = + pour tout t > 1,
1 − t 2(1 − t)
µ t2
∃µ ∈ R, f (t) = + pour tout t < 1.
1 − t 2(1 − t)
Prenons t = 1 dans l’équation, on obtient : 0 × f 0 (1) − f (1) = 1, soit f (1) = −1.
La fonction f doit être continue et dérivable en 1 :
λ t 2 +∞ si λ > −1/2
lim + = −∞ si λ < −1/2 .
x→1+ 1 − t 2(1 − t)
−1 si λ = −1/2
Ainsi on a nécessairement λ = −1/2. De même, µ = −1/2. On obtient alors f (t) = − 21 − 2t sur R. En particulier
f est bien dérivable sur R.
Finalement, on obtient qu’il y a une unique solutions de (E) sur R, qui est :
1 t
t 7→ − − .
2 2
8
PCSI5 Lycée Saint Louis, Paris
( (
(E) (E)
Le problème de Cauchy n’a pas de solution dans ce cas. Le problème de Cauchy a
y(1) = 0 y(1) = −1
pour sa part une unique solution.
2.1 Généralités
Définition.
Soient I un intervalle de R, a, b, c ∈ K (a 6= 0), et f : I → K une fonction continue. On appelle équation
différentielle linéaire du second ordre à coefficients constants :
ay 00 + by 0 + cy = f (t) (E)
La fonction f s’appelle le second membre. On appelle solution sur I de (E) toute fonction y : I → K
deux fois dérivable sur I telle que :
La courbe représentative d’une solution y sur I est alors appelée courbe intégrale sur I.
Remarque. Soit y une solution de (E) sur I, alors f est de classe C 2 sur I. En effet, elle est supposée deux
fois dérivable sur I. De plus f 00 est continue sur I puisque y 00 (t) = a1 f (t) − ab y 0 (t) − ac y(t).
L’intensité I dans un circuit LC constitué d’un condensateur de capacité R et d’une bobine d’inductance L
vérifie l’équation différentielle :
1
I 00 + I = 0.
LC
ay 00 + by 0 + cy = f (t) (E)
On appelle équation différentielle homogène (sans second membre) associée à (E) l’équation :
ay 00 + by 0 + cy = 0. (E0 )
Propriété 10
L’ensemble des solutions S0 de (E0 ) est stable par combinaison linéaire, c’est à dire :
∀λ, µ ∈ K, ∀f, g ∈ S0 , λf + µg ∈ S0 .
Preuve.
9
PCSI5 Lycée Saint Louis, Paris
Propriété 11
Définition.
On appelle équation caractéristique associée à (E) l’équation ar2 + br + c = 0 d’inconnue r ∈ C.
ay 00 + by 0 + c = 0 (E0 )
n o
S0 = t ∈ R 7→ (λ + µt)ert ; (λ, µ) ∈ C2 .
Preuve. Soit r une solution de l’équation caractéristique (qui existe toujours sur C). La fonction t 7→ ert est
solution de (E0 ) sur R.
Posons maintenant y(t) = u(t)ert avec u une fonction de classe C 2 sur R. On montre alors que y est solution de
(E0 ) sur R si et seulement si u est solution de l’équation :
10
PCSI5 Lycée Saint Louis, Paris
ay 00 + by 0 + c = 0 (E0 )
Si ∆ < 0 : l’équation caractéristique admet deux solutions complexes non réelles α ± iβ, (avec
(α, β) ∈ R2 ). La solution générale de (E) est x 7→ eαx (λ cos(βx) + µ sin(βx)), (λ, µ) ∈ R2 .
Remarque. Comme vu dans le chapitre précédente, le troisième point peut se réécrire sous la forme x 7→
eαx (A cos(βx + φ)), avec (A, φ) ∈ R2 .
Preuve. Notons d’abord que les solutions à valeurs réelles de cette équation sont des solutions à valeurs
complexes caractérisées par le fait qu’elles sont égales à leur conjugué). On a donc selon le signe du discriminant
:
si ∆ > 0, l’équation caractéristique a deux racines distinctes r < r0 . Les solutions complexes sont les
fonctions de la forme : 0
t 7→ λert + µer t
avec (λ, µ) ∈ C2 . Ces solutions sont réelles si pour tout t ∈ R
0 0
λert + µer t = λ̄ert + µ̄er t
t 7→ λe(α+iβ)t + µe(α−iβ)t
Remarque. On dit que l’ensemble des solutions de (E0 ) forme un plan vectoriel, puisque celles-ci sont combi-
naisons linéaires de deux fonctions non proportionnelles (que ce soit dans le cas réel ou complexe).
11
PCSI5 Lycée Saint Louis, Paris
ay 00 + by 0 + cy = f (t) (E)
Alors toute solution de (E) s’obtient en ajoutant à une solution particulière yP de (E) une solution
de l’équation homogène associée (E0 ). Ainsi l’ensemble S des solutions de (E) satisfait :
S = yp + S0 .
Preuve.
Remarque.
Pour obtenir l’ensemble des solutions de (E), il faudra donc déterminer une solution particulière de (E),
et lui ajouter l’ensemble des solutions de (E0 ) (dont on a vu la forme).
L’ensemble S est donc la somme d’un point yp et de combinaisons linéaires de deux fonctions non propor-
tionnelles. On dit alors que S est un plan affine.
Si y1 est solution sur I de (E1 ), y2 est solution sur I de (E2 ), alors pour tout λ, µ ∈ K, λy1 + µy2
est solution sur I de :
ay 00 + by 0 + cy = λf1 + µf2 .
Preuve.
On considère des nombres a, b, c, λ ∈ K (où a 6= 0), une fonction polynomiale P à coefficients dans
K et l’équation différentielle linéaire :
Alors l’équation (E) a au moins une solution particulière de la forme t 7→ eλt tm Q(t) où Q est une
fonction polynomiales à coefficients dans K telle que deg(Q) = deg(P ) et :
m = 0 si λ n’est pas racine du polynôme caractéristique,
12
PCSI5 Lycée Saint Louis, Paris
L’équation caractéristique de (E) est r2 − 3r + 2 = 0, dont les solutions sont 1 et 2. Ainsi la solution générale
de l’équation homogène associée à (E) est x 7→ λe2x + µex , (λ, µ) ∈ R2 .
Comme 2 est solution simple de l’équation caractéristique, on cherche une solution particulière sous la forme
y : x 7→ xQ(x)e2x , avec Q fonction polynomiale de degré 1, donc il existe (a, b) ∈ R2 telles que Q(x) = (ax+b). y
est deux fois dérivable comme produit de fonctions qui le sont, et y 0 (x) = (2ax+b)e2x +2(ax2 +bx)e2x = (2ax2 +
2bx + 2ax + b)e2x , puis y 00 (x) = (4ax + 2b + 2a)e2x + 2(2ax2 + 2bx + 2ax + b)e2x = (4ax2 + 4bx + 8ax + 4b + 2a)e2x .
En reportant dans (E), il vient
xex = (4ax2 + 4bx + 8ax + 4b + 2a)e2x − 3(2ax2 + 2bx + 2ax + b)e2x + 2(ax2 + bx)e2x = 2axe2x + (b + 2a)e2x
1
donc 2a = 1 et b + 2a = 0 par unicité des coefficients d’un polynôme, i.e. a = 2 et b = −1.
Une solution particulière est donc x 7→ ( 21 x2 − x)e 2x
et la solution générale de (E) est x 7→ ( 21 x2 − x)e2x +
λe2x + µex , (λ, µ) ∈ R2 .
ay 00 + by 0 + cy = f (t). (E)
Le problème de Cauchy associé à un triplet (t0 , y0 , y00 ) ∈ I × K2 est la recherche des solutions de (E)
vérifiant les deux conditions initiales y(t0 ) = y0 , y 0 (t0 ) = y00 .
Propriété 17
Pour (t0 , y0 , y00 ) ∈ I × K2 , le problème de Cauchy associé aux conditions initiales y(t0 ) = y0 , y 0 (t0 ) =
y00 admet une solution unique sur I.
Propriété 18
ay 00 + by 0 + cy = f (t). (E)
(1) Une seule courbe intégrale définie sur I passe par le point M0 (t0 , y0 ) de I × R avec une tangente
de pente donnée y00 ∈ R en ce point.
(2) Deux telles courbes intégrales distinctes ne peuvent se couper dans I × R avec une tangente
commune.
Preuve. Le premier point est une reformulation du résultat précédent. Montrons le deuxième point : supposons
que deux courbes intégrales se coupent en M0 (t0 , y0 ) avec une tangente commune de pente y00 ∈ R. Notons y1
et y2 les solutions de (E) sur I correspondantes. Elles sont donc toutes deux solutions du même problème de
Cauchy :
(E)
y(t0 ) = y0
0
y (t0 ) = y00
Remarque. La seule solution de (E0 ) : ay 00 + by 0 + cy = 0 telle que y(t0 ) = y 0 (t0 ) = 0 est la fonction nulle. En
effet, la fonction nulle est solution de (E0 ). Et par la proposition précédente, c’est l’unique solution telle que
y(t0 ) = y 0 (t0 ) = 0.
13
PCSI5 Lycée Saint Louis, Paris
1. En utilisant le principe fondamental de la dynamique, montrer que l’équation du mouvement est régie par
l’équation différentielle :
α
(
00 0 2
2λ = m le coefficient d’amortissement de l’oscillateur
z + 2λz + ω0 z = 0 avec q
k
ω0 = m qui représente sa pulsation propre
Si l’origine est en haut du ressort et l’axe vertical orienté vers le bas, la loi fondamentale de la dynamique
donne :
mz 00 = −αz 0 − K(z − L) + mg ou mz 00 + αz 0 + Kz = KL + mg.
Cette équation équation différentielle a un second membre constant, et on vérifie que z0 = L + mg/K est
solution constante (correspondant à l’équilibre). En portant l’origine en z0 , i.e. en posant Z = z − z0 , on
obtient :
mZ 00 + αZ 0 + KZ = 0 ou Z 00 + 2λZ 0 + ω02 Z = 0.
2. On précise les conditions initiales: z(0) = 10 et z 0 (0) = 0 (c’est à dire qu’on suppose l’oscillateur lâché
sans vitesse initiale après avoir tiré d’une longueur 10 à partir de sa position d’équilibre).
(a) Avec λ = 0, ω0 = 1, déterminer l’expression explicite du mouvement.
(b) Avec λ = ω0 = 1, déterminer l’expression explicite du mouvement.
√
(c) Avec λ = 1, ω0 = 5, déterminer l’expression explicite du mouvement.
Interprétation géométrique. Dans ce dernier cas, on observe un comportement transitoire qui présente
quelques oscillations périodiques avant de retrouver un état stable.
√
3. Supposons toujours λ = 1, ω0 = 5. On impose à présent à l’extrémité supérieure du ressort un mou-
vement verticale sinusoı̈dale de sorte que son ordonnée à l’instant t est f (t) = 5 cos(t). Déterminer une
expression explicite du mouvement. Déterminer le comportement asymptotique des solutions (état stable).
Vocabulaire. Le second membre d’une équation différentielle ay 00 + by 0 + cy = f (t) est également appelé signal
d’entrée et les solutions de l’équation définiront la réponse du système étudié. On parle de régime libre
si f est la fonction nulle, et de régime forcé sinon. Dans le cas d’un oscillateur harmonique amorti en régime
libre, on parle de :
régime apériodique quand λ > ω0 ,
14