段淑珍,郭 英,韩晓冬,王 琦,毕京学

(山东科技大学测绘科学与工程学院,山东 青岛 266590)

0 引 言

全球卫星定位系统的迅速发展,引起了各国军事部门和广大民用部门的普遍关注，同时,也对GPS卫星星历的精度和实时性提出了越来越高的要求。星历按照精度可分为精密星历和广播星历，由于精密星历只能在卫星观测的11天后获得,无法为实时定位、精确导航等提供有效服务,所以要实时计算卫星的位置只能利用GPS卫星广播星历,其中含有开谱勒轨道参数、轨道摄动参数等信息[1]。由于用广播星历每一个小时发布一次,每2 h每颗卫星有一组独立的参数,用相邻两组星历计算同一时刻的卫星位置,将会有不同的结果,因此，计算卫星轨道是一项繁琐的过程[2-4]。

本文主要根据Bernese软件叙述了利用FORTRAN语言编译的卫星轨道程序处理广播星历的步骤,利用自编程序处理广播星历并与Bernese GPS 5.0软件处理的广播星历和精密星历的轨道坐标进行精度比较，并得出结论。

1 卫星轨道的概述

由瑞士伯尔尼大学天文研究所研究开发的GPS数据处理软件——Bernese GPS 4.2软件,是目前国际上比较成熟的高端GPS精密应用软件之一，升级后的 5.0 版本既可以进行 GPS 精密定位定轨、电离层和对流层监测应用, 也可以应用到低轨卫星定轨和卫星钟差的计算,以及高精度时频比对等[5]，由于Bernese GPS软件的源代码是开放的,也便于为大学、研究机构和高精度的国家测绘机构等单位进行研究以及二次开发。

其中轨道部分的源代码与其他部分相对独立,便于进行研究。……