Equation différentielle linéaire à coefficients constants

Introduction

On cherche les solutions aux équations du type :

a 1 dy dx ( x ) + a 0 y ( x ) = b ( x ) a_1 dy over dx (x) + a_0 y(x) = b(x)

a 0 a_0 et a 1 a_1 sont des constantes réelles et a 1 0 a_1 <> 0

Nous débuterons par l'étude des équations homogènes (i.e. avec un second membre nul). Puis en utilisant les résultats, nous verrons comment résoudre les équations avec second membre.

Equation homogène

Soit l'équation à résoudre :

a 1 dy dx ( x ) + a 0 y ( x ) = 0 a_1 dy over dx (x) + a_0 y(x) = 0

a 0 a_0 et a 1 a_1 sont des constantes réelles et a 1 0 a_1 <> 0 .

Supposons que y ( x ) 0 y(x) <> 0 , x forall x alors, nous pouvons écrire :

y ' ( x ) y ( x ) + a 0 a 1 = 0 y ' ( x ) y ( x ) = a 0 a 1 {y'(x)}over {y(x)} + a_0 over a_1 = 0 ~dlrarrow~ {y'(x)}over {y(x)} =- {a_0 over a_1 }

Intégrons par rapport à x x . La primitive du membre de gauche est : ln | y ( x ) | + A ln abs {y(x)} + A , donc nous obtenons :

ln | y ( x ) | = a 0 a 1 x + B ln abs {y(x)} = -{{a_0}over{a_1}}x + B

B B est une constante réelle qui tient compte des constantes d'intégration des deux membres.

Prenons l'exponentielle de la relation précédente, on obtient :

| y ( x ) | = e a 0 a 1 x + B abs {y(x)} = e^{-{{a_0}over{a_1}}x + B }

Or e α + β = e α e β e^{%alpha + %beta} = e^{%alpha} e^{%beta} , d'où :

| y ( x ) | = e B e a 0 a 1 x abs {y(x)} = e^B e^{-{{a_0}over{a_1}}x}

e B e^B est une constante réelle positive. Soit K K cette constante. On a alors :

y ( x ) = ± K e a 0 a 1 x y(x) = +-K e^{-{{a_0}over{a_1}}x}

L'expression ± K +-K représente une constante positive ou négative, notons C C cette constante réelle.

On a alors :

y ( x ) = C e a 0 a 1 x y(x) = C e^{-{{a_0}over{a_1}}x}

Nous obtenons donc une famille de fonctions. La constante C C pourra être calculée pour répondre à un cas précis à partir, par exemple, d'une condition particulière du type : y ( x 1 ) = D y(x_1) = D .

Exemple

Résoudre l'équation

y ' ( x ) + y ( x ) = 0 y'(x) + y(x) = 0

  1. en prenant la condition y ( 0 ) = 1 y(0) = 1

  2. en prenant la condition y ( 2 ) = 3 2 y(2) = 3 over 2

Solution

Partons de l'équation différentielle :

y ' ( x ) + y ( x ) = 0 y'(x) + y(x) = 0

en supposant que y ( x ) 0 y(x) <> 0 , x forall x , nous pouvons diviser cette équation par y ( x ) y(x) , soit :

y ' ( x ) y ( x ) + 1 = 0 {y'(x)} over {y(x)} +1 = 0

soit :

y ' ( x ) y ( x ) = 1 {y'(x)} over {y(x)} = -1

soit :

| y ( x ) | = e x + A abs {y(x)} = e^{-x+A}

soit :

y ( x ) = C e x y(x) = C e^{-x}

La constante C C dépend de la condition initiale :

  1. pour y ( 0 ) = 1 y(0) = 1 , on obtient C = 1 C=1 , soit y ( x ) = e x y(x) = e^{-x}

  2. pour y ( 2 ) = 3 2 y(2) = 3 over 2 , on obtient C = 3 2 e 2 C= 3 over 2 e^2 , soit y ( x ) = 3 2 e 2 x y(x) = 3 over 2 e^{2-x}

Equation avec second membre

L'équation différentielle du premier ordre à coefficients constants avec second membre s'écrit :

a 1 dy dx ( x ) + a 0 y ( x ) = b ( x ) a_1 dy over dx (x) + a_0 y(x) = b(x)

a 0 a_0 et a 1 a_1 sont des constantes réelles et a 1 0 a_1 <> 0 .

Propriété

(Méthode de variation de la constante) La méthode de résolution se décompose en trois points :

  1. on résout l'équation homogène, la solution générale est de la forme :

    y ( x ) = C e a 0 a 1 x y(x) = C e^{-{{a_0}over{a_1}}x} , a 1 0 a_1 <> 0 , C C in setR

  2. On considère que la constante C C dépend de x x (d'où le nom de la méthode : méthode de la variation de la constante), on pose :

    y ( x ) = C ( x ) e a 0 a 1 x y(x) = C(x) e^{-{{a_0}over{a_1}}x}

    on dérive, on reporte dans l'équation avec second membre et on résout l'équation différentielle d'inconnue C C .

  3. la solution générale de l'équation complète est obtenue en remplaçant C C par la famille de primitives trouvée au point précédent.

Application

Appliquons cette méthode sur le cas général. Calculons la dérivée de : y ( x ) = C ( x ) e a 0 a 1 x y(x) = C(x) e^{-{{a_0}over{a_1}}x} .

On obtient :

y ' ( x ) = C ' ( x ) e a 0 a 1 x a 0 a 1 C ( x ) e a 0 a 1 x y'(x) = C'(x) e^{-{{a_0}over{a_1}}x} - {a_0 over a_1} C(x) e^{-{{a_0}over{a_1}}x}

En injectant cette expression dans l'équation différentielle à résoudre, nous obtenons :

a 1 ( C ' ( x ) e a 0 a 1 x a 0 a 1 C ( x ) e a 0 a 1 x ) + a 0 ( C ( x ) e a 0 a 1 x ) = b ( x ) a_1 left ( C'(x) e^{-{{a_0}over{a_1}}x} - {a_0 over a_1} C(x) e^{-{{a_0}over{a_1}}x} right ) + a_0 left ( C(x) e^{-{{a_0}over{a_1}}x} right ) = b(x)

soit :

a 1 C ' ( x ) e a 0 a 1 x = b ( x ) a_1 C'(x) e^{-{{a_0}over{a_1}}x} = b(x)

soit :

C ' ( x ) = b ( x ) a 1 e a 0 a 1 x C'(x)= {b(x)} over a_1 e^{{{a_0}over{a_1}}x}

En intégrant, on obtient :

C ( x ) = b ( t ) a 1 e a 0 a 1 t dt + D C(x)= int {{b(t)} over a_1 e^{{{a_0}over{a_1}}t}} dt + D

Si on peut calculer explicitement cette intégrale, on en déduit la forme générale de la solution :

y ( x ) = ( b ( t ) a 1 e a 0 a 1 t dt + D ) e a 0 a 1 x y(x)= left (int {{b(t)} over a_1 e^{{{a_0}over{a_1}}t}} dt + D right ) e^{{{a_0}over{a_1}}x}

La valeur de D D est obtenue à partir d'une condition aux limites.

Exemple

Exemple

y ' ( x ) + y ( x ) = 1 y'(x) + y(x) = 1

avec la condition initiale y ( 0 ) = 2 y(0) = 2

L'équation homogène associée est :

y ' ( x ) + y ( x ) = 0 y'(x) + y(x) = 0

Sa solution a été déterminée dans l'exercice précédent, nous avons donc :

y ( x ) = B e x y(x) = B e^{-x}

Supposons à présent que B B dépend de la variable x x . On a alors : y ( x ) = B ( x ) e x y(x) = B(x) e^{-x} .

y ' ( x ) = B ' ( x ) e x B ( x ) e x y'(x) = B'(x) e^{-x}-B(x) e^{-x}

En injectant ces expressions dans l'équation différentielle, nous obtenons :

B ' ( x ) e x B ( x ) e x + B ( x ) e x = 1 B'(x) e^{-x}-B(x) e^{-x}+ B(x) e^{-x} = 1

Soit :

B ' ( x ) e x = 1 B ' ( x ) = e x B'(x) e^{-x} = 1 ~drarrow~ B'(x) = e^{x}

Une primitive particulière de e x e^x peut s'écrire :

e t dt = e x int e^t dt = e^{x}

d'où

B ( x ) = e x + C B(x) = e^{x} + C

L'expression générale de la solution de l'équation différentielle est donc :

y ( x ) = [ e x + C ] e x = 1 + C e x y(x) = [e^x + C]e^{-x} = 1+C e^-x

La constante C C peut être calculée à partir de la condition initiale y ( 0 ) = 2 y(0) = 2 . On obtient alors :

y ( 0 ) = 1 + C e 0 = 1 + C y(0) =1+C e^0 =1+C

d'où C = 1 C=1 . La solution est donc :

y ( x ) = 1 + e x y(x) = 1+ e^-x