En mathématiques, et plus précisément en analyse, la méthode de variation des constantes (ou méthode de Lagrange) est une méthode de résolution des équations différentielles. Elle permet en particulier de déterminer les solutions d'une équation différentielle avec second membre, connaissant les solutions de l'équation homogène (c'est-à-dire sans second membre) associée.
La méthode a été inventée par le mathématicien et physicien Pierre-Simon de Laplace, pour la résolution des équations différentielles linéaires. Elle tire son nom de ce que, pour l'essentiel, elle consiste à chercher les solutions sous une forme analogue à celle déjà trouvée pour une équation associée plus simple, mais en remplaçant la ou les constantes de cette solution par de nouvelles fonctions inconnues.
Cas du premier ordre
Pour une équation différentielle linéaire d'ordre 1, si la solution générale de l'équation homogène
![{\displaystyle az'+bz=0}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95a478165b1a666e4197283bb58fb69a2f7a431c)
est
![{\displaystyle z_{K}(x)=Kz_{1}(x),K\in \mathbb {R} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/11a4dba31dd1167f77360a7f01cbe60f85229efa)
on cherche celle de
![{\displaystyle ay'+by=c}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e645b259dd199a965fe39603f277d505b4a8d89b)
sous la forme
![{\displaystyle y(x)=k(x)z_{1}(x).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c349e6352a13ef5bb8feb3223468dd0a70d35be9)
En reportant dans l'équation initiale, on obtient une équation équivalente à l'équation initiale mais portant sur k :
![{\displaystyle k'={\frac {c}{az_{1}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8dce8171f27f2217b122a9e186382aa8d210a9b2)
En notant k0 une primitive de la fonction c/(az1), la solution générale k s'exprime sous la forme
ce qui permet de remonter à l'expression de la solution générale yK = y0 + zK :
Pour expliciter z1 puis k0, il faut réaliser deux calculs de primitives. De ce fait, la solution ne s'exprime le plus souvent pas à l'aide des fonctions usuelles (voir à ce sujet le théorème de Liouville).
Cas du second ordre
Pour une équation différentielle linéaire d'ordre deux, mise sous la forme
:
Notons
et
deux solutions formant une base des solutions de l'équation homogène. On cherchera une solution particulière y sous la forme
(d'où le nom de « méthode de variation des constantes »).
On montre que si les fonctions
et
vérifient le système suivant
![{\displaystyle {\begin{cases}\lambda '(x)\cdot y_{1}(x)+\mu '(x)\cdot y_{2}(x)=0\\\lambda '(x)\cdot y'_{1}(x)+\mu '(x)\cdot y'_{2}(x)=d(x)\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aff2f6700e5b6e518c5017dc454dfaa941cc410e)
alors la fonction y ci-dessus est une solution particulière.
Remarque:
1) Puisque le Wronskien ne s'annule pas,
,
et
s'obtiennent en utilisant la règle de Cramer.
Cas général
Pour une équation différentielle linéaire d'ordre n avec second membre, on cherchera une solution particulière
combinaison linéaire d'un système fondamental de solutions
, c.-à-d. d'une base de l'espace vectoriel des solutions de l'équation homogène. Les coefficients
de la combinaison linéaire sont maintenant des fonctions que l'on cherche à déterminer. C'est une simple généralisation du cas n=2, cependant il existe une reformulation matricielle.
Une équation différentielle ordinaire linéaire non homogène s'écrit de manière générale
où
est la dérivée k-ième de
. On suppose au préalable disposer de n solutions linéairement indépendantes
de l'équation homogène
. En regroupant dans une colonne les dérivées successives de chaque solution
, on forme la matrice suivante
. L'indépendance des n solutions peut se vérifier pour chaque x fixé, en calculant le déterminant de cette matrice qui ne doit pas s'annuler, cf. wronskien.
La méthode de variation de la constante consiste à chercher une solution particulière de (1) sous la forme
avec des fonctions
au moins
. En réalité, on introduit à cette étape beaucoup trop d'inconnues et il faut imposer une égalité similaire sur les dérivées supérieures:
ce qui revient, par récurrence, à imposer
. Cependant, la dérivée n-ième est laissée telle quelle ; on obtient en dérivant (5) avec k := n – 1 :
. Remarque : on comprend mieux l'origine des conditions (6) et (7) dans la formulation matricielle, équation (16).
En insérant maintenant (4), (5) et (7) dans l'équation non homogène (1), on obtient
. Or puisque les
sont solution de (2), la première somme disparaît, laissant
. Théorème : Si les
vérifient les équations différentielles (6) et (8) alors l'expression (4) est solution de (1).
(6) et (8) peuvent se réécrire sous forme matricielle
ce qui se résout comme dans le cas du second ordre avec la règle de Cramer. Notons
le wronskien, que l'on a supposé non nul, alors
Soit
. Notons
le vecteur colonne
Une fonction
est solution de (1) si et seulement si
satisfait
Notons
la matrice du membre de droite. On a transformé une équation d'ordre n en une équation d'ordre 1.
La résolvante de l'équation homogène associée est l'application
qui envoie la valeur
de toute solution au point x1 à sa valeur
au point x2 (justification de l'existence, cf. Théorème de Cauchy-Lipschitz), c.-à-d. pour toute solution
de l'équation homogène associée à (11)
. Si l'on regroupe n solutions
linéairement indépendantes de l'équation homogène dans la matrice
de (3), on a naturellement
ainsi que
. Remarque:
- (13) et (14) restent valables pour toute famille
de solutions, linéairement indépendantes ou pas. - Si l'on connait explicitement un système fondamental de solutions
, alors on peut expliciter la résolvante à partir de (14). En effet, par la formule de Liouville, si
est inversible pour un certain x0, alors elle l'est pour tout x ; on peut ainsi écrire
.
Attention : Nous allons utiliser, sans la démontrer, l'invariance de la résolvante par translation, en particulier que
(abus de notation) et que vu comme fonction à une variable, c'est un « groupe à 1 paramètre » partout où il est défini, c.-à-d.
et
lorsque cela est défini. Par ailleurs,
satisfait aussi (13).
Dans cette formulation, la méthode de variation des constantes consiste à faire le "changement de variable"
(11) se réécrit donc
. Or puisque
satisfait aussi (13), il ne reste que
ce qui est équivalent à (9), avec ici la dépendance explicite aux valeurs des solutions fondamentales en x0.
En remarquant que
(matrice identité) et avec (15), on intègre le vecteur composante par composante
. En réutilisant (15) et
partout où cela est défini, on obtient finalement
. Exemple d'application à la physique
L'équation différentielle du second ordre à coefficients constants
intervient en physique dans l'étude des systèmes oscillants à un degré de liberté, lorsque l'excitation (force, courant…) appliquée au système oscillant est nulle.
La méthode de l'équation caractéristique (découverte par Euler) donne la solution de cette équation différentielle homogène, qui est une combinaison linéaire de fonctions exponentielles (complexes).
Lorsque l'on applique une excitation
, l'équation devient :
. La méthode de variation de la constante permet d'en trouver la solution générale.
Liens externes
Portail de l'analyse