0 Daumen
910 Aufrufe

Aufgabe:

Lösungen von folgender Differentialgleichung:

y'(x) = y2xy2xy3+xy+x2 -\frac{y^{2}-xy}{2xy^{3}+xy+x^{2}}


Problem/Ansatz:

Ich hätte es mit Substitution probiert, die Gleichung versucht umzuformen, aber konnte keine Lösung finden.

Ich würde um Hilfe bitten wie ich diese Differentialgleichung lösen könnte.

Avatar von

2 Antworten

0 Daumen
 
Beste Antwort

Hallo,

um ehrlich zu sein, wußte ich nicht (s.u.) wie man so eine DGL löst. In Zeiten von Wolfram Alpha & Co kann man sich aber damit behelfen eine Herleitung zu 'konstruieren'. Die (oder eine) Lösung wärey2xy+ln(y)+ln(x)=Cy^{2}-\frac{x}{y}+\ln\left(y\right)+\ln\left(x\right)=Cleitet man diese Lösung nach xx ab und formt man sie etwas um, so kommt man am Ende auf obige DGL - probier das mal selber. Und wenn man diesen Weg rückwärts beschreitet, so kommt dabei folgendes heraus:y=y2xy2xy3+xy+x22xy3y+xyy+x2y=y2+xy÷xy22yy+yy+xyy2=1x+1y1y2yy+yy+xyy21y=1x2yy+yy+xyy2yy2=1x2yy+yyyxyy2=1xdxy2+ln(y)xy=ln(x)+C+ln(x)y2+ln(y)xy+ln(x)=C\begin{aligned}y'& =-\frac{y^{2}-xy}{2xy^{3}+xy+x^{2}} \\ 2xy^3y' +xyy' +x^2y' &= -y^2 + xy &&|\, \div xy^2 \\ 2yy' + \frac{y'}{y} + \frac{xy'}{y^2} &= - \frac{1}{x} +\frac{1}{y} &&|\,-\frac{1}{y}\\ 2yy' + \frac{y'}{y} + \frac{xy'}{y^2} -\frac{1}{y} &= - \frac{1}{x}\\ 2yy' + \frac{y'}{y} + \frac{xy'}{y^2} -\frac{y}{y^2} &= - \frac{1}{x}\\ 2yy' + \frac{y'}{y} - \frac{y-xy'}{y^2} &= - \frac{1}{x} &&|\,\int\text{d}x\\ y^2 + \ln(y) - \frac{x}{y} &= -\ln(x) + C &&|\, +\ln(x) \\ y^2 + \ln(y) - \frac{x}{y} +\ln(x)&= C \\\end{aligned}D.h. die Terme werden solange 'geknetet' bis sich nach xx integrierbare Ausdrücke ergeben.

Das ganze in graphischer Form:


Den Punkt kann man mit der Maus verschieben. Es wird immer die Richtung der Ableitung - also y=f(x,y)y'=f(x,y) - angezeigt.

Gruß Werner

Avatar von 49 k

Vielen Dank!

Ich habe mich dank des Tipps von Grosserloewe schlau gemacht.

Die Umwandlung in die im Wiki-Artikel beschrieben Form der exakten DGL sieht so aus:y=y2xy2xy3+xy+x2(y2xy)=p(x,y)+y(2xy3+xy+x2)=q(x,y)=0y'=-\frac{y^{2}-xy}{2xy^{3}+xy+x^{2}} \\ \underbrace{(y^{2}-xy)}_{=p(x,y)} + y'\underbrace{(2xy^{3}+xy+x^{2})}_{=q(x,y)} = 0

Da die Integrabilitätsbedingung p(x,y)y=?q(x,y)x\frac{\partial p(x,y)}{\partial y} \stackrel{?}{=} \frac{\partial q(x,y)}{\partial x} hier nicht erfüllt ist, muss man zunächst nach einer Funktion μ(x,y)\mu(x,y) suchen, mit der man die DGL multipliziert. Dazu multipliziere ich pp und qq mit μ(x,y)\mu(x,y) und setze es in die Integrabilitätsbedingung ein:μ(x,y)p(x,y)y=μ(x,y)q(x,y)xμ(x,y)y(y2xy)+μ(x,y)(2yx)=μ(x,y)x(2xy3+xy+x2)+μ(x,y)(2y3+y+2x)μ(x,y)(2yx)μ(x,y)(2y3+y+2x)=μ(x,y)x(2xy3+xy+x2)+μ(x,y)y(xyy2)μ(x,y)(y3x2y3)=μ(x,y)x(2xy3+xy+x2)+μ(x,y)y(xyy2)(xy23x2y2xy4)(2xy4+xy2+x2y)+(x2yxy2)\begin{aligned} \frac{\partial \mu(x,y)p(x,y)}{\partial y} &= \frac{\partial \mu(x,y)q(x,y)}{\partial x} \\ \frac{\partial \mu(x,y)}{\partial y}(y^2 - xy) + \mu(x,y)(2y - x) &= \frac{\partial \mu(x,y)}{\partial x} \left(2xy^3 +xy+x^2\right)+\mu(x,y)\left( 2y^3+y+2x \right)\\ \mu(x,y)(2y - x) -\mu(x,y)\left( 2y^3+y+2x \right)&= \frac{\partial \mu(x,y)}{\partial x} \left(2xy^3 +xy+x^2\right)+\frac{\partial \mu(x,y)}{\partial y}(xy -y^2) \\ \mu(x,y)\left(y -3x- 2y^3 \right)&= \frac{\partial \mu(x,y)}{\partial x} \left(2xy^3 +xy+x^2\right)+\frac{\partial \mu(x,y)}{\partial y}(xy -y^2) \\ \left({\color{red}xy^2} {\color{blue}-3x^2y}- 2xy^4 \right)&\leftrightarrow \left(2xy^4 {\color{red}+xy^2}+{\color{blue}x^2y}\right)+({\color{blue}x^2y} {\color{red}-xy^2}) \\ \end{aligned}die letzte Zeile mit den Farben ist keine formale Gleichung, sondern nur eine 'Denkhilfe'! Da alle Terme Potenzen von xx und yy sind, kann man wohl davon ausgehen, dass es bei μ(x,y)\mu(x,y) auch so ist. Und wenn man nach xx bzw. yy ableitet, dann vermindert sich der Exponent jeweils um 1.

Also füge ich beim nicht abgeleiteten Term ein xyxy hinzu und bei der Ableitung nach xx nehme ich ein xx weg und füge nur ein yy hinzu und beim nach yy abgeleiteten Term bleibt dann ein xx. Und dann vergleiche ich die Terme mit identischen Exponenten ... die habe ich farblich markiert. Daraus folgt dann z.B.:μ(x,y)(2xy4)=μ(x,y)x(2xy4)\mu(x,y) \cdot (-2xy^4) = \frac{\partial \mu(x,y)}{\partial x} \cdot (2xy^4)und dies geht wegen des Vorzeichenwechsels nur, wenn u(x,y)u(x,y) den Faktor x1=1/xx^{-1}=1/x enthält.

Daselbe mache ich nun für die rot markierten Terme. Natürlich ohe die hinzugefügten xx bzw. yy:μ(x,y)=ykx    μ(x,y)x=ykx2,μ(x,y)y=kyk1xykxy=ykx2xy+kyk1x(y2)yk=ykkyk    k=2    μ(x,y)=1xy2\mu(x,y)=\frac{y^k}{x} \implies \frac{\partial \mu(x,y)}{\partial x} = -\frac{y^k}{x^2}, \quad \frac{\partial \mu(x,y)}{\partial y}=k\frac{y^{k-1}}{x} \\ \begin{aligned} \frac{y^k}{x} {\color{red}y} &= -\frac{y^k}{x^2} {\color{red}xy} +k\frac{y^{k-1}}{x} ({\color{red}-y^2})\\ y^k &= -y^k - ky^k \\ \implies k &= -2\\ \implies \mu(x,y)&= \frac{1}{xy^2} \end{aligned}Für die blauen kannst Du es ja noch mal kontrollieren.

Jetzt kann man dann schreiben:μ(x,y)=1xy2(y2xy)+y(2xy3+xy+x2)=01xy2(1x1y)=p(x,y)+y(2y+1y+xy2)=q(x,y)=0\mu(x,y)= \frac{1}{xy^2} \\ \begin{aligned} \left(y^2-xy\right) + y'(2xy^3+xy+x^2) &= 0 &&|\,\cdot \frac{1}{xy^2} \\ \underbrace{\left(\frac{1}{x}-\frac{1}{y}\right)}_{=p(x,y)} + y'\underbrace{\left(2y+ \frac{1}{y} + \frac{x}{y^2}\right)}_{=q(x,y)} &= 0 \end{aligned}Und das weiter Vorgehen ist im Wiki-Artikel beschrieben. Hier noch meine Rechnung:Φ(x,y)=p(x,y)dx+φ(x)=ln(x)xy+φ(y)dφ(y)dy=q(x,y)yp(x,y)dx=2y+1y+xy2y(ln(x)xy)=2y+1y+xy2xy2=2y+1ydyφ(y)=y2+ln(x)+C    Φ(x,y)=ln(x)xy+y2+ln(x)+C\begin{aligned} \varPhi(x,y) &= \int p(x,y)\,\text{d}x + \varphi(x)\\ &= \ln(x) - \frac{x}{y} + \varphi(y) \\\end{aligned}\\\begin{aligned} \frac{\text{d} \varphi(y)}{\text{d} y} &= q(x,y) - \frac{\partial}{\partial y}\int p(x,y)\,\text{d}x\\ &= 2y+ \frac{1}{y} + \frac{x}{y^2} - \frac{\partial}{\partial y}\left(\ln(x) - \frac{x}{y}\right)\\ &= 2y + \frac{1}{y} + \frac{x}{y^2} - \frac{x}{y^2}\\ &= 2y + \frac{1}{y}&&|\,\int \text{d}y\\ \varphi(y) &= y^2 + \ln(x) + C \\ \end{aligned}\\ \implies \varPhi(x,y)=\ln(x) - \frac{x}{y}+y^2 + \ln(x) + C kommt tatsächlich das gleiche heraus. Wieder was dazu gelernt!

+1 Daumen

Hallo,

diese DGL kannst Du nach dem Verfahren "Exakte DGL" lösen.

Prüfe zuerst , ob Py=Qx ist .(Integrabilitätsbedingung)

Wenn nein, mußt Du einen "Integrierenden Faktor" berechnen und die gesamte DGL damit multiplizieren.

Dann kannst Du die DGL berechnen.


y=y2xy2xy3+xy+x2dydx=y2xy2xy3+xy+x2(2xy3+xy+x2Q)dy+(y2xyP)dx=0 \begin{array}{l}y^{\prime}=-\frac{y^{2}-x y}{2 x y^{3}+x y+x^{2}} \\ \frac{d y}{d x}=-\frac{y^{2}-x y}{2 x y^{3}+x y+x^{2}} \\ (\underbrace{2 x y^{3}+x y+x^{2}}_{Q}) d y+(\underbrace{y^{2}-x y}_{P}) d x=0\end{array}

Py=2yxQx=2y3+y+2xPyQx nicht exakt  integrieˋrender  Faktor  Ansatz :  m(x,y)=xaybPdm(x,y)dyQdm(x,y)dx+m(x,y)(dPdydQdx)=0 \begin{array}{l}\begin{array}{l}P_{y}=2 y-x \\ Q_{x}=2 y^{3}+y+2 x \Rightarrow P_{y} \neq Q_{x}\end{array} \\ \text { nicht exakt } \\ \Rightarrow \text { integrièrender } \\ \text { Faktor } \\ \text { Ansatz: } m(x, y)=x^{a} y^{b} \\ P \frac{d m(x, y)}{d y}-Q \frac{d m(x, y)}{d x}+m(x, y)\left(\frac{d P}{d y}-\frac{d Q}{d x}\right)=0 \\\end{array}


blob.png

blob.png

Avatar von 121 k 🚀

Dankeschön, werde es später probieren.

diese DGL kannst Du nach dem Verfahren "Exakte DGL" lösen

Der Begriff "Exakte DGL" war mir gänzlich fremd. Danke für den Tipp ... es funktioniert (s Kommentar unter meiner Antwort)

@ Danie

Viel Spaß mit dieser  Sonntagslektüre :-) , siehe oben, habe meinen Beitrag ergänzt

Ein anderes Problem?

Stell deine Frage