为了精确计算某个时刻月球的位置,需要在月球的经度、纬度和距离中考虑数百个周期项。如果考虑最重要的周期项上海厚成股票配资案,可以让月球经度达到的精度,纬度达到的精度。
这里计算以地心为原点的月球中心的黄道坐标,黄经,黄纬,及其与地球中心的距离。由月球的赤道面视差得到
\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在中额外加上(金星、地球扁率、木星的影响)在中额外加上
所以月球的黄道坐标是
由此可以转换为赤道坐标
图片上海厚成股票配资案
本站仅提供存储服务,所有内容均由用户发布,如发现有害或侵权内容,请点击举报。启泰网提示:文章来自网络,不代表本站观点。