行星位置计算.doc

日月、行星位置计算 高健北京师范大学 天文系关键词：、 1.1 计算步骤概要 1、输入观测时间（格林尼治UTC），并给定测站在WGS-84下的大地经纬度L、B和大地高H。利用表，在站心J2000.0平赤道坐标系下的在站心J2000.0平赤道坐标系下的并转换得到观测瞬时在站心地平坐标系下的方位角和地平高度。输出观测瞬时，从测站观测到的的位置坐标（站心J2000.0平赤道坐标系下的赤经α赤纬δ，站心地平坐标系下的方位角和地平高度）.2 日月、行星轨道运动特性及位置来源 太阳系是由太阳、8颗大行星和一些矮行星（如冥王星）以及它们的卫星、众多的小天体（小行星、彗星和流星体）以及行星际物质组成的天体系统。地球是其中的一颗大行星，而月球是地球唯一的自然卫星。太阳是太阳系中最重要的天体，其他天体都在它的引力作用下运动。大行星在太阳的引力作用下，以椭圆轨道绕太阳公转，而太阳位于椭圆轨道的一个焦点上。虽然行星绕太阳运动的规律已经研究的比较清楚，但由于行星之间也存在着引力作用，与其它太阳系天体之间也有相互的摄动作用，因此行星绕太阳公转运动的轨道是在不断变化的，简单的二体运动轨道和行星平根数不能满足高精度研究的需要。而要精确地描述行星绕太阳的运动，需要建立完善的运动方程。虽然建立这样的多体运动方程并不太难，但要求解这一运动方程是需要花费许多人力和物力的。美国喷气推进实验室(JPL)[1]。自1984年以来，JPL的行星/月球历表已经成为历算工作中广泛使用的历表，同时也广泛应用于空间导航、行星和月球探测中。 我们采用JPL的行星和月球历表DE405，其位置精度完全满足项目所要求的[1]。可以提供在历表所包含时间范围内任意时刻的日、月和行星的位置和速度。DE405是JPL所提供的精度最高的行星和月球历表，包含了从公元前1599年12月9日到公元2201年2月20日范围内的日、月和大行星位置和速度，并且采用了IAU所推荐的国际天球参考系ICRS[2]，可以为本项目提供精确的日、月及大行星的位置和速度数据。有关JPL-DE405历表的介绍和使用方法，在多个书籍和文献中都有介绍，这里不再详述，请相关网站及文献，如/、/ [1][3][4]。 图一 日月和行星站心位置计算流程图 1.3 日月和行星的位置计算 图一中给出了本方案计算日月和行星的站心位置的流程图。具体步骤参考Kaplan等人[5][6]所给出的天体位置计算方法，以及参考文献[4][8]。从JPL-DE405行星和月球历表中能得到任意时刻日月和行星的太阳系质心坐标(BCRS下的位置和速度)。要得到从地面测站观测时，日月和行星的J2000.0站心平赤道坐标，需要将日月和行星的位置从BCRS变换到J2000.0平赤道坐标系[7]，并且将坐标系原点从太阳系质心平移到地面测站。.3.1测站的J2000.0质心平赤道位置和速度 观测瞬时()，观测者位于地理经纬度L、B大地高H的测站，测站相对于太阳系质心的位置和速度可由下列步骤计算[6]： I．计算观测瞬时(TT)测站的地方恒星时。 从§A.2可知，从UTC时间容易转换得到TT时间和TDB时间。已知测站的地理经纬度L、B和大地高H据§A.2，易得观测瞬时()对应的测站地方恒星时II．计算观测瞬时(TT)，测站在观测瞬时地心真赤道坐标系相对于地球质心的位置矢量和速度矢量。 已知测站的地理经纬度L、B和大地高H，以及地方恒星时，有： ， a是地球赤道半径，为赤道处的地球自转，f为参考椭球椭率，为z轴方向的单位矢量。这一步并未考虑地球极移带来的影响。极移引起的地面测站地心位置的变化约为10m（0″.3），考虑离地球最近的月球，极移带来的月球地心位置的变化也仅为几个毫角秒 III．计算测站的J2000.0地心平赤道位置矢量和速度矢量。观测瞬时(TT)地心真赤道的位置矢量和速度矢量测站的J2000.0地心平赤道位置矢量和速度矢量 为星表参考历元，即J2000.0，这里为时间尺度为从J2000.0到观测瞬时()的岁差转换矩阵，为观测瞬时()的章动转换矩阵。的具体形式请参看附录§A.3及相关参考文献[8][9]。 IV. 地心相对于太阳系 利用JPL-DE405历表可直接获得(TDB)，地球质心相对于太阳系质心的位置和速度．计算测站的J2000.0平赤道位置和速度。(TDB)，地球质心相对于太阳系质心的位置和速度 其中，，为观测瞬时(TDB)，测站的J2000.0地心平赤道位置和速度 至此，得到观测瞬时(TDB)，测站相对于太阳系质心的J2000.0平赤道和表示。 1.3.2 日月和行星地心位置的计算(J2000.0地心平赤道) 日月和行星在BCRS下的位置和速度利用常矩阵B[7]变换得到在J2000.0平赤道惯性系下的位置和速度。以行星为例，介绍日月和行星地心位置计算过