Bernoulli-DGL

Bernoulli-DGL

Studium
y′+a(x) y=b(x) yn  →  z = y1−n    z′+(1−n) a(x) z=(1−n) b(x)y' + a(x)\,y = b(x)\,y^n \;\xrightarrow{\;z \,=\, y^{1-n}\;}\; z' + (1-n)\,a(x)\,z = (1-n)\,b(x)

Voraussetzung: n∈R∖{0,1},n \in \mathbb{R} \setminus \{0, 1\}, y>0y > 0

yy
gesuchte Funktion y(x)
y′y'
Ableitung von y nach x
a(x),b(x)a(x), b(x)
stetige Koeffizientenfunktionen
nn
Exponent der Nichtlinearität
zz
neue Funktion z = y^(1−n), erfüllt eine lineare DGL
z′z'
Ableitung von z nach x

Die Substitution z = y^(1−n) macht aus der nichtlinearen Bernoulli-DGL eine lineare DGL erster Ordnung für z.

Erklärung

Eine Bernoulli-Differentialgleichung sieht aus wie eine lineare DGL erster Ordnung, hat rechts aber den Faktor yny^n:

y′+a(x) y=b(x) yny' + a(x)\,y = b(x)\,y^n

Für n=0n = 0 ist sie linear, für n=1n = 1 linear und homogen; interessant sind alle anderen nn. Man teilt durch yny^n und substituiert z=y1−nz = y^{1-n}. Wegen z′=(1−n) y−n y′z' = (1-n)\,y^{-n}\,y' entsteht die lineare Gleichung

z′+(1−n) a(x) z=(1−n) b(x)z' + (1-n)\,a(x)\,z = (1-n)\,b(x)

Diese löst man mit der Formel für lineare DGL oder mit Variation der Konstanten und setzt am Ende y=z11−ny = z^{\frac{1}{1-n}} zurück.

Die Bedingung y>0y > 0 sorgt dafür, dass yny^n und y1−ny^{1-n} für jedes reelle nn definiert sind. Bei ganzzahligem nn darf yy auch negativ sein, solange y≠0y \neq 0 ist, wo durch yy geteilt wird. Für n>0n > 0 ist außerdem y≡0y \equiv 0 eine Lösung, die bei der Substitution verloren geht.

Beispiele

Die logistische Gleichung y′=y−y2y' = y - y^2 ist eine Bernoulli-DGL mit a=−1a = -1, b=−1b = -1, n=2n = 2. Mit z=y−1z = y^{-1} wird daraus

z′+(1−2)(−1) z=(1−2)(−1)⇒z′+z=1z' + (1-2)(-1)\,z = (1-2)(-1) \quad\Rightarrow\quad z' + z = 1

mit der Lösung z=1+C e−xz = 1 + C\,e^{-x}. Also ist y=11+C e−xy = \frac{1}{1 + C\,e^{-x}}. Für y(0)=0,5y(0) = 0{,}5 ist C=1C = 1, und es gilt y(2)=11+e−2≈0,8808y(2) = \frac{1}{1 + e^{-2}} \approx 0{,}8808.

y′+y=y3y' + y = y^3 mit n=3n = 3. Die Substitution z=y−2z = y^{-2} ergibt

z′−2z=−2⇒z=1+C e2x⇒y=11+C e2xz' - 2z = -2 \quad\Rightarrow\quad z = 1 + C\,e^{2x} \quad\Rightarrow\quad y = \frac{1}{\sqrt{1 + C\,e^{2x}}}

für positive Lösungen. Probe: y′=−C e2x (1+C e2x)−3/2y' = -C\,e^{2x}\,(1 + C\,e^{2x})^{-3/2}, und y′+yy' + y ergibt genau (1+C e2x)−3/2=y3(1 + C\,e^{2x})^{-3/2} = y^3.

Typische Fehler

  • Den Faktor (1−n)(1-n) nur auf einer Seite anbringen. Er steht vor a(x)a(x) und vor b(x)b(x).
  • z=yn−1z = y^{n-1} statt z=y1−nz = y^{1-n} substituieren. Dann entsteht keine lineare Gleichung.
  • Am Ende vergessen, von zz auf yy zurückzurechnen.
  • Die Lösung y≡0y \equiv 0 übersehen, wenn n>0n > 0 ist.

Verwandte Formeln