CN108037521A - 一种基于北斗超宽巷约束的bds/gps宽巷模糊度单历元固定方法 - Google Patents
一种基于北斗超宽巷约束的bds/gps宽巷模糊度单历元固定方法 Download PDFInfo
- Publication number
- CN108037521A CN108037521A CN201810082409.0A CN201810082409A CN108037521A CN 108037521 A CN108037521 A CN 108037521A CN 201810082409 A CN201810082409 A CN 201810082409A CN 108037521 A CN108037521 A CN 108037521A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- wide
- bds
- gps
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
- G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种基于北斗超宽巷约束的BDS/GPS宽巷模糊度单历元固定方法,包括如下步骤:利用GNSS接收机接收的BDS三频伪距、载波观测数据组合单历元解算BDS的超宽巷(0,‑1,1)和(1,4,‑5)的组合模糊度;组合两个固定后的超宽巷模糊度,计算观测噪声和电离层延迟较小的北斗(1,‑1,0)的宽巷组合模糊度;当BDS宽巷模糊度固定之后,将其作为高精度的约束条件,与接收机接收的GPS双频观测数据组合得到的GPS宽巷观测方程相联立;对上一步得到的观测模型进行卡尔曼滤波解算,得到GPS双差宽巷模糊度的浮点解和方差协方差矩阵,采用LAMBDA算法进行整数固定,得到GPS双差宽巷模糊度整数解。本发明实现BDS/GPS宽巷模糊度的单历元固定,有效提升GNSS基线解算的时效性。
Description
技术领域
本发明涉及全球卫星导航系统定位技术领域,尤其是一种基于北斗超宽巷约束的BDS/GPS宽巷模糊度单历元固定方法。
背景技术
2012年底,北斗二代卫星导航定位系统正式向亚太地区提供服务,作为国际上首个提供三频服务的卫星导航定位系统,北斗开启了GNSS三频定位时代。此外,随着GPS、GLONASS的现代化以及Galileo的不断完善,多频逐渐取代双频成为GNSS的发展趋势。
利用多频GNSS信号可以构造出许多具有长波长、小噪声、弱电离层影响的线性组合观测值,极大的丰富了GNSS的解算方案。多频GNSS模糊度解算可避开传统搜索方法的复杂计算,根据不同组合观测值的波长及误差特点,采用简单的舍入取整法逐级固定各组合模糊度,目前GPS、GLONASS系统播发三频信号相对较少,Galileo系统还不够完善,在一些高遮挡环境下很难利用单BDS系统实现多频解算,且由于GEO卫星带来的影响,在某些地区图形结构弱,无法实现高精度定位。宽巷模糊度的快速解算是目前基线解算中必不可少的环节,目前常用的宽巷模糊度解算方法有双频P码和相位观测值线性组合法(MW组合),但是MW组合受伪距噪声影响较大,需要一段时间的平滑才能实现宽巷模糊度的正确固定。充分利用BDS三频超宽巷、宽巷模糊度易于固定的优势,辅助GPS等其余星座宽巷模糊度的快速固定,是推进GNSS高精度应用的重要内容。
发明内容
本发明所要解决的技术问题在于,提供一种基于北斗超宽巷约束的BDS/GPS宽巷模糊度单历元固定方法,能够实现BDS/GPS宽巷模糊度的单历元固定,有效提升了GNSS基线解算的时效性。
为解决上述技术问题,本发明提供一种基于北斗超宽巷约束的BDS/GPS宽巷模糊度单历元固定方法,包括如下步骤:
(1)利用GNSS接收机接收的BDS三频伪距、载波观测数据组合单历元解算BDS的超宽巷(0,-1,1)和(1,4,-5)的组合模糊度;
(2)组合两个固定后的超宽巷模糊度,计算观测噪声和电离层延迟较小的北斗(1,-1,0)的宽巷组合模糊度;
(3)当BDS宽巷模糊度固定之后,将其作为高精度的约束条件,与接收机接收的GPS双频观测数据组合得到的GPS宽巷观测方程相联立;
(4)对步骤(3)得到的观测模型进行卡尔曼滤波解算,得到GPS双差宽巷模糊度的浮点解和方差协方差矩阵,采用LAMBDA算法进行整数固定,得到GPS双差宽巷模糊度整数解。
优选的,利用GNSS接收机接收的BDS三频伪距、载波观测数据组合单历元解算BDS的超宽巷(0,-1,1)和(1,4,-5)的组合模糊度具体为:
利用载波和伪距组合单历元解算BDS的超宽巷(0,-1,1)组合模糊度,如式(1)所示:
式中,为双差算子,为(0,-1,1)组合模糊度,[·]四舍五入算子,为以米为单位的(0,-1,1)组合的载波观测值,为以米为单位的(0,1,1)组合的伪距观测值,λ(0,-1,1)为(0,-1,1)组合载波观测值波长;
利用载波、伪距组合单历元解算BDS的超宽巷(1,4,-5)组合模糊度,如式(2)所示:
式中,为双差算子,为(1,4,-5)组合模糊度,[·]四舍五入算子,为以米为单位的(1,4,-5)组合的载波观测值,为以米为单位的(0,1,1)组合的伪距观测值,λ(1,4,-5)为(1,4,-5)组合载波观测值波长。
优选的,步骤(2)中,组合两个固定后的超宽巷模糊度,计算观测噪声和电离层延迟较小的北斗(1,-1,0)的宽巷组合模糊度具体为:
组合两个固定后的超宽巷模糊度,计算观测噪声和电离层延迟较小的北斗(1,-1,0)宽巷组合模糊度,如式(3)所示。
式中,为双差算子,为(1,-1,0)组合模糊度,为(0,-1,1)组合模糊度,为(1,4,-5)组合模糊度。
优选的,步骤(3)中,当BDS宽巷模糊度固定之后,将其作为高精度的约束条件,与接收机接收的GPS双频观测数据组合得到的GPS宽巷观测方程相联立具体为:
当BDS宽巷模糊度固定之后,将其作为高精度约束条件,与GPS宽巷观测方程相联立,公式如(4)式所示:
其中:
式中AC,AG分别代表BDS和GPS双差基线向量系数,λG表示GPS宽巷模糊度波长,ΔX表示基线向量,表示GPS(1,-1)组合宽巷模糊度,表示已经固定的BDS(1,-1,0)组合宽巷模糊度,分别代表BDS和GPS观测值矩阵,代表BDS宽巷组合观测值,λC(1,-1,0)代表BDS宽巷组合波长,代表GPS宽巷组合观测值,分别代表双差站星距和双差对流层延迟,代表B1上的电离层延迟值,βC(1,-1,0)、βG(1,-1,0)分别代表BDS和GPS的宽巷电离层延迟因子。
本发明的有益效果为:本发明首先利用北斗三频易于组成长波长、弱电离层、低噪声的组合的优势,单历元可靠固定BDS(0,-1,1)组合和(1,4,-5)组合超宽模糊度,然后利用模糊度“两个超宽巷模糊度固定之后,即可通过整数组合得到任一组合系数为0的宽巷模糊度”这一特性,组合两个超宽巷模糊度得到北斗(1,-1,0)宽巷模糊度,其次利用北斗宽巷模糊度约束GPS宽巷模糊度,得到GPS宽巷模糊度的浮点解和其方差协方差阵,最后采用LAMBDA算法实现了GPS宽巷模糊度的固定;相比于传统MW组合法计算宽巷模糊度时需要多历元平滑,本发明可实现BDS/GPS宽巷模糊度的单历元固定,有效提升了GNSS基线解算的时效性。
附图说明
图1为本发明的方法流程示意图。
图2为本发明基线JNNF-LHNF上北斗超宽巷组合(0,-1,1)模糊度单历元解算偏差示意图。
图3为本发明基线JNNF-LHNF上北斗超宽巷组合(1,4,-5)模糊度单历元解算偏差示意图。
图4为本发明基线LHNF-MAQN上北斗超宽巷组合(0,-1,1)模糊度单历元解算偏差示意图。
图5为本发明基线LHNF-MQAN上北斗超宽巷组合(1,4,-5)模糊度单历元解算偏差示意图。
图6为本发明基线JNNF-LHNF上GPS宽巷模糊度固定Ratio值示意图。
图7为本发明基线LHNF-MAQN上GPS宽巷模糊度固定Ratio值示意图。
具体实施方式
如图1所示,一种基于北斗超宽巷约束的BDS/GPS宽巷模糊度单历元固定方法,包括如下步骤:
(1)利用GNSS接收机接收的BDS三频伪距、载波观测数据组合单历元解算BDS的超宽巷(0,-1,1)和(1,4,-5)的组合模糊度;
(2)组合两个固定后的超宽巷模糊度,计算观测噪声和电离层延迟较小的北斗(1,-1,0)的宽巷组合模糊度;
(3)当BDS宽巷模糊度固定之后,将其作为高精度的约束条件,与接收机接收的GPS双频观测数据组合得到的GPS宽巷观测方程相联立;
(4)对步骤(3)得到的观测模型进行卡尔曼滤波解算,得到GPS双差宽巷模糊度的浮点解和方差协方差矩阵,采用LAMBDA算法进行整数固定,得到GPS双差宽巷模糊度整数解。
一种北斗超宽巷约束的BDS/GPS宽巷模糊度单历元解算方法,包括如下具体步骤:
步骤1),利用载波和伪距组合单历元解算BDS的超宽巷(0,-1,1)组合模糊度,如式(1)所示:
式中,为双差算子,为(0,-1,1)组合模糊度,[·]四舍五入算子,为以米为单位的(0,-1,1)组合的载波观测值,为以米为单位的(0,1,1)组合的伪距观测值,λ(0,-1,1)为(0,-1,1)组合载波观测值波长;
步骤2),利用载波、伪距组合单历元解算BDS的超宽巷(1,4,-5)组合模糊度,如式(2)所示:
式中,为双差算子,为(1,4,-5)组合模糊度,[·]四舍五入算子,为以米为单位的(1,4,-5)组合的载波观测值,为以米为单位的(1,1,0)组合的伪距观测值,λ(1,4,-5)为(1,4,-5)组合载波观测值波长;
步骤3),组合两个固定后的超宽巷模糊度,计算观测噪声和电离层延迟较小的北斗(1,-1,0)宽巷组合模糊度,如式(3)所示。
式中,为双差算子,为(1,-1,0)组合模糊度,为(0,-1,1)组合模糊度,为(1,4,-5)组合模糊度。
步骤4),当BDS宽巷模糊度固定之后,将其作为高精度约束条件,与GPS宽巷观测方程相联系,公式如下:
式中AC,AG分别代表BDS和GPS双差基线向量系数,λG表示GPS宽巷模糊度波长,ΔX表示基线向量,表示GPS(1,-1)组合宽巷模糊度,表示已经固定的BDS(1,-1,0)组合宽巷模糊度,分别代表BDS和GPS观测值矩阵。代表BDS宽巷组合观测值,λC(1,-1,0)代表BDS宽巷组合波长,代表GPS宽巷组合观测值,分别代表双差站星距和双差对流层延迟。代表B1上的电离层延迟值,βC(1,-1,0)、βG(1,-1,0)分别代表BDS和GPS的宽巷电离层延迟因子。
对中长基线而言,利用改正模型可消除绝大部分的双差对流层延迟,式(4)中GPS宽巷模糊度浮点解解算精度主要受载波噪声和电离层延迟影响,GPS宽巷模糊度浮点解精度为:
式中表示GPS宽巷浮点解精度,λG表示GPS宽巷模糊度波长,f1和f2分别表示载波L1和L2上的频率,表示电离层延迟,表示载波噪声。假定载波噪声为0.5cm,双差电离层延迟为0.1m,0.15m,0.2m,则GPS宽巷浮点解精度如下表所示:
表1电离层延迟和载波噪声对GPS宽巷模糊度精度影响
可以看出,当电离层延迟变大的时候,若采用单历元取整GPS模糊度存在一定偏差,采用精确的BDS宽巷模型进行约束,可以准确的计算出基线向量,从而减弱GPS宽巷模糊度之间的相关性,将GPS宽巷模糊度搜索空间压缩到少数几种组合,从而利用LAMBDA算法可以搜索出正确的模糊度。
步骤5),对所述式(4)的观测模型进行卡尔曼滤解算,得到GPS双差宽巷模糊度的浮点解和方差协方差矩阵;然后,采用LAMBDA算法进行整数固定,得到GPS双差宽巷模糊度整数解。
实施例选取:采用南京计量院CORS长度分别为26km和41km两条基线进行模糊度解算实验,数据采样间隔为1s,数据采样时间为2017年7月15号UTC时0:00-1:00;共计3600个历元。
图2-5中,C代表BDS卫星,C后数字代表卫星PRN编号,从图2和图4可以看出,两组基线(0,1,-1)组合模糊度偏差均很小,均在±0.2周之内,因此按照单历元四舍五入直接取整可获取到模糊度的固定解。从图3和图5可以看出两组基线(1,4,-5)组合的模糊度偏差均大于(0,1,-1)组合模糊度偏差,但其偏差仍然在±0.4周之内,因此在解算区间内按照四舍五入原则也能直接取整得到模糊度的固定解。
图6-7代表两组基线采用LAMBDA算法搜索GPS宽巷模糊度时的Ratio值,下表2代表Ratio值的统计结果。
表2 BDS约束GPS宽巷模糊度固定Ratio值统计结果
从图6-7和表2可以看出,对两条中长基线而言,采用BDS宽巷模糊度辅助固定GPS宽巷模糊度,其Ratio值均大于2,大于5的超过97%,因此,采用这种方式实现GPS宽巷模糊度的单历元固定是可行的。
本发明首先利用北斗三频易于组成长波长、弱电离层、低噪声的组合的优势,单历元可靠固定BDS(0,-1,1)组合和(1,4,-5)组合超宽模糊度,然后利用模糊度“两个超宽巷模糊度固定之后,即可通过整数组合得到任一组合系数为0的宽巷模糊度”这一特性,组合两个超宽巷模糊度得到北斗(1,-1,0)宽巷模糊度,其次利用北斗宽巷模糊度约束GPS宽巷模糊度,得到GPS宽巷模糊度的浮点解和其方差协方差阵,最后采用LAMBDA算法实现了GPS宽巷模糊度的固定;相比于传统MW组合法计算宽巷模糊度时需要多历元平滑,本发明可实现BDS/GPS宽巷模糊度的单历元固定,有效提升了GNSS基线解算的时效性。
Claims (4)
1.一种基于北斗超宽巷约束的BDS/GPS宽巷模糊度单历元固定方法,其特征在于,包括如下步骤:
(1)利用GNSS接收机接收的BDS三频伪距、载波观测数据组合单历元解算BDS的超宽巷(0,-1,1)和(1,4,-5)的组合模糊度;
(2)组合两个固定后的超宽巷模糊度,计算观测噪声和电离层延迟较小的北斗(1,-1,0)的宽巷组合模糊度;
(3)当BDS宽巷模糊度固定之后,将其作为高精度的约束条件,与接收机接收的GPS双频观测数据组合得到的GPS宽巷观测方程相联立;
(4)对步骤(3)得到的观测模型进行卡尔曼滤波解算,得到GPS双差宽巷模糊度的浮点解和方差协方差矩阵,采用LAMBDA算法进行整数固定,得到GPS双差宽巷模糊度整数解。
2.如权利要求1所述的基于北斗超宽巷约束的BDS/GPS宽巷模糊度单历元固定方法,其特征在于,步骤(1)中,利用GNSS接收机接收的BDS三频伪距、载波观测数据组合单历元解算BDS的超宽巷(0,-1,1)和(1,4,-5)的组合模糊度具体为:
利用载波和伪距组合单历元解算BDS的超宽巷(0,-1,1)组合模糊度,如式(1)所示:
<mrow>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msub>
</mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,为双差算子,为(0,-1,1)组合模糊度,[·]四舍五入算子,为以米为单位的(0,-1,1)组合的载波观测值,为以米为单位的(0,1,1)组合的伪距观测值,λ(0,-1,1)为(0,-1,1)组合载波观测值波长;
利用载波、伪距组合单历元解算BDS的超宽巷(1,4,-5)组合模糊度,如式(2)所示:
<mrow>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
<mo>,</mo>
<mo>-</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
<mo>,</mo>
<mo>-</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>P</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msub>
</mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
<mo>,</mo>
<mo>-</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</msub>
</mfrac>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,为双差算子,为(1,4,-5)组合模糊度,[·]四舍五入算子,为以米为单位的(1,4,-5)组合的载波观测值,为以米为单位的(1,1,0)组合的伪距观测值,λ(1,4,-5)为(1,4,-5)组合载波观测值波长。
3.如权利要求1所述的基于北斗超宽巷约束的BDS/GPS宽巷模糊度单历元固定方法,其特征在于,步骤(2)中,组合两个固定后的超宽巷模糊度,计算观测噪声和电离层延迟较小的北斗(1,-1,0)的宽巷组合模糊度具体为:
组合两个固定后的超宽巷模糊度,计算观测噪声和电离层延迟较小的北斗(1,-1,0)宽巷组合模糊度,如式(3)所示。
<mrow>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<mn>5</mn>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</msub>
<mo>+</mo>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>N</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mn>4</mn>
<mo>,</mo>
<mo>-</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,为双差算子,为(1,-1,0)组合模糊度,为(0,-1,1)组合模糊度,为(1,4,-5)组合模糊度。
4.如权利要求1所述的基于北斗超宽巷约束的BDS/GPS宽巷模糊度单历元固定方法,其特征在于,步骤(3)中,当BDS宽巷模糊度固定之后,将其作为高精度的约束条件,与接收机接收的GPS双频观测数据组合得到的GPS宽巷观测方程相联立具体为:
当BDS宽巷模糊度固定之后,将其作为高精度约束条件,与GPS宽巷观测方程相联立,公式如(4)式所示:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mi>C</mi>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mi>G</mi>
</msub>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mi>G</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&CenterDot;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mi>X</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>N</mi>
<mi>G</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>L</mi>
<mi>C</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>L</mi>
<mi>G</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>L</mi>
<mi>C</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>L</mi>
<mi>G</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<mi>T</mi>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<mi>I</mi>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>N</mi>
<mi>C</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<mi>T</mi>
<mo>+</mo>
<msub>
<mi>&beta;</mi>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mi>&Delta;</mi>
<mo>&dtri;</mo>
<mi>I</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中AC,AG分别代表BDS和GPS双差基线向量系数,λG表示GPS宽巷模糊度波长,ΔX表示基线向量,表示GPS(1,-1)组合宽巷模糊度,表示已经固定的BDS(1,-1,0)组合宽巷模糊度,分别代表BDS和GPS观测值矩阵,代表BDS宽巷组合观测值,λC(1,-1,0)代表BDS宽巷组合波长,代表GPS宽巷组合观测值, 分别代表双差站星距和双差对流层延迟,代表B1上的电离层延迟值,βC(1,-1,0)、βG(1,-1,0)分别代表BDS和GPS的宽巷电离层延迟因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810082409.0A CN108037521A (zh) | 2018-01-29 | 2018-01-29 | 一种基于北斗超宽巷约束的bds/gps宽巷模糊度单历元固定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810082409.0A CN108037521A (zh) | 2018-01-29 | 2018-01-29 | 一种基于北斗超宽巷约束的bds/gps宽巷模糊度单历元固定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108037521A true CN108037521A (zh) | 2018-05-15 |
Family
ID=62096839
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810082409.0A Pending CN108037521A (zh) | 2018-01-29 | 2018-01-29 | 一种基于北斗超宽巷约束的bds/gps宽巷模糊度单历元固定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108037521A (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085628A (zh) * | 2018-08-27 | 2018-12-25 | 桂林电子科技大学 | 一种整周模糊度的固定方法及系统 |
CN109581455A (zh) * | 2019-01-24 | 2019-04-05 | 东南大学 | 一种bds和gps融合的三频宽巷紧组合定位方法 |
CN109977499A (zh) * | 2019-03-11 | 2019-07-05 | 东南大学 | 基于位置约束的北斗三频/静力水准仪斜拉桥监测方法 |
CN110109166A (zh) * | 2019-04-30 | 2019-08-09 | 东南大学 | 一种快速得到高可靠性卫星定位整数解的方法 |
CN110501735A (zh) * | 2019-07-19 | 2019-11-26 | 中国人民解放军92942部队 | 一种基于北斗三频信号的超窄巷rtk定位方法 |
CN110780323A (zh) * | 2019-11-08 | 2020-02-11 | 厦门理工学院 | 一种长距离下基于北斗三频信号的实时分米级定位方法 |
CN110954932A (zh) * | 2019-12-19 | 2020-04-03 | 中国测绘科学研究院 | 一种三频rtk单历元模糊度快速搜索方法 |
CN111751853A (zh) * | 2020-06-20 | 2020-10-09 | 北京华龙通科技有限公司 | 一种gnss双频载波相位整周模糊度解算方法 |
CN112285745A (zh) * | 2020-11-18 | 2021-01-29 | 青岛杰瑞自动化有限公司 | 基于北斗三号卫星导航系统的三频模糊度固定方法及系统 |
CN112526569A (zh) * | 2021-02-18 | 2021-03-19 | 中国人民解放军国防科技大学 | 一种惯导辅助卫导相对定位多历元逐级模糊度求解方法 |
CN113655510A (zh) * | 2021-08-25 | 2021-11-16 | 广州市城市规划勘测设计研究院 | 一种bds模糊度解算方法、装置、设备和存储介质 |
CN113970772A (zh) * | 2021-10-27 | 2022-01-25 | 东南大学 | 一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法 |
US20220058322A1 (en) * | 2018-09-21 | 2022-02-24 | Trimble Inc. | Correction information integrity monitoring in navigation satellite system positioning methods, systems, and devices |
CN115494535A (zh) * | 2022-09-27 | 2022-12-20 | 中国测绘科学研究院 | 一种北斗三频非差相位模糊度的解算方法及装置 |
CN118011445A (zh) * | 2024-04-08 | 2024-05-10 | 高速铁路建造技术国家工程研究中心 | 一种多gnss天线的模糊度固定方法及系统 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353969A (zh) * | 2011-09-02 | 2012-02-15 | 东南大学 | 精密单点定位技术中相位偏差的估计方法 |
CN103675835A (zh) * | 2013-12-04 | 2014-03-26 | 航天恒星科技有限公司 | 一种北斗三频信号载波相位整周模糊度单历元确定方法 |
CN103760582A (zh) * | 2014-01-03 | 2014-04-30 | 东南大学 | 一种遮挡环境下卫星双差观测结构的优化方法 |
CN104680008A (zh) * | 2015-02-12 | 2015-06-03 | 东南大学 | 一种基于多参考站的网络rtk区域大气误差建模方法 |
CN105158782A (zh) * | 2015-05-29 | 2015-12-16 | 东南大学 | 一种bds和gps观测信息融合的宽巷模糊度解算方法 |
CN105676250A (zh) * | 2016-01-15 | 2016-06-15 | 北京航空航天大学 | 一种基于gnss的单历元三频模糊度解算方法 |
US20170269228A1 (en) * | 2016-03-18 | 2017-09-21 | Deere & Company | Satellite Navigation Receiver With Improved Ambiguity Resolution |
CN107193029A (zh) * | 2017-05-17 | 2017-09-22 | 武汉大学 | 北斗三频信号的网络rtk基准站间模糊度快速确定方法 |
CN107193028A (zh) * | 2017-03-29 | 2017-09-22 | 中国航空无线电电子研究所 | 基于GNSS的Kalman相对定位方法 |
CN107462904A (zh) * | 2017-07-28 | 2017-12-12 | 东南大学 | Gnss和ins融合的高精度gnss终端动态定位检测车及检测方法 |
-
2018
- 2018-01-29 CN CN201810082409.0A patent/CN108037521A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102353969A (zh) * | 2011-09-02 | 2012-02-15 | 东南大学 | 精密单点定位技术中相位偏差的估计方法 |
CN103675835A (zh) * | 2013-12-04 | 2014-03-26 | 航天恒星科技有限公司 | 一种北斗三频信号载波相位整周模糊度单历元确定方法 |
CN103760582A (zh) * | 2014-01-03 | 2014-04-30 | 东南大学 | 一种遮挡环境下卫星双差观测结构的优化方法 |
CN104680008A (zh) * | 2015-02-12 | 2015-06-03 | 东南大学 | 一种基于多参考站的网络rtk区域大气误差建模方法 |
CN105158782A (zh) * | 2015-05-29 | 2015-12-16 | 东南大学 | 一种bds和gps观测信息融合的宽巷模糊度解算方法 |
CN105676250A (zh) * | 2016-01-15 | 2016-06-15 | 北京航空航天大学 | 一种基于gnss的单历元三频模糊度解算方法 |
US20170269228A1 (en) * | 2016-03-18 | 2017-09-21 | Deere & Company | Satellite Navigation Receiver With Improved Ambiguity Resolution |
CN107193028A (zh) * | 2017-03-29 | 2017-09-22 | 中国航空无线电电子研究所 | 基于GNSS的Kalman相对定位方法 |
CN107193029A (zh) * | 2017-05-17 | 2017-09-22 | 武汉大学 | 北斗三频信号的网络rtk基准站间模糊度快速确定方法 |
CN107462904A (zh) * | 2017-07-28 | 2017-12-12 | 东南大学 | Gnss和ins融合的高精度gnss终端动态定位检测车及检测方法 |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085628A (zh) * | 2018-08-27 | 2018-12-25 | 桂林电子科技大学 | 一种整周模糊度的固定方法及系统 |
CN109085628B (zh) * | 2018-08-27 | 2022-09-30 | 桂林电子科技大学 | 一种整周模糊度的固定方法及系统 |
US11709280B2 (en) * | 2018-09-21 | 2023-07-25 | Trimble Inc. | Correction information integrity monitoring in navigation satellite system positioning methods, systems, and devices |
US20220058322A1 (en) * | 2018-09-21 | 2022-02-24 | Trimble Inc. | Correction information integrity monitoring in navigation satellite system positioning methods, systems, and devices |
CN109581455A (zh) * | 2019-01-24 | 2019-04-05 | 东南大学 | 一种bds和gps融合的三频宽巷紧组合定位方法 |
CN109977499A (zh) * | 2019-03-11 | 2019-07-05 | 东南大学 | 基于位置约束的北斗三频/静力水准仪斜拉桥监测方法 |
CN109977499B (zh) * | 2019-03-11 | 2023-04-07 | 东南大学 | 基于位置约束的北斗三频/静力水准仪斜拉桥监测方法 |
CN110109166A (zh) * | 2019-04-30 | 2019-08-09 | 东南大学 | 一种快速得到高可靠性卫星定位整数解的方法 |
US11686860B2 (en) | 2019-04-30 | 2023-06-27 | Southeast University | Method for quickly acquiring highly reliable integer solution for satellite positioning |
CN110109166B (zh) * | 2019-04-30 | 2020-06-09 | 东南大学 | 一种快速得到高可靠性卫星定位整数解的方法 |
WO2020220579A1 (zh) * | 2019-04-30 | 2020-11-05 | 东南大学 | 一种快速得到高可靠性卫星定位整数解的方法 |
CN110501735A (zh) * | 2019-07-19 | 2019-11-26 | 中国人民解放军92942部队 | 一种基于北斗三频信号的超窄巷rtk定位方法 |
CN110780323B (zh) * | 2019-11-08 | 2021-10-15 | 厦门理工学院 | 一种长距离下基于北斗三频信号的实时分米级定位方法 |
CN110780323A (zh) * | 2019-11-08 | 2020-02-11 | 厦门理工学院 | 一种长距离下基于北斗三频信号的实时分米级定位方法 |
CN110954932A (zh) * | 2019-12-19 | 2020-04-03 | 中国测绘科学研究院 | 一种三频rtk单历元模糊度快速搜索方法 |
CN110954932B (zh) * | 2019-12-19 | 2021-08-03 | 中国测绘科学研究院 | 一种三频rtk单历元模糊度快速搜索方法 |
CN111751853A (zh) * | 2020-06-20 | 2020-10-09 | 北京华龙通科技有限公司 | 一种gnss双频载波相位整周模糊度解算方法 |
CN111751853B (zh) * | 2020-06-20 | 2023-10-03 | 北京华龙通科技有限公司 | 一种gnss双频载波相位整周模糊度解算方法 |
CN112285745B (zh) * | 2020-11-18 | 2022-09-23 | 青岛杰瑞自动化有限公司 | 基于北斗三号卫星导航系统的三频模糊度固定方法及系统 |
CN112285745A (zh) * | 2020-11-18 | 2021-01-29 | 青岛杰瑞自动化有限公司 | 基于北斗三号卫星导航系统的三频模糊度固定方法及系统 |
CN112526569A (zh) * | 2021-02-18 | 2021-03-19 | 中国人民解放军国防科技大学 | 一种惯导辅助卫导相对定位多历元逐级模糊度求解方法 |
CN112526569B (zh) * | 2021-02-18 | 2021-05-07 | 中国人民解放军国防科技大学 | 一种惯导辅助卫导相对定位多历元逐级模糊度求解方法 |
CN113655510B (zh) * | 2021-08-25 | 2023-02-24 | 广州市城市规划勘测设计研究院 | 一种bds模糊度解算方法、装置、设备和存储介质 |
CN113655510A (zh) * | 2021-08-25 | 2021-11-16 | 广州市城市规划勘测设计研究院 | 一种bds模糊度解算方法、装置、设备和存储介质 |
CN113970772A (zh) * | 2021-10-27 | 2022-01-25 | 东南大学 | 一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法 |
CN113970772B (zh) * | 2021-10-27 | 2024-04-12 | 东南大学 | 一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法 |
CN115494535A (zh) * | 2022-09-27 | 2022-12-20 | 中国测绘科学研究院 | 一种北斗三频非差相位模糊度的解算方法及装置 |
CN115494535B (zh) * | 2022-09-27 | 2023-05-23 | 中国测绘科学研究院 | 一种北斗三频非差相位模糊度的解算方法及装置 |
CN118011445A (zh) * | 2024-04-08 | 2024-05-10 | 高速铁路建造技术国家工程研究中心 | 一种多gnss天线的模糊度固定方法及系统 |
CN118011445B (zh) * | 2024-04-08 | 2024-06-25 | 高速铁路建造技术国家工程研究中心 | 一种多gnss天线的模糊度固定方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108037521A (zh) | 一种基于北斗超宽巷约束的bds/gps宽巷模糊度单历元固定方法 | |
CN109581452B (zh) | 一种gnss参考站载波相位整周模糊度解算方法 | |
CN108445518B (zh) | 一种基于双差模糊度固定解约束的gnss精密时间传递方法 | |
WO2019174113A1 (zh) | 一种gps/bds紧组合载波差分定位方法 | |
CN105607089B (zh) | 一种基于伪距辅助载波相位的北斗共视时间传递的方法 | |
CN105158782B (zh) | 一种bds和gps观测信息融合的宽巷模糊度解算方法 | |
CN108828642B (zh) | 一种ins辅助bds单频接收机的模糊度快速解算方法 | |
CN109765589B (zh) | 一种基于无电离层组合的三频gnss实时周跳固定技术 | |
CN108363084A (zh) | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 | |
Li et al. | Benefits of BDS-3 B1C/B1I/B2a triple-frequency signals on precise positioning and ambiguity resolution | |
Li et al. | Principle and performance of multi-frequency and multi-GNSS PPP-RTK | |
CN111638535B (zh) | 一种用于gnss实时精密单点定位的混合模糊度固定方法 | |
CN109581455B (zh) | 一种bds和gps融合的三频宽巷紧组合定位方法 | |
CN106569242A (zh) | 固定参考星的gnss单差处理方法 | |
CN101089650A (zh) | 基于地基增强系统的误差分离方法及地基增强系统 | |
CN105929430B (zh) | 一种gnss零基线参考站间模糊度快速固定方法 | |
CN105974440B (zh) | 一种获取观测噪声的方法和装置 | |
CN110806594B (zh) | 一种基于5g通信的车辆间协同导航定位方法 | |
CN106896386A (zh) | Glonass频间偏差精确估计方法 | |
CN116299615B (zh) | 一种实现单北斗实时ppp模糊固定的相位偏差估计方法 | |
CN116148909B (zh) | 多频多模非组合精密单点定位瞬时窄巷模糊度固定方法 | |
CN106125113A (zh) | 一种利用多系统gnss观测值的高精度基线解算方法 | |
Gao et al. | Method and assessment of BDS triple-frequency ambiguity resolution for long-baseline network RTK | |
CN116540303A (zh) | 一种历元间差分观测方程建立方法和地震位移计算方法 | |
CN107121689A (zh) | Glonass频间偏差单历元快速估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180515 |