Trennung der Veränderlichen

Proportionale Differentialgleichung Erster Ordnung lösen[1] durch Trennung der Veränderlichen.[2]
Lineare Differentialgleichung lösen[3] durch Trennung der Veränderlichen.[2]

Die Methode der Trennung der Veränderlichen, Trennung der Variablen, Separationsmethode oder Separation der Variablen ist ein Verfahren aus der Theorie der gewöhnlichen Differentialgleichungen. Mit ihr lassen sich separierbare Differentialgleichungen erster Ordnung lösen. Das sind Differentialgleichungen, bei denen die erste Ableitung ein Produkt aus einer nur von x {\displaystyle x} und einer nur von y {\displaystyle y} abhängigen Funktion ist: y = f ( y ) g ( x ) . {\displaystyle y'=f(y)g(x).} Der Begriff „Trennung der Veränderlichen“ geht auf Johann I Bernoulli zurück, der ihn 1694 in einem Brief an Gottfried Wilhelm Leibniz verwendete.[4]

Ein ähnliches Verfahren für bestimmte partielle Differentialgleichungen ist der Separationsansatz. Ein wichtiger Anwendungsfall in Koordinatensystemen, deren Koordinatenflächen konfokale Quadriken sind, ist die Laplace- und die Helmholtz-Gleichung, siehe dort.

Lösung des Anfangswertproblems

Wir untersuchen das Anfangswertproblem

y ( x ) = f ( y ( x ) ) g ( x )   ,   y ( x 0 ) = y 0 {\displaystyle y'(x)=f(y(x))g(x)\ ,\ y(x_{0})=y_{0}}

für stetige (reelle) Funktionen f {\displaystyle f} und g {\displaystyle g} . Falls f ( y 0 ) = 0 {\displaystyle f(y_{0})=0} , so wird dieses Anfangswertproblem durch die konstante Funktion y ( x ) := y 0 {\displaystyle y(x):=y_{0}} gelöst. Diese Lösung muss unter den angegebenen Bedingungen nicht eindeutig sein.

Formulierung des Satzes

Voraussetzungen

U {\displaystyle U} sei ein offenes Intervall, y 0 U {\displaystyle y_{0}\in U} und f : U R {\displaystyle f\colon U\to \mathbb {R} } eine stetige Funktion mit f ( x ) 0 {\displaystyle f(x)\neq 0} für alle x U {\displaystyle x\in U} . Dann gilt nach dem Zwischenwertsatz entweder f ( x ) > 0 {\displaystyle f(x)>0} für alle x U {\displaystyle x\in U} , oder f ( x ) < 0 {\displaystyle f(x)<0} für alle x U {\displaystyle x\in U} . Also ist die Funktion

Φ : U R y y 0 y 1 f ( s ) d s {\displaystyle {\begin{aligned}\Phi \colon U&\to \mathbb {R} \\y&\mapsto \int _{y_{0}}^{y}{\frac {1}{f(s)}}{\rm {d}}s\end{aligned}}}

streng monoton (das folgt aus dem Hauptsatz der Differential- und Integralrechnung und dem Mittelwertsatz). Das heißt, Φ {\displaystyle \Phi } ist injektiv und es gibt die Umkehrfunktion Φ 1 : Φ ( U ) U {\displaystyle \Phi ^{-1}\colon \Phi (U)\to U} .
Ferner sei V {\displaystyle V} ein offenes Intervall, x 0 V {\displaystyle x_{0}\in V} und g : V R {\displaystyle g\colon V\to \mathbb {R} } eine stetige Funktion. Dann ist die Funktion

h : V R x x 0 x g ( s ) d s {\displaystyle {\begin{aligned}h\colon V&\to \mathbb {R} \\x&\mapsto \int _{x_{0}}^{x}g(s){\rm {d}}s\end{aligned}}}

wohldefiniert und differenzierbar. Wir wollen die Lösungsmenge M {\displaystyle M} des Anfangswertproblems bestimmen:

M := { u C 1 ( V , U ) : u ( x 0 ) = y 0  und  u = ( f u ) g } {\displaystyle M:=\{u\in C^{1}(V,U):u(x_{0})=y_{0}{\text{ und }}u'=(f\circ u)\cdot g\}}

Der Satz

Unter den oben genannten Voraussetzungen gilt:

M = { u C 1 ( V , U ) : Φ u = h } {\displaystyle M=\{u\in C^{1}(V,U):\Phi \circ u=h\}}

Das heißt, im Fall h ( V ) Φ ( U ) {\displaystyle h(V)\subset \Phi (U)} hat das Anfangswertproblem genau eine Lösung – nämlich die Funktion Φ 1 h {\displaystyle \Phi ^{-1}\circ h} – und andernfalls ist M {\displaystyle M} leer.

Beweis

Sei N := { u C 1 ( V , U ) : Φ u = h } {\displaystyle N:=\{u\in C^{1}(V,U):\Phi \circ u=h\}} . Wir beweisen zuerst M N {\displaystyle M\subset N} und dann N M {\displaystyle N\subset M} :

1.

Sei u M {\displaystyle u\in M} , dann gilt nach der Substitutions-Regel

h ( x ) = x 0 x g ( s ) d s = x 0 x u ( s ) f ( u ( s ) ) d s = u ( x 0 ) u ( x ) 1 f ( y ) d y = y 0 u ( x ) 1 f ( y ) d y = Φ ( u ( x ) ) Φ ( U ) {\displaystyle h(x)=\int _{x_{0}}^{x}g(s){\rm {d}}s=\int _{x_{0}}^{x}{\frac {u'(s)}{f(u(s))}}{\rm {d}}s=\int _{u(x_{0})}^{u(x)}{\frac {1}{f(y)}}{\rm {d}}y=\int _{y_{0}}^{u(x)}{\frac {1}{f(y)}}{\rm {d}}y=\Phi (u(x))\in \Phi (U)}

für alle x V {\displaystyle x\in V} , also Φ u = h {\displaystyle \Phi \circ u=h} .

2.

Nun bleibt zu zeigen, dass für den Fall h ( V ) Φ ( U ) {\displaystyle h(V)\subset \Phi (U)} das einzige Element von N {\displaystyle N} – die Funktion Φ 1 h {\displaystyle \Phi ^{-1}\circ h} – eine Lösung des Anfangswertproblems ist, also Φ 1 h M {\displaystyle \Phi ^{-1}\circ h\in M} gilt: Nach der Kettenregel, der Umkehrregel und dem Hauptsatz der Differential- und Integralrechnung gilt

( Φ 1 h ) ( x ) = ( Φ 1 ) ( h ( x ) ) h ( x ) = 1 Φ ( Φ 1 ( h ( x ) ) ) h ( x ) = 1 Φ ( u ( x ) ) h ( x ) = f ( u ( x ) ) g ( x ) {\displaystyle (\Phi ^{-1}\circ h)'(x)=(\Phi ^{-1})'(h(x))h'(x)={\frac {1}{\Phi '(\Phi ^{-1}(h(x)))}}h'(x)={\frac {1}{\Phi '(u(x))}}h'(x)=f(u(x))g(x)}

für alle x V {\displaystyle x\in V} . Natürlich ist u ( x 0 ) = y 0 {\displaystyle u(x_{0})=y_{0}} .

Bemerkung

D {\displaystyle D} und D {\displaystyle D'} seien Teilmengen der reellen Zahlen, f : D R {\displaystyle f\colon D\to \mathbb {R} } und g : D R {\displaystyle g\colon D'\to \mathbb {R} } stetige Funktionen, x 0 {\displaystyle x_{0}} sei ein innerer Punkt von D {\displaystyle D} , y 0 {\displaystyle y_{0}} ein innerer Punkt von D {\displaystyle D'} und f ( y 0 ) 0 {\displaystyle f(y_{0})\neq 0} . Dann gilt:

Ist f ( y 0 ) > 0 {\displaystyle f(y_{0})>0} ( f ( y 0 ) < 0 ) {\displaystyle (f(y_{0})<0)} , dann gibt es wegen der Stetigkeit von f {\displaystyle f} ein y 0 {\displaystyle y_{0}} umfassendes offenes Intervall U R {\displaystyle U\subset \mathbb {R} } mit f ( y ) > 0 {\displaystyle f(y)>0} ( f ( y ) < 0 ) {\displaystyle (f(y)<0)} für alle y U {\displaystyle y\in U} . Weil Φ {\displaystyle \Phi } auf U {\displaystyle U} stetig ist, ist Φ ( U ) {\displaystyle \Phi (U)} nach dem Zwischenwertsatz ein Intervall und es gilt h ( x 0 ) = 0 = Φ ( y 0 ) {\displaystyle h(x_{0})=0=\Phi (y_{0})} . Deswegen gibt es ein x 0 {\displaystyle x_{0}} umfassendes offenes Intervall V R {\displaystyle V\subset \mathbb {R} } , sodass die Abbildung

x x 0 x g ( s ) d s {\displaystyle x\mapsto \int _{x_{0}}^{x}g(s){\rm {d}}s}

für alle x V {\displaystyle x\in V} Werte in Φ ( U ) {\displaystyle \Phi (U)} hat. Das heißt, die Restriktionen f | U {\displaystyle f|_{U}} und g | V {\displaystyle g|_{V}} erfüllen die Bedingungen des oben formulierten Satzes.

Beispiel

Gesucht sei die Lösung y {\displaystyle y} des Anfangswertproblems

y = x y 2 + x   ,   y ( 0 ) = 1 {\displaystyle y'=xy^{2}+x\ ,\ y(0)=1} .

Hierbei handelt es sich um eine Differentialgleichung mit getrennten Variablen:

y = x ( y 2 + 1 ) {\displaystyle y'=x(y^{2}+1)} .

Setze also

Φ ( y ) := 1 y 1 1 + s 2 d s = arctan y arctan 1 = arctan y π 4 {\displaystyle \Phi (y):=\int _{1}^{y}{\frac {1}{1+s^{2}}}{\rm {d}}s=\arctan y-\arctan 1=\arctan y-{\frac {\pi }{4}}} .

Die Umkehrfunktion lautet

Φ 1 ( y ) = tan ( y + π 4 ) {\displaystyle \Phi ^{-1}(y)=\tan \left(y+{\frac {\pi }{4}}\right)} .

Also ist die Lösung des Anfangswertproblems gegeben durch

y ( x ) = tan ( 0 x s d s + π 4 ) = tan ( x 2 2 + π 4 ) {\displaystyle y(x)=\tan \left(\int _{0}^{x}s{\rm {d}}s+{\frac {\pi }{4}}\right)=\tan \left({\frac {x^{2}}{2}}+{\frac {\pi }{4}}\right)} .

Differentiale als anschauliche Rechenhilfe

Anschaulich besagt der Satz von der Trennung der Veränderlichen, dass das folgende Vorgehen erlaubt ist, d. h. zu richtigen Ergebnissen führt (obwohl die Differentiale d x {\displaystyle \mathrm {d} x} und d y {\displaystyle \mathrm {d} y} eigentlich nur Symbole sind, mit denen man streng genommen nicht rechnen kann):

  • Schreibe die Ableitung konsequent als d y d x {\displaystyle {\tfrac {\mathrm {d} y}{\mathrm {d} x}}} .
  • Bringe alle Terme, in denen ein x {\displaystyle x} vorkommt – einschließlich des d x {\displaystyle \mathrm {d} x} – auf die rechte, und alle anderen – einschließlich des d y {\displaystyle \mathrm {d} y} – auf die linke Seite, unter Anwendung gewöhnlicher Bruchrechnung.
  • Es sollte dann links im Zähler ein d y {\displaystyle \mathrm {d} y} und rechts im Zähler ein d x {\displaystyle \mathrm {d} x} stehen.
  • Setze einfach vor beide Seiten ein Integralsymbol und integriere.
  • Löse die Gleichung gegebenenfalls nach y {\displaystyle y} auf.
  • Ermittle die Integrationskonstante C {\displaystyle C} mithilfe der Anfangsbedingung.

Die Rechnung für das obige Beispiel würde dann auf folgende Weise ablaufen:

d y d x = x y 2 + x d y 1 + y 2 = x d x arctan ( y ) = x 2 2 + C y = tan ( x 2 2 + C ) {\displaystyle {\frac {\mathrm {d} y}{\mathrm {d} x}}=xy^{2}+x\;\Longrightarrow \,\int {\frac {\mathrm {d} y}{1+y^{2}}}=\int x\,\mathrm {d} x\;\Longrightarrow \,\arctan(y)={\frac {x^{2}}{2}}+C\;\Longrightarrow \,y=\tan \left({\frac {x^{2}}{2}}+C\right)}

mit 1 = y ( 0 ) = tan ( C ) {\displaystyle 1=y(0)=\tan(C)} , also C = arctan ( 1 ) = π 4 {\displaystyle C=\arctan(1)={\frac {\pi }{4}}} .

Computerprogramm

Die CAS-Software Xcas kann Trennung der Veränderlichen mit folgendem Befehl[5] durchführen: split((x+1)*(y-2),[x,y]) = [x+1,y-2]

Literatur

  • Wolfgang Walter: Gewöhnliche Differentialgleichungen. 4. überarbeitete Auflage. Springer, 1990, ISBN 3-540-52017-1, S. 13–20. 
  • Kurt Endl, Wolfgang Luh: Analysis I. 9. Auflage. Aula-Verlag, Wiesbaden 1989, ISBN 3-89104-498-4, S. 316–333. 
  • Harro Heuser: Gewöhnliche Differentialgleichungen. Einführung in Lehre und Gebrauch. 6. aktualisierte Auflage. Vieweg+Teubner, 2009, ISBN 978-3-8348-0705-2, S. 102–122 (eingeschränkte Vorschau in der Google-Buchsuche). 
  • Jochen Merker: Differentialgleichungen. (pdf) In: Skript, Sommersemester. Uni Rostock, 1. April 2011, S. 12–14, abgerufen am 14. März 2024 (englisch). 
  • Eric W. Weisstein: Separation of Variables. In: MathWorld (englisch).
  • Paul Dawkins: Separation of Variables. In: Paul’s Online Math Notes. Lamar University, 16. November 2022, abgerufen am 1. Mai 2024 (englisch). 
  • Ron Larson: Separation of Variables. (pdf) In: Calculus: Applied approach. Abgerufen am 1. Mai 2024 (englisch). 

Einzelnachweise

  1. How do you solve this differential equation using the separation of variables dy/dx= (y-2)/x? Abgerufen am 27. Januar 2022 (englisch). 
  2. a b Trennung der Variablen: Erklärung und Beispiel. Abgerufen am 18. September 2021. 
  3. How do you solve this differential equation using the separation of variables dy/dx= (y-2)/x? Abgerufen am 27. Januar 2022 (englisch). 
  4. Harro Heuser: Gewöhnliche Differentialgleichungen. Einführung in Lehre und Gebrauch. 2. Auflage. Teubner, Stuttgart 1991, ISBN 3-519-12227-8, S. 128
  5. Bernard Parisse: Symbolic algebra and Mathematics with Xcas. Abgerufen am 23. August 2021.