摘要:天体视位置是进行天文定姿的必备数据,目前视位置计算的研究多基于地面和海上,天基平台的发展需进一步研究。基于新天文参考框架提出一种太阳系内天体视位置计算模型,可计算任意时刻天体的地心视位置和站心视位置。与天文年历地心视位置相比,外层行星类天体的赤经差值在10ms以内,赤纬差值在220mas以内;与STK仿真的站心视位置相比,赤经差值在0.15s以内,赤纬差值在2''以内。实验表明,该模型可以正确计算测瞬太阳系内天体的视位置,其精度满足天文定姿需求。
关 键 词:天文定姿;视位置;天文年历;STK
星敏感器利用天体视位置和星空图像矢量信息进行天文定姿。天文定姿利用恒星、日月以及行星等自然天体作为导航信标确定航天器在轨道坐标系中的三轴姿态,天文定姿包括初始模式和跟踪模式两种方式[1]。首先,星敏感器对视场内天体成像,经过星图处理和星点提取等技术手段得到天体的图像矢量信息,对星图上的星点数据进行星图匹配,识别出天球坐标系中对应的天体,将这些天体的位置归算到观测历元下测站对应的视位置,天体视位置和天体的图像矢量信息可以确定航天器姿态矩阵,解算出航天器三轴姿态信息。因此,天体视位置是天文定姿的必备数据,为提高定姿精度,确保天文定姿的实时性要求,需对恒星以及太阳系内天体的实时视位置计算进行研究。
目前为止,国内外关于天体视位置计算的研究主要应用于地基和海基观测。文献[2]基于DE405历表数据,采用天文矢量处理方法,提出了一种高精度行星视位置计算的数学模型,主要是运用于舰船平台。文献[3]归纳并建立了常数变易法和坐标摄动法两种计算行星视位置的方法,在船舶天文导航系统中使用。文献[4][5]介绍了基于地面观测的新型天文测量系统中行星视位置的计算方法。文献[6]通过理论分析得出一种行星视位置长时段更新算法,但未考虑载体运动产生的影响。文献[7][8]对恒星视位置的计算做出相关研究,对于太阳系内天体视位置的计算有一定的参考作用。
天文定姿时星敏感器位于高速运动的航天器上,与地基海基的运用情况最大的差别就是测站本身产生的光行差量级较大,不可忽略(实验结果显示影响达角秒量级)。本文研究了新天文参考框架下天文定姿时太阳系内天体视位置计算,并建立了一种实用解算模型。与2011年、2013年、2014年和2015年的天文年历上天体视位置进行比对,验证地心视位置计算的精确度,利用STK仿真软件模拟卫星的在轨运行,通过得到的位置和速度信息计算天体的序列站心视位置,与STK中的太阳系内天体视位置数据比对分析站心视位置计算的准确度。
1.太阳系内天体视位置转换
1.1新天文参考框架
新天文参考框架规定了新的坐标系和岁差章动模型,天体视位置计算需做相应的改变。国际天文协会(International Astronomical Union,IAU)1991年引入以广义相对论为基础的国际天球参考系(International Celestial Reference System,ICRS),定义了两套运动学无旋转时空坐标系太阳系质心天球参考系 (Barycentric Celestial Reference System,BCRS)和地心天球参考系(Geocentric Celestial Reference System,GCRS)。03年开始采用新的运动瞬时参考框架,该框架同动力学参考架没有关系,由新的岁差章动模型所确定。IAU2000岁差章动模型描述了天球中间极CIP在GCRS中的运动[9][10]。
(1)P03岁差模型
天文学中采用过几种岁差模型,1896年以来一直使用的是纽康岁差系统,1986年后开始采用IAU1976岁差模型,2000年后采用IAU2000岁差模型,2006年的IAU大会决定2009年1月1日之后采用最新的P03模型。P03岁差模型的岁差量如式(1)
式中T为距里历元J2000.0的儒略世纪数。
(2)IAU2000章动模型
1984年章动的计算采用的是IAU1980章动理论,2005年其正式采用IAU2000章动模型。其中,IAU2000章动模型分为2000A和2000B模型,A模型精度最高,优于0.2mas,它包含了678个日月项和687个行星项,总共1365项,数量上是过去模型的十多倍。为了简化计算,提高解算效率减少计算时间,对A模型进行了简化取其中的78项组成模型B,它的精度在1995-2050年之间可以保持在1mas。
式中各参数详细见文献[11]。
1.2视位置转换
DE/LE历表里包含了太阳系内所有天体的基本信息,天文定姿时太阳系内天体的测瞬站心视位置计算需要充分考虑各种因素的影响。1)天体的自行(proper motion):从岁首历元时刻到测瞬时刻时间长度的天体自行,一般投影到星空中自行的影响小于1(''/年);2)视差(parallax):两个不同地点看到的物体的视方向之差,视差通常用来换算天体之间的距离,它的量级小于1'';3)光行差(aberration):由观测者与观测对象之间的相对运动和光速有限引起的位置视觉差异,在地面进行天体的观测时光行差的影响可以达到21'',而航天器在轨道上的运行速度对视位置的计算影响更大,是研究星敏上恒星视位置计算与地面观测之间的主要差别;4)还有光引力偏折、岁差和章动等。图1给出了天文定姿中天体各种位置之间的转换关系。
图1 各种视位置关系
Fig.1 Relationship between variety apparent position
2.视位置计算模型
同恒星的视位置计算相比,太阳系内天体的视位置计算最大的区别处在于需要考虑由光行差引起的位置影响。因为测站同被观测天体距离近相对运动速度大,故需要考虑周日视差和周日光行差改正。具体计算步骤如下:
(1)计算观测瞬间的地球时TT对应的儒略日和TDB时间
鉴于早时所使用的历书时存在不连续和不一致等问题以及原子频标技术的高速发展,1976年起采用力学时以广义相对论作为时间工作的理论基础。目前采用两种力学时,太阳系质心参考星BCRS中用质心力学时TDB,地心视位置的解算用地球时TT。TT是建立在国际原子时TAI的基础上,两者关系如式:
LS为跳秒,可根据IERS上面发布的公报B参数数据确定跳秒的值。
(2) 通过DE/LE历表内插得到J2000历元下太阳系内各天体质心的位置矢量和地球在TDB时刻对应的BCRS中的质心位置和速度,是地球绕太阳的公转速度,单位为AU/日。
(3)计算航天器TT时刻在真赤道真春分点坐标体系下相对于地心的位置和速度,以km为单位,以km/s为单位,从而得到航天器在BCRS中的位置和速度
(4)计算太阳系天体的光行时。光行时的近似值为:
根据DE405历表获取时太阳系内天体的质心位置矢量,则光行时改进值为
依此循环迭代,直到光行时满足我们的精度要求。
(5)计算时太阳系内天体相对于T时航天器的坐标方向的单位矢量:
(6)进行光线引力弯曲改正、周年光行差改正和岁差和章动改正:
(7)将位置矢量转换为球面坐标:
(8)加入周日视差改正:
为行星的周日地平视差,可根据历表获取;为测站天文纬度,t为天体的时角。