gravitation_-_planetenbahnen
Unterschiede
Hier werden die Unterschiede zwischen zwei Versionen angezeigt.
Beide Seiten der vorigen RevisionVorhergehende ÜberarbeitungNächste Überarbeitung | Vorhergehende Überarbeitung | ||
gravitation_-_planetenbahnen [2024/01/13 16:25] – [Algorithmus] torsten.roehl | gravitation_-_planetenbahnen [2024/01/13 17:37] (aktuell) – torsten.roehl | ||
---|---|---|---|
Zeile 36: | Zeile 36: | ||
| ..... | ..... | ... weitere simulierte Datenpunkte| | | ..... | ..... | ... weitere simulierte Datenpunkte| | ||
|x(tn)| y(tn)| tn | | |x(tn)| y(tn)| tn | | ||
- | + | |//Tabelle 1: Fiktive Tabelle zur Darstellung der Planetenbahn// | |
- | //Tabelle 1: Fiktive Tabelle zur Darstellung der Planetenbahn// | + | |
<WRAP center round help 100%> | <WRAP center round help 100%> | ||
Es stellt sich also die Frage, wie wir die benötigten Koordinaten erzeugen können. | Es stellt sich also die Frage, wie wir die benötigten Koordinaten erzeugen können. | ||
Zeile 74: | Zeile 73: | ||
** | ** | ||
Heuristische Herleitung** | Heuristische Herleitung** | ||
- | {{ : | + | |{{ : |
- | //Die Diskretisierung entspricht einem Übergang vom Differentialquotienten zum Differenzenquotienten. Die Schrittweite Δt darf dabei nicht zu groß gewählt werden.// | + | |//Die Diskretisierung entspricht einem Übergang vom Differentialquotienten zum Differenzenquotienten. Die Schrittweite Δt darf dabei nicht zu groß gewählt werden.//| |
- | ---- | ||
Die letzten Gleichungen (Gl. 8 und Gl. 9) liefern uns die Positionen (x und y) in Abhängigkeit | Die letzten Gleichungen (Gl. 8 und Gl. 9) liefern uns die Positionen (x und y) in Abhängigkeit | ||
Zeile 106: | Zeile 104: | ||
|vy0|Y-Komponente der Geschwindigkeit zum Startzeitpunkt|vy0|3.43066 <color # | |vy0|Y-Komponente der Geschwindigkeit zum Startzeitpunkt|vy0|3.43066 <color # | ||
|Δt|Schrittweite|timeStep|0.005 < | |Δt|Schrittweite|timeStep|0.005 < | ||
- | // | + | |// |
{{: | {{: | ||
Zeile 116: | Zeile 114: | ||
{{: | {{: | ||
- | |1 AE|= $150\cdot 10^9} m|→ 1 m = $.\frac{1}{150\cdot 10^9}$ AE| | + | |1 AE|= $150\cdot 10^9$ m|→ 1 m = 1150⋅109 AE| |
|1 M|= 2⋅1030 kg|→ 1 kg = 12⋅1030 M| | |1 M|= 2⋅1030 kg|→ 1 kg = 12⋅1030 M| | ||
- | |1d|= 86400 s = 0.005 Δt|→ 1 s =$\frac{0.005}{86400} Δt| | + | |1d|= 86400 s = 0.005 Δt|→ 1 s =$\frac{0.005}{86400}$ Δt| |
Den neuen Zahlenwert erhalten wir durch Einsetzen: | Den neuen Zahlenwert erhalten wir durch Einsetzen: | ||
- | $G = 6.67 \cdot 10^{-11}\cdot \frac{ m^3}{ kg \cdot s^{2}} \rightarrow \\ 6.67\cdot10^{-11} | + | $G = 6.67 \cdot 10^{-11}\cdot \frac{ m^3}{ kg \cdot s^{2}} \rightarrow \\ 6.67\cdot10^{-11} |
+ | |||
+ | |||
+ | {{: | ||
+ | Welchen Zahlenwert | ||
+ | \\ | ||
+ | <color # | ||
+ | \\ | ||
+ | \\ | ||
+ | \\ | ||
+ | ^Gleichung^Bedeutung^ | ||
+ | |r=√x2+y2 |Abstand Sonne-Planet| | ||
+ | |vx(t+Δt)=vx(t)−(GMxr3)Δt |x-Komponente der Geschwindigkeit| | ||
+ | |vy(t+Δt)=vy(t)−(GMyr3)Δt |y-Komponente der Geschwindigkeit| | ||
+ | |x(t+Δt)=x(t)+vx(t+Δt)Δt |x-Komponente des Ortes| | ||
+ | |y(t+Δt)=y(t)+vy(t+Δt)Δt |y-Komponente des Ortes| | ||
+ | |Tabelle: Gleichungen für die Simulation|| | ||
+ | == Algorithmus == | ||
+ | - Wähle geeignete Startbedingungen und initialisiere alle Parameter (insbesondere die Schrittweite Δt) | ||
+ | - Zeichne die Position x(0) und y(0) in ein geeignetes Koordinatensystem. | ||
+ | - Berechne den Abstand Sonne-Planet. | ||
+ | - Berechne die Geschwindigkeiten vx und vy für den nächsten Zeitschritt. | ||
+ | - Berechne mit dem **Euler-Cromer-Algorithmus** die Positionen von x und y zum nächsten Zeitschritt. | ||
+ | - Zeichne die ermittelten Positionen x und y in ein geeignetes Koordinatensystem | ||
+ | - Wiederhole ab Schritt 3. bis das Programm abgebrochen wird oder die gewünschte Zahl an Iterationen erreicht wurde. | ||
gravitation_-_planetenbahnen.1705163148.txt.gz · Zuletzt geändert: 2024/01/13 16:25 von torsten.roehl