-
近年来,随着北斗卫星导航系统(BeiDou satellite navigation system, BDS)的迅猛发展,北斗系统已具备在亚太大部分地区提供卫星导航服务的能力,结合现有的GPS等全球卫星导航定位系统,GNSS组合导航定位已成为趋势。多模、多频信号的使用将会提高RTK(real-time kinematic)用户的精度和可靠性,但其关键技术是快速准确固定载波信号的模糊度,然而过去的研究表明,随着模糊度维数的陡增,模糊度正确固定的概率会降低[1]。对于多模、多频组合定位,由于卫星数较多,以及观测噪声、大气残余误差等影响,很难在短时间内成功固定所有卫星的模糊度。基于此,国内外学者提出一系列关于部分模糊度固定的方法[1, 4-8]来提高模糊度成功固定的概率,这些方法归纳起来有以下3种:① 根据信噪比[7]、高度角[8]、最优解与次最优解是否相等[4]等标准选取噪声较小的卫星频段的模糊度参与固定,去除噪声较大、残余误差较大的卫星;② 将各频段模糊度浮点解及方差阵进行组合形成宽窄巷,先固定宽巷模糊度,在宽巷模糊度固定后,通过提高窄巷模糊度方差阵的精度进而促进窄巷模糊度的固定[3];③ 直接选取模糊度的子集进行固定,在降相关后根据成功率[1, 5, 7]等准则在投影空间选取模糊度的最优子集,提高模糊度固定的概率。这些方法各具特点,在不同的环境下表现出的效果也不尽相同,例如在双频环境下,模糊度难以固定时,可以采取先固定宽巷的方式加快坐标收敛的速度[3],在卫星数非常多时直接采取高截止高度角的方式得到较好的定位结果[14]。已有文献中大多都是针对GPS单系统或GPS和Galileo模拟数据组合[4, 5, 9],本文研究在BDS/GPS双系统组合模式下的模糊度固定效果以及坐标的收敛速度情况,为当前的GNSS组合导航定位用户提供参考。
-
模糊度是载波信号首次跟踪时未知周数的整数部分,模糊度参数的求解以简化的方程表示为:
$$ \mathit{\boldsymbol{y}} = \mathit{\boldsymbol{Aa}} + \mathit{\boldsymbol{Bb}} + \mathit{\boldsymbol{\varepsilon }} $$ (1) 式中,a是双差整周模糊度参数向量;b是实参数向量包括基线分量及可能的电离层、对流层参数等;A、B是与之对应的系数矩阵;y是观测值向量,包括伪距和载波相位观测值; ε是随机噪声向量。根据式(1) 可以解出模糊度参数和其他参数的浮点解X以及方差协方差矩阵Q:
$$ \mathit{\boldsymbol{X}} = \left[ {\begin{array}{*{20}{c}} {\mathit{\boldsymbol{\hat a}}}\\ {\mathit{\boldsymbol{\hat b}}} \end{array}} \right];\mathit{\boldsymbol{Q}} = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{Q}}_{\hat a}}}&{{\mathit{\boldsymbol{Q}}_{\hat a\hat b}}}\\ {{\mathit{\boldsymbol{Q}}_{\hat b\hat a}}}&{{\mathit{\boldsymbol{Q}}_{\hat b}}} \end{array}} \right] $$ (2) -
目前,使用最广泛的是Lambda方法[15]:
$$ \mathit{\boldsymbol{\overset{\lower0.5em\hbox{$\smash{\scriptscriptstyle\smile}$}}{z} }} = \mathop {\min }\limits_{z \in \mathbb{Z}} {\left( {\mathit{\boldsymbol{\hat z}} - \mathit{\boldsymbol{z}}} \right)^{\rm{T}}}\mathit{\boldsymbol{Q}}_{\hat z}^{ - 1}\left( {\mathit{\boldsymbol{\hat z}} - \mathit{\boldsymbol{z}}} \right) $$ (3) 式中, ${\mathit{\boldsymbol{\hat{z}}}} $和$ {\mathit{\boldsymbol{Q}}_{\hat z}} $的表达式为:
$$ \left\{ \begin{array}{l} \mathit{\boldsymbol{\hat z}} = \mathit{\boldsymbol{Z\hat a}}\\ {\mathit{\boldsymbol{Q}}_{\hat z}} = \mathit{\boldsymbol{Z}}{\mathit{\boldsymbol{Q}}_{\hat a}}{\mathit{\boldsymbol{Z}}^{\rm{T}}} \end{array} \right. $$ (4) 式中, Z为降相关矩阵; ℤ代表整数集。在搜索结束后结合Z将固定后的模糊度度$ {\mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over z} }}} $恢复到原始空间,然后通过模糊度检验指标确定是否接受搜索模糊度的最佳备选组,常用的模糊度检验方法有Ratio检验、FFR[10]检验等,这些方法主要计算次最优备选组和最优备选组残差二次型的比值:
$$ RT = \frac{{{{\left\| {\mathit{\boldsymbol{\hat a}} - \mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over a} }},} \right\|}_{{\mathit{\boldsymbol{Q}}_{\hat a}}}}}}{{{{\left\| {\mathit{\boldsymbol{\hat a}} - \mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over a} }}} \right\|}_{{\mathit{\boldsymbol{Q}}_{\hat a}}}}}} \ge t $$ (5) 式中,$ {\mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over a} }}} $是次最优整周模糊度备选组;t是阈值,一般情况设t=2或3,对于FFRT,t根据先验的失败率计算[10],是动态变化的。
若最优模糊度备选组通过检验,可以用模糊度的固定解根据条件平差得到基线的固定解:
$$ \mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over b} = \hat b}} - {\mathit{\boldsymbol{Q}}_{\hat b\hat a}}\mathit{\boldsymbol{Q}}_{\mathit{\boldsymbol{\hat a}}}^{ - 1}\left( {\mathit{\boldsymbol{\hat a}} - \mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over a} }}} \right) $$ (6) -
当所有卫星的模糊度用Lambda方法难以固定时可以考虑用部分模糊度固定的方法。先固定那些较为容易固定的模糊度,假设将模糊度向量$ {\hat a} $分解为两部分,它和对应的协方差矩阵Q为:
$$ \mathit{\boldsymbol{\hat a}} = \left[ {\begin{array}{*{20}{c}} {{{\hat a}_1}}\\ {{{\hat a}_2}} \end{array}} \right];\mathit{\boldsymbol{Q}} = \left[ {\begin{array}{*{20}{c}} {{\mathit{\boldsymbol{Q}}_{{{\hat a}_1}}}}&{{\mathit{\boldsymbol{Q}}_{{{\hat a}_1}{{\hat a}_2}}}}\\ {{\mathit{\boldsymbol{Q}}_{{{\mathit{\boldsymbol{\hat a}}}_2}{{\mathit{\boldsymbol{\hat a}}}_1}}}}&{{\mathit{\boldsymbol{Q}}_{{{\hat a}_2}}}} \end{array}} \right] $$ (7) 式中,$ {{\mathit{\boldsymbol{\hat a}}}_1} $代表在原始空间按某些准则选定的容易固定部分卫星频段或者双频组合形成的宽巷模糊度;$ {{\mathit{\boldsymbol{\hat a}}}_2} $代表噪声较大的卫星或者窄巷模糊度部分。若采用Lambda方法固定$ {{\mathit{\boldsymbol{\hat a}}}_1} $后,类似于坐标分量的修正过程,用${{\mathit{\boldsymbol{\hat a}}}_1} $的固定解修正$ {{\mathit{\boldsymbol{\hat a}}}_2} $及其方差协方差阵:
$$ \left\{ \begin{array}{l} {{\mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over a} }}}_2} = {{\mathit{\boldsymbol{\hat a}}}_2} - {\mathit{\boldsymbol{Q}}_{{{\mathit{\boldsymbol{\hat a}}}_2}{{\mathit{\boldsymbol{\hat a}}}_1}}}\mathit{\boldsymbol{Q}}_{{{\mathit{\boldsymbol{\hat a}}}_1}}^1\left( {\mathit{\boldsymbol{\hat a}} - \mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over a} }}} \right)\\ {\mathit{\boldsymbol{Q}}_{{{\mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over a} }}}_2}}} = {\mathit{\boldsymbol{Q}}_{{{\hat a}_2}}} - {\mathit{\boldsymbol{Q}}_{{{\mathit{\boldsymbol{\hat a}}}_2}{{\mathit{\boldsymbol{\hat a}}}_1}}}\mathit{\boldsymbol{Q}}_{{{\mathit{\boldsymbol{\hat a}}}_1}}^1{\mathit{\boldsymbol{Q}}_{\mathit{\boldsymbol{\hat a\hat b}}}} \end{array} \right. $$ (8) 然后再用Lambda方法固定$ {{\mathit{\boldsymbol{\hat a}}}_2} $部分模糊度,若只有${{\mathit{\boldsymbol{\hat a}}}_1} $固定则可以用式(6) 更新坐标参数得到固定解,若$ {{\mathit{\boldsymbol{\hat a}}}_2}$固定则可用其再次修正坐标参数解。
对于引言中的方法3,首先用bootstrapping方法计算成功率,根据设定的成功率阈值P0挑选降相关后模糊度的最优子集[13]:
$$ {P_{s,{B_i}}} = \prod\limits_{j = i}^n {\left( {2\mathit{\Phi }\left( {\frac{1}{{2{\sigma _{{{\boldsymbol{\hat z}}_{j|J}}}}}}} \right) - 1} \right) > {P_0}} $$ (9) 式中,Φ(x)是正态分布函数;$ {\sigma _{{{\mathit{\boldsymbol{\hat z}}}_{j|J}}}} $代表条件方差,其计算方法如下,式(4) 中降相关后的方差-协方差矩阵$ {\mathit{\boldsymbol{Q}}_{\hat z}} $经过三角分解得:
$$ {\mathit{\boldsymbol{Q}}_{\mathit{\boldsymbol{\hat z}}}} = {\mathit{\boldsymbol{L}}^{\rm{T}}}\mathit{\boldsymbol{DL}} $$ (10) 其中L和D满足:
$$ \left\{ \begin{array}{l} {l_{j,i}} = {\sigma _{{{\hat z}_i}{{\hat z}_{j|J}}}}\sigma _{{{\hat z}_{j|J}}}^{ - 2}\\ {d_j} = {\sigma _{{{\mathit{\boldsymbol{\hat z}}}_{j|J}}}} \end{array} \right. $$ (11) 式中,${\sigma _{{{\mathit{\boldsymbol{\hat z}}}_{j|J}}}} $即为$ {{\mathit{\boldsymbol{\hat z}}}_j} $的条件方差,其条件是$ {{\hat z}_J} $得到固定解,其中J={j+1, j+2, …, n};$ {\sigma _{{{\hat z}_{n|n}}}} = {\sigma _{{{\hat z}_n}}} $是精度最高的浮点解模糊度。D是对角阵,矩阵中的元素是按顺序排列的,满足$ {\sigma _{{{\hat z}_{j|J}}}} \ge {\sigma _{{{\hat z}_{k|J}}}} $,i < k。
根据式(9) 即可从降相关后的$ {\mathit{\boldsymbol{\hat z}}} $向量中挑选出n-i+1个元素的模糊度子集,因为Lambda算法中,经过Z变换后${{\mathit{\boldsymbol{\hat z}}}_j} $的排列顺序与条件方差$ {\sigma _{{{\hat z}_{j|J}}}} $是对应的,所以选取的子集为$ {{\mathit{\boldsymbol{\hat z}}}_P} = \{ {{\hat z}_i}, {{\hat z}_{i + 1}}, \cdots {{\hat z}_n}\} $,其对应的Z变换矩阵为ZP, ZP即为Z矩阵的第i到n列元素。当$ {{\mathit{\boldsymbol{\hat z}}}_P} $固定后,可以根据式(8) 更新另一部分模糊度的浮点解,然后根据Z矩阵将其恢复到原始空间,由于$ {\mathit{\boldsymbol{\hat z}}} $中部分是浮点数,恢复后所有的模糊度仍是浮点数,也可以不用恢复,直接根据式(12) 得到基线的“固定解”(原始模糊度并未固定为整数):
$$ {{\mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over b} }}}_P} = \mathit{\boldsymbol{\hat b}} - {\mathit{\boldsymbol{Q}}_{\mathit{\boldsymbol{\hat b}}{{\mathit{\boldsymbol{\hat z}}}_P}}}\mathit{\boldsymbol{Q}}_{{{\mathit{\boldsymbol{\hat z}}}_P}}^{ - 1}\left( {{{\mathit{\boldsymbol{\hat z}}}_P} - {{\mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over z} }}}_P}} \right) $$ (12) 其中,$ {\mathit{\boldsymbol{Q}}_{\hat b{{\hat z}_P}}} = {\mathit{\boldsymbol{Q}}_{\hat b\hat a}}{\mathit{\boldsymbol{Z}}_P}, {\mathit{\boldsymbol{Q}}_{{{\hat z}_P}}} = \mathit{\boldsymbol{Z}}_P^{\rm{T}}{Q_{\hat a}}{\mathit{\boldsymbol{Z}}_P} $。
-
对于BDS/GPS双系统组合定位时需采用统一的时间和坐标基准。WGS-84与CGCS2000的系统差异理论上在0~0.105 mm范围内,对短距离的相对定位来说完全可以忽略[2]。对于时间基准同步,BDT和GPST除了相差1 356周14 s外还存在微小的同步误差,为了消除不同卫星系统对定位的影响,本文在双差组合时,在各自系统内部分别选取参考卫星。实时定位时采用伪距和载波观测值,双差观测方程为:
$$ \left\{ \begin{array}{l} \Delta \nabla {P^S} = \Delta \nabla {R^S} + \mathit{\boldsymbol{A}}{\rm{d}}X + \Delta \nabla {T^S} + \\ \;\;\;\Delta \nabla {I^S} + \Delta \nabla e_P^S\\ {\lambda ^S}\Delta \nabla {\varphi ^S} = \Delta \nabla {R^S}\mathit{\boldsymbol{ + A}}{\rm{d}}X + {\lambda ^S} \cdot \Delta \nabla {N^S} + \\ \;\;\;\Delta \nabla {T^S} - \Delta \nabla {I^S} + \Delta \nabla e_\varphi^S \end{array} \right. $$ (13) 式中,S代表系统类型;λ为波长;P为伪距观测值;φ为载波观测值;R为卫星到接收机相位中心的近似几何距离;N为模糊度;T为对流层延迟误差;I为电离层延迟误差;e为未模型化的误差和观测噪声;dX表示相对坐标改正向量,是双系统共有的未知参数;A为坐标改正向量的系数矩阵。若基线较短,可不考虑电离层误差和对流层误差。
-
为了分析部分模糊度对于BDS/GPS组合实时定位的效果,本文从总结的3种类型选取了相应典型的部分模糊度固定方法:① 从Lambda算法中得到的最优与次最优模糊度备选组挑选相等的模糊度重新用Lambda方法固定直到其通过Ratio检验,若挑选的子集中模糊度个数太少则结果采用浮点解(FS方法);② 通过频率间的组合,先选取波长较长的宽巷模糊度参与固定,宽巷固定后再固定窄巷(WN方法);③ 基于成功率准则在投影空间选取成功率较高的模糊度子集(SRC方法)。本文主要从模糊度子集的成功率、固定率(通过模糊度检验的历元所占的比例)、基线解算精度、收敛时间等几方面进行评价部分模糊度固定的效果。实验数据采用上海司南卫星导航技术有限公司官网上下载的两条基线数据,分别是位于广州市天河区的长约8 km的多系统观测数据(基线一)和位于上海闵行区的长约17 km的多系统观测数据(基线二),两条基线的采样率都为10 s,本文采用其中的BDS/GPS双系统L1L2B1B2频段数据用于实验分析。
-
成功率是模糊度是否能正确固定的重要参考指标之一[1],同时也反应了浮点解的精度。现有的成功率计算方法很多,其中Lambda最小二乘准则搜索模糊度的成功率是最高的,但是其难以直接计算,而bootstrapping的成功率与最小二乘成功率是相当接近的[16],因此一般用其作为最小二乘成功率的下界。模糊度固定的另一个重要指标是Ratio值,在Lambda算法中一般以其作为模糊度是否固定的最终指标。本文首先处理了基线一的BDS/GPS组合数据,处理方式为动态单历元解算。图 1给出了全部模糊度参与固定(FULL方法)及上述3种部分模糊度固定方法的成功率及Ratio值,其中蓝色线条表示Ratio值,青色代表模糊度子集的成功率,红色虚线是指参考的Ratio阈值,大小为2。图 2是3种部分模糊度固定方法选取子集的模糊度以及全部模糊的个数,FULL代表所有模糊度的个数。
本文中SRC方法设置的成功率阈值P0=0.99,FS方法中挑选最优子集的截止条件是Ratio值大于给定的阈值2且模糊度个数不少于10(若小于10则该历元是浮点解),WN方法由于宽巷固定后再固定窄巷,因此前半段的成功率是宽相模糊度的成功率(窄巷未固定),后半段是宽巷固定后窄巷模糊度参与固定时的成功率,Ratio值与之对应。从图 1可以看出WN方法的成功率也有了一定的提高,但不明显,主要原因是该时段一些观测值误差较大,通过增加波长的方法并不能排除误差较大的卫星频段的模糊度,FS和SRC方法部分模糊度选取子集的成功率都显著提高了,成功率接近于1,结合图 2可知,这两种方法的模糊度个数都有所减少,自动排除了噪声较大的模糊度,而FS方法的模糊度个数减少相对较多,接近一半,原因可能是利用该准则选取最优子集时可能过于苛刻,误删了一些精度较高的模糊度。从Ratio值方面来看,WN方法在宽巷模糊度固定后窄巷模糊度的Ratio值有了显著的提高,FS和SRC方法在排除精度较差的模糊度后Ratio值也有了显著提高,固定率大大提升,将Ratio阈值设为2时,FS方法的固定率接近100%。
由图 1可以看出,SRC方法虽然成功率很高,当Ratio阈值为2时,仍有很多历元不能通过模糊度度检验,当阈值设为3时,固定率会更低,因此Ratio阈值的设置比较关键。特别对于多系统而言,模糊度个数较多,Ratio值一般较低[10],用某一确定的阈值模糊度固定率较低,Teunissen等提出用基于失败率的FFRT方法得到Ratio阈值(动态阈值),该阈值与成功率和模糊度的个数有关[10],是动态变化的,有更好的适应性。表 1是基线一在Ratio阈值设置为2、3以及失败率为0.01的FFRT阈值时得到的4种模糊度固定方法的固定率情况。
表 1 四种方法在不同的Ratio阈值下的固定率/%
Table 1. The Proportion of Successfully Fixed Epochs Under Different Ratio Threshold/%
Ratio阈值 FULL FS WN SRC 3.0 32.9 63.8 61.8 41.4 2.0 50.4 99.8 69.0 67.2 动态阈值 69.1 99.9 69.2 100 由表 1结合图 1可以看出,Ratio阈值对固定率的影响很大,当Ratio阈值从2变为3时, FS和SRC方法固定率都下降了20%以上,同时也可以看出当Ratio阈值为定值时FS方法的固定率最高,采用FFRT方法设置阈值时各种方法固定率有了明显的提高,SRC方法固定率提高到100%。由此可以看出BDS/GPS组合时,Ratio值一般较小,采用动态的Ratio阈值更适合多系统组合时模糊度个数较多的情况; 利用FS方法时,由于其选取子集时以Ratio值作为选取条件之一,该种方法的固定率受Ratio阈值的影响相对较小。
-
通过§3.1的分析可知, BDS/GPS组合时,不同的Ratio阈值设置方法对模糊度固定率影响较大,因此, 只用固定率评判部分模糊度固定的效果是不充分的。事实上,模糊度通过Ratio检验也可能出现固定错误的情况。如果选取的模糊度子集(子集一)固定错误,较难固定的模糊度子集(子集二)根据式(8) 更新相应的模糊度浮点解及方差协方差阵必然也会受到影响,因此,子集二在采用Lambda方法固定时,一般不会通过Ratio检验,即使通过,子集二错误固定的可能性也较大。为了进一步检验其效果,需要分析模糊度固定后基线解算的精度。本文用公式$ \left\| {\mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over b}-b}}} \right\| = c \times \left\| {\mathit{\boldsymbol{\hat b-b}}} \right\| $评定部分模糊度基线解算的精度,其中, b代表真实坐标向量; ${\mathit{\boldsymbol{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\smile$}} \over b} }}} $代表部分模糊度方法解算的坐标向量; $ {\mathit{\boldsymbol{\hat b}}} $代表全部模糊度参与固定的解算结果; c是部分模糊度方法的三维坐标误差与FULL方法坐标误差的比值。图 3是基线一单历元动态解算时部分模糊度基线解算效果情况,其中图 3(a)是全部模糊度参与固定(FULL)时基线误差情况,图 3(b)是相应的部分模糊度方法解算的误差的比值。当c=1时代表该历元的解算结果都为浮点解或者都固定且结果一致,当c < 1时表示部分模糊度方法解算的结果精度有所提高。从图 3中可以看出, FS和SRC方法对于FULL方法未固定的模糊度其比值大多都在0.2以下,基线解算精度有明显的提升,WN方法在某些历元精度也有所提高,由此表明采用部分模糊度固定的方法不仅可以提高固定率,且在模糊度固定后,基线解算精度有大幅的提升。
图 3 3种部分模糊度与FULL方法基线解算精度比较
Figure 3. Contrast of Three Partial Ambiguities Resolution Methods with "FULL" Method about the Precision of Baseline Solution
图 4是采用滤波的方式动态解算基线二时的NEU坐标误差随着时间的变化结果,图中蓝色点代表固定解,红色点代表浮点解。图 5是4种不同的方法所对应的模糊度固定的个数,表 2是该条基线用4种方法解算时首次固定时所需要的时间以及每种方法解算坐标误差NEU方向的STD值。由图 4和表 2可知, 3种部分模糊度固定方法都明显缩短了初始化时间,FS方法只需用20 s即可得到固定解,WN方法虽有很多历元是浮点解,然而其坐标误差较小(0.012, 0.016, 0.042),结合图 5可知其原因是红点部分对应的窄巷模糊度虽没有固定,但从第3个历元开始宽相模糊度已固定,可以得到较为准确的坐标解;图中FULL方法出现固定后又浮点的情况,主要原因是视野中出现的新升起的卫星,其高度角较低,噪声较大,影响模糊度的固定,而用部分模糊度的方法可以较好地解决此问题。
图 4 4种模糊度固定方法解算的坐标误差结果
Figure 4. Performance of Four Ambiguity > Resolution Methods for "Baseline 2"
表 2 4种解算方法的初始化时间及坐标误差结果
Table 2. The Initialization Time and Baseline Solution Error of Four Methods
模糊度固定方法 FULL FS WN SRC 初始化时间/s 910 20 500 640 N/m 0.040 0.012 0.014 0.030 E/m 0.063 0.016 0.016 0.036 U/m 0.111 0.030 0.042 0.082 -
本文针对BDS/GPS组合动态定位模糊度难以固定的问题,归纳总结了现有的部分模糊度固定方法,分析了每种方法的特点,最后用实测数据分析了3种部分模糊度固定的效果。实验结果表明, 部分模糊度固定方法对于BDS/GPS组合定位可以显著提高模糊度固定的成功率及固定率,并可以缩短RTK定位的初始化时间,得到更高精度的基线解算结果。其中宽窄巷(WN)方法通过先固定波长较长的宽巷模糊度再固定窄巷,宽巷固定后不仅可以得到较高精度的基线解,而且又能促进窄巷的固定,缩短初始化时间。另外两种方法通过一些准则挑选最优模糊度子集,可以自动排除噪声较大的卫星或浮点解精度较低的模糊度,从而提高模糊度固定的成功率,加快模糊度的固定,缩短坐标的收敛时间。在后续的研究中若能结合各种部分模糊度固定方法的优势,则可能会得到更好的结果。
Performance Assessment of Partial Ambiguity Resolution Based on BDS/GPS Combined Positioning
-
摘要: 多系统组合有利于提高卫星导航定位的精度及可靠性,然而对于载波差分定位由于模糊度维数的陡增、观测噪声、大气残余误差等原因用传统的Lambda方法很难得到所有模糊度的固定解,采用部分模糊度方法固定最优的模糊度子集则相对容易。总结了现有的部分模糊度固定方法,分析了不同方法的特点,并用实测数据分析了BDS/GPS组合动态定位时部分模糊度固定的效果。实验结果表明,部分模糊度方法可以显著提高模糊度固定时的成功率及Ratio值,并且可以缩短RTK定位时的初始化时间,加快坐标的收敛速度,提高组合系统动态定位结果的精度。Abstract: Multi-GNSS combination is helpful to improve the accuracy and relibality of satellite navigation and positioning. However, in case of BDS/GPS combined carrier phase positioning systems, it is difficult to fix all the ambiguities due to the increasing number of ambiguities, high measurement noises or residual atmosphere delays with the traditional Lambda method. But it is of greater probability to fix a subset of ambiguities. In this paper, we divided current partial ambiguity fixing methods into three categories and analyzed the characteristics of every method. Finally, the effect of the three partial resolution methods were tested with the measured BDS/GPS data. The results show that the success rate and ratio value is obviously improved when using partial ambiguity fixing, at the same time, the initialization time of RTK is shortened, the precision of kinematic positioning is also improved.
-
Key words:
- partial ambiguity resolution /
- BDS/GPS /
- success rate /
- RTK /
- Lambda
-
表 1 四种方法在不同的Ratio阈值下的固定率/%
Table 1. The Proportion of Successfully Fixed Epochs Under Different Ratio Threshold/%
Ratio阈值 FULL FS WN SRC 3.0 32.9 63.8 61.8 41.4 2.0 50.4 99.8 69.0 67.2 动态阈值 69.1 99.9 69.2 100 表 2 4种解算方法的初始化时间及坐标误差结果
Table 2. The Initialization Time and Baseline Solution Error of Four Methods
模糊度固定方法 FULL FS WN SRC 初始化时间/s 910 20 500 640 N/m 0.040 0.012 0.014 0.030 E/m 0.063 0.016 0.016 0.036 U/m 0.111 0.030 0.042 0.082 -
[1] Teunissen P J G, et al. Geometry-free Ambiguity Success Rates in Case of Partial Fixing[C]. Proceedings of ION-NTM, San Diego, CA, 1999 http://www.researchgate.net/publication/266016106_Geometry-free_Ambiguity_Success_Rates_in_Case_of_Partial_Fixing [2] 高星伟, 过静君, 程鹏飞, 等.基于时空系统统一的北斗与GPS融合定位[J].测绘学报, 2012, 41(5):743-748 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=chxb201205021&dbname=CJFD&dbcode=CJFQ Gao Xingwei, Guo Jingjun, Cheng Pengfei, et al. Fusion Positioning of BeiDou/GPS Based on Spatio Temporal System Unification[J].Acta Geodaetica et Cartographica Sinica, 2012, 41(5):743-748 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=chxb201205021&dbname=CJFD&dbcode=CJFQ [3] Li Jinlong, Yang Yuanxi, Xu Junyi, et al. GNSS Multi-carrier Fast Partial Ambiguity Resolution Strategy Tested with Real BDS/GPS Dual-and Triple-frequency Observations[J]. GPS Solutions, 2015, 19(1):5-13 doi: 10.1007/s10291-013-0360-6 [4] Dai Liwen, Eslinger D, Sharpe T. Innovative Algorithms to Improve Long Range RTK Reliability and Availability[C]. Proceedings of the 2007 National Technical Meeting of the Institute of Navigation, San Diego, CA, 2007 http://www.researchgate.net/publication/242267642_Innovative_Algorithms_to_Improve_Long_Range_RTK_Reliability_and_Availability [5] Wang Jun, Feng Yanming. Reliability of Partial Ambiguity Fixing with Multiple GNSS Constellations[J]. Journal of Geodesy, 2013, 87(1):1-14 doi: 10.1007/s00190-012-0573-4 [6] Li Pan, Zhang Xiaohong. Precise Point Positioning with Partiatial Ambiguity Fixing[J]. Sensors, 2015, 15(6):13627-13643 doi: 10.3390/s150613627 [7] Parkins A. Increasing GNSS RTK Availability with a New Single-epoch Batch Partial Ambiguity Resolution Algorithm[J]. GPS Solutions, 2011, 15(4):391-402 doi: 10.1007/s10291-010-0198-0 [8] Takasu T, Yasuda A. Kalman-filter-based Integer Ambiguity Resolution Strategy for Long-baseline RTK with Ionosphere and Troposphere Estimation[C]. Proceedings of the 23rd International Technical Meeting of the Satellite Division of the Institute of Navigation, Portland, OR, 2010 https://www.researchgate.net/publication/260387204_Kalman-Filter-Based_Integer_Ambiguity_Resolution_Strategy_for_Long-Baseline_RTK_with_Ionosphere_and_Troposphere_Estimation [9] Cao Wei, O'Keefe K, Cannon M. Partial Ambiguity Fixing within Multiple Frequencies and Systems[C]. Proceedings of ION GNSS, Fort Worth, TX, 2007 http://www.researchgate.net/publication/250305744_Partial_Ambiguity_Fixing_within_Multiple_Frequencies_and_Systems [10] Verhagen S, Teunissen P J G. The Ratio Test for Future GNSS Ambiguity Resolution[J]. GPS Solutions, 2013, 17(4):535-548 doi: 10.1007/s10291-012-0299-z [11] Verhagen S, Teunissen P J G, van der Marel H, et al. GNSS Ambiguity Resolution:Which Subset to Fix[C]. Proceedings of IGNSS Symposium, Sydney, Australia, 2011 http://www.researchgate.net/publication/236213381_GNSS_ambiguity_resolution_which_subset_to_fix [12] Teunissen P J G. GNSS Ambiguity Bootstrapping:Theory and Application[C]. Proceedings of International Symposium on Kinematic Systems in Geodesy, Geomatics and Navigation, Banff, Canada, 2001 http://www.researchgate.net/publication/228870630_GNSS_ambiguity_bootstrapping_theory_and_application [13] Teunissen P J G. Success Probability of Integer GPS Ambiguity Rounding and Bootstrapping[J]. Journal of Geodesy, 1998, 72(10):606-612 doi: 10.1007/s001900050199 [14] Teunissen P J G. Instantaneous BeiDou+ GPS RTK Positioning with High Cut-off Elevation Angles[J]. Journal of Geodesy, 2014, 88(4):335-350 doi: 10.1007/s00190-013-0686-4 [15] Teunissen P J G. The Least-squares Ambiguity Decorrelation Adjustment:A Method for Fast GPS Integer Ambiguity Estimation[J]. Journal of Geodesy, 1995, 70(1/2):65-82 https://www.researchgate.net/publication/224969472_The_least-squares_ambiguity_decorrelation_adjustment_A_method_for_fast_GPS_integer_ambiguity_estimation [16] Verhagen S. On the Approximation of the Integerleast-squares Success Rate:Which Lower or Upperbound to Use[J]. Journal of Global Positioning Systems, 2003, 2(2):117-124 doi: 10.5081/jgps -