gps星历的轨道误差分析.docxVIP

gps星历的轨道误差分析

一、卫星星历的精度和实时性

gps（全球定位系统）是美国开发的下一代卫星网络系统。由于它具有整个领域、全球范围和测站之间不需要通视、观测速度快、观测精度高等优点，因此在测量、导航、气候等领域得到了广泛应用。同时,也对GPS卫星星历的精度和实时性提出了越来越高的要求。星历按照精度可分为精密星历和广播星历。精密星历是由国际GPS服务中心(IGS)通过Internet(互联网)发布,它的轨道精度可达到±10cm左右,足以满足精密定位的需要,但是精密星历只能在卫星观测的11天后获得,无法为实时定位、导航、气象等实时性要求很强的应用提供有效的服务。广播星历是通过接收机接收卫星发射的含有轨道信息的导航电文,经过解码获得的卫星星历推算得到卫星位置,精度只在±100m左右,无法满足许多GPS应用的需求。本文将尝试对广播星历的轨道误差进行探讨,这对精确确定卫星的位置,提高定位精度是很有必要的。

二、开普勒轨道参数定义

GPS卫星的广播星历提供16个星历参数,其中包括1个参考时刻、6个相应参考时刻的开普勒轨道参数和9个轨道摄动修正参数。这些参数的定义见表1。

用广播星历参数计算观测时刻t的卫星位置的步骤如下。

1.偏近点角t-toe

计算平均运动角速度n和平近点角M

n=GM/a3??????√?M=M0+(n+Δn)(t?toe)n=GΜ/a3?Μ=Μ0+(n+Δn)(t-toe)

利用开普勒方程迭代公式计算偏近点角,当|Ei+1-Ei|ξ=10-12时停止迭代

Ei+1=M+esinEi,E0=M

计算真近点角u和升交点角距φ

u=arctan1?e2√sinEcosE?eu=arctan1-e2sinEcosE-e

φ=u+ω

2.计算摄动纠正项目

δr=Crccos2φ+Crssin2φ

δi=Ciccos2φ+Cissin2φ

δu=Cuccos2φ+Cussin2φ

3.计算出的方向、真近角、轨道倾角和上升点的高度

4.计算卫星在轨道平面上的位置

x=rcosu,y=rsinu

5.计算卫星在地固系wgs84中的位置

X=xcosΩ-ycosisinΩ

Y=xsinΩ+ycosicosΩ

Z=ysini

三、广播星历和精密星历的坐标矢量模的差值

GPS接收机接收到的广播星历一般每隔2h为一组,为了保持卫星坐标的精度,一般采用限制星历外推时间间隔的方法。为此本文选择外推的时间间隔最大不超过1h。卫星的坐标是按照观测时刻最接近的一组广播星历数据来计算的,每次计算时都要进行时间比较,利用不同时刻的广播星历参数,最后根据上面的轨道坐标计算公式求出各个时刻卫星的位置。

本文对2002年7月11日～15日的GPS广播星历进行计算,将获得的卫星轨道坐标与对应时刻的GPS精密星历的轨道坐标进行比较。在假设精密星历为精确值的前提下,先读取广播星历中的开谱勒轨道参数,用上面的轨道坐标计算公式进行计算求得每颗卫星每15min的轨道坐标;再从精密星历中读取相对应的轨道坐标;最后计算两者矢量模的差值

Δr=(Δx)2+(Δy)2+(Δz)2??????????????????√Δr=(Δx)2+(Δy)2+(Δz)2

现以一号卫星为例,横轴以15min为间隔,纵轴以米为单位,根据求得结果绘制5日的卫星坐标矢量模的差值图(如图1所示)。

作者从现在正在使用的28颗卫星中选取6颗,这6颗卫星分别位于GPS星座的6个轨道面内。分析发现广播星历与精密星历的坐标矢量的模的差值一般达到米级,但由于各种原因如:电离层、太阳黑子、月球引力、气候等,基本上每颗卫星每天在某个时段的差值都会有一些跳跃,其跳跃值可达几十米,甚至可达几百米。现把7月11日～15日计算的6颗GPS卫星坐标矢量模的差值的均值列于表2。

分析表2可以看出,7月11日～13日这3天的广播星历与精密星历的坐标矢量模的差值达到十几米以上,而14日和15日两天的差值则在几米的范围内。这是由于在11日～13日3天间爆发了一次小规模的太阳耀斑,太阳耀斑爆发的直接效应是电离层的突然扰动。这样,异常的电离层活动使得卫星轨道的位置发生偏离,但是在按照轨道模型确定未来的广播星历时不能很好的预见太阳耀斑事件的发生,使得轨道预测发生了一定的误差。所以,卫星轨道的位置会随着电离层的变化而变化。在电离层活动比较稳定的情况下,卫星会沿着轨道运行,广播星历与精密星历的坐标矢量模的差值也就比较小;相反,在电离层活动异常情况下,卫星的位置将会偏离轨道位置,使得差值变的很大。

图2是分布在GPS星座的6个轨道面上的6颗卫星的广播星历与精密星历的坐标矢量模的差值图。从图2中可以看出每天各个卫星的坐标矢量模的差值Δr明显地表现出周期变化趋势,并且周期与卫星的运行周期11h58