Название говорит об этом, правда. Насколько я понимаю, наше суммарное пекулярное движение относительно реликтового излучения составляет 631 (+/-20) км/ч, в направлении 10,5ч RA; -24 градус декабря (точнее, далеко от 22,5h RA; +24 градус декабря, верно? Это близко к направлению Лямбда Пегаса).
Итак, если мы проведем линию от Луны до лямбда Пегаса, пересечет ли ее когда-нибудь Земля?
Я думаю, что ответ отрицательный, но мне не хватает модели, чтобы доказать это самому себе.
Наклон лунной орбиты составляет 5,14 градуса, а Земля простирается примерно на 1,2 градуса от Луны, но траектория Луны проходит на 24 градуса выше эклиптики. Значит, Земля не пересечет его, верно?
Все это предполагает, что приведенная выше векторная информация является общей и уже учитывает движение солнца вокруг Млечного пути. Если это не так, пожалуйста, дайте мне знать.
По крайней мере, не раньше, чем через плюс/минус несколько сотен лет! Плоскость орбиты Луны вокруг Земли (рассматриваемая в инерциальной системе отсчета) немного качается взад и вперед ( прецессия лунных узлов с периодом около 18,6 лет), но здесь это не поможет.
Вы можете играть с числами в этой программе.
import numpy as np
import matplotlib.pyplot as plt
from skyfield.api import Loader
degs = 180./np.pi
load = Loader('~/Documents/yournamehere/SkyData')
data = load('de405.bsp')
ts = load.timescale()
timez = []
dayz = np.arange(1, 31, 0.1)
dayz = np.arange(1, 366, 0.1)
#for year in range(1700, 2200):
for year in range(2017, 2028):
#for year in range(2000, 2030):
#time = ts.utc(year, 7, dayz) # July
time = ts.utc(year, 1, dayz) # January
timez.append(time)
earth = data['earth']
moon = data['moon']
RAdegs, Decdegs = [], []
for time in timez:
RA, Dec, dist = earth.at(time).observe(moon).radec()
RAdeg, Decdeg = degs*RA.radians, degs*Dec.radians
# lazy way to block the wrapparound
dRAdeg = RAdeg[1:] - RAdeg[:-1]
dropit = np.abs(dRAdeg) > 100
RAdeg[:-1][dropit] = np.nan
Decdeg[:-1][dropit] = np.nan
RAdegs.append(RAdeg)
Decdegs.append(Decdeg)
RAdegpts = [157.5, 157.5+180] # coordinates from the question
Decdegpts = [-24.0, +24.0]
if 1 == 1:
plt.figure()
for RAdeg, Decdeg in zip(RAdegs, Decdegs): #[::5]: # plot every 5th year
plt.plot(RAdeg, Decdeg)
plt.scatter(RAdegpts, Decdegpts)
plt.title("2017 throuth 2027 continuously", fontsize=20)
plt.show()
ПрофРоб
Дэвид Кертис