Équation différentielle de Bernoulli

{{#invoke:Bandeau|ébauche}}

Une équation différentielle de Bernoulli est une équation différentielle du premier ordre de la forme <math> y'(x) + a(x)y(x)= b(x)y^m(x)~</math>.

Description

On considère donc l'équation :

<math> y'(x) + a(x)y(x)= b(x)y^m(x)~</math>

où Modèle:Mvar est un réel différent de 0 et 1 et où Modèle:Mvar et Modèle:Mvar sont des applications définies sur un intervalle ouvert Modèle:Mvar de <math>\mathbb R</math> et à valeurs réelles. En général, <math>m</math> est un entier naturel, mais on peut prendre Modèle:Mvar réel à condition de chercher Modèle:Mvar à valeurs strictement positives. En général, Modèle:Mvar et Modèle:Mvar sont des fonctions continues.

Cette forme d'équation a été proposée par Jacques Bernoulli en 1695 et résolue un an plus tard par Leibniz grâce à un changement de fonction qui ramène à une équation différentielle linéaire.

En supposant que la fonction Modèle:Mvar est à valeurs strictement positives sur l'intervalle Modèle:Mvar, on peut diviser l'équation par Modèle:Math et on obtient

<math> \frac{y'(x)}{y^m(x)} + a(x)\frac{1}{y^{m-1}(x)} = b(x)</math>

On pose

<math> u(x)=\frac1{y^{m-1}(x)}=y^{1-m}(x)</math>

donc

<math>u'(x)=(1-m)y^{-m}(x)y'(x)=\frac{(1-m)y'(x)}{y^m(x)}.</math>

L'équation de Bernoulli sur Modèle:Mvar équivaut donc à l'équation différentielle linéaire d'ordre un sur Modèle:Mvar :

<math>\frac1{1-m}u'(x) + a(x)u(x) = b(x)</math>

dont la solution générale est

<math> u(x) = \exp\left(-(1-m)\int a(t)\mathrm dt\right)\left(C + (1-m)\int b(t)\exp\left((1-m)\int a(s)\mathrm ds\right)~\mathrm dt\right),</math>

ce qui donne pour la fonction <math>y = u^{1/(1-m)}</math> :

<math> y(x) = \exp\left(- \int a(t)\mathrm dt\right)\left(C + (1-m)\int b(t)\exp\left((1-m)\int a(s)\mathrm ds\right)~\mathrm dt\right)^{1\over1-m}</math>

La solution de cette équation qui passe par le point Modèle:Math est la fonction Modèle:Mvar définie par :

<math>y(x)=y_0\exp\left(-\int_{x_0}^xa(t)\mathrm dt\right)\left(1+(1-m)y_0^{m-1}\int_{x_0}^xb(t)\left[\exp\left(-\int_{x_0}^ta(s)\mathrm ds\right)\right]^{m-1}~\mathrm dt\right)^{\frac1{1-m}}</math>.

Des solutions peuvent être cherchées parmi les fonctions qui ne sont pas partout positives sur leur domaine de définition, mais alors de nombreuses précautions doivent être prises quant aux domaines de validité des solutions.

Voir aussi

Sources et bibliographie

• {{#invoke:Langue|indicationDeLangue}} John E. Sasser, Modèle:Pdf - History of ordinary differential equations - The first hundred year
• J. Lelong-Ferrand et J.-M. Arnaudiès, Cours de mathématiques, tome 4

Article connexe

• Jacques Bernoulli