Ein Kostenvoranschlag ist in Ordnung. Aber wenn möglich, suche ich nach einer Referenz, die ich zitieren kann.
Dies ist der Anfang einer Teilantwort, ich freue mich, wenn jemand sie übernimmt oder ergänzt.
Warum nicht alle sechs Missionen statt nur einer?
Unter Verwendung von Skripten aus diesen 1 , 2 habe ich die Entfernung zwischen den Erd- und Mondzentren als schwarze Linie gezeichnet, dann sind die dicken Farbbalken diese Entfernungen minus dem Erdradius plus dem Radius des Mondes plus einer Höhe von 183 auf der anderen Seite km.
Die Längen der Balken sind die ungefähren Gesamtmissionslängen und daher müssen wir die Zeiten jedes Ein- und Ausstiegs beim Passieren hinter dem Mond hinzufügen.
Ich bin mir nicht sicher, wo ich diese finden kann, aber wenn mich jemand in die richtige Richtung weisen kann, kann ich die Nummern kratzen und sie hier hinzufügen.
from skyfield.api import Topos
from skyfield.api import Loader
import numpy as np
import matplotlib.pyplot as plt
from skyfield.api import load
loaddata = Loader('~/Documents/fishing/SkyData') # avoids multiple copies of large files
ts = loaddata.timescale() # include builtin=True if you want to use older files (you may miss some leap-seconds)
eph = loaddata('de421.bsp')
earth, moon, venus = [eph[x] for x in ('earth', 'moon', 'venus')]
apollos = [(10, 1969, 5, 18, 26), (11, 1969, 7, 16, 18),
(12, 1969, 11, 14, 24), (13, 1970, 4, 11, 17),
(14, 1971, 1, 31, 40), (15, 1971, 7, 26, 38),
(16, 1972, 4, 16, 27), (17, 1972, 12, 7, 19)]
# https://en.wikipedia.org/wiki/Apollo_program
timez_apollo = []
for n, year, month, d_start, d_stop in apollos:
times = ts.utc(year, month, range(d_start, d_stop+1))
timez_apollo.append(times)
days = 1 + np.arange(5*365.2564+1)
times = ts.utc(1969, 1, days)
years = days/365.2564
t_1969 = times.tt[0]
earthpos, moonpos = [x.at(times).position.km for x in (earth, moon)]
r = np.sqrt(((moonpos - earthpos)**2).sum(axis=0))
Rmoon, Rearth = 1737., 6378.
alt = 183.
fig = plt.figure()
ax1 = fig.add_subplot(1, 1, 1)
ax1.plot(years, r, '-k', linewidth=0.5)
for timez in timez_apollo:
yearz = (timez.tt - t_1969) / 365.2564
earthpoz, moonpoz = [x.at(timez).position.km for x in (earth, moon)]
r = np.sqrt(((moonpoz - earthpoz)**2).sum(axis=0))
d = r + Rmoon + alt - Rearth
ax1.plot(yearz, d, linewidth=2.5)
ax1.set_xlim(0.2, 4.0)
ax1.set_ylim(357500, 405500)
ax1.set_ylabel('km')
ax1.set_title('Earth Moon and Apollo far-side orbit to Earth surface distances')
ax1.set_xlabel('years since 1969-01-01')
plt.show()
Jo
äh
äh
SE - hör auf, die Guten zu feuern
äh
1969-07-16 00:00 UTC
und endet am1969-07-18 00:00 UTC
. Mit einer horizontalen Skala von 6 Monaten pro Zentimeter ist das Diagramm nicht zum Abrufen präziser numerischer Größen gedacht. Das heißt, Ihre Nummer liegt weit unter der gesamten Linie von Apollo 11 (orange), also werde ich nachforschen. Bleiben Sie dran...äh
SE - hör auf, die Guten zu feuern