-
整周模糊度的快速解算是动态GNSS及GNSS/INS精密定位的重点和难点问题。目前常用的动态模糊度解算方法包括模糊度函数法[1, 2]、最小二乘搜索法[3, 4]、快速模糊求解法(fast ambiguity resolution approach,FARA)[5]、快速模糊值搜索滤波器 (fast ambiguity search filter,FASF)[6]和LAMBDA[7]等。这些方法一般需要利用多个历元的观测数据,若发生周跳,需重新求解模糊度,否则会导致解算失败[8]。基于此,许多学者提出单历元解算方法,例如,双频P码法[9]需要使用P码,普通用户难以获得高精度的伪距观测值;而单频单历元模糊度解算一般借助约束条件来实现[10, 11];文献[12]在变形监测背景下,提出了模糊度直接解法 (direct calculation,DC),要求监测点位移小于0.16 m;文献[13]使用多种组合观测量,监测点误差可以扩大到0.5 m,但(-3,4) 和(4,-5) 组合对电离层延迟误差的放大超过两倍,降低了方法的可靠性。与变形监测不同,在动态定位中,由于无法获得较为可信的坐标初值,难以使用直接法,对此,文献[14]先基于双频相关法(dual frequency correlation method,DUFCOM)固定部分一级卫星的模糊度并得到坐标参数局部解,再基于局部解对其它卫星使用直接法固定模糊度。由于仍需要对一级卫星进行搜索,因此该方法不完全属于直接法。
本文选择4颗高度角较高的卫星作为基本卫星,使用INS位置作为坐标初值,估计这些卫星各个频点整周模糊度参数值。对可能的GNSS失锁造成的INS位置误差积累现象,构建波长更长的(1,-1) 和(-3,4) 组合模糊度。针对(-3,4) 组合放大电离层误差的问题,对其进行多项式拟合,得到电离层延迟误差估值,对组合模糊度进行电离层补偿,削弱电离层误差等系统性影响后,直接取整固定组合模糊度。利用组合模糊度确定各频点模糊度固定解后,回代求解接收机坐标,基于此坐标直接固定其他卫星各频点的模糊度,避免低高度角卫星因数据连续性不足导致电离层误差难以拟合的问题。最后采用实测GNSS/INS数据对算法进行了验证。
-
在精密定位中,模糊度的解算是关键问题。GNSS/INS组合下,INS位置提供了额外的几何信息,使得动态情况下估计模糊度成为可能。
设组合系统接收机为j,则对应于卫星q的载波相位观测方程为:
(1) 式中,λ为波长;Φjq为以周为单位的载波相位值;ρjq为几何距离;Njq为整周模糊度;δt为钟差;eΦ,jq为测量噪声,式中略去了大气折射延迟dI和dT。
设参考星为p,则可得星间单差观测方程为:
(2) 若有基准站i,且在此站同步观测到卫星p、q,则可得站星双差观测方程为:
(3) 将此时INS计算的位置参数XI代入到式(3) 中的∇ρjp,q中,略去上下标符号,得到:
(5) 由式(4) 可对双差模糊度按式(6) 进行估计:
(6) 估值∇Δ${\hat{N}}$受INS位置误差、残余大气折射延迟误差、多路径效应和观测噪声影响,在随后的论述中使用表示双差模糊度∇Δ${\hat{N}}$。
若要采用直接取整法固定模糊度,则${\hat{N}}$必须有足够的精度,进而要求位置参数具有较高的精度。文献[12]从中误差的角度推导了变形监测点坐标初值的误差限值,且未考虑载波相位观测噪声影响,所以并不严密。这里,直接从真误差的角度进行分析。
${\hat{N}}$要取整确定,其误差d${\hat{N}}$应满足:
(7) 基准站位置固定,误差视为零,则要求:
(8) 暂不考虑残余大气折射延迟和多路径效应,设载波相位中误差σΦ=0.01,则有:
(9) 取99.97%置信区间的情况下,∇ΔΦ误差的取值范围为:
(10) 在式(8) 中顾及这部分噪声影响,则有:
(11) 对于频点L1和L2,分别有:
(13) 即INS导出的等效距离误差必须小于0.084 m,才能通过取整正确固定模糊度。易知:
(14) 式中,l为方向余弦;dXI=[dxIdyIdzI]T为INS位置误差。(lq-lp)的模为:
(15) 式中,β为lkp和lkq的夹角,即接收机到卫星p和参考星q的夹角,夹角越小,误差放大因子越小;夹角越大,误差放大因子越大。若β∈[90,270],误差会放大,否则误差缩小。
若假设|lq-lp|=1且INS各方向误差相同,即dxI=dyI=dzI,则有:
(16) (17) 即INS三维坐标误差应小于0.145 m。
在GNSS/INS组合精密定位系统中,INS的位置参数在GNSS更新点上被GNSS载波相位观测值更新,对于k时刻的INS位置XI,设其之前某时刻m有固定解,则其误差可表示为:
(18) 式中,σXI(m)为固定模糊度后的载波相位观测值与INS的组合解,可称GNSS/INS固定解,精度高;而σINS(m,k)为m到k时刻INS单独导航的累积误差。
若卫星几何结构良好,各历元模糊度正确固定,则m=k-1。m与k之间的时间间隔越大,INS误差积累越大。若长时间失锁或观测条件不佳导致模糊度一直未固定,则导致INS误差无法得到修正,此时对模糊度估值直接取整会造成错误。且在前面论述中,没有考虑大气折射延迟的影响,如果双差残余大气误差较大,会进一步降低模糊度估值精度。
图 1显示的是某一组GPS/INS观测数据(PRN6) 的 L1频点模糊度N1的直接估计误差。大部分时间观测条件较好,可直接取整固定单频模糊度,而在267 579 s~267 659 s之间,几何机构不佳,INS误差积累,模糊度估计误差较大,无法通过取整固定。
某型INS (标称陀螺漂移1 deg/h)在模拟的3次失锁(每次均为20 s)情况下的误差积累结果如表 1所示,其中第2次失锁的X轴误差情况如图 2所示。可以看出,对于一般的战术级INS,当失锁导致INS误差积累到一定程度时,难以取整固定单频模糊度。
表 1 模拟失锁(20 s)时INS误差结果/m
Table 1. INS Errors Under the Simulated Losslock (20 s)/m
失锁序号 1 2 3 X 0.219 0.368 -0.102 Y 0.172 0.103 -0.170 Z -0.025 0.114 0.206 -
单一频点模糊度估值对INS误差的要求较高,限制了单历元直接取整固定方法的应用范围。由式(12) 、式(13) 易知,如果能增加载波信号的波长,便能降低直接解法对INS精度的要求。
能够增加波长的双频宽巷组合主要有(1,-1) 、(-3,4) 、(4,-5) 、(-7,9) 等,对应的模糊度组合分别用NWL11、NWL34、NWL45、NWL79表示。
L1频点INS辅助模糊度估计误差可表示为:
(19) (1,-1) 组合模糊度估计误差为:
(20) 因此,(1,-1) 组合大大削弱了INS误差影响,结合式(17) 、式(19) 、式(20) 可得INS点位误差需满足:
(21) 由于(1,-1) 组合的噪声大于L1,实际的点位误差限值应比0.657更小些,对于一般的战术级惯性测量单元 (inertial measurement unit,IMU),在GNSS完全失锁10~15 s以内均能使误差低于该水平。
上述推导过程未考虑残余电离层延迟误差影响。设L1频点双差残余电离层延迟误差为∇ΔdI,1=K/f12,以周为单位表示为:
(22) 则dNWL11中电离层误差为:
(23) 可以看出,使用(1,-1) 组合,虽然观测噪声略微放大,但是电离层误差影响减小,可容许的INS误差限值增大。
对于(-3,4) 模糊度组合,有:
(24) 其电离层延迟误差为:
(25) 相对单频 L1模糊度,(-3,4) 组合模糊度能将INS误差影响削弱到 L1频点的0.12倍,但是电离层误差影响是 L1的2倍以上。图 3所示为PRN6的(-3,4) 组合模糊度结果,实线为真值,可以看出,受电离层误差的影响,直接对(-3,4) 组合模糊度取整会产生错误的结果。
(4,-5) 和(-7,9) 对INS误差的削弱作用更强,而电离层误差和噪声的放大倍数也更高。使用后三种组合,均需对电离层误差进行补偿,下面以(-3,4) 组合为例说明。
这里采用多项式拟合的方式对电离层延迟误差进行拟合并外推到当前历元。设在k-n到k-m历元,有(-3,4) 组合固定解序列[NWL34,k-n,…,NWL34,k-m],结合估值序列[${\hat{N}}$WL34,k-n,…,${\hat{N}}$WL34,k-m],则可得误差序列:
(26) 由于组合模糊度受电离层影响远大于噪声,使用误差序列进行多项式拟合后,拟合结果基本上能反映电离层误差的趋势,将拟合结果外推到当前历元k,得到当前误差估值:
(27) 即可对模糊度估值${\hat{N}}$WL34,k进行补偿:
(28) 实际计算得到的模糊度组合中的误差项${\hat{e}}$k,I不完全由电离层误差构成,还可能包括对流层残余误差和多路径效应等,只是电离层误差影响最大,这里的补偿实际上也削弱了对流层延迟和多路径效应影响。若观测质量较好,一般m=k-1,n可取15,即采用15个历元的数据进行拟合。图 4所示为PRN6电离层补偿后的(-3,4) 组合模糊度结果,拟合中采用的是最简单的一阶多项式,最初部分历元数据短,无法进行拟合,为图形显示方便,这部分历元所示的是模糊度真值。电离层拟合并补偿后,组合模糊度估值的系统误差部分得到了明显的控制,模糊度估计误差均在0.5周以内,直接取整便能确定组合模糊度值。若要得到更好的结果,还可以选择阶数更高的多项式。
-
使用模糊度组合可以削弱INS误差的影响,利用多项式拟合能抑制模糊度组合中电离层误差的影响。多项式拟合需要数据具有较好的连续性,对于低高度角的卫星,难以保证持续锁定,不利于电离层误差的拟合和补偿。4颗观测卫星能唯一确定接收机位置,因此,在所有卫星中,若能固定其中4颗卫星的模糊度,便能得到比INS更为精确的位置坐标,以该坐标可直接估计各个频点的模糊度并取整固定。INS辅助模糊度单历元直接解法可归纳为以下5个步骤。
(1) 使用INS位置坐标估计高度角最高的4颗基本卫星各频点模糊度参数值;
(2) 用单频模糊度来构建组合模糊度,综合噪声和波长特性并顾及式(21) 的需求,这里选(1,-1) 和(-3,4) 组合;
(3) 对(-3,4) 组合使用多项式拟合进行电离层补偿,获得去除系统误差项的组合模糊度估值;
(4) 组合模糊度取整,反算各频点模糊度:
N1=4NWL11+NWL34
N2=3NWL11+NWL34
(5) 利用已固定模糊度的4颗卫星估计接收机坐标,使用此坐标取代INS位置,估计其它卫星的模糊度参数并直接取整固定。
固定所有模糊度后,即可回代所有载波相位观测方程,对滤波器进行观测更新,获得高精度的导航参数。
-
使用采集于武汉的动态车载GPS/INS组合系统数据验证INS辅助单历元直接解法的有效性,车载天线和基站之间全程平均相距15 km。INS采样率200 Hz,标称陀螺漂移1 deg/h,标称加速度计偏置1e-4 g,GPS为双频接收机,采样率1 Hz。267 579 s~267 659 s历元内,可视卫星数较少,其他历元观测质量较好,最大可视卫星数达9颗。分别采用LAMBDA方法和本文直接解法解算模糊度,具体分为3个方案进行比较。
方案1 LAMBDA方法搜索模糊度,采用的是先固定宽巷再固定 L1的策略。
方案2 对全部卫星采用INS辅助直接解法确定模糊度。
方案3 本文INS辅助直接解算方案,先解算4颗基本卫星的模糊度,确定接收机坐标后回代求解其它卫星各频点模糊度。
三个方案均采用单历元瞬时固定模式,其中方案2和方案3的区别仅仅在于是否将卫星分为两部分分别求解。固定后的模糊度均回代更新INS位置参数,模糊度真值由顾及前后历元关系的后处理方法得到。观测高度角最高的4颗卫星星号为PRN3、6、19、27;作为基本卫星,其中PRN3作为组差参考星;PRN6为高度角较低的卫星(特别是在初期),这里进行重点评价。图 5和图 6所示为方案1采用LAMBDA方法搜索的 L1频点双差模糊度最优结果的误差,方案2和方案3的错误固定结果较少,不再以图形显示。表 2显示的是3种方案 L1频点模糊度解的错误率,历元总数为770。上述卫星中,PRN3观测历元769,PRN6观测历元729,其余均为770。PRN6的 (-3,4) 组合及电离层补偿结果如图 3、4所示,其它基本卫星的电离层误差更小,这里不再显示。
表 2 模糊度解算的错误率
Table 2. Error Rates of Ambiguity Resolution
方案 1 方案 2 方案 3 PRN3 14/769 0/769 0/769 PRN6 6/729 4/729 0/729 PRN19 12/770 0/770 0/770 PRN27 9/770 0/770 0/770 从以上结果可以得出以下结论。
1) 本文INS辅助模糊度单历元直接解法正确固定了所有的模糊度,即使在卫星观测条件不佳的情况下,仍能通过合理的解算方案固定模糊度;
2) 若对全部卫星采用INS辅助直接解法确定模糊度,则低高度角的卫星会因观测数据的不连续行而难以准确拟合出(-3,4) 模糊度组合的电离层误差,造成取整错误;
3) 在观测质量较优的情况下,LAMBDA方法能够准确固定模糊度;当观测质量较差时,LAMBDA方法无法搜索到最优结果。在没有其它约束条件的情况下,使用LAMBDA方法难以实现单历元模糊度固定。
-
LAMBDA等搜索算法难以实现单历元模糊度固定;使用INS发挥几何信息的作用,可直接辅助模糊度固定;单频模糊度估计结果受INS误差影响较大,使用波长较长的模糊度组合能削弱INS误差影响,提高对INS精度水平的容忍度。长波长的(-3,4) 、(4,-5) 、(-7,9) 组合在很大程度上放大了电离层延迟误差等系统误差项的值,直接取整存在较大风险,而采用多项式拟合方法可以对电离层误差进行补偿,去除模糊度组合中的系统误差部分,便于直接取整固定。先固定基本卫星的模糊度后,其它卫星每个频点的模糊度均可直接取整固定,效率较高。
INS Aided Direct Calculation Algorithm for Single-Epoch Ambiguity Resolution
-
摘要: 提出了一种无须进行搜索的单历元模糊度解法。首先,选择4颗基本卫星,使用INS位置坐标估计这些卫星各频点整周模糊度参数值;构建(1,-1)和(-3,4)组合模糊度,对(-3,4)组合模糊度进行多项式拟合,得到电离层延迟误差估值并进行误差补偿,削弱电离层误差等系统性影响后,直接取整固定组合模糊度;利用组合模糊度确定各频点模糊度后,回代求解接收机坐标,基于此坐标直接固定其他卫星各频点的模糊度。采用实测GNSS/INS数据对算法实现了有效性验证。Abstract: A single-epoch ambiguity resolution algorithm without a search process is proposed. Four basic satellites are selected, whose integer ambiguities are estimated by using the INS position, then (1,-1) and (-3,4) ambiguity combinations are formed. The ionospheric error of the (-3,4) combination is polynomially fitted. After systematic influences such as the ionospheric error are compensated, the ambiguity combination can be fixed by simple rounding. The ambiguities of each frequency can be determined by the combination and returned to resolve the receiver position. Based on the receiver position, the ambiguities of all other satellites can be calculated directly and fixed by rounding. The algorithm was validated by field test GNSS/INS data.
-
Key words:
- INS /
- integer ambiguity /
- single-epoch /
- direct calculation /
- polynomial fitting /
- basic satellites
-
表 1 模拟失锁(20 s)时INS误差结果/m
Table 1. INS Errors Under the Simulated Losslock (20 s)/m
失锁序号 1 2 3 X 0.219 0.368 -0.102 Y 0.172 0.103 -0.170 Z -0.025 0.114 0.206 表 2 模糊度解算的错误率
Table 2. Error Rates of Ambiguity Resolution
方案 1 方案 2 方案 3 PRN3 14/769 0/769 0/769 PRN6 6/729 4/729 0/729 PRN19 12/770 0/770 0/770 PRN27 9/770 0/770 0/770 -
[1] Han S, Rizos C. Improving the Computational Efficiency of the Ambiguity Function Algorithm[J]. Journal of Geodesy, 1996, 70(6):330-341 doi: 10.1007/BF00868185 [2] Corvett S J, Cross P A. GPS Single Epoch Ambiguity Resolution[J]. Survey Review, 1995, 33(2):149-160 [3] Hatch R. Ambiguity Resolution in the Fast Lane[C]. ION GPS-89, Second International Technical Meeting of the Satellite Division of the Institute of Navigation, Colorado, 1989 [4] Euler H J, Hatch R. Comparison of Several AROTF Kinematic Techniques[C]. ION GPS-94, Salt Lake City, Utah, 1994 [5] Frei E, Beutler G. Rapid Static Positioning Based on the Fast Ambiguity Resolution Approach:Theory and First Results[J]. Manuscript Geodetic, 1990, 15(6):325-356 http://cn.bing.com/academic/profile?id=2466282545&encoded=0&v=paper_preview&mkt=zh-cn [6] Chen D. Fast Ambiguity Search Filter (FASF):A Novel Concept for GPS Ambiguity Resolution[C]. ION GPS-93, Salt Lake City, Utah, 1993 [7] Teunissen P J G. The Least-Square Ambiguity Decorrelation Adjustment:A Method for Fast GPS Integer Ambiguity Estimation[J]. Journal of Geodesy, 1995, 70:65-82 doi: 10.1007/BF00863419 [8] 刘宁,熊永良,王德军,等.一种新的GPS整周模糊度单历元求解算法[J].武汉大学学报·信息科学版,2013,38(3):291-294 http://ch.whu.edu.cn/CN/abstract/abstract57.shtml Liu Ning, Xiong Yongliang, Wang Dejun, et al. A New Algorithm for Solving GPS Integer Ambiguity Using Single Epoch Data[J]. Geomatics and Information Science of Wuhan University, 2013, 38(3):291-294 http://ch.whu.edu.cn/CN/abstract/abstract57.shtml [9] Sjoberg L E. Unbiased vs Biased Estimation of GPS Phase Ambiguities Form Dual-Frequency Code and Phase Observable[J]. Journal of Geodesy, 1993, 73:118-124 http://cn.bing.com/academic/profile?id=1987928167&encoded=0&v=paper_preview&mkt=zh-cn [10] 唐卫明,孙红星,刘经南.附有基线长度约束的单频数据单历元LAMBDA方法整周模糊度确定[J].武汉大学学报·信息科学版,2005,30(5):444-446 http://ch.whu.edu.cn/CN/abstract/abstract2189.shtml Tang Weiming, Sun Hongxing, Liu Jingnan. Ambiguity Resolution of Single Epoch Single Frequency Data with Baseline Length Constraint Using LAMBDA Algorithm[J]. Geomatics and Information Science of Wuhan University, 2005, 30(5):444-446 http://ch.whu.edu.cn/CN/abstract/abstract2189.shtml [11] 李博峰, 沈云中.顾及基线先验信息的GPS模糊度快速解算[J].测绘学报,2008,37(4):423-428 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB200804005.htm Li Bofeng, Shen Yunzhong. Fast GPS Ambiguity Resolution Constraint to Available Conditions[J]. Acta Geodaetica Et Cartographica Sinica, 2008, 37(4):423-428 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB200804005.htm [12] 王新洲,花向红,邱蕾.GPS变形监测中整周模糊度解算的新方法[J].武汉大学学报·信息科学版,2007,32(1):24-26 http://ch.whu.edu.cn/CN/abstract/abstract1795.shtml Wang Xinzhou, Hua Xianghong, Qiu Lei. A New Method for Integer Ambiguity Resolution in GPS Deformation Monitoring[J]. Geomatics and Information Science of Wuhan University, 2007,32(1):24-26 http://ch.whu.edu.cn/CN/abstract/abstract1795.shtml [13] 薛志宏,楼益栋, 易文婷,等.一种适用于动态变形测量的双频模糊度实时解算方法[J].武汉大学学报·信息科学版,2012,37(6):683-687 http://ch.whu.edu.cn/CN/abstract/abstract245.shtml Xue Zhihoong, Lou Yidong, Yi Wenting, et al. An Instantaneous Ambiguity Resolution for Kinematic Deformation Surveying Using GPS[J]. Geomatics and Information Science of Wuhan University, 2012, 37(6):683-687 http://ch.whu.edu.cn/CN/abstract/abstract245.shtml [14] 郭际明,周命端,谢翔,等.利用DUFCOM和DC算法的GPS单历元双差整周模糊度快速确定算法[J].武汉大学学报·信息科学版,2013,38(7):813-817 http://ch.whu.edu.cn/CN/abstract/abstract2701.shtml Guo Jiming, Zhou Mingduan, Xie Xiang, et al. A Fast Fixed Algorithm of GPS Single Epoch Ambiguity Resolution Based on DUFCOM and DC Algorithms[J]. Geomatics and Information Science of Wuhan University, 201237(6):813-817 http://ch.whu.edu.cn/CN/abstract/abstract2701.shtml -