基于多系统组合的卫星定位系统.docxVIP

基于多系统组合的卫星定位系统 0 多系统组合组合观测系统 全球卫星网络系统（gnss）已广泛应用于许多领域，如大地测量、地壳运动监测、精确测量和导航。因此,包括我国在内的许多国家都在建设自己的全球定位系统,当今运营和在建的全球卫星导航系统除GPS和GLONASS外,还包括GALILEO及我国的北斗全球导航系统。从目前GNSS工程应用来看,单一系统存在以下问题,即在城市高楼密集区和位于深山峡谷的水库、电站、矿山等环境下采用单卫星系统进行定位时,受到地形及周围环境的影响,用户接收机可见的卫星数目少且分布不佳,导致系统的可用性降低,无法满足定位的最低要求。卫星精密测量应用于形变监测领域时,其三维定位的精度和可用性等方面还存在一些缺陷,表现为垂直方向上的定位误差通常为水平方向的2~3倍,不能满足大坝、桥梁和高层建筑等安全监测的要求。 为解决单一导航系统存在的以上缺陷,可采用地面增强技术。但是,地面增强系统受远近效应、对流层延迟、多路径效应、时间同步、成本高和维护不便等一系列问题的影响,不能很好的解决上述问题。近年来,随着GPS的现代化、GLONASS的逐渐完善和GALILEO计划的启动,可用卫星数量将大大增加,因此,采用多系统组合是解决以上问题的一条有效途径。由于GPS、GLONASS和GALILEO采用不同的轨道参数,对其进行有针对性的组合,可有效提高覆盖区域的一致性,使组合系统在全球范围内有更好的覆盖;多系统组合可以增加可观测卫星的数量;增强观测卫星的几何图形强度;提高整个卫星定位系统的可靠性和可用性。此外,采用多系统组合使位于两极的用户同样可以接收到卫星的信号,这样高纬度甚至两极用户也可用之定位。多系统有更多的可用频率,用来消除电离层延迟影响;并可以利用不同系统的数据消除或减弱系统本身误差以及传播误差的影响;不同系统间频率的重合可以使用差分的方法消除和减弱与频率有关误差的影响;采用多频组合能够大大提高相位模糊度的固定效率。针对以上问题,利用GPS、GLONASS、GALILEO的轨道参数模拟卫星星座,通过计算地面格网上观测到的卫星站心坐标及其高度角;引入精度因子作为GNSS系统可用性的衡量指标,对GNSS及其组合系统的可用性加以分析。 1 卫星站中心坐标研究和可用性评估指标 1.1 星历参数仿真计算 根据卫星的轨道参数,设to为参考时刻,根据开普勒轨道参数A√A、e、i0、ω、Ω0、M0分别为参考时间的轨道半长轴的平方根、轨道偏心率、轨道倾角、轨道近地点角距、轨道升交点赤经、卫星平近点角。l、Ω、Δn分别为轨道倾角、轨道升交点赤经、卫星平近点角随时间的变化率。Cuc、Cus、Crc、Crs、Cic、Cis分别为卫星相对于升交点角距、轨道半长轴和归到角调和改正项的振幅。坐标系统采用WGS-84坐标系,GM为地球引力常数,a为椭球的半长轴,e为轨道偏心率,ωe为地球自转角速度。 用星历参数仿真计算t时刻卫星位置的步骤如下: 1) 计算t时刻卫星的真近点角u和卫星相对升交点角距: 计算平均运动角速度n和平近点角m: n=GM/a3??????√(1)n=GΜ/a3(1) m=m0+n·Δt(2) 式中Δt=t-t0,t为当前时间,t0为历元时间。 利用开普勒方程迭代方式计算偏近点角,当|Ei+1-Ei|=10-12时停止迭代。 Ei+1=Ei?(Ei?e?sinEi?m)1?e?cosEi,E0=m(3)Ei+1=Ei-(Ei-e?sinEi-m)1-e?cosEi,E0=m(3) 计算真近点角u和升交点角距φ: u=arctan(1?e2√sinE1?e?cosE,cosE?e1?e?cosE)(4)u=arctan(1-e2sinE1-e?cosE,cosE-e1-e?cosE)(4) φ=u+ω. (5) 2) 计算摄动改正项: δr=Crccos2φ+Crssin2φ(6) δi=Ciccos2φ+Cissin2φ(7) δu=Cuccos2φ+Cussin2φ(8) 3) 计算改正后的向径、真近点角、轨道倾角和升交点经度: r=a(1-ecosE)+δr(9) u=u+δu(10) i=i0+i·Δt+δi(11) Ω=Ω0+(Ω-ωe)Δt-ωetoe(12) 式中ωe为地球自转角速度。 4) 计算卫星在轨道平面中的位置 x=rcosU,y=rsinU(13) 5) 计算卫星在地心坐标系中的位置 ???XYZ???=(?Ω)(?i)???xy0???=[xcosΩ?ycosisinΩxsinΩ+ycosicosΩ](14)[XYΖ]=(-Ω)(-i)[xy0]=[xcosΩ-ycosisinΩxsinΩ+ycosicosΩ](14) 1.2 gnss及其组合系统的发挥 可用性是系统在某一指定覆盖区域内提供