1时间坐标系统转换方法研究1.1不同时间类型研究内容中涉及到7种不同时间类型，分别是协调世界时（UTC）、地球动力学时（TT）、国际原子时（TAI）、太阳系质心动力学时（TDB）、地心坐标时（TCG）、GPS时（GPST）和北斗时（BDT）。UTC是协调世界时，协调世界时的秒长严格等于原子时的秒长，而协调世界时与世界时UT间的时刻差规定需要保持在0.9s以内，否则将采取闰秒的方式进行调整。闰秒一般发生在6月30日及12月31日。地球动力学时（TDT）是建立在国际原子时TAI的基础上的，其秒长与国际原子时相等。1991年，第21届IAU大会决定将地球动力学时（TDT）改称为地球时（TT）。地球时（TT）和国际原子时（TAI）之间的关系式可以表示为：32.184TTTAIs（1-1）国际原子时间（TAI），是地球上的时间基准，它由国际时间局从多个国家的原子钟分析得出，被定义为：32.184()TAITTsUTC跳秒（1-2）太阳系质心动力学时有时也被简称为质心动力学时。这是一种用以解算坐标原点位于太阳系质心的运动方程（如行星运动方程）并编制其星表时所用的时间系统。质心动力学时（TDB）和地球时的（TT）之间没有长期漂移只有周期项变化，即0.001658sinsTDBTTM0e20.000014sin2()sMVXXc（1-3）其中M为地球绕日公转的平近点角；eV为地球质心在太阳系质心坐标系中的公转速度矢量；0X为地心在太阳系质心坐标系中的位置矢量；X为地面钟在太阳系质心坐标系中的位置矢量；0XX实际上就是地面钟在地心坐标系中的位置矢量；c为真空中的光速。地心坐标时（TCG）是原点位于地心的天球坐标系中所使用的第四维坐标—时间坐标，用于讨论绕地球运行的卫星等天体的运动规律、编制相应的星历。GPST是由GPS星载原子钟和地面监测站原子钟组成的一种原子时基准，与国际原子时保持19秒的常数差，起始历元时间是1980年1月6日0时0分0秒（UTC）。北斗卫星导航系统的系统时间称北斗时间BDT，北斗时属于原子时，起始历元时间是2006年1月1日0时0分0秒（UTC）。其中这7大类时间之间的转换原理如如图1-1所示：TAI+ΔUT-ΔATUTCUT1TT+ΔUT1+ΔT+ΔTT+32.184sTCGTDBGPSTBDT图1-1各个时间类型转换原理图儒略日（JDate）是一种不涉及年、月等概念的长期连续记日法。根据公历年（Y）、月（M）、日（D）来计算对应的儒略日。79275*JD=1721013.5+367*Y-int{*[Y+int()]}+d+int()412249MhM（1-4）式中常数1721013.5为公元1年1月1日0时的儒略日；Y、M、D分别为公历中的年月日h为世界时的小时数。简化儒略日（ModJDate）是IAU在1973年采用的一种更为简便的连续计时法，定义为：MJD=JD-2400000.5（1-4）1.2不同坐标系五种不同的地球坐标系Fixed、MeanofDate、TrueofDate、ICRF和J2000。Fixed是地心地固坐标系；MeanofDate是瞬时平天球坐标系，TrueofDate为瞬时真天球坐标系，ICRF是国际天球参考架，它是天球位置参照的理想质心坐标系，相对于遥远河外天体是运动学不转动的。J2000是平赤道地心坐标系，原点在地球质心，xy平面为J2000时刻的地球平赤道面，x轴指向J2000时刻的平春分点。地心地固定系转换到J2000需要经过四个步骤：Fixed经过极移修正后得到真地球坐标系，真地球坐标系绕Z轴旋转格林尼治真恒星时（GAST）后可转换到TrueofDate；TrueofDate经过章动改正可得到MeanofDate；MeanofDate经过岁差改正得到GCRS，GCRS与J2000间相差一个常值偏差矩阵B。坐标系转换如图1-2所示：CoordinateTransformationOfEarthCoordinateSystemFixedTrueofDateMeanofDateGCRSJ2000PolarMotionAndGASTNutationPrecessionArrayOfDeviation图1-2Fixed转换J2000步骤以IAU2000岁差模型和IAU2000章动模型为例，在某历元UTC时刻，ITRS到GCRS的转换矩阵可写成：()()()()GCRSITRSITRSTrQtRtWtrHGtr（1-5）其中，ITRSr、GCRSr分别对应同一向量在ITRS和GCRS坐标系中的坐标，上式中，()Wt、()Rt和()Qt分别对应极移，自转和岁差章动转换矩阵。在计算()Rt和的时候，会有两种计算方法，我们分别称为CEO-based转换方法和equinox-based转换方法。其中前者为IAU提出的新的计算方法。a.极移矩阵()Wt'()()()()zypxpWtRsRxRy（1-6）上式中，s’为：-0.047mas.t。极移量的求解为：(,)(,)(,)(,)ppIE...
