CN103675835B - 一种北斗三频信号载波相位整周模糊度单历元确定方法 - Google Patents

Abstract

一种北斗三频信号载波相位整周模糊度单历元确定方法，基于几何无关模型由确定北斗三频组合模型、三频超宽项组合模糊度单历元的解算、最优宽项模糊度的解算、最优窄巷模糊度的解算这些步骤配合，逐级分步确定各组合模糊度，完成三频信号模糊度单历元的解算；本发明在三频模糊度解算方面，采用了求解模糊度成功率最高的LAMBDA方法，并进行了改进，进一步提高了其单历元搜索成功率,该方法固定模糊度效率更高，更能适用于实时动态的环境中。

Description

一种北斗三频信号载波相位整周模糊度单历元确定方法

技术领域

本发明涉及一种模糊度单历元固定方法，特别是北斗三频信号载波相位整周模糊度单历元确定方法,可用于卫星导航系统定位。

背景技术

RTK系统能为终端用户提供实时厘米级定位，实时快速解算参考站的整周模糊度是保证该系统实现高精度动态核心关键问题之一。目前，RTK技术已经推广应用，随着全球卫星导航系统逐步进入多频多系统的联合定位新时期，单一的GPSRTK模式逐步被多频多模的GNSSRTK所替代。该发明实现北斗三频单历元固定整周模糊度的方法。根据不同应用场合，观测方程模型可表示为几何相关模型或几何无关模型。几何相关模型在用于求解模糊度参数时观测时间较长才能保证求救成功率较高，且单历元内为秩亏方程，增加历元虽然能获得满秩，但是对微小扰动非常敏感；几何无关模型是一种线性模型，不以基线分量为未知数，几个历元便可以解算出模糊度，特别是三频观测，甚至单历元即可成功固定模糊度。目前用于采用几何相关模型应用较为广泛，但是无法实现单历元解算。多个频率的观测模式更加有利于形成多种特性较优的线性组合，从而给模糊度快速解算带来了机遇。多频模糊度解算方法最著名的是TCAR法、CIR法，两种算法的主要思想都是根据不同最优组合观测值的波长及其误差特点，采用简单的舍入取整法逐级地依次固定超宽巷、宽巷及窄巷或者中巷模糊度，最终确定基础载波的整周模糊度。这类方法虽然降低了运算复杂度，但是需要长时间对观测进行平滑，并且可靠性不高。整周模糊度的单历元固定，不仅需要较高的组合整周模糊度单历元解算成功率，更需要提高窄巷及基础各单频的整周模糊度的单历元解算成功率。本发明基于常规三频模糊度解算方法，结合北斗系统特点，结合LAMBDA算法，采用新的逐级固定方法，对解算的每个步骤进行了相应的优化，保证单历元情况下，测站可以得到较高可靠性的整周模糊度。

发明内容

本发明的技术解决问题：克服现有技术中双频动态确定模糊度效率低、基线短和可靠性不高的不足，提出一种北斗三频信号载波相位整周模糊度单历元确定方法，采用新的逐级固定方法，对解算的每个步骤进行了相应的优化，保证单历元情况下，测站可以得到较高可靠性的整周模糊度。

本发明的技术解决方案：一种北斗三频信号载波相位整周模糊度单历元确定方法，基于几何无关模型，逐级分步确定各组合模糊度，完成三频信号载波相位整周模糊度单历元的解算，所述逐级分布步骤如下：

步骤1：依据基线长度和波长、电离层放大倍数，确定北斗三频组合载波相位双差观测模型，北斗三频组合载波相位双差观测模型写为：

$\▿\mathrm{\Δ}{\mathrm{\Φ}}_C=\mathrm{\α}\▿\mathrm{\Δ}{\mathrm{\Φ}}_1+\mathrm{\β}\▿\mathrm{\Δ}{\mathrm{\Φ}}_2+\mathrm{\γ}\▿\mathrm{\Δ}{\mathrm{\Φ}}_3---\left(1\right)$

其中，为组合载波相位双差值，为B1载波相位双差值、为B2载波相位双差值、为B3载波相位双差值；

为使几何距离不受观测值组合不同而变化，将组合观测模型表示为一般形式为：

$\▿\mathrm{\Δ}{\mathrm{\Φ}}_C=\frac{\▿\mathrm{\Δ\ρ}}{{\mathrm{\λ}}_C}+\▿\mathrm{\Δ}{N}_C-\frac{{\mathrm{\μ}}_I(\mathrm{\α},\mathrm{\β},\mathrm{\γ})\▿\mathrm{\ΔI}}{{\mathrm{\λ}}_C}+\frac{\▿\mathrm{\Δ}{\mathrm{\ϵ}}_C}{{\mathrm{\λ}}_C}---\left(2\right)$

式中，为双差算子，为卫星与测站之间的几何距离双差值（m），并包含所有与频率无关的误差残差项（如对流层残差、钟差残差等），为B1载波频率上的电离层延迟双差值，为组合观测噪声；

组合电离层比例因子μ_I(α,β,γ)：

${\mathrm{\μ}}_I(\mathrm{\α},\mathrm{\β},\mathrm{\γ})=\frac{{f_1}^2(\mathrm{\α}/f_1+\mathrm{\β}/f_2+\mathrm{\γ}/f_3)}{\mathrm{\α}{f}_1+\mathrm{\β}{f}_2+\mathrm{\γ}{f}_3}---\left(3\right)$

f₁,f₂和f₃为北斗三频B1、B2和B3的中心频率；

组合后的波长λ_C及模糊度分别为：

$\left\{\begin{array}{c}{\mathrm{\λ}}_C=\frac{{\mathrm{\λ}}_1{\mathrm{\λ}}_2{\mathrm{\λ}}_3}{\mathrm{\α}{\mathrm{\λ}}_2{\mathrm{\λ}}_3+{\mathrm{\β\λ}}_1{\mathrm{\λ}}_3+\mathrm{\γ}{\mathrm{\λ}}_1{\mathrm{\λ}}_2}\\ \▿\mathrm{\Δ}{N}_C=\mathrm{\α}\▿\mathrm{\Δ}{N}_1+\mathrm{\β}\▿\mathrm{\Δ}{N}_2+\mathrm{\γ}\▿\mathrm{\Δ}{N}_3\end{array}\right.---\left(4\right)$

根据不同α,β和γ的取值，构成超宽巷、宽巷和窄巷的不同组合，根据不同基线长度、波长、电离层放大倍数确定长短基线组合的模型；确定长短基线组合的模型依据是：短基线采用（-4,1,4）超宽巷组合，波长8.1403m，电离层放大2.21倍；长基线时可采用最优超宽巷（1,4,-5）组合，波长6.37m，电离层放大误差0.019倍；另外（0,-1,1）超宽巷组合，波长4.88m，电离层误差放大0.06倍。

步骤2：依据伪距观测量和超宽巷载波观测量线性组合方法，确定单历元三频组合超宽巷载波相位整周模糊度，步骤如下：

步骤2.1：采集新的历元，短基线时采用载波观测组合Φ(-4,1,4)作为第一超宽巷和Φ(0,-1,1)作为第二超宽巷和B1、B2和B3的伪距观测量以及电离层延迟观测量联立方程，采用Neill模型削弱对流层延迟误差的影响，对三个伪距观测量以及两个载波相位观测量进行处理，得到几何无关模型方程，如下

$\left[\begin{array}{c}\▿{\mathrm{\Δp}}_1\\ \▿\mathrm{\Δ}{p}_2\\ \▿\mathrm{\Δ}{p}_3\\ {\mathrm{\λ}}_{\mathrm{EWL}}^1\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{EWL}}^1\\ {\mathrm{\λ}}_{\mathrm{EWL}}^2\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{EWL}}^2\\ \▿\mathrm{\Δ}{I}_{\mathrm{\υ}}\end{array}\right]=\left[\begin{array}{c}0\\ 0\\ 0\\ {-\mathrm{\λ}}_{\mathrm{EWL}}^1\\ -{\mathrm{\λ}}_{\mathrm{EWL}}^2\\ 0\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^1\\ \▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^2\end{array}\right]+\left[\begin{array}{cc}1& 1\\ 1& {f_1}^2/{f_2}^2\\ 1& {f_1}^2/{f_3}^2\\ 1& {\mathrm{\μ}}_I(-\mathrm{4,1,4})\\ 1& {\mathrm{\μ}}_I\left(\mathrm{0,1,1}\right)\\ 0& 1\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ\ρ}\\ \▿\mathrm{\ΔI}\end{array}\right]---\left(5\right)$

长基线时采用载波观测组合Φ(1,4,-5)作为第一超宽巷和Φ(0,-1,1)作为第二超宽巷和B1、B2和B3的伪距观测量以及电离层延迟观测量联立方程，采用Neill模型削弱对流层延迟误差的影响，对三个伪距观测量以及两个载波观测量进行处理，得到几何无关模型方程，如下：

$\left[\begin{array}{c}\▿{\mathrm{\Δp}}_1\\ \▿\mathrm{\Δ}{p}_2\\ \▿\mathrm{\Δ}{p}_3\\ {\mathrm{\λ}}_{\mathrm{EWL}}^1\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{EWL}}^1\\ {\mathrm{\λ}}_{\mathrm{EWL}}^2\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{EWL}}^2\\ \▿\mathrm{\Δ}{I}_{\mathrm{\υ}}\end{array}\right]=\left[\begin{array}{c}0\\ 0\\ 0\\ {-\mathrm{\λ}}_{\mathrm{EWL}}^1\\ -{\mathrm{\λ}}_{\mathrm{EWL}}^2\\ 0\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^1\\ \▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^2\end{array}\right]+\left[\begin{array}{cc}1& 1\\ 1& {f_1}^2/{f_2}^2\\ 1& {f_1}^2/{f_3}^2\\ 1& {\mathrm{\μ}}_I(-\mathrm{1,4},-5)\\ 1& {\mathrm{\μ}}_I\left(\mathrm{0,1,1}\right)\\ 0& 1\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ\ρ}\\ \▿\mathrm{\ΔI}\end{array}\right]---\left(6\right)$

式中，和为B1、B2和B3伪距观测量双差值（m），和为第一、第二超宽巷波长（m/cycle），和为第一、二超宽巷载波相位观测量双差值（cycle），和为第一、二超宽巷整周模糊度双差值，为电离层伪距观测量双差值（m），采用最小二乘平差算法处理方程（5）或（6），获得超宽巷双差模糊度和的浮点解；

步骤2.2：采用LAMBDA算法搜索出超宽巷模糊度，根据步骤2.1中的和的浮点解确定其固定解，在搜索得到的超宽巷模糊度N组可能解的基础上增加20%的范围，得到M组可能解，带入最小二乘平差方程，计算出M组可能解对应的M组残差平方和V^TPV，从M组残差平方和V^TPV中选出最小值和次小值，和的固定解应满足残差平方和V^TPV最小，且次小残差平方和V^TPV和最小残差平方和VTPV比值大于预设阈值，如下式：

其中，T为一个正数，V为残差，V^T为V的转置矩阵，P为观测权矩阵，min为最小值。在实际解算过程中，不固定T的值，通过固定模糊度解算中能容忍的最小失败率，反算出ratio检验中的阈值T代入（7）式，和式中的Ratio值进行比较，如果满足（7）式的判断条件，则认为判断成功。

采用超宽巷双差模糊度和的固定解和超宽巷双差模糊度和的浮点解的单位权中误差的一致性检验方法，进一步判断满足（7）式的超宽巷双差模糊度和固定解是否合理，若满足（8）式即认定超宽巷双差模糊度和固定解合理；

${\mathrm{\ξ}}_{{\mathrm{\χ}}^2}(f,\mathrm{\α}/2)\≤\frac{{\mathrm{\σ}}_A}{{\mathrm{\σ}}_0}\≤{\mathrm{\ξ}}_{{\mathrm{\χ}}^2}(f,1-\mathrm{\α}/2)---\left(8\right)$

式中，整数解的单位权中误差为σ_A，初始解的单位权中误差为σ₀，置信度为1-α取0.995，f为参数估计中的自由度，上式检验也称方差因子的χ²检验。

步骤3：根据宽巷载波相位观测量与步骤2所得的超宽巷载波相位观测量的线性关系，确定宽巷载波模糊度固定解，步骤如下

步骤3.1：将宽巷观测量定义为与超宽巷观测量Φ(-4,1,4)和Φ(0,-1,1)或Φ(1,4,-5)和Φ(0,-1,1)线性相关的组合观测量。短基线时表示为：

Φ_WL(i,j,k)=κΦ(-4,1,4)+ηΦ(0,-1,1)(9)

或者长基线时表示为：

Φ_WL(i,j,k)=κΦ(1,4,-5)+ηΦ(0,-1,1)(10)

式中，κ,η为不同时取零的任意整数，从满足条件的κ,η组合中经测试选出一组宽巷观测量，短基线时选出的宽巷载波相位观测组合为长基线时选出的宽巷载波相位观测组合为由于宽巷载波相位观测量与超宽巷载波相位观测量有线性关系，可以直接由步骤2所得超宽巷双差模糊度和固定解获得宽巷双差模糊度固定解。

步骤4：根据步骤3所得的宽巷双差模糊度固定解，与宽巷组合载波相位观测量以及窄巷组合载波相位观测量，确定窄巷载波模糊度；

窄巷组合载波相位观测量Φ_NL(m,n,l)与任意两组超宽巷组合载波相位观测量Φ_EWL(i,j,k)或者宽巷组合载波相位观测量Φ_WL(u,v,w)线性无关，并且窄巷组合载波相位观测量中观测噪声影响应小于0.5周，窄巷组合载波相位观测量应该具有最小的电离层折射误差，并且考虑到其他误差，窄巷组合载波相位观测量波长不宜过短，综上获得窄巷组合载波相位观测量Φ_NL(2,0,-2)。

短基线时，窄巷组合载波相位观测量：Φ_NL(2,0,-2)与步骤3中的组合联立得：

$\left[\begin{array}{c}{\mathrm{\λ}}_{\mathrm{WL}}^1\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{WL}}^1\\ {\mathrm{\λ}}_{\mathrm{NL}}\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{NL}}\end{array}\right]=\left[\begin{array}{c}-{\mathrm{\λ}}_{\mathrm{WL}}^1\\ {\mathrm{\λ}}_{\mathrm{NL}}\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{\mathrm{WL}}^1\\ \▿\mathrm{\Δ}{N}_{\mathrm{NL}}\end{array}\right]+\left[\begin{array}{cc}1& {\mathrm{\μ}}_I(-4,-\mathrm{3,8})\\ 1& {\mathrm{\μ}}_I(\mathrm{2,0},-2)\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ\ρ}\\ \▿\mathrm{\ΔI}\end{array}\right]---\left(11\right)$

长基线时，窄巷组合载波观测量：Φ_NL(2,0,-2)与步骤3中的组合联立得：

$\left[\begin{array}{c}{\mathrm{\λ}}_{\mathrm{WL}}^2\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{WL}}^2\\ {\mathrm{\λ}}_{\mathrm{NL}}\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{NL}}\end{array}\right]=\left[\begin{array}{c}-{\mathrm{\λ}}_{\mathrm{WL}}^2\\ {\mathrm{\λ}}_{\mathrm{NL}}\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{\mathrm{WL}}^2\\ \▿\mathrm{\Δ}{N}_{\mathrm{NL}}\end{array}\right]+\left[\begin{array}{cc}1& {\mathrm{\μ}}_I(-1,-\mathrm{1,0})\\ 1& {\mathrm{\μ}}_I(\mathrm{2,0},-2)\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ\ρ}\\ \▿\mathrm{\ΔI}\end{array}\right]---\left(12\right)$

式中，和由步骤2固定后的和代入(5)式或(6)式重新解算得到。

据上式最小二乘平差可求得窄巷载波模糊度浮点解，然后根据步骤2.3获得窄巷组合模糊度的固定解。

步骤5：依据步骤2所得超宽巷双差模糊度和的固定解及步骤4所得窄巷组合模糊度的固定解，确定基础载波模糊度固定，方程如下：

短基线情况，

$\left[\begin{array}{ccc}2& 0& -2\\ -4& 1& 4\\ 0& -1& 1\end{array}\right]\left[\begin{array}{c}{\▿\mathrm{\ΔN}}_1\\ {\▿\mathrm{\ΔN}}_2\\ {\▿\mathrm{\ΔN}}_3\end{array}\right]=\left[\begin{array}{c}{\▿\mathrm{\ΔN}}_{\mathrm{NL}}\\ \▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^1\\ {\▿\mathrm{\ΔN}}_{\mathrm{EWL}}^2\end{array}\right]---\left(13\right)$

长基线情况，

$\left[\begin{array}{ccc}2& 0& -2\\ 1& 4& -5\\ 0& -1& 1\end{array}\right]\left[\begin{array}{c}{\▿\mathrm{\ΔN}}_1\\ {\▿\mathrm{\ΔN}}_2\\ {\▿\mathrm{\ΔN}}_3\end{array}\right]=\left[\begin{array}{c}{\▿\mathrm{\ΔN}}_{\mathrm{NL}}\\ \▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^1\\ {\▿\mathrm{\ΔN}}_{\mathrm{EWL}}^2\end{array}\right]---\left(14\right)$

求解方程即得最终获得基础载波模糊度和从而完成三频模糊度的单历元固定。

本发明与现有技术相比具有如下优点：

（1）本发明基于常规三频模糊度解算方法，结合北斗系统特点，结合LAMBDA算法，采用新的逐级固定方法，保证单历元情况下，测站可以得到较高可靠性的整周模糊度；

（2）本发明能够提供长距离高精度的实时相位定位结果；

（3）本发明确定基础载波模糊度效率更高，能适应于实时动态的环境中。

附图说明

图1为本发明的具体实施例工作流程的说明。

具体实施方式

一种北斗三频信号载波相位整周模糊度单历元确定方法，基于几何无关模型，逐级分步确定各组合模糊度，完成三频信号载波相位整周模糊度单历元的解算，所述逐级分布步骤如下：

步骤1：依据基线长度和波长、电离层放大倍数，确定北斗三频组合载波相位双差观测模型，北斗三频组合载波相位双差观测模型写为：

$\▿\mathrm{\Δ}{\mathrm{\Φ}}_C=\mathrm{\α}\▿\mathrm{\Δ}{\mathrm{\Φ}}_1+\mathrm{\β}\▿\mathrm{\Δ}{\mathrm{\Φ}}_2+\mathrm{\γ}\▿\mathrm{\Δ}{\mathrm{\Φ}}_3---\left(1\right)$

其中，为组合载波相位双差值，为B1载波相位双差值、为B2载波相位双差值、为B3载波相位双差值；

为使几何距离不受观测值组合不同而变化，将组合观测模型表示为一般形式为：

$\▿\mathrm{\Δ}{\mathrm{\Φ}}_C=\frac{\▿\mathrm{\Δ\ρ}}{{\mathrm{\λ}}_C}+\▿\mathrm{\Δ}{N}_C-\frac{{\mathrm{\μ}}_I(\mathrm{\α},\mathrm{\β},\mathrm{\γ})\▿\mathrm{\ΔI}}{{\mathrm{\λ}}_C}+\frac{\▿\mathrm{\Δ}{\mathrm{\ϵ}}_C}{{\mathrm{\λ}}_C}---\left(2\right)$

式中，为双差算子，为卫星与测站之间的几何距离双差值（m），并包含所有与频率无关的误差残差项（如对流层残差、钟差残差等），为B1载波频率上的电离层延迟双差值，为组合观测噪声。

组合电离层比例因子μ_I(α,β,γ)：

${\mathrm{\μ}}_I(\mathrm{\α},\mathrm{\β},\mathrm{\γ})=\frac{{f_1}^2(\mathrm{\α}/f_1+\mathrm{\β}/f_2+\mathrm{\γ}/f_3)}{\mathrm{\α}{f}_1+\mathrm{\β}{f}_2+\mathrm{\γ}{f}_3}---\left(3\right)$

f₁,f₂和f₃为北斗三频B1、B2和B3的中心频率。

组合后的波长λ_C及模糊度分别为：

$\left\{\begin{array}{c}{\mathrm{\λ}}_C=\frac{{\mathrm{\λ}}_1{\mathrm{\λ}}_2{\mathrm{\λ}}_3}{\mathrm{\α}{\mathrm{\λ}}_2{\mathrm{\λ}}_3+{\mathrm{\β\λ}}_1{\mathrm{\λ}}_3+\mathrm{\γ}{\mathrm{\λ}}_1{\mathrm{\λ}}_2}\\ \▿\mathrm{\Δ}{N}_C=\mathrm{\α}\▿\mathrm{\Δ}{N}_1+\mathrm{\β}\▿\mathrm{\Δ}{N}_2+\mathrm{\γ}\▿\mathrm{\Δ}{N}_3\end{array}\right.---\left(4\right)$

根据不同α,β和γ的取值，构成超宽巷、宽巷和窄巷的不同组合，根据不同基线长度、波长、电离层放大倍数确定长短基线组合的模型；确定长短基线组合的模型依据是：短基线采用（-4,1,4）超宽巷组合，波长8.1403m，电离层放大2.21倍；长基线时可采用最优超宽巷（1,4,-5）组合，波长6.37m，电离层放大误差0.019倍；另外（0,-1,1）超宽巷组合，波长4.88m，电离层误差放大0.06倍。

步骤2：依据伪距观测量和超宽巷载波观测量线性组合方法，确定单历元三频组合超宽巷载波相位整周模糊度，步骤如下：

步骤2.1：采集新的历元，短基线时采用载波观测组合Φ(-4,1,4)作为第一超宽巷和Φ(0,-1,1)作为第二超宽巷和B1、B2和B3的伪距观测量以及电离层延迟观测量联立方程，采用Neill模型削弱对流层延迟误差的影响，对三个伪距观测量以及两个载波相位观测量进行处理，得到几何无关模型方程，如下

$\left[\begin{array}{c}\▿{\mathrm{\Δp}}_1\\ \▿\mathrm{\Δ}{p}_2\\ \▿\mathrm{\Δ}{p}_3\\ {\mathrm{\λ}}_{\mathrm{EWL}}^1\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{EWL}}^1\\ {\mathrm{\λ}}_{\mathrm{EWL}}^2\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{EWL}}^2\\ \▿\mathrm{\Δ}{I}_{\mathrm{\υ}}\end{array}\right]=\left[\begin{array}{c}0\\ 0\\ 0\\ {-\mathrm{\λ}}_{\mathrm{EWL}}^1\\ -{\mathrm{\λ}}_{\mathrm{EWL}}^2\\ 0\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^1\\ \▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^2\end{array}\right]+\left[\begin{array}{cc}1& 1\\ 1& {f_1}^2/{f_2}^2\\ 1& {f_1}^2/{f_3}^2\\ 1& {\mathrm{\μ}}_I(-\mathrm{4,1,4})\\ 1& {\mathrm{\μ}}_I\left(\mathrm{0,1,1}\right)\\ 0& 1\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ\ρ}\\ \▿\mathrm{\ΔI}\end{array}\right]---\left(5\right)$

长基线时采用载波观测组合Φ(1,4,-5)作为第一超宽巷和Φ(0,-1,1)作为第二超宽巷和B1、B2和B3的伪距观测量以及电离层延迟观测量联立方程，采用Neill模型削弱对流层延迟误差的影响，对三个伪距观测量以及两个载波观测量进行处理，得到几何无关模型方程，如下：

$\left[\begin{array}{c}\▿{\mathrm{\Δp}}_1\\ \▿\mathrm{\Δ}{p}_2\\ \▿\mathrm{\Δ}{p}_3\\ {\mathrm{\λ}}_{\mathrm{EWL}}^1\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{EWL}}^1\\ {\mathrm{\λ}}_{\mathrm{EWL}}^2\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{EWL}}^2\\ \▿\mathrm{\Δ}{I}_{\mathrm{\υ}}\end{array}\right]=\left[\begin{array}{c}0\\ 0\\ 0\\ {-\mathrm{\λ}}_{\mathrm{EWL}}^1\\ -{\mathrm{\λ}}_{\mathrm{EWL}}^2\\ 0\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^1\\ \▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^2\end{array}\right]+\left[\begin{array}{cc}1& 1\\ 1& {f_1}^2/{f_2}^2\\ 1& {f_1}^2/{f_3}^2\\ 1& {\mathrm{\μ}}_I(-\mathrm{1,4},-5)\\ 1& {\mathrm{\μ}}_I\left(\mathrm{0,1,1}\right)\\ 0& 1\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ\ρ}\\ \▿\mathrm{\ΔI}\end{array}\right]---\left(6\right)$

式中，和为B1、B2和B3伪距观测量双差值（m），和为第一、第二超宽巷波长（m/cycle），和为第一、二超宽巷载波相位观测量双差值（cycle），和为第一、二超宽巷整周模糊度双差值，为电离层伪距观测量双差值（m），采用最小二乘平差算法处理方程（5）或（6），获得超宽巷双差模糊度和的浮点解；

步骤2.2：采用LAMBDA算法搜索出超宽巷模糊度，根据步骤2.1中的和的浮点解确定其固定解，在搜索得到的超宽巷模糊度N组可能解的基础上增加20%的范围，得到M组可能解，带入最小二乘平差方程，计算出M组可能解对应的M组残差平方和V^TPV，从M组残差平方和V^TPV中选出最小值和次小值，和的固定解应满足残差平方和V^TPV最小，且次小残差平方和V^TPV和最小残差平方和V^TPV比值大于预设阈值，如下式：

其中，T为一个正数，一般大于2.0，V为残差，V^T为V的转置矩阵，P为观测权矩阵，min为最小值。在实际解算过程中，不固定T的值，通过固定模糊度解算中能容忍的最小失败率（如0.005、0.001等），反算出ratio检验中的阈值T代入（7）式，和式中的Ratio值进行比较，如果满足（7）式的判断条件，则认为判断成功。

采用超宽巷双差模糊度和的固定解和超宽巷双差模糊度和的浮点解的单位权中误差的一致性检验方法，进一步判断满足（7）式的超宽巷双差模糊度和固定解是否合理，若满足（8）式即认定超宽巷双差模糊度和固定解合理；

${\mathrm{\ξ}}_{{\mathrm{\χ}}^2}(f,\mathrm{\α}/2)\≤\frac{{\mathrm{\σ}}_A}{{\mathrm{\σ}}_0}\≤{\mathrm{\ξ}}_{{\mathrm{\χ}}^2}(f,1-\mathrm{\α}/2)---\left(8\right)$

式中，整数解的单位权中误差为σ_A，初始解的单位权中误差为σ₀，置信度为1-α取0.995，f为参数估计中的自由度，上式检验也称方差因子的χ²检验。

步骤3：根据宽巷载波相位观测量与步骤2所得的超宽巷载波相位观测量的线性关系，确定宽巷载波模糊度固定解，步骤如下

步骤3.1：将宽巷观测量定义为与超宽巷观测量Φ(-4,1,4)和Φ(0,-1,1)或Φ(1,4,-5)和Φ(0,-1,1)线性相关的组合观测量。短基线时表示为：

Φ_WL(i,j,k)=κΦ(-4,1,4)+ηΦ(0,-1,1)(9)

或者长基线时表示为：

Φ_WL(i,j,k)=κΦ(1,4,-5)+ηΦ(0,-1,1)(10)

式中，κ,η为不同时取零的任意整数，从满足条件的κ,η组合中经测试选出一组宽巷观测量，短基线时选出的宽巷载波相位观测组合为长基线时选出的宽巷载波相位观测组合为由于宽巷载波相位观测量与超宽巷载波相位观测量有线性关系，可以直接由步骤2所得超宽巷双差模糊度和固定解获得宽巷双差模糊度固定解。

步骤4：根据步骤3所得的宽巷双差模糊度固定解，与宽巷组合载波相位观测量以及窄巷组合载波相位观测量，确定窄巷载波模糊度。

窄巷组合载波相位观测量Φ_NL(m,n,l)与任意两组超宽巷组合载波相位观测量Φ_EWL(i,j,k)或者宽巷组合载波相位观测量Φ_WL(u,v,w)线性无关，并且窄巷组合载波相位观测量中观测噪声影响应小于0.5周，窄巷组合载波相位观测量应该具有最小的电离层折射误差，并且考虑到其他误差，窄巷组合载波相位观测量波长不宜过短，综上获得窄巷组合载波相位观测量Φ_NL(2,0,-2)。

短基线时，窄巷组合载波相位观测量：Φ_NL(2,0,-2)与步骤3中的组合联立得：

$\left[\begin{array}{c}{\mathrm{\λ}}_{\mathrm{WL}}^1\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{WL}}^1\\ {\mathrm{\λ}}_{\mathrm{NL}}\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{NL}}\end{array}\right]=\left[\begin{array}{c}-{\mathrm{\λ}}_{\mathrm{WL}}^1\\ {\mathrm{\λ}}_{\mathrm{NL}}\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{\mathrm{WL}}^1\\ \▿\mathrm{\Δ}{N}_{\mathrm{NL}}\end{array}\right]+\left[\begin{array}{cc}1& {\mathrm{\μ}}_I(-4,-\mathrm{3,8})\\ 1& {\mathrm{\μ}}_I(\mathrm{2,0},-2)\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ\ρ}\\ \▿\mathrm{\ΔI}\end{array}\right]---\left(11\right)$

长基线时，窄巷组合载波观测量：Φ_NL(2,0,-2)与步骤3中的组合联立得：

$\left[\begin{array}{c}{\mathrm{\λ}}_{\mathrm{WL}}^2\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{WL}}^2\\ {\mathrm{\λ}}_{\mathrm{NL}}\▿\mathrm{\Δ}{\mathrm{\Φ}}_{\mathrm{NL}}\end{array}\right]=\left[\begin{array}{c}-{\mathrm{\λ}}_{\mathrm{WL}}^2\\ {\mathrm{\λ}}_{\mathrm{NL}}\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{\mathrm{WL}}^2\\ \▿\mathrm{\Δ}{N}_{\mathrm{NL}}\end{array}\right]+\left[\begin{array}{cc}1& {\mathrm{\μ}}_I(-1,-\mathrm{1,0})\\ 1& {\mathrm{\μ}}_I(\mathrm{2,0},-2)\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ\ρ}\\ \▿\mathrm{\ΔI}\end{array}\right]---\left(12\right)$

式中，和由步骤2固定后的和代入(5)式或(6)式重新解算得到。

据上式最小二乘平差可求得窄巷载波模糊度浮点解，然后根据步骤2.3获得窄巷组合模糊度的固定解。

步骤5：依据步骤2所得超宽巷双差模糊度和的固定解及步骤4所得窄巷组合模糊度的固定解，确定基础载波模糊度固定，方程如下：

短基线情况，

$\left[\begin{array}{ccc}2& 0& -2\\ -4& 1& 4\\ 0& -1& 1\end{array}\right]\left[\begin{array}{c}{\▿\mathrm{\ΔN}}_1\\ {\▿\mathrm{\ΔN}}_2\\ {\▿\mathrm{\ΔN}}_3\end{array}\right]=\left[\begin{array}{c}{\▿\mathrm{\ΔN}}_{\mathrm{NL}}\\ \▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^1\\ {\▿\mathrm{\ΔN}}_{\mathrm{EWL}}^2\end{array}\right]---\left(13\right)$

长基线情况，

$\left[\begin{array}{ccc}2& 0& -2\\ 1& 4& -5\\ 0& -1& 1\end{array}\right]\left[\begin{array}{c}{\▿\mathrm{\ΔN}}_1\\ {\▿\mathrm{\ΔN}}_2\\ {\▿\mathrm{\ΔN}}_3\end{array}\right]=\left[\begin{array}{c}{\▿\mathrm{\ΔN}}_{\mathrm{NL}}\\ \▿\mathrm{\Δ}{N}_{\mathrm{EWL}}^1\\ {\▿\mathrm{\ΔN}}_{\mathrm{EWL}}^2\end{array}\right]---\left(14\right)$

求解方程即得最终获得基础载波模糊度和从而完成三频模糊度的单历元固定。

Claims (2)

1.一种北斗三频信号载波相位整周模糊度单历元确定方法，其特征在于基于几何无关模型，逐级分步确定各组合模糊度，完成三频信号载波相位整周模糊度单历元的解算，所述逐级分布步骤如下：

步骤1：依据基线长度和波长、电离层放大倍数，确定北斗三频组合载波相位双差观测模型，北斗三频组合载波相位双差观测模型写为：

$\▿{\mathrm{\Δ\Φ}}_C=\mathrm{\α}\▿{\mathrm{\Δ\Φ}}_1+\mathrm{\β}\▿{\mathrm{\Δ\Φ}}_2+\mathrm{\γ}\▿{\mathrm{\Δ\Φ}}_3---\left(1\right)$

其中，为组合载波相位双差值，为B1载波相位双差值、为B2载波相位双差值、为B3载波相位双差值；

为使几何距离不受观测值组合不同而变化，将组合观测模型表示为一般形式为：

$\▿{\mathrm{\Δ\Φ}}_C=\frac{\▿\mathrm{\Δ}\mathrm{\ρ}}{{\mathrm{\λ}}_C}+\▿{\mathrm{\ΔN}}_C-\frac{{\mathrm{\μ}}_I(\mathrm{\α},\mathrm{\β},\mathrm{\γ})\▿\mathrm{\Δ}I}{{\mathrm{\λ}}_C}+\frac{\▿{\mathrm{\Δ\ϵ}}_C}{{\mathrm{\λ}}_C}---\left(2\right)$

式中，为双差算子，为卫星与测站之间的几何距离双差值m，并包含所有与频率无关的误差残差项，为B1载波频率上的电离层延迟双差值，为组合观测噪声；

组合电离层比例因子μ_I(α,β,γ)：

${\mathrm{\μ}}_I(\mathrm{\α},\mathrm{\β},\mathrm{\γ})=\frac{f_1^2(\mathrm{\α}/f_1+\mathrm{\β}/f_2+\mathrm{\γ}/f_3)}{{\mathrm{\αf}}_1+{\mathrm{\βf}}_2+{\mathrm{\γf}}_3}---\left(3\right)$

f₁,f₂和f₃为北斗三频B1、B2和B3的中心频率；

组合后的波长λ_C及模糊度分别为：

$\left\{\begin{array}{c}{\mathrm{\λ}}_C=\frac{{\mathrm{\λ}}_1{\mathrm{\λ}}_2{\mathrm{\λ}}_3}{{\mathrm{\α\λ}}_2{\mathrm{\λ}}_3+{\mathrm{\β\λ}}_1{\mathrm{\λ}}_3+{\mathrm{\γ\λ}}_1{\mathrm{\λ}}_2}\\ \▿{\mathrm{\ΔN}}_C=\mathrm{\α}\▿{\mathrm{\ΔN}}_1+\mathrm{\β}\▿{\mathrm{\ΔN}}_2+\mathrm{\γ}\▿{\mathrm{\ΔN}}_3\end{array}\right.---\left(4\right)$

根据不同α,β和γ的取值，构成超宽巷、宽巷和窄巷的不同组合，根据不同基线长度、波长、电离层放大倍数确定长短基线组合的模型；

步骤2：依据伪距观测量和超宽巷载波观测量线性组合方法，确定单历元三频组合超宽巷载波相位整周模糊度，步骤如下：

步骤2.1：采集新的历元，短基线时采用载波观测组合Φ(-4,1,4)作为第一超宽巷和Φ(0,-1,1)作为第二超宽巷和B1、B2和B3的伪距观测量以及电离层延迟观测量联立方程，采用Neill模型削弱对流层延迟误差的影响，对三个伪距观测量以及两个载波相位观测量进行处理，得到几何无关模型方程，如下：

$\left[\begin{array}{c}\▿{\mathrm{\Δp}}_1\\ \▿{\mathrm{\Δp}}_2\\ \▿{\mathrm{\Δp}}_3\\ {\mathrm{\λ}}_{EWL}^1\▿{\mathrm{\Δ\Φ}}_{EWL}^1\\ {\mathrm{\λ}}_{EWL}^2\▿{\mathrm{\Δ\Φ}}_{EWL}^2\\ \▿{\mathrm{\ΔI}}_{\mathrm{\υ}}\end{array}\right]=\left[\begin{array}{c}0\\ 0\\ 0\\ -{\mathrm{\λ}}_{EWL}^1\\ -{\mathrm{\λ}}_{EWL}^2\\ 0\end{array}\right]=\left[\begin{array}{c}\▿{\mathrm{\ΔN}}_{EWL}^1\\ \▿{\mathrm{\ΔN}}_{EWL}^2\end{array}\right]+\left[\begin{array}{cc}1& 1\\ 1& f_1^2/f_2^2\\ 1& f_1^2/f_3^2\\ 1& {\mathrm{\μ}}_I(-4,1,4)\\ 1& {\mathrm{\μ}}_I(0,-1,1)\\ 0& 1\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}\mathrm{\ρ}\\ \▿\mathrm{\Δ}I\end{array}\right]---\left(5\right)$

长基线时采用载波观测组合Φ(1,4,-5)作为第一超宽巷和Φ(0,-1,1)作为第二超宽巷和B1、B2和B3的伪距观测量以及电离层延迟观测量联立方程，采用Neill模型削弱对流层延迟误差的影响，对三个伪距观测量以及两个载波观测量进行处理，得到几何无关模型方程，如下：

$\left[\begin{array}{c}\▿{\mathrm{\Δp}}_1\\ \▿{\mathrm{\Δp}}_2\\ \▿{\mathrm{\Δp}}_3\\ {\mathrm{\λ}}_{EWL}^1\▿{\mathrm{\Δ\Φ}}_{EWL}^1\\ {\mathrm{\λ}}_{EWL}^2\▿{\mathrm{\Δ\Φ}}_{EWL}^2\\ \▿{\mathrm{\ΔI}}_{\mathrm{\υ}}\end{array}\right]=\left[\begin{array}{c}0\\ 0\\ 0\\ -{\mathrm{\λ}}_{EWL}^1\\ -{\mathrm{\λ}}_{EWL}^2\\ 0\end{array}\right]=\left[\begin{array}{c}\▿{\mathrm{\ΔN}}_{EWL}^1\\ \▿{\mathrm{\ΔN}}_{EWL}^2\end{array}\right]+\left[\begin{array}{cc}1& 1\\ 1& f_1^2/f_2^2\\ 1& f_1^2/f_3^2\\ 1& {\mathrm{\μ}}_I(1,4,-5)\\ 1& {\mathrm{\μ}}_I(0,-1,1)\\ 0& 1\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}\mathrm{\ρ}\\ \▿\mathrm{\Δ}I\end{array}\right]---\left(6\right)$

式中，和为B1、B2和B3伪距观测量双差值(m)，和为第一、第二超宽巷波长(m/cycle)，和为第一、二超宽巷载波相位观测量双差值(cycle)，和为第一、二超宽巷整周模糊度双差值，为电离层伪距观测量双差值(m)，采用最小二乘平差算法处理方程(5)或(6)，获得超宽巷载波相位双差模糊度和的浮点解及其协方差阵或

步骤2.2：采用LAMBDA算法搜索出超宽巷载波相位双差模糊度，将步骤2.1中得到的或的浮点解及其协方差阵或带入LAMBDA算法确定其固定解，算法中将得到的超宽巷载波相位双差模糊度N组可能解增加20％的范围，得到M组可能解，计算出M组可能解对应的M组残差平方和V^TPV，从M组残差平方和V^TPV中选出最小值和次小值，和的固定解应满足残差平方和V^TPV最小，且次小残差平方和V^TPV和最小残差平方和V^TPV比值大于预设阈值，如下式：

其中，T为一个正数，V为残差，V^T为V的转置矩阵，P为观测权矩阵，min为最小值；在实际解算过程中，不固定T的值，通过固定模糊度解算中能容忍的最小失败率，反算出ratio检验中的阈值T代入(7)式，和式中的Ratio值进行比较，如果满足(7)式的判断条件，则认为判断成功；

采用超宽巷双差模糊度和的固定解和超宽巷双差模糊度和的浮点解的单位权中误差的一致性检验方法，进一步判断满足(7)式的超宽巷双差模糊度和固定解是否合理，若满足(8)式即认定超宽巷双差模糊度和固定解合理；

${\mathrm{\ξ}}_{{\mathrm{\χ}}^2}(f,\mathrm{\α}/2)\≤\frac{{\mathrm{\σ}}_A}{{\mathrm{\σ}}_0}\≤{\mathrm{\ξ}}_{{\mathrm{\χ}}^2}(f,1-\mathrm{\α}/2)---\left(8\right)$

式中，整数解的单位权中误差为σ_A，初始解的单位权中误差为σ₀，置信度为1-α，f为参数估计中的自由度，上式检验也称方差因子的χ²检验；

步骤3：根据宽巷载波相位观测量与步骤2所得的超宽巷载波相位观测量的线性关系，确定宽巷载波模糊度固定解，步骤如下：

步骤3.1：将宽巷观测量定义为与超宽巷观测量Φ(-4,1,4)和Φ(0,-1,1)或Φ(1,4,-5)和Φ(0,-1,1)线性相关的组合观测量，短基线时表示为：

Φ_WL(i,j,k)＝κΦ(-4,1,4)+ηΦ(0,-1,1)(9)

或者长基线时表示为：

Φ_WL(i,j,k)＝κΦ(1,4,-5)+ηΦ(0,-1,1)(10)

式中，κ,η为不同时取零的任意整数，从满足条件的κ,η组合中经测试选出一组宽巷观测量，短基线时选出的宽巷载波相位观测组合为长基线时选出的宽巷载波相位观测组合为由于宽巷载波相位观测量与超宽巷载波相位观测量有线性关系，可以直接由步骤2所得超宽巷双差模糊度和固定解获得宽巷双差模糊度固定解；

步骤4：根据步骤3所得的宽巷双差模糊度固定解，与宽巷组合载波相位观测量以及窄巷组合载波相位观测量，确定窄巷载波模糊度；窄巷组合载波相位观测量Φ_NL(m,n,l)与任意两组超宽巷组合载波相位观测量Φ_EWL(i,j,k)或者宽巷组合载波相位观测量Φ_WL(u,v,w)线性无关，并且窄巷组合载波相位观测量中观测噪声影响应小于0.5周，窄巷组合载波相位观测量应该具有最小的电离层折射误差，并且考虑到其他误差，窄巷组合载波相位观测量波长不宜过短，综上获得窄巷组合载波相位观测量Φ_NL(2,0,-2)；

短基线时，窄巷组合载波相位观测量：Φ_NL(2,0,-2)与步骤3中的组合联立得：

$\left[\begin{array}{c}{\mathrm{\λ}}_{WL}^1\▿\mathrm{\Δ}{\mathrm{\Φ}}_{WL}^1\\ {\mathrm{\λ}}_{NL}\▿{\mathrm{\Δ\Φ}}_{NL}\end{array}\right]=\left[\begin{array}{c}-{\mathrm{\λ}}_{WL}^1\\ {\mathrm{\λ}}_{NL}\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{WL}^1\\ \▿{\mathrm{\ΔN}}_{NL}\end{array}\right]+\left[\begin{array}{cc}1& {\mathrm{\μ}}_I(-4,-3,8)\\ 1& {\mathrm{\μ}}_I(2,0,-2)\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}\mathrm{\ρ}\\ \▿\mathrm{\Δ}I\end{array}\right]---\left(11\right)$

长基线时，窄巷组合载波观测量：Φ_NL(2,0,-2)与步骤3中的组合联立得：

$\left[\begin{array}{c}{\mathrm{\λ}}_{WL}^2\▿\mathrm{\Δ}{\mathrm{\Φ}}_{WL}^2\\ {\mathrm{\λ}}_{NL}\▿{\mathrm{\Δ\Φ}}_{NL}\end{array}\right]=\left[\begin{array}{c}-{\mathrm{\λ}}_{WL}^2\\ {\mathrm{\λ}}_{NL}\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{WL}^2\\ \▿{\mathrm{\ΔN}}_{NL}\end{array}\right]+\left[\begin{array}{cc}1& {\mathrm{\μ}}_I(1,-1,0)\\ 1& {\mathrm{\μ}}_I(2,0,-2)\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}\mathrm{\ρ}\\ \▿\mathrm{\Δ}I\end{array}\right]---\left(12\right)$

式中，和由步骤2固定后的和代入(5)式或(6)式重新解算得到；

据上式最小二乘平差可求得窄巷载波模糊度浮点解，然后根据步骤2.3获得窄巷组合模糊度的固定解；

步骤5：依据步骤2所得超宽巷双差模糊度和的固定解及步骤4所得窄巷组合模糊度的固定解，确定基础载波模糊度固定，方程如下：

短基线情况：

$\left[\begin{array}{ccc}2& 0& -2\\ -4& 1& 4\\ 0& -1& 1\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_1\\ \▿{\mathrm{\ΔN}}_2\\ \▿{\mathrm{\ΔN}}_3\end{array}\right]=\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{NL}\\ \▿{\mathrm{\ΔN}}_{EWL}^1\\ \▿{\mathrm{\ΔN}}_{EWL}^2\end{array}\right]---\left(13\right)$

长基线情况：

$\left[\begin{array}{ccc}2& 0& -2\\ 1& 4& -5\\ 0& -1& 1\end{array}\right]\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_1\\ \▿{\mathrm{\ΔN}}_2\\ \▿{\mathrm{\ΔN}}_3\end{array}\right]=\left[\begin{array}{c}\▿\mathrm{\Δ}{N}_{NL}\\ \▿{\mathrm{\ΔN}}_{EWL}^1\\ \▿{\mathrm{\ΔN}}_{EWL}^2\end{array}\right]---\left(14\right)$

求解方程即得最终获得基础载波模糊度和从而完成三频模糊度的单历元固定。

2.根据权利要求1所述的一种北斗三频信号载波相位整周模糊度单历元确定方法，其特征在于：所述步骤1中确定长短基线组合的模型依据是：短基线采用(-4,1,4)超宽巷组合，波长8.1403m，电离层放大2.21倍；长基线时可采用最优超宽巷(1,4,-5)组合，波长6.37m，电离层放大误差0.019倍；另外(0,-1,1)超宽巷组合，波长4.88m，电离层误差放大0.06倍。