ChronoMath, une chronologie des MATHÉMATIQUES
à l'usage des professeurs de mathématiques, des étudiants et des élèves des lycées & collèges

Équations différentielles linéaires du 1er ordre           Second ordre
    
Méthode de la variation de la constante       Exercices , équation différentielle homogène

Au lycée, les équations différentielles du 1er ordre sont entrevues en Terminale en liaison avec les sciences physiques. Deux cas sont étudiés :

y' = ay et y' = ay + b    où a et b sont des constantes réelles données.

Ce sont des cas particuliers d'équations linéaires du 1er ordre. Les solutions sont respectivement :

y = keax   et   y = keax -b/a

  Par linéaire, (du latin linea = ligne droite), on veut signifier, depuis le 18è siècle avec d'Alembert dans son Encyclopédie (1752), que les inconnues (ici y et y') n'apparaissent qu'au premier degré : pas d'exposant entier ou fractionnaire. En algèbre les équations ax + b = 0 ou ax + by + c = 0 sont des équations linéaires respectivement à 1 et 2 variables. En termes de fonction, on distingue aujourd'hui les appellations f(x) = ax : fonction linéaire et f(x) = ax + b : fonction affine (terme moderne, du latin affinis = voisin, allié), dont les représentations graphiques sont des droites (des lignes).

On étudie ici le cas général. La fonction exponentielle correspondant à y' = y et y(0) = 1 étant supposée connue, par exemple en tant que fonction réciproque de la fonction logarithme népérien ou en tant que fonction f, non nulle, dérivable sur R, telle que f(x + y) = f(x)f(y) et f(0) = (0) :

Différentes approches de la fonction exponentielle :

Les fonctions xa(x), xb(x) et xc(x) étant continues sur un intervalle J de R, l'équation différentielle linéaire (par rapport à y et y') du premier ordre s'écrit sous la forme :

a(x).y' + b(x).y = c(x) , xJ         (e)

et se résout facilement si l'on en connaît une solution particulière. En effet, soit yo = φ(x) une solution de (e). Par différence et en posant Y = y - yo, on obtient :

a(x).Y' + b(x).Y = 0       (h)

Une telle équation est dite sans second membre ou, parfois, homogène : le second membre est nul.

A noter que le qualificatif homogène n'est pas très heureux car il entretient une confusion avec les équations différentielles homogènes. De plus, ne pas oublier que yo est une fonction et non une constante, raison pour laquelle on a posé yo = φ(x) pour éviter toute ambiguïté dans la suite de l'exposé.

  Un point important :  on peut supposer maintenant que x a(x) ne s'annule pas sur tout un sous-intervalle J' de J : l'équation ne serait plus différentielle sur J'... Les cas où a(x) = 0 pour des valeurs isolées de x se traitent en revenant à l'équation initiale : b(x)y = c(x). Donc, quitte à décomposer J en une réunion finie ou dénombrable d'intervalles sur lesquels a(x) est non nul, on peut supposer a(x) non nul sur J.

L'équation sans second membre (h) peut se résoudre par quadrature car l'on peut séparer les variables x et Y en écrivant, sous la condition Y(x) non nul en tout point de J :

Y'/Y = -b(x)/a(x)

Si f(x) désigne une primitive de -b(x)/a(x), on a successivement ln |Y| = f(x) + C, ln désignant le logarithme népérien . Puis |Y| = ef(x) + C = eCef(x) avec C une constante arbitraire.

La solution générale de (h) peut donc s'écrire :

Y = k.ef(x)

où f(x) désigne une primitive de -b(x)/a(x) et k une constante arbitraire. La solution de (e) s'obtient en ajoutant yo= φ(x) :

y = φ(x) + k.ef(x)

La valeur de k sera calculée compte tenu d'un condition initiale : valeur connue de y ou de y' en un point donné.

La solution Y précédemment trouvée ne s'annule effectivement pas sur J (comme supposé) si la constante k est choisie non nulle. Mais la question est de savoir si cette résolution fournit toutes les solutions de (h) ?

En remarquant, avec les notations précédentes, que xef(x) est effectivement une solution de (h) ne prenant jamais la valeur 0, on est en droit de poser g(x) = Y/ef(x), c'est à dire Y = g(x)ef(x) où g désigne une fonction dérivable sur J. En reportant dans (h), et en simplifiant par ef(x), on obtient :

a(x)g'(x) + a(x)g(x)(x) + b(x)g(x) = 0 pour tout x de J

Mais f '(x) = -b(x)/a(x). Il nous reste donc a(x)g'(x) = 0 pour tout x de J. C'est dire que g'(x) = 0 sur J, donc que g est constante sur J : la solution de (h) est donc bien Y = k.ef(x)

  Un exemple de résolution , autres exercices

Recherche d'une solution particulière de l'équation complète :     

On a vu ci-dessus que la résolution de (e) réside dans :

Concernant ce second point, On doit à Laplace une méthode astucieuse si le premier est résolu, dite de la variation de la constante :

Si Y est solution de (h) et ne s'annule pas sur J, on peut écrire qu'une solution particulière de (e) est y = k(x).Y(x) où k = y/Y est une fonction à déterminer. Or, k est dérivable et on a : y' = k'Y + kY'. En reportant dans (e) on obtient :

[a(x).Y' + b(x).Y].k + a(x).k'Y = c(x)

Le crochet est identiquement nul, donc :

et on obtiendra k par une simple quadrature (si tout se passe bien...).

 Remarquer que l'ensemble des solutions de (h), donc de la forme Y = k.ef(x), constitue un espace vectoriel de dimension 1. Par suite, si y1 et y2 sont des solutions particulières linéairement indépendantes de l'équation complète (e), alors y1 - y2 est solution de (h) et par conséquent, k(y1 - y2), où k est un nombre réel arbitraire, est la solution générale de cette équation.



Résoudre les équations différentielles :
1.
y' + y + 1 = ex       2. y' + y = cos x       3. xy' + 2y = ex , y(1) =0     4. y'' - 8y' + 15y = e3xcos2x     (on posera y = ze3x)

5. Soit l'équation différentielle du 1er ordre (x2 - x)y' - (x - 2)y = x4ex

Rep : y = kx2/(x - 1), k réel quelconque    décomposition en éléments simples d'une fraction rationnelle.
b) Donner toutes les solutions par la méthode de la variation de la constante.
Rep : k'(x) = xex, et finalement y = x2[ex + k/(x - 1)]

Autres exercices avec solution :

Résolution de l'équation2y' - y =  x et étude d'une courbe intégrale.

Résolution & étude des courbes intégrales de x2y' - y =  x2 - x + 1

Résolution & étude des courbes intégrales de x3y' + 3x2y = 1

Résolution par changement de variable et décomposition en éléments simples  de xy2 + x2(x + 1)yy' = 5 - 3x

Radioactivité période du radium #1 , dm = -km.dt

Radioactivité période du radium #2

Refroidissement d'un réacteur   1er ordre linéaire : y + 10y' = 16

Carbone 14  datation , premier ordre linéaire

Circuit RC  décharge d'un condensateur , premier ordre linéaire

 

Équation différentielle linéaire du second ordre :


© Serge Mehl - www.chronomath.com