2016, Vol. 37
依附于全球卫星导航系统(global navigation satellite systems，GNSS)的增强系统，能够减少飞行成本，降低恶劣天气对定位精度的影响，提高整个空中运输系统的效率.这要求基于卫星导航的增强辅助定位服务必须满足定位精度、完好性、连续性和可用性等要求.然而，利用惯性导航系统(inertial navigation system，INS)实现对卫星导航系统的辅助导航，还存在很多问题，不能够达到精密进近所需导航性能需求(required navigation performance，RNP)[1]，尤其是精密进近(precision approach)阶段.
为提高飞行阶段的各方面的性能，引入了卫星导航地基增强系统(ground-based augmentation systems，GBAS)提高机载设备的导航性能，包括导航定位精度[2].地基增强系统采用相位平滑伪距差分修正辅助定位算法[3 - 4]，这种定位方法相对于伪距定位提高了定位精度，同时不存在整周模糊度的解算问题[5].传统的相位平滑伪距算法假定了电离层延迟保持不变，然而这一点不总正确.首先如果电离层发生较快、较大的变化，如出现电离层闪烁[6]，就会出现两倍的电离层误差.其次，如果平滑滤波器长时间不重新设置，就会使得电离层的误差累计，出现较大的电离层的残差，导致系统的定位精度不高、载波相位平滑伪距的差分修正量误差较大[7 - 8].由于电离层误差受电离层刺穿点和天气状况等影响，动态的预测和滤波成为研究难点.
本文针对以上问题，以提高卫星导航定位精度为目标，提出了GBAS相位平滑伪距差分修正改进算法.
1 GBAS载波相位平滑伪距差分修正改进算法地基增强就是利用相位平滑伪距差分修正实现对星基导航的辅助定位，差分修正过程如图 1所示.
|
图 1 载波相位平滑伪距差分修正过程 Fig.1 Carrier phase smoothing pseudorange |
根据卫星播发报文中的信息，解调出码伪距测量值和载波相位测量值，则伪距方程为
|
(1) |
载波相位方程为
|
(2) |
其中：
利用载波相位差对每个接收机接收到的伪距作平滑处理[9].先得到t1时刻的伪距值和载波相位值：
|
(3) |
|
(4) |
得到相位观测方程后，对卫星信号连续观测，则t1，t2两时刻的相位测量之差为
|
(5) |
根据伪距方程，在t2时刻的伪距观测量为
|
(6) |
将式(3)，式(6)带入到式(5)中，得
|
(7) |
从式(7)的推导可知，t1时刻的伪距值等t2时刻的伪距值与两个时刻的载波相位差的差值和两倍的电离层误差的和.但是因为电离层对码伪距测量误差具有滞后作用，而对于载波相位测量具有超前作用，导致利用载波相位在伪距平滑过程中加剧了电离层误差，当电离层出现严重的抖动的时候，电离层延迟会严重影响定位精度[10].利用无码载偏离平滑(divergence free smoothing，DFS)算法，对电离层的双倍叠加做修正，修正流程如图 2所示.
|
图 2 无码载偏离平滑(DFS) Fig.2 Divergence free smoothing |
其中:
|
(8) |
式中：ηρ1，ηφ1分别为码相位和载波相位的滤波衰减因子；ηD为平滑伪距滤波器的衰减因子；F是一个低通滤波器.首先对信号传播中的误差值做预测模型，针对其前一时刻的误差模型对当前时刻做移步预测值.根据当前时刻预测到的误差值，给出相位平滑伪距中滤波器的衰减因子.然后对伪距和载波相位线性组合后的值相减，通过滤波器F，消除电离层偏差影响得到χ，将χ再输入滤波器消噪，消噪后得到无码载偏离影响的载波相位平滑伪距值
|
(9) |
利用迭代法求出k个历元时刻观测伪距值和k个历元与第一个历元的载波相位差值，带入式(9)能够得到k个关于t1时刻伪距的方程式.利用载波相位差值以及观测的各个时刻的载波相位差值求出k个t1时刻的载波相位平滑伪距值.对这k个载波相位平滑伪距值做最优化处理，设定一个阈值ρmi(t2)，ρmi(t2)为卫星到接收机的实际距离.设定θ为误差偏离度，其值如下：
|
(10) |
|
(11) |
得到的相位平滑伪距值应满足：
|
(12) |
如果不满足上述条件，则将其剔除，将剩下的相位平滑伪距值利用最小二乘法取平均，即
|
(13) |
通过对卫星导航定位误差的研究，首先根据北斗的报文播发参数以及北斗发布的修正公式，对每个接收机对应的每颗卫星的钟差参数及星历参数做修正.
通过载波相位平滑滤波对钟差和星历的修正后，得到的伪距信息如下：
|
(14) |
相位平滑后的伪距值减弱了在导航信号传播中的多径效应以及噪声的影响.为了辅助用户进行定位，对得到的相位平滑伪距差分修正，得到差分修正量如下：
|
(15) |
假设有M颗卫星，则有4M个差分修正量.根据上文得出的数据，对差分修正量求差分修正量误差的标准差σpr_gnd.为了排除非理想误差的干扰(如多径误差和电离层风暴等)，增大差分修正量的可用性，保证完好性，地面主控站除了通过报文广播差分修正量、伪距差分修正变化率外，还要广播每颗卫星的相位平滑伪距的差分修正量误差的标准差.以此来提高地基增强系统辅助精密进近的完好性水平和可用性水平.
将膨胀后的差分修正量、差分修正量的变化率以及差分修正量误差的标准差通过VDB发送给用户.用户利用收到信息修正用户的伪距信息：
|
(16) |
其中：Pcorrected为用户伪距修正信息；Pn为用户伪距；RRC*(t－tzcount)为平滑伪距修正率(t为接收到的伪距时间；tzcount为用户使用伪距修正量的时间)；Δtsv为用户卫星钟差修正.
电离层残差的多少直接影响定位精度，由图 3可知，平滑滤波后电离层残差已经控制在0～0.1m之间，降低了电离层误差对定位的影响.对总的误差(见图 4)分析可知，改进的相位平滑差分修正算法在总的定位误差上比传统的算法的总的误差要低，总体上提高了定位的精度.
|
图 3 电离层残差 Fig.3 Ionospheric residual |
|
图 4 平滑伪距差分修正误差 Fig.4 Smoothing pseudorange differential |
RMS即一组数据的均方根值，对于地基增强系统，主要通过播发差分修正信息来对飞机实现辅助定位，对差分修正后数据的RMS值分析，可以改进航空机载用户的精度、完好性、可用性.
由图 5～图 7可知，在X，Y，Z轴上，RMS值都有相应的降低，而且随着历元的增加，RMS值较改进前更加平稳.修正误差小于2m，达到Ⅱ类盲降标准，有利于差分修正值的完好性监测，使得差分修正量的可用性更高，定位更加稳定、精确.利用差分修正值的标准差进行地面段的完好性计算，对误差滤波修正后，由于差分修正量的变化稳定在2m之内，也就是说4台接收机接收到的差分修正值比较接近，差分修正量的标准差比较小，不容易超出完好性要求的阈值水平.将原来差分修正量标准差达到完好性阈值的时间提高了5%，提高了差分修正量的标准差对于达到完好性阈值的可用性.
|
图 5 X轴的RMS随历元的变化 Fig.5 RMS with the epoch change on the X axis |
|
图 6 Y轴的RMS随历元的变化 Fig.6 RMS with the epoch change on the Y axis |
|
图 7 Z轴的RMS随历元的变化 Fig.7 RMS with the epoch change on the Z axis |
通过对以上数据的分析，在相位平滑伪距差分定位中，误差修正滤波后的电离层误差显著降低，定位精度有所提升.同时载波相位平滑伪距的差分修正值在X，Y，Z方向上的RMS有所减小，并且随着历元的增加趋于平稳.滤波状态稳定，没有误差积累.平滑伪距差分修正量误差小于2m，精度稳定在2m之内，所以垂直精度小于2.5m，达到Ⅱ类精密进近着陆引导精度标准，提高了地基增强辅助定位系统的完好水平性和差分修正量以及差分修正量标准差的可用性水平.
| [1] |
Cassell R,Smith A.Development of required navigation performance (RNP) requirements for airport surface movement guidance and control[C]// Digital Avionics Systems Conference.Alexandria:IEEE,1995:57-64.
( 0)
|
| [2] |
Guenter D,Dennis J.Initial operational experience with CAT I ground based augmentation system (GBAS)[C]// IEEE Conference Publications.Herdon:IEEE,2015:1-14.
(0)
|
| [3] |
Wang S,Zhu Y,Li R.The effect of navigation frequency on smoothed pseudorange variance in SBAS[C]// 2011 IEEE 5th International Conference on Cybernetics and Intelligent Systems (CIS).Qingdao:IEEE,2011:283-287.
(0)
|
| [4] |
Reddy A S,Sriharsha M R,Sarma A D,et al.Performance evaluation of single and dual frequency carrier smoothing techniques for LAAS[C]// Applied Electromagnetics Conference (AEMC).Kolkata:IEEE,2011:1-4.
(0)
|
| [5] |
Xue R, Zhang J, Zhu Y B.
Cascade dual frequency smoothing for local area augmentation system[J]. Chinese Journal of Aeronautics , 2009, 22 (1) : 49–55.
DOI:10.1016/S1000-9361(08)60068-0 (0)
|
| [6] |
Malla R,Fries R.Multi-layer modeling and simulation of the effects of ionospheric scintillation on service availability of the GPS augmentation systems [C]//Position Location and Navigation Symposium (PLANS).Indian Wells:IEEE,2010:688-693.
(0)
|
| [7] |
刘玥.支持在航的载波相位差分定位方法研究[D].哈尔滨:哈尔滨工程大学,2013.
( Liu Yue.Research on carrier phase difference method to support on-the-fly positioning [D].Harbin:Harbin Engineering University,2013. ) ( 0)
|
| [8] |
Ma L,Yang X,Li Z.A method for determination of clock on board satellite combining L-band phase smoothed code and C-band TWTRS measurements[C]// IEEE International Frequency Control Symposium.Baltimore:IEEE,2012:1-4.
(0)
|
| [9] |
Konno H.Design of an aircraft landing system using dual-frequency GNSS [D].Palo Alto:Stanford University,2007.
(0)
|
| [10] |
甘兴利.GPS局域增强系统的完善性监测技术研究[D].哈尔滨:哈尔滨工程大学，2008.
( Gan Xing-li.Research on the integrity monitoring of GPS local area augmentation system [D].Harbin:Harbin Engineering University,2008. ) ( 0)
|