Ist die Halo-Umlaufzeit des JWST mit der elliptischen Umlaufbahn der Erde um die Sonne synchronisiert?

Ist die Halo-Umlaufzeit des JWST mit der elliptischen Umlaufbahn der Erde um die Sonne synchronisiert?

Wow, es ist überraschend nah!

Zuerst dachte ich, dass es aufgrund früherer Informationen und Fragen und Antworten, die auf einer früheren JWST-Flugbahn basierten, näher an 6,5 oder 7 Monaten liegen würde, aber wow, beim Versuch, die Periode einer komplizierten Umlaufbahn mit 3 oder 4 Körpern (einschließlich des Mondeffekts) zu definieren. Auf einfache Weise bin ich auf 180 Tage gekommen!

Ich habe die neue vorhergesagte Umlaufbahn für JWST in Horizons genommen

Revised: Jan 28, 2022
2Y_SCHEDULE_2022027000000_01U.OEM.V0.1   2022-Jan-27 00:01  2024-Jan-27 00:01

lud die heliozentrischen Positionen von JWST und Erde herunter, subtrahierte das Erd-/Mond-Schwerpunktzentrum von JWST, drehte es dann um den Erdwinkel, um in einem pseudo-synodischen Rahmen zu sein (es dreht sich mit dem EM/Bary-Winkel, nicht stetig), und erhielt Folgendes (Daten beginnen mit 2022-Feb-08, endet 2024-Jan-27):

Der Versatz in X ist die durchschnittliche Entfernung von der Erde, nahe den 1,5 Millionen km der klassischen L2-Entfernung.

Wenn ich die Maxima herausgreife und einen Zeitraum für die X-, Y- und Z-Bewegung berechne, erhalte ich 179,3, 180,0 und 180,7 Tage!

Da dies eine Ad-hoc -Methode ist, sage ich, dass diese Zeiträume *derzeit nicht von einem halben Jahr zu unterscheiden sind , daher ist meine Antwort ein eingeschränktes, aber ziemlich überraschtes "Ja!"

import numpy as np
import matplotlib.pyplot as plt

fnames = ('Earth orbit heliocentric horizons_results.txt',
          'JWST orbit heliocentric horizons_results.txt')

JDs, datas = [], []
for fname in fnames:
    n_offset = 10
    with open(fname, 'r') as infile:
        lines = infile.readlines()
    a = [i for (i, line) in enumerate(lines) if 'SOE' in line][0]
    b = [i for (i, line) in enumerate(lines) if 'EOE' in line][0]
    lines = lines[a+1+n_offset: b]
    info = [line.split(',') for line in lines]
    JD = np.array([float(line[0]) for line in info])
    data = np.array([ [float(thing) for thing in line[2:8]] for line in info])
    JDs.append(JD)
    datas.append(data)

JD = JDs[0]
Earth, JWST = [thing.T.copy() for thing in datas]

dJWST = JWST - Earth

z = dJWST[2]

th = np.arctan2(Earth[1], Earth[0])

s, c = [f(-th) for f in (np.sin, np.cos)]

x = dJWST[0] * c - dJWST[1] * s
y = dJWST[1] * c + dJWST[0] * s

fig, axes = plt.subplots(3, 1)
for thing, name, ax in zip([x, y, z], 'XYZ', axes):
    ax.plot((JD - JD[0]) / 365.2564, thing)
    ax.set_ylabel(name + ' (km)')
plt.show()

maxima = [np.where((p[2:] < p[1:-1]) * (p[1:-1] >= p[:-2])) for p in (x, y, z)]

maxima = [thing[0] for thing in maxima]

periods = [(m[-1] - m[0]) / (len(m) - 1) for m in maxima]

print(periods)