DGL mit Maple < Maple < Mathe-Software < Mathe < Vorhilfe
|
| Aufgabe | | Meine Aufgabe besteht darin einen Raketenstart zu simulieren. Die Simulation mit Vensim steht schon. Jetzt möchte ich mein Ergebniss mit Maple überprüfen und hoffentlich auf die selben Ergebnisse kommen! |
Hallo,
ich bin gerade dabei eine DGL zu lösen welche die Höhe (s(t)) einer Rakete beschreibt. Dabei habe ich Luftreibungskraft und Gravitationskraft berücksichtigt. Beide Kräfte ändern sich mit der Höhe.
Die Differentailgleichung sieht folgendermaßen aus:
dgl := diff(s(t), t, t) = [mm] (85*10^4+8.8*10^6)*(1/795600)-(6.673*10^{-11}*5.677)*10^24/(636800+s(t))^2-(.1*(19.96+9.58)*(1.2*exp(-(1/7990)*s(t))))*(diff(s(t), t))^2/(2*795600)
[/mm]
Mein Problem liegt jetzt speziell in der Maplegerechten Formulierung der Anfangswerte und eine spezielle Lösung zu erhalten. Anschließend möchte ich einen Plot von s(t) erhalten.
Die allgemeine Lösung mit dem Befehl : "dsolve(dgl, s(t))" berechnet mir Maple ohne Probleme. Nun möchte ich die Konstanten C1 und C2 belegen, sodas gilt s(0)=0 sowie v(0)=0 (diff(s(t);t)=0)
Anfangshöhe sowie Anfangsgeschwindigkeit sollen Null sein.
Vielen Danke schonmal für eure Hilfe!
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
|
|
| |
|
Hallo,
es ist eigentlich eben so einfach, wie es in der Hilfe zu Maple beschrieben wird:
soln:=dsolve({dgl,s(0)=0,D(s)(0)=0},numeric,s(t));
Wahlweise kannst Du auch noch range=0..sonstwohin oder eine der anderen ausführlich beschriebenen Optionen angeben.
Zum Plotten das Paket "plots" per with(plots) hinzuladen und z.B. mit odeplot(soln,0..60) die erste Minute des Fluges darstellen.
Ansonsten empfehle ich ![[]](/images/popup.gif) RTFM ![[buchlesen]](/images/smileys/buchlesen.gif "[buchlesen]")
Gruß,
Peter
edit: ich wollte ja noch ein ![[a]](/images/paperclip.gif "Dateianhänge") Mini-Worksheet anhängen...
Dateianhänge:
Anhang Nr. 1 (Typ: mw) [nicht öffentlich]
|
|
|
|
