| |
| Autor |
  DGL mittels Substitution lösen |
|
maxmustermann9991
Aktiv 
Dabei seit: 28.02.2016
Mitteilungen: 271
 |
Ich habe eine DGL, die ich mittels Substitution lösen soll.
  
y'(x)+4*y^2/x^2-y/x=0 mit Anfangsbedingung: y(1)=2 gesucht: Werte für x=1,1 und x=1,4
Mein Vorgehen: erst einmal y' auf eine Seite bringen:
y'(x)=y/x-4*y^2/x^2 Subst. u=y/x y'(x)=u-4u^2 du/dx=u-4u^2 1/(u-4u^2) du=dx ln(u)-ln(4u-1)=x+C e-Funktion auf beiden Seiten anwenden: u-(4u-1)=e^(x+C) u=(e^(x+C)-1)/3 Rücksubstitution: u=y/x y/x=(e^(x+C)-1)/3 y=(e^(x+C)-1)/3*x Anfangsbedingung: y(1)=2 (e^(1+C)-1)/3*1=2 C=ln(7)-1 y(1,1)=(e^(1,1+ln(7)-1)-1)/3*1,1=2,46994 y(1,4)=4,40663
Wo mache ich hierbei meinen Fehler, da die Lösungen nicht korrekt sind?
|
 Notiz  Profil
 Quote
 Link |
MontyPythagoras
Senior 
Dabei seit: 13.05.2014
Mitteilungen: 2290
Aus: Hattingen
 |      Beitrag No.1, eingetragen 2019-11-14
|
Hallo mm9991,
den Fehler machst Du schon ganz am Anfang. Du kannst nicht $u=\frac yx$ substituieren, und dann ganz selbstverständlich von einem Schritt zum nächsten $y'$ durch $u'$ ersetzen. Nach der Substitution ist doch wohl $y=u\cdot x$. Was folgt daraus für $y'$?
Ciao,
Thomas
|
Notiz Profil
Quote
Link |
Diophant
Senior
Dabei seit: 18.01.2019
Mitteilungen: 4238
Aus: Rosenfeld, BW
 | Beitrag No.2, eingetragen 2019-11-14
|
\(\begingroup\)\(\newcommand{\ba}{\begin{aligned}}
\newcommand{\ea}{\end{aligned}}
\newcommand{\bpm}{\begin{pmatrix}}
\newcommand{\epm}{\end{pmatrix}}
\newcommand{\bc}{\begin{cases}}
\newcommand{\ec}{\end{cases}}\)
Hallo,
2019-11-14 11:02 - maxmustermann9991 im Themenstart schreibt:
Ich habe eine DGL, die ich mittels Substitution lösen soll.
y'(x)+4*y^2/x^2-y/x=0 mit Anfangsbedingung: y(1)=2 gesucht: Werte für x=1,1 und x=1,4
Mein Vorgehen: erst einmal y' auf eine Seite bringen:
y'(x)=y/x-4*y^2/x^2 Subst. u=y/x y'(x)=u-4u^2 du/dx=u-4u^2 1/(u-4u^2) du=dx ln(u)-ln(4u-1)=x+C e-Funktion auf beiden Seiten anwenden: u-(4u-1)=e^(x+C) u=(e^(x+C)-1)/3 Rücksubstitution: u=y/x y/x=(e^(x+C)-1)/3 y=(e^(x+C)-1)/3*x Anfangsbedingung: y(1)=2 (e^(1+C)-1)/3*1=2 C=ln(7)-1 y(1,1)=(e^(1,1+ln(7)-1)-1)/3*1,1=2,46994 y(1,4)=4,40663
Wo mache ich hierbei meinen Fehler, da die Lösungen nicht korrekt sind?
Gleich zu Beginn, wo du \(y'=\frac{du}{dx}\) setzt. Das geht so natürlich nicht.
Du musst zunächst die Substitutionsgleichung \(u=\frac{y}{x}\) geeignet* ableiten, um das Differential \(dy\) irgendwie durch \(du\) ersetzt zu bekommen.
Gruß, Diophant
* Siehe hierzu den vorherigen Beitrag von MontyPythagoras.
[Die Antwort wurde vor Beitrag No.1 begonnen.]
[Verschoben aus Forum 'Differentialgleichungen' in Forum 'DGLen 1. Ordnung' von Diophant]\(\endgroup\)
|
Notiz Profil
Quote
Link |
maxmustermann9991
Aktiv
Dabei seit: 28.02.2016
Mitteilungen: 271
| Beitrag No.3, vom Themenstarter, eingetragen 2019-11-14
|
Notiz Profil
Quote
Link |
Diophant
Senior
Dabei seit: 18.01.2019
Mitteilungen: 4238
Aus: Rosenfeld, BW
| Beitrag No.4, eingetragen 2019-11-14
|
\(\begingroup\)\(\newcommand{\ba}{\begin{aligned}}
\newcommand{\ea}{\end{aligned}}
\newcommand{\bpm}{\begin{pmatrix}}
\newcommand{\epm}{\end{pmatrix}}
\newcommand{\bc}{\begin{cases}}
\newcommand{\ec}{\end{cases}}\)
Hallo,
die erste Version ist im Prinzip richtig. Nur macht das alles keinen Sinn, so lange da \(u'\) steht. Da nach x abgeleitet wird, bekommst du zunächst
\[\frac{dy}{dx}=\frac{du}{dx}x+u\]
Damit kannst du nun den Differentialquotenten \(\frac{dy}{dx}\) substituieren.
Gruß, Diophant
\(\endgroup\)
|
Notiz Profil
Quote
Link |
maxmustermann9991
Aktiv
Dabei seit: 28.02.2016
Mitteilungen: 271
| Beitrag No.5, vom Themenstarter, eingetragen 2019-11-14
|
Sieht das Folgende soweit richtig aus? Bevor ich hier falsch weiterrechne.
dy/dx=du/dx*x+u=u-4u^2 du/dx*x=-4u^2 du*x=-4u^2 dx 1/(-4u^2) du=1/x dx Integrieren: 1/4u=ln(x)+C Rücksubstitution: u=y/x 1/(4*y/x)=ln(x)+C x/4y=ln(x)+C 4y/x=1/(ln(x)+C) y=x/(4*(ln(x)+C))
|
Notiz Profil
Quote
Link |
Diophant
Senior
Dabei seit: 18.01.2019
Mitteilungen: 4238
Aus: Rosenfeld, BW
| Beitrag No.6, eingetragen 2019-11-14
|
\(\begingroup\)\(\newcommand{\ba}{\begin{aligned}}
\newcommand{\ea}{\end{aligned}}
\newcommand{\bpm}{\begin{pmatrix}}
\newcommand{\epm}{\end{pmatrix}}
\newcommand{\bc}{\begin{cases}}
\newcommand{\ec}{\end{cases}}\)
Hallo,
es ist noch nicht ganz richtig. Das Integral der Hyperbelfunktion ist ja
\[\int{\frac{dx}{x}}=\ln\left|x\right|+C\]
Außerdem würde ich zu deiner Rücksubstitution anmerken wollen: quer durch die Brust ins Auge... 😉
Warum löst du nicht zunächst nach u auf, substituierst zurück und löst dann vollends nach y?
Außerdem könnte man - so man möchte - die Integrationskonstante noch in den Logarithmus hereinziehen. Das ist aber Geschmacksache.
Bis auf die fehlenden Betragsklammern ist es jedoch richtig gerechnet.
Gruß, Diophant\(\endgroup\)
|
Notiz Profil
Quote
Link |
maxmustermann9991
Aktiv
Dabei seit: 28.02.2016
Mitteilungen: 271
| Beitrag No.7, vom Themenstarter, eingetragen 2019-11-14
|
2019-11-14 11:58 - Diophant in Beitrag No. 6 schreibt:
Außerdem würde ich zu deiner Rücksubtsitution anmerken wollen: quer durch die Brust ins Auge... 😉
Warumk löst du nicht zunächst nach u auf, substitierst zurück und löst dann vollends nach y? Normalerweise, mit Stift und Papier, würde ich es genauso tun wie du es hier sagst. Weshalb ich diesmal den umständlicheren Weg bei der Eingabe verwendet habe, das weiß nur der liebe Gott.
Außerdem könnte man so man möchte die Integrationskonstante noch in den Logarithmus hereinziehen. Das ist aber Geschmacksache.
Stimmt.
Bis auf die fehlenden Betragsklammern ist es jedoch richtig gerechnet.
Die sind tatsächlich der Faulheit bei der Eingabe mit dem Editor geschuldet. 😁
Nun weiter mit der Aufgabe:
y=x/(4*(ln abs(x)+C)) Anfangsbedingung: y(1)=2 habe im Folgnden ln(1) einfach mitgeschleppt, mir ist bewusst, dass ln(1)=0 ist 1/(4*(ln(1)+C))=2 1/(ln(1)+C)=8 ln(1)+C=1/8 C=1/8-ln(1) C=1/8 Lösung mit Anfangswertbedingung: y=x/(4*(ln abs(x)+1/8)) y(1,1)=1,1/(4*(ln(1,1)+1/8))=1,24824 y(1,4)=1,4/(4*(ln(1,4)+1/8))=0,78544
Sind die Lösungen nun korrekt? 😄
|
Notiz Profil
Quote
Link |
Diophant
Senior
Dabei seit: 18.01.2019
Mitteilungen: 4238
Aus: Rosenfeld, BW
| Beitrag No.8, eingetragen 2019-11-14
|
\(\begingroup\)\(\newcommand{\ba}{\begin{aligned}}
\newcommand{\ea}{\end{aligned}}
\newcommand{\bpm}{\begin{pmatrix}}
\newcommand{\epm}{\end{pmatrix}}
\newcommand{\bc}{\begin{cases}}
\newcommand{\ec}{\end{cases}}\)
Hallo,
die Konstante ist richtig, ebenso der Funktionswert \(y(1.1)\). Beim zweiten Wert hast du im Ergebnis noch einen Zahlendreher drin, das sollte \(y(1.4)\approx 0.758\) heißen.
Gruß, Diophant\(\endgroup\)
|
Notiz Profil
Quote
Link |
maxmustermann9991
Aktiv
Dabei seit: 28.02.2016
Mitteilungen: 271
| Beitrag No.9, vom Themenstarter, eingetragen 2019-11-14
|
Habs gerade nochmal eingegeben, jawohl war ein Zahlendreher.
Vielen Dank nochmal für die Hilfen! :-)
|
Notiz Profil
Quote
Link |
maxmustermann9991 hat die Antworten auf ihre/seine Frage gesehen.
maxmustermann9991 hat selbst das Ok-Häkchen gesetzt. |  [Neues Thema]  [Druckversion] |
|
All logos and trademarks in this site are property of their respective owner. The comments are property of their posters, all the rest © 2001-2020 by Matroids Matheplanet
This web site was originally made with PHP-Nuke, a former web portal system written in PHP that seems no longer to be maintained nor supported. PHP-Nuke is Free Software released under the GNU/GPL license.
Ich distanziere mich von rechtswidrigen oder anstößigen Inhalten, die sich trotz aufmerksamer Prüfung hinter hier verwendeten Links verbergen mögen.
Lesen Sie die
Nutzungsbedingungen,
die Distanzierung,
die Datenschutzerklärung und das Impressum.
[Seitenanfang]
|
Home / Home ohne Frame
Wie unterstütze ich den Matheplaneten mit Geld?
