天文算法笔记（6）月球位置的计算

2025-10-26 05:48

为了精确计算某个时刻月球的位置，需要在月球的经度、纬度和距离中考虑数百个周期项。如果考虑最重要的周期项，可以让月球经度达到的精度，纬度达到的精度。

这里计算以地心为原点的月球中心的黄道坐标，黄经，黄纬，及其与地球中心的距离。由月球的赤道面视差得到

\sin\pi=\frac{6378.14}{\Delta}

首先计算时刻（相对于J2000.0，以世纪为单位）

T=\frac{{\rm JD}-2451545}{36525}.

然后计算一些中间量，这些量都以度为单位（如果值太大，可以除以360取余数，也就是对360取模），月球的平均经度

L'=218.3164477+481267.88123421 T-0.0015786 T^2 +T^3/538841-T^4/65194000

月球的平均角距（mean elongation）为

D=297.8501921+445267.1114034 T-0.0018819 T^2+ T^3/545868-T^4/113065000

太阳的平近点角（mean anomaly）

M=357.5291092 + 35999.0502909 T-0.0001536 T^2 +T^3/24490000

月球的平近点角

M'=134.9633964+477198.8675055 T+0.0087414 T^2 +T^3/69699-T^4/14712000

月球和升交点的平均距离

F=93.2720950 + 483202.0175233 T-0.0036539 T^2-T^3/3526000 +T^4/863310000

此外还需要三项

A_1=119^{\circ}.75+131^{\circ}.849 TA_2=53^{\circ}.09+479264^{\circ}.290 TA_3=313^{\circ}.45+481266^{\circ}.484 T

然后根据下面的表计算周期项、和。表中给出了的系数，举个例子，如果表中的系数组是，对应的的系数是，那么需要计算的对的贡献就是。

然而，含有角的项依赖于地球绕太阳轨道的偏心率，这个值目前随时间减小，为此，含有或的项的幅度需要乘以一个改正因子，含有或的项的幅度要乘以，

上面提到的那一项应该变为

表1

的系数加到的sin项的系数加到的sin项的系数6288774-209053551274027-3699111658314-2955968213618-569925-18511648888-114332-31495879324615857066-15213853322-17073345758-204586-40923-129620-34720108743-303831047551532710321-12528109807966110675-3478210034-232108548-21636-788824208-676630824-5163-83794987-166754036-128313994-104453861-11650366514403-2689-7003-2602239010056-234863222236-9884-21205751-20692048-4950-17734130-15951215-3958-1110-8923258-8102616759-1897-713-2117-7002354691596549-1423537-1117520-1571-487-1739-399-381-4421351-340330327-32311652992948752

表2

的系数加到的sin项的系数51281222806022776931732375541346271325731719892668822821643244200-3359246322112065-18701828-1794-1749-1565-1491-1475-1410-1344-133511071021833777671607596491-451439422421-366-351331315302-283-229223223-220-220-185181-177176166-164132-119115107在中额外加上（金星、地球扁率、木星的影响）

在中额外加上

所以月球的黄道坐标是

由此可以转换为赤道坐标

图片

本站仅提供存储服务，所有内容均由用户发布，如发现有害或侵权内容，请点击举报。