Lösen eines DGLsystems < Mathematica < Mathe-Software < Mathe < Vorhilfe
|
Hi @ all,
ich habe folgendes Problem:
Ich möchte gerne ein elektro-mechanisches System, dass aus 8 Differentialgleichungen
erster Ordnung besteht, in Mathematica simulieren bzw. das DGLsystem lösen.
Dazu bin ich hingegangen und habe das System mit den 8 teilweise nichtlinearen DGLs
und den Anfangswerte wie folgt definiert:
System=
[mm] \qquad[/mm] [mm] \begin{array}{ll}
\{ x^{'}_{1}(t)= x_{2}(t), & x_{1}(0)=-\bruch{\pi}{3},
\\ x^{'}_{2}(t)= f ( x_{1}(t),x_{2}(t),x_{3}(t),x_{6}(t)), & x_{2}(0)=0,
\\x^{'}_{3}(t)= x_{4}(t),& x_{3}(0)=0,
\\ x^{'}_{4}(t)=f(x_{3}(t),x_{4}(t)),& x_{4}(0)=0,
\\ x^{'}_{5}(t)= x_{6}(t),\,& x_{5}(0)=-\bruch{\pi}{2},
\\ x^{'}_{6}(t)=f( x_{3}(t),x_{4}(t),x_{6}(t),x_{7}(t)),& x_{6}(0)=0,
\\x^{'}_{7}(t)= x_{8}(t),\, & x_{7}(0)=0,
\\ x^{'}_{8}(t)=f(x_{7}(t),x_{8}(t)),& x_{8}(0)=0 \}
\end{array} [/mm]
Dann habe ich versucht das System mit
DSolve= [mm] \left[System, \left\{ x_{1}(t), x_{2}(t), x_{3}(t), x_{4}(t), x_{5}(t), x_{6}(t), x_{7}(t), x_{8}(t) \right\},t \right]
[/mm]
bzw.
NDSolve= [mm] \left[System, \left\{ x_{1}(t), x_{2}(t), x_{3}(t), x_{4}(t), x_{5}(t), x_{6}(t), x_{7}(t), x_{8}(t) \right\},\{t,0,10\} \right]
[/mm]
zulösen aber ohne Erfolg! Mathematica löst das DGLsystem einfach nicht!?!
Was habe ich falsch gemacht?
Grüße und Danke im Voraus
BlackComet
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
Dateianhänge:
Anhang Nr. 1 (Typ: nb) [nicht öffentlich]
|
|
| |
|
Kannst du vielleicht mal das $f(x)$ angeben... Denn ansonsten ist das ein Problem, das zu allgemein ist, um es zu lösen.
MfG Sunny
|
|
|
| |
|
Hi Leute,
ich habe zu Vollständigkeitszwecken den Mathematica-File als Anhang dem Original-Post beigefügt.
Vielleicht kann man hier schneller den Fehler entdecken.
Vielen Dank
BlackComet
|
|
|
| |
|
Hallo black-comet,
Mit DSolve durchzukommen ist natürlich klar hoff-
nungslos ...
NDSolve scheitert offenbar, weil eine der Ableitungen
schon zu Beginn (t=0) keinen Zahlenwert liefert.
Mir ist noch aufgefallen, dass du zwar [mm] u_1=1 [/mm] als Konstante
definiert hast, dann aber in einer (der längsten) Formel
ein [mm] u_1[/mm] [t] stehen hast. Dies kann wohl nicht gut gehen !
LG Al-Chwarizmi
|
|
|
|
