导航卫星的广播星历拟合是卫星导航系统中一个重要的技术环节,影响到用户的定位精度。本文推导了星历参数的拟合算法,并且通过算例得到了较高精度的参数拟合结果。参数的拟合结果中大数量级的拟合参数达到了1.0e+06及以上的拟合精度,小量级的参数拟合精度受到计算机浮点数的精度影响较大,有待进一步提高。文中分析了不同拟合历元数对参数拟合结果的影响,通过探测突变的拟合历元,有助于寻找最佳的拟合时长。 Broadcast ephemeris fitting of navigation satellites is an important technical aspect in the satellite navigation system that affects the user’s location accuracy. Ephemeris parameters fitting algorithm is derived, and the higher precision of parameter fitting with examples is obtained. Fitting arguments of large magnitude in parameters fitting result achieve the fitting accuracy of 1.0e+06 and above, and minor-level parameter fitting precision greater affected by computer floating- point precision needs to be further improved. This paper analyzes the different fitting epoch number effects on the results of parameter fitting, by detecting mutation of fitting epoch, helpful to find the best fitting duration.
郭宋1,2,苏晋展1,2
1同济大学测绘与地理信息学院,上海
2现代工程测量国家测绘地理信息局重点实验室,上海
收稿日期:2016年7月7日;录用日期:2016年7月23日;发布日期:2016年7月29日
导航卫星的广播星历拟合是卫星导航系统中一个重要的技术环节,影响到用户的定位精度。本文推导了星历参数的拟合算法,并且通过算例得到了较高精度的参数拟合结果。参数的拟合结果中大数量级的拟合参数达到了1.0e+06及以上的拟合精度,小量级的参数拟合精度受到计算机浮点数的精度影响较大,有待进一步提高。文中分析了不同拟合历元数对参数拟合结果的影响,通过探测突变的拟合历元,有助于寻找最佳的拟合时长。
关键词 :拟合广播星历,最佳拟合时长,浮点数位数,小量级参数拟合
在GPS定位中,为了确定地面点的坐标,必须先确定卫星的坐标。卫星坐标的计算主要依靠卫星的星历。GPS星历分为精密星历和广播星历两类。它们各有优缺点,前者精度较高,但是存在延时,主要用于事后处理,例如IGS精密星历,通常采用SP3格式,以卫星在地固坐标系中离散点的坐标(和速度)提供给用户使用;后者的经度较低,但是能够实时转发给用户,用于实时的导航和定位。广播星历由GPS的地面运行控制中心站计算得到,以广播星历参数的形式提供给用户使用 [
目前,许多国家都在建设自己的卫星导航系统,导航卫星的广播星历拟合是卫星导航系统中一个重要的技术环节 [
GPS广播星历具有16个参数,包括星历参考时刻
GPS星历表中的时间和坐标分别属于GPS时间系统和WGS-84坐标系。在GPS数据处理中,常利用卫星的星历表参数来计算卫星天线相位中心在WGS-84坐标系中的位置,其具体计算过程及相应的表达式如下 [
1) 半长轴:
2) 卫星平均角速度:
3) 瞬时历元到参考历元的时间差:
4) 平角速度:
5) 平近点角:
6) 偏近点角
7) 真近点角
8) 维度参数:
9) 周期改正项:
10) 改正后的纬度参数:
11) 改正后的向径:
12) 改正后的倾角:
13) 升交点的经度:
14) 卫星在轨道平面内的坐标:
15) 卫星在WGS-84坐标系中的坐标:
其中
将上节中卫星的坐标表示成位置向量:
因此广播星历计算卫星坐标可以表示为如下形式 [
其中
由(3)式可得,单历元的误差方程为:
其中
令
其中
式中
根据(5)、(6)、(7)式进行间接平差迭代求解。迭代的收敛条件为:
其中,
B阵为(3)中观测方程系数的矩阵形式,可以通过解析法求偏导计算,但解析法较为复杂不利于计算机运算,本文采用数值导数的方式计算。
由偏导数的定义:
当
结合本文系数阵B阵的实际,则有
因此误差方程系数阵
其他偏导系数同理可依次求出,此处略。
本文采用GPS星期第1765周第一天即2013年11月03日的IGS精密星历作为观测值,拟合当日格林尼治时间12 h参考时刻的广播星历参数。IGS精密星历提供的卫星坐标时间间隔为15分钟。由计算模型可知,只要历元数大于5便能根据最小二乘原理进行间接平差。本文用不同的历元数进行拟合计算,并比较不同拟合历元数所得广播星历的精度。将拟合的广播星历参数计算参考历元时刻的卫星坐标与精密星历对应时刻的坐标差异作为广播星历精度的评价指标。最后将最优的拟合参数与当日卫星所发送的广播星历进行比较,说明实际的应用意义。
15个广播星历参数中9个摄动参数的初值设为0,6个轨道根数的初值如表1所示。
本文数值导数的增量为1.0e−08。
参数 | ||||||
---|---|---|---|---|---|---|
初值 | 5154 | 0.9 | 0.01 | 2 | 2.5 | 2.5 |
表1. 六个轨道根数的初值
不同拟合历元,各参数的变化情况略有不同,但在某一特定历元大部分参数出现了曲线的突变,这样的变化有助于我们找出最佳的拟合历元。下面以卫星11号,2013年11月03日格林尼治12:00的广播参数拟合为例,说明各参数不同拟合历元的参数变化,各参数的变化曲线见图1。
图1. (a) 广播星历参数不同拟合历元数的参数结果变化曲线(1~9参数)。其中第1到第9参数分别是:
从各参数的变化情况可以看出,参数变化的出现在两个历元出。其一为9个拟合历元处,另一个为17个拟合历元处。其中在9个拟合历元处出现转折的参数为:
为了比较拟合广播星历的实际应用价值,下面将拟合计算所得的参数与相同参考时刻相同卫星发布的广播星历进行比较,比较结果见表2。
表2中给出了卫星11号2013年11月03日格林尼治12:00为参考时刻的,15个拟合历元数过拟合的结果与当日卫星广播的结果的比较。从表中可以看出,拟合参数均达到了较高的精度,证明算法的正确性以及使用性。同时结果中存在不足的是,对于极小量级的参数的拟合效果不佳。这个主要原因为计算机浮点数的精度有限,因此为了提高参数的拟合参数因该加大浮点数的位数,以获得高精度的拟合结果。
本文推导了拟合广播星历的方法,并且在求取误差方程系数阵时引入数值导数的方法。最后通过算例证明了方法是正确的。
通过算例的分析可知,不同的拟合历元数所得参数的结果具有差异,并且参数值在特定的拟合历元数出现明显的突变与转折。该转折的具体原因有待进一步的研究。不过笔者断定转折有助于我们决定最
参数 | 拟合值 | 广播星历值 | 差值 | 相对误差 |
---|---|---|---|---|
1 | 5.1536774269E+03 | 5.1536772938E+03 | −1.3307419249E−04 | −2.58E−08 |
2 | 7.0354957738E−09 | 7.0088633759E−09 | −2.6632397903E−11 | −3.80E−03 |
3 | −3.7920565558E−01 | −3.7920578950E−01 | −1.3391738113E−07 | 3.53E−07 |
4 | 1.4751157084E−02 | 1.4751165989E−02 | 8.9049194511E−09 | 6.04E−07 |
5 | 1.2355322795E+00 | 1.2355323861E+00 | 1.0657108063E−07 | 8.63E−08 |
6 | 3.6510600404E−06 | 3.6619603634E−06 | 1.0900322987E−08 | 2.98E−03 |
7 | −1.2675470388E+01 | −1.2625000000E+01 | 5.0470387678E−02 | −4.00E−03 |
8 | 1.5255678019E−07 | 1.6391277313E−07 | 1.1355992938E−08 | 6.93E−02 |
9 | −5.1516177681E−07 | −5.0105154514E−07 | 1.4110231665E−08 | −2.82E−02 |
10 | 2.6970820145E+02 | 2.6965625000E+02 | −5.1951445255E−02 | −1.93E−04 |
11 | −9.4441570740E−08 | −1.0803341866E−07 | −1.3591847916E−08 | 1.26E−01 |
12 | 8.8906911228E−01 | 8.8906914338E−01 | 3.1097491582E−08 | 3.50E−08 |
13 | −1.7800478874E−10 | −1.8679349498E−10 | −8.7887062406E−12 | 4.71E−02 |
14 | 1.6303709456E+00 | 1.6303710645E+00 | 1.1895197471E−07 | 7.30E−08 |
15 | −9.6087944750E−09 | −9.5975426333E−09 | 1.1251841732E−11 | −1.17E−03 |
表2. 拟合参数与卫星广播参数的比较
佳的拟合历元数。最后本文将拟合的结果与实际的广播参数进行了比较。从比较的结果中可以看出,应用本文所述的方法得到了较高精度的广播星历参数,具有很高的实际运用价值。但是极小量级的参数收到了计算机浮点数精度的影响,未能达到很高的精度。所以应该增加浮点数的位数,以获得高精度小量级参数。
国家自然科学基金,编号:41201425。
郭宋,苏晋展. GPS广播星历的拟合与拟合历元数的影响分析The Influence Analysis of GPS Broadcast Ephemeris Fitting and Fitting Epoch Number[J]. 测绘科学技术, 2016, 04(03): 100-107. http://dx.doi.org/10.12677/GST.2016.43012