Substitutionen zur Berechnung eines Integrals |
04.10.2015, 23:09 | DrummerS | Auf diesen Beitrag antworten » | ||||||
Substitutionen zur Berechnung eines Integrals Hallo, ausgehend von einer physikalischen Betrachtung habe ich folgende Gleichung mit den Konstanten sowie den Variablen hergeleitet: Ich suche eine Lösung der Form Die Größen bezeichnen Längen/Strecken, bezeichnet die Gravitationsbeschleunigung in Abhängigkeit des Abstandes zur Erde und Zeitpunkte. Meine Ideen: Durch Substitutionen mithilfe von , mit als Geschwindigkeit, käme ich zu: mit Damit wäre . Ich vermute die Ursache bei den gewählten Substitutionen. Könnte mir da jemand auf die Sprünge helfen ? Es fehlt natürlich noch der Schritt von |
||||||||
05.10.2015, 10:24 | klarsoweit | Auf diesen Beitrag antworten » | ||||||
RE: Substitutionen zur Berechnung eines Integrals Abgesehen von dieser Beziehung:
die mir irgendwie suspekt ist, ist mir die physikalische Betrachtung noch nicht so eingängig, vor allem da die Gravitationsbeschleunigung eine weg- und damit auch zeitabhängige Größe ist. Anfreunden könnte ich mich noch mit: oder eher noch mit dieser DGL: |
||||||||
05.10.2015, 10:59 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Ja, ich hab auch schon gerätselt, ob da nicht eigentlich das Integral gemeint ist. Ziemlich undeutlich dargestellt, da sollte der Threadersteller mal genau klären, ob das so oder anders gemeint ist: Bei allen anderen (von unabhängigen) Zeitpunkten kann man das ja schließlich als Konstante aus dem Integral herausziehen: P.S.: Zum physikalischen Verständnis: Laut Gravitationgesetz gilt . Wie kommst du von da zu deiner (mir suspekten) Integralgleichung , (wenn es denn so gemeint ist) ? |
||||||||
05.10.2015, 13:52 | DrummerS | Auf diesen Beitrag antworten » | ||||||
@klarsoweit und @Hal 9000 Hallo, erstmal vielen Dank für eure Kommentare! Leider neige ich manchmal dazu, Probleme zu verkomplizieren, indem ich versuche, Bilanzen aufzustellen, die Integrale enthalten. Die Darstellungen sind dann leider durch meine Gedankenfehler falsch...Ich sollte mir daher vielleicht mal antrainieren, mit einer DGL beginnen Problemstellung: Eine Kugel befindet sich in einem bestimmten Abstand zum Planeten Erde in Ruhe . Die Kugel beginnt dann aufgrund der Gravitation sich auf die Erde zuzubewegen. Am Zeitpunkt schlägt die Kugel auf der Erde auf. Wie lautet die Funktion ? Die Eigenrotation der Erde, die Erdatmosphäre und die Bewegung der Erde um die Sonne können für die Berechnung vernachlässigt werden. @klarsoweit Hatte mir zur Sinnhaftigkeit von auch Gedanken gemacht. Bin dann aber zu dem vermutlich naiven Schluss gekommen, dass bei die Ableitungen - mathematisch gesehen - durchrotiert werden können, auch wenn es physikalisch wenig Sinn macht. Habe mal die DGL h'(t) = -g(h(t))t berechnet... Die Diskrininate (Link) ergibt , womit genau eine reelle Lösung existiert. Mit p = 0 komme ich auf: Da , stimmt da noch etwas nicht. @Hal 9000 Sorry, hatte tatsächlich das im Integral vergessen, sodass . Hier nun meine Schritte: An der rechtsseitigen eckigen Klammer ist unten und oben hinzuzufügen. Jedenfalls sieht es so aus, dass die Darstellung leider komliziert aussehen wird. Bin ab jetzt erstmal "afk", aber nachher wieder online. VG, DrummerS |
||||||||
05.10.2015, 14:05 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Mit der Zeile gehe ich noch mit, mit und (also ) kann man damit auch berechnen. Aber wie kommst du von dieser DGL auf
das ist mir ein Rätsel. Durch Trennung der Variablen stände für diese DGL m.E. zunächst , da, was es nun zu integrieren gilt. Das "-" steht da übrigens, weil die Kugel ja fallen wird. Wenn ich dies richtig integriert habe, kommt da heraus. mit eingesetzt ergibt sich . Die Endzeit mit kann man damit dann gemäß berechnen. Eine explizite Umstellung scheint dagegen nur schwer möglich. |
||||||||
05.10.2015, 18:15 | DrummerS | Auf diesen Beitrag antworten » | ||||||
Ich hatte dieses Tool (Link) verwendet. Leider kann ich mein Ergebnis nicht mehr reproduzieren, was mir gerade wirklich peinlich ist. Ich verstehe, das so berechnet werden kann, wie du es beschrieben hast. Da hätte ich auch drauf kommen können....
Habe eben auch versucht, das Integral zu berechnen. Meine Vorgehensweise: Ab hier steige ich aus. Habe verschiedene Substitutions-Versuche durchgeführt, die zu keinem integrierbaren Ausdruck führten. Sogar die Konkurrenz-Internetseite (Link2) gibt eine Lösung mit Imaginärteil aus. Habe daher großen Respekt vor deinen mathematischen Fähigkeiten ! Da das Resultat so unschön ist, werde ich nun von der Aufgabe ablassen. Vielen Dank nochmals für eure Zeit! |
||||||||
Anzeige | ||||||||
|
||||||||
05.10.2015, 23:46 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Könnte man die Frage nicht mit einer einfacheren DGL behandeln? Die Kugel (als Massepunkt betrachtet) gehorcht im Schwerefeld der Erde doch der DGL Jetzt kann man den bekannten "Energiesatz-Trick" anwenden: Und beidseitig integrieren: Es folgt (negativ, weil die Kugel fallen soll). Mit der Methode der getrennen Variablen erhält man nun die Gleichung Hieraus sollte sich die Zeit bereits ermittleln lassen. Daraus folgen außerdem auch die Lösungen der DGL: Ob der andere Ansatz auch funktioniert habe ich ehrlich gesagt nicht verstanden . Allein dieses Integral ... |
||||||||
06.10.2015, 00:01 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Hmm beim Energiesatz habe ich wohl Integrationskonstanten ignoriert. Diese legen wahrscheinlich aber nur die Anfangsenergie fest, die man problemlos auf Null setzen kann. Oder? |
||||||||
06.10.2015, 07:09 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Hier fehlt jetzt eine Integrationskonstante, die sich aus der Anfangsbedingung ergibt. Wenn du die einfach 0 setzt, vereinfacht sich zwar der weitere Weg, aber du löst damit nicht das Problem. Ansonsten wird oben nichts anderes gemacht, nur dass dein da heißt (mit h als Variable). Von "anderem" Ansatz zu sprechen ist also klar verfehlt - es ist derselbe, nur mit einer konstanten R-Verschiebung (was den Aufschrieb etwas verlängert). Mit statt sowie statt geschrieben lautet die richtige DGL nach der ersten Integration mit eben dann der Lösung .
Hochmut kommt vor dem Fall. |
||||||||
06.10.2015, 08:09 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Hmm ok ... Diese Konstante macht also alles sehr viel komplizierter, das hatte ich übersehen. Wäre ja auch etwa zu einfach alles sonst. Oder die Realität ist etwas zu schwierig, wie man es nimmt! Schade, ich fände explizite Lösungen schön! (Ist ja ein relativ grundlegendes Problem aus der Newtonschen Gravitationsphysik.) Wahrscheinlich keine Chance, weitere Lösungen aus dem Spezialfall mit Integrationskonstante=0 zu konstruieren, oder? Ich habe schon verschiedenes probiert aber nichts gefunden. Bei linearen DGLs 1. Ordnung könnte man z.B. Variation der Konstanten benutzen, Gibt es Ähnliches eigentlich auch für (bestimmte, z.B. vorliegende) nichtlineare DGLs? |
||||||||
06.10.2015, 08:14 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Und welche DGL beschreibt das? |
||||||||
06.10.2015, 08:25 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Das Originalproblem hat die Anfangsbedingungen sowie . Bezogen auf
geht das nur mit , also gar nicht. |
||||||||
06.10.2015, 08:37 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Ich meinte: Lässt sich unter Umständen z.B. eine (wie auch immer geartete) Transformation des Spezialfalls in die ursprüngliche Bwgl 2. Ordnung einsetzen, um weitere Lösungen (mit anderen Anfangsbedingungen) zu erhalten? Ich habe z.B. schon den Ansatz versucht, aber wieder eine komplizierte nichtlineare DGL erhalten. Oder meintest du, dass auch diese Transformationen die Anfangsbedingungen unmöglich erfüllen können? |
||||||||
06.10.2015, 09:00 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Wenn du eine Anfangsbedingungen und hast mit , dann gilt deine Lösung. In allen Fällen mit lässt es sich auf meine Lösung oben zurückführen (mit einer "Zeitverschiebung"). Bleibt noch der Fall , da läuft die Integration wohl etwas anders (ich vermute mal, mit oder statt o.ä., hab es jetzt aber nicht durchgeführt). Jedenfalls entspricht der Grenzfall der Fluchtgeschwindigkeit (zweite kosmische Geschwindigkeit). |
||||||||
06.10.2015, 09:11 | Ehos | Auf diesen Beitrag antworten » | ||||||
Das ursprüngliche Problem von @DrummerS ist das sogenannte Zweikörperproblem (oder Kepler-Problem), dessen Lösung man z.B. im Buch "Mechanik" von L.D.Landau findet. Die Lösung ist nicht trivial und lässt sich als Integral darstellen. Im genannten Buch wird sogar der allgemeine Fall behandelt, wo beide Körper einen Drehimpuls haben dürfen und sich auf einer Ellipse umkreisen (also mit variablem Abstand). Im hier vorliegenden Fall verschwindet der Drehimpuls natürlich, was die Sache speziell macht. |
||||||||
06.10.2015, 09:16 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Klar ist es das Zweikörperproblem, und zwar der Spezialfall, dass und kollinear sind, womit das ganze in Dimension 1 abhandelbar wird. |
||||||||
06.10.2015, 12:17 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Was ich mich noch frage: Müsste das Ergebnis nicht mit dem Keplerschen Gesetz für die Umlaufzeiten kompatibel sein? |
||||||||
06.10.2015, 12:20 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Hmm wohl nicht da unklar ist, welcher Anteil der Umlaufzeit bis zum Aufprall benötigt wird. |
||||||||
06.10.2015, 13:34 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Die Frage ist, was du hier als Umlauf siehst: Wenn wir die Erde als "durchlässig" ansehen würden und nicht homogen, sondern als Massepunkt in der Mitte, dann stürzt die Kugel startend bei und Geschwindigkeit 0 mit zunehmender Geschwindigkeit bis zum Erdmittelpunkt, und dann mit abnehmender Geschwindigkeit darüber hinaus bis zur Tiefe und Geschwindigkeit 0 - und dann wieder zurück usw. periodisch. Sie beschreibt also im Kepler-Sinne eine entartete Ellipse (die eine Halbachse ist gleich Null), die Zeit für eine vollständige Periode ist dann gleich der vierfachen Zeit für in , das wäre dann , also tatsächlich wie im 3.Keplerschen Gesetz vermerkt. |
||||||||
06.10.2015, 14:08 | DrummerS | Auf diesen Beitrag antworten » | ||||||
Einen weiteren interessanten Aspekt an der anfangs von mir erwähnten Aufgabe fände ich, wenn der Betrag von zu einer sehr großen Relativgeschwindigkeit führt, bei der Relativitätseffekte (hier die relativistische Masse der Kugel) beginnen bedeutend zu werden. Im Abschnitt "Relativistische Masse" (Quelle) wird jedoch erwähnt, dass diese bei Berücksichtigung im newtonschen Gravitationsgesetz zu falschen Ergebnissen führen. Das ist dann aber wohl eher ein Thema für das Physikerboard. |
||||||||
06.10.2015, 14:16 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Das ist richtig: Die Rechnung hier ist nur gültig im Rahmen der Newtonschen, nichtrelativistischen Mechanik, d.h. für vergleichsweise kleine Geschwindigkeiten. Und wo Zeit noch die normale gemütliche Zeit ist, nicht irgendwie gedehnt... |
||||||||
06.10.2015, 14:21 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Klasse! |
||||||||
06.10.2015, 14:28 | klarsoweit | Auf diesen Beitrag antworten » | ||||||
Da kann man mal sehen, wo wir von der Frage nach der Berechnung eines Integrals nun gelandet sind. |
||||||||
06.10.2015, 15:47 | Mathema | Auf diesen Beitrag antworten » | ||||||
Da ich kein Physiker bin und mich die Berechnung des Integrals interessiert, gehe ich mal dort hin zurück. Man sollte ja nicht gleich die Flinte ins Korn werfen... Also - es geht doch um dieses Integral, oder? Das habe ich umgeformt zu: So: und Die Faktoren vor den Integralen lasse ich erstmal weg und bastel das später zusammen: und Resub.: Zusammenbasteln liefert mir: Da ich auch großen Respekt vor HAL`s mathematischen Fähigkeiten habe, vermute ich einen Fehler in meiner Rechnung?! |
||||||||
06.10.2015, 15:51 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Warum soll da ein Fehler sein? Für alle gilt . Ich hab nur deswegen statt gewählt, damit die Integrationskonstante für die vorliegende Anfangsbedingung "elegant" zu Null wird. |
||||||||
06.10.2015, 15:55 | Mathema | Auf diesen Beitrag antworten » | ||||||
Mist - da war ja was Diese "Eleganz" geht mir noch etwas ab. Danke dir! edit:
Das passt doch nicht meiner Meinung nach. Es müsste doch lauten: Und somit: Oder irre ich mich nun? |
||||||||
06.10.2015, 16:55 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Was sagt denn die Probe, also Differentation dieses Ausdrucks? Ich hab das oben (mit Ausnahme des arsin/arccos-Schlenkers) alles vom CAS machen lassen. ------------------------------------ Ein Anhaltspunkt liefern auch die physikalischen Einheiten - nicht für die Richtigkeit, wohl aber für die Falschheit einer Formel: hat die Einheit , natürlich die Zeiteinheit , und alles andere die Längeneinheit . Und in der Hinsicht sieht es dann nicht gut bestellt aus um deine letzte Formel... |
||||||||
06.10.2015, 17:01 | Mathema | Auf diesen Beitrag antworten » | ||||||
Ich geh noch mal auf Fehlersuche... |
||||||||
07.10.2015, 10:32 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Wikipedia sagt, die Lösung des allgemeinen Kepler-Problems wird numerisch berechnet, weil:
Kann jemand vielleicht mehr dazu sagen, woran dies liegt und/oder wie man das beweisen kann? |
||||||||
07.10.2015, 10:34 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Ups, da ist mir das Zitat missglückt! Kann das vielleicht jemand beheben? Vielen Dank! EDIT: erledigt (klarsoweit) |
||||||||
07.10.2015, 12:00 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Du redest jetzt über das -Körperproblem mit ? Denn für gibt es meines Wissens nach eine geschlossene Lösung: Bei nicht allzu großer kinetischer Energie umkreisen beide Körper in Ellipsenbahnen, die sich in derselben Ebene befinden, den gemeinsamen Masse-Schwerpunkt. An einem "Umkipppunkt" ist es dann eine Parabelbahn, und für noch größere Energie dann eine Hyperbelbahn. |
||||||||
07.10.2015, 12:49 | Huggy | Auf diesen Beitrag antworten » | ||||||
Die geschlossene Lösung bei gibt es für die Bahnkurve und die Zeit als Funktion des Ortes. Die Umkehrfunktion Ort als Funktion der Zeit ist auch in diesem Fall nur numerisch ermittelbar. |
||||||||
07.10.2015, 13:55 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Ja, das ist jetzt Definitionsfrage, was man als "geschlossene Lösung" akzeptiert: Zumindest auf trifft das ja zu, wie ja auch bei der eindimensionalen Vereinfachung hier im Thread zu sehen ist. |
||||||||
07.10.2015, 14:33 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Hmm physikalisch müsste es ja eine Lösung in allen Fällen geben, auch für n > 2. Kann man die Mathematik da vllt noch irgendwie verbessern? Viele physikalische Probleme sind ja auch erst durch DGLs behandelbar geworden. |
||||||||
07.10.2015, 14:41 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Ja klar gibt es eine Lösung, und numerisch kann man die ja auch bestimmen: Ich würde mal behaupten, der übergroße Anteil der Rechenkapazitäten der diversen Supercomputer bzw. Computer-Cluster beschäftigt sich mit der numerischen Lösung von Differentialgleichungen. Der eine oder andere davon vielleicht auch mit dem N-Körperproblem. Ich selbst hab mal vor ca. 20 Jahren ein diesbezügliches PC-Programm geschrieben, natürlich nur für kleines N. Bei gewissen Konstellationen der Anfangsbedingungen ergeben sich dann ganz lustige, chaotische Trajektorien. |
||||||||
07.10.2015, 14:50 | DieLösung | Auf diesen Beitrag antworten » | ||||||
Cool, das muss ich auch mal probieren. Ich frage mich aber trotzdem, ob in dreihundert Jahren eine einfache explizite Formel z.B. für r(t) im ZKP existieren wird. |
||||||||
07.10.2015, 15:28 | HAL 9000 | Auf diesen Beitrag antworten » | ||||||
Heutzutage kann man mit der beachtlichen Parallel-Rechenleistung moderner Grafikkarten schon auf einem normalen PC das etwas weiter nach oben schrauben: http://http.developer.nvidia.com/GPUGems3/gpugems3_ch31.html (Bitte jetzt nicht als Werbung für NVidia verstehen - mit den Grafikkarten des anderen großen Herstellers AMD/ATI geht's so ähnlich auch.) |
||||||||
09.10.2015, 17:04 | Mathema | Auf diesen Beitrag antworten » | ||||||
Nur noch mal der Vollständigkeit halber:
Hab an dieser Stelle den roten Teil verschlampt. Dann kürzt sich zum Schluss der Nenner raus und es bleibt übrig:
Und das passt denn zur von HAL präsentierten Lösung. Manchmal ist man doch blind und es braucht etwas Abstand um den Fehler zu finden. |
|
Verwandte Themen
Die Beliebtesten » |
Die Größten » |
|
Die Neuesten » |
|