月球位置计算.doc

月球位置计算 月球位置计算 月球的位置计算是所有天体计算里比较复杂的部分。往往为了提高一点计算的精确度，必须使用数百个扰动修正项，目前这些修正项在许多专业的天文书籍里都找得到。由於月球每天在天空移动的数度很快，因此计算时间的正确性要求的比较多。之前计算太阳与行星时，都采用地方平均太阳时，但是计算月球位置必须采用力学时来减低误差。这里提供的方法，月球位置的误差约小於1/15度。虽然在许多观测上，这个精密度可能稍嫌不足，但是却能够简单地推算月球的位置也算是实用的方法。 虽然月球绕行地球公转，但是严格说起来，月球应该看成绕太阳运行。由於受到地球的扰动，所以月球以蛇行的方式环绕太阳。月球的轨道面会产生很大的位移，这也是他与其他行星最大的不同点。因为迅速移动的轨道面，造成实际计算非常不易。大致来说，月球相对於黄道面的轨道进动周期为18.6年（如同地球的岁差，我们又称为升交点进动）。月球的近日点也会周期性的运动，主要的项次约为8.85年。另外月球绕地运行的轨道半径每年约增加数公分之多。 计算时间 时间 JD1= JD－JD2000.0 说明：这里的JD是使用力学时(ET)来计算，相关的资料请参考「时间」的网页。 距离2000年的儒略世纪 T＝（JD－JD2000）/ 36525 太阳的轨道资料 轨道周期 P＝1.00004回归年 距离春分点平均离角 ε＝280.466457＋0.985647358 JD1＋0.000304 T2 近地点离角 ω＝282.937348＋0.00004707624 JD1＋0.0004569 T2 轨道偏心率 E=00T 轨道半径 A＝1AU 计算结果 近日点平均离角 M＝ε－ω 利用级数展开的克卜勒方程式（Kepler’s Equation）求得真实的近日点离角 ν＝M＋360 E sin(M)/π＋900 E2 sin(2M)/4π－180E3sin(M)/ 4π 太阳的黄道经度 λ＝ν＋ω 月球的轨道资料（以黄道面来定义） 月球的平均经度 L＝218.316646＋13JD1－0.0014664 T2 月球的平均近日点离角 Mm＝L－（83.353243＋0.11140352394 JD1－0.0103217 T2） 升交点引数 N＝125.044555－0.05295376277 JD1＋0.0020756 T2 月面轨道倾角 I＝5 扰动修正项 Ev＝1.2739 sin( 2L-2λ－Mm ) Ae＝0.1858 sin( M ) A3＝0.37 sin( M ) 修正月球的平均近日点离角 Mm＝Mm＋Ev－Ae－A3 Ec＝6.2886 sin( Mm ) A4＝0.214 sin( 2Mm ) 修正月球的平均经度 L＝L＋Ev＋Ec－Ae＋A4 V＝0.6583 sin( 2L－2λ) 再次修正月球的平均经度 L＝L＋V 修正升交点引数 N＝N－0.16 sin( M ) 计算月球的黄道经度与纬度 sin(θ)＝sin( L－N ) cos( I ) cos(θ)＝cos( L－N ) 判断θ的象限后，求得θ。 月球的黄道经度 λm＝θ＋N 月球的黄道纬度 βm＝sin-1[ sin( L－N ) sin( I ) ] 若需要月球的赤道经度与纬度，可参考座标系统的变换部分。将黄道座标转成赤道座标。或者更进一步转成地平方为角。 范例 2000年9月1日 AM10:30 月球的黄道经纬度 、赤道经纬度 与地平方为角 查得当年力学时与平均太阳时的差值 ΔT＝64秒 ET=UT+64秒 计算ET儒略日得到 JD＝2451788.60492 JD1＝243.60492 儒略世纪 T＝0.00667 太阳轨道资料 ε＝520.492123 ω＝282.948812 E＝0.01671 太阳黄道位置计算结果 M＝237.543311 λ＝158.8947 β＝0 月球的轨道资料（以黄道面来定义） 月球的平均经度 L＝188.3759761 月球的平均近日点离角 Mm＝77 升交点引数 N＝112 月面轨道倾角 I＝5 扰动修正项 Ev＝－0 Ae＝－0 A3＝－0 修正月球的平均近日点离角 Mm＝77 Ec＝6 A4＝0 修正月球的平均经度 L＝194 V＝0 再次