Hallo Byanus, hallo Ralf,
der Optimierungsalgorithmus für die beste Bahn hatte sich festgefressen. Ich habe mittels einen näherungsweisen Korrelationsmatrix die Optimierung fortgesetzt und habe eine Lösung mit deutlich kleinerem Residuum gefunden.
Statt ein Residuum von 0.1912 Bogensekunden zu den Beobachtungen ist das Residuum der neuen Lösung nur 0.1388 Bogensekunden. Leider hat die beste Lösung eine Exzentrizität von mehr als 1. Deswegen habe sie offenbar die Exzentrizität beim MPC festgelegt.
Auch wird die Monte Carlo Simulation schwieriger, da in einen größeren Abstand von der besten Lösung gesucht werden muss und weniger gebundene Lösungen vorkommen. Daher habe ich diesmal 55 Millionen Bahnen für die Monte-Carlo-Simulation verwendet, von denen nur wenige Tausende gebunden waren.
Nach der Wichtung waren 4% der (gewichteten) Bahnen waren gebunden, 96% über Fluchtgeschwindigkeit.
Die Ergebnisse der Monte-Carlo-Simulation sind aber nahezu die gleichen wie zuvor.
Für die Inklinationen gilt nun
i: Mittelwert = 39°, Median = 23°, 95%CI = [12 - 149]°
Nur im Fall, wenn ich die Verteilung der Inklinationen berücksichtige, dass es nur wenige mit hoher Inklination von mehr als 30° gibt ändert sich das Ganze. Bahnen mit großer Inklination, dist, kleiner Geschwindigkeit und kleiner Periheldistanz werden unwahrscheinlich. Damit wird 2018 VG18, falls es nicht ein seltenes Objekt mit hoher Inklination ist, wahrscheinlich sein Perihel nicht innerhalb der Neptunbahn haben.
Für q<30 AE ist die Wahrscheinlichkeit nun nur noch 3%, für q<50 AE 16%.
Die große Halbachse a und die Apheldistanz Q lässt sich nur schlecht bestimmen. Allerdings sollten sich Objekte mit sehr großem Q nur selten nahe genug an der Sonnen aufhalten um entdeckt zu werden, wodurch die Wahrscheinlichkeit für ein großes Q geringer sein sollte, als berechnet.
Etwa 14% der (gewichteten) Bahnen waren gebunden 86% über Fluchtgeschwindigkeit.
Bester Orbit mit Residuum 0.1388 Bogensekunden zu Beobachtungen:
dist = 119.0778 AE, v = 4.2430 km/s, a = -288.1 AE, e = 1.1943, q = 55.982 AE, Q = -632.2 AE, i = 15.3981°
Jetzt Mittelwert, Median und 95% Konfidenzintervall (0.025 - 0.975) für die Parameter aus der Monte-Carlo Simulation.
dist: Mittelwert = 119.1 AE, Median = 119.0 AE, 95%CI = [114.6 - 124.2] AE
v: Mittelwert = 3.17 km/s, Median = 3.25 km/s, 95%CI = [2.03 - 3.88] km/s
a: Mittelwert = ca 1900 AE, Median = 204 AE, 95%CI = [87 - ca 6800] AE
e: Mittelwert = 0.61, Median = 0.64, 95%CI = [0.17 - 0.981] AE
q: Mittelwert = 78 AE, Median = 80 AE, 95%CI = [28 - 116] AE
Q: Mittelwert = ca 3800 AE, Median = 325 AE, 95%CI = [123 - ca 13500] AE
i: Mittelwert = 16.3°, Median = 15.4°, 95%CI = [11.0 - 27.0]°
Allen ein frohes Weihnachtsfest!
UMa