行星位置计算.doc

日月、行星位置计算 高健 北京师范大学 天文系 摘要：本文给出从地面站观测时，日、月和行星的位置计算过程。 关键词：历表：JPL-DE405；行星：视位置计算 ============================================================= 一、日月、行星位置的计算 1.1 计算步骤概要 1、输入观测时间（格林尼治UTC），并给定测站在WGS-84下的大地经纬度L、B和大地高H。 2、利用JPL-DE405月球/行星历表，得到观测瞬时日月和行星在太阳系质心坐标系（ICRS）下的位置和速度，并计算得到日月和行星在站心J2000.0平赤道坐标系下的位置和速度。同时，也计算日月和行星到测站的距离。 3、利用日月和行星在站心J2000.0平赤道坐标系下的位置和速度转换得到日月和行星在站心真赤道坐标系下的位置和速度，并转换得到观测瞬时日月和行星在站心地平坐标系下的方位角A和地平高度h。 4、输出观测瞬时，从测站观测到的日月和行星的位置坐标（站心J2000.0平赤道坐标系下的赤经α和赤纬δ，站心地平坐标系下的方位角A和地平高度h）。 1.2 日月、行星轨道运动特性及位置来源 太阳系是由太阳、8颗大行星和一些矮行星（如冥王星）以及它们的卫星、众多的小天体（小行星、彗星和流星体）以及行星际物质组成的天体系统。地球是其中的一颗大行星，而月球是地球唯一的自然卫星。太阳是太阳系中最重要的天体，其他天体都在它的引力作用下运动。大行星在太阳的引力作用下，以椭圆轨道绕太阳公转，而太阳位于椭圆轨道的一个焦点上。虽然行星绕太阳运动的规律已经研究的比较清楚，但由于行星之间也存在着引力作用，与其它太阳系天体之间也有相互的摄动作用，因此行星绕太阳公转运动的轨道是在不断变化的，简单的二体运动轨道和行星平根数不能满足高精度研究的需要。而要精确地描述行星绕太阳的运动，需要建立完善的运动方程。虽然建立这样的多体运动方程并不太难，但要求解这一运动方程是需要花费许多人力和物力的。美国喷气推进实验室(JPL)使用最新确定的天文常数和天体运动理论计算了太阳、月球、地月质心、水星、金星、火星、木星、土星、天王星、海王星和冥王星的位置，以及地球黄经章动和黄赤交角章动的数值 REF _Ref262590231 \r \h [1]。自1984年以来，JPL的行星/月球历表已经成为历算工作中广泛使用的历表，同时也广泛应用于空间导航、行星和月球探测中。 我们采用JPL的行星和月球历表DE405来获得行星和月球的位置信息，其位置精度完全满足项目所要求的 REF _Ref262590231 \r \h [1]。它可以提供在历表所包含的时间范围内任意时刻的日、月和行星的位置和速度。DE405是目前JPL所提供的精度最高的行星和月球历表，包含了从公元前1599年12月9日到公元2201年2月20日范围内的日、月和大行星位置和速度，并且采用了IAU所推荐的国际天球参考系ICRS REF _Ref262590274 \r \h [2]，可以为本项目提供精确的日、月及大行星的位置和速度数据。 有关JPL-DE405历表的详细介绍和使用方法，在多个书籍和文献中都有介绍，这里不再详述，请查阅相关网站及文献，如/、/ ，参考文献 REF _Ref262590231 \r \h [1] REF _Ref262590336 \r \h [3] REF _Ref262590454 \r \h [4]。 JPL DE405 JPL DE405 日、月和行星ICRS下位置、速度 输入观测时刻UTC 计算对应的UT1，TT，TDB 修正周年视差 修正光行差 观测瞬时日、月和行星J2000.0站心平赤道位置、速度矢量 日、月和行星J2000质心平赤道位置、速度 修正周日视差 修正光行时 输出观测瞬时日月和行星在J2000.0站心平赤道坐标系下赤经α和赤纬δ，站心地平坐标系下的方位角A和地平高h。 测站WGS-84经纬度和大地高 观测瞬时日、月和行星站心真赤道位置、速度矢量 观测瞬时日、月和行星的站心地平坐标 极移 恒星时 章动转换 岁差转换 岁差转换 章动转换 观测瞬时测站的J2000.0地心平赤道位置、速度矢量 观测瞬时测站的J2000.0质心平赤道位置、速度矢量 修正光线引力弯曲 图一 日月和行星站心位置计算流程图 1.3 日月和行星的位置计算 图一中给出了本方案计算日月和行星的站心位置的流程图。具体步骤参考Kaplan等人 REF _Ref262590624 \r \h [5] REF _Ref262590631 \r \h [6]所给出的天体位置计算方法，以及参考文献 REF _Ref262590454 \r \h [4] REF _Ref262597247 \