CN110095796B - 一种跨系统mw紧组合实时动态精密导航定位方法 - Google Patents
一种跨系统mw紧组合实时动态精密导航定位方法 Download PDFInfo
- Publication number
- CN110095796B CN110095796B CN201910387854.2A CN201910387854A CN110095796B CN 110095796 B CN110095796 B CN 110095796B CN 201910387854 A CN201910387854 A CN 201910387854A CN 110095796 B CN110095796 B CN 110095796B
- Authority
- CN
- China
- Prior art keywords
- satellite
- observed quantity
- wide lane
- difference
- frequency
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 239000005433 ionosphere Substances 0.000 claims description 5
- 238000009795 derivation Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000010521 absorption reaction Methods 0.000 claims 1
- 230000010354 integration Effects 0.000 claims 1
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 239000005436 troposphere Substances 0.000 description 2
- 125000004432 carbon atom Chemical group C* 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 125000001153 fluoro group Chemical group F* 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
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
本发明涉及一种跨系统MW紧组合实时动态精密导航定位方法,属于卫星导航与定位领域。包括获取基站与移动站接收机输出的原始伪距与载波相位观测量;构造MW紧组合,得到MW组合下观测量LMW;求卫星q和卫星s宽巷单差观测量、宽巷双差观测量、跨系统双差观测量;计算MW宽巷ISBMW;将MW宽巷组合模型代替双频伪距观测量,得到跨系统固定双频模糊度方程。本发明在双模双频紧组合的基础上引入宽巷模型的形式,保留了双模双频紧组合在收星受限情况下实现模糊度固定的优势,并且能获得更准确的载波相位ISB估计值,同时,跨系统MW紧组合利用伪距和载波相位组合构成的MW宽巷观测量替代双频伪距观测量,有利于降低直接引用伪距带来的误差风险,提高紧组合定位性能。
Description
技术领域
本发明涉及一种跨系统MW紧组合实时动态精密导航定位方法,属于卫星导航与定位领域。
背景技术
模糊度固定是基于载波相位差分的实时动态方法获得高精度与高可靠性定位的关键技术之一。双模MW松组合是根据载波相位和伪距观测量特点,构造宽巷组合模型,获得更长波长和弱电离层特点的新观测量,有利于辅助提高模糊度固定率的方法。该方法需要满足收星数至少5颗,或者2颗北斗卫星加上2颗GPS卫星的基本条件。在恶劣的环境下,可用卫星数相对较少,如何在收星受限的情况下固定模糊度是提升模糊度固定成功率的有效途径之一。
为了解决在收星受限情况下难以固定模糊度的问题,双模双频紧组合是一种可行的解决方案。该组合方式可以在只有4颗卫星甚至3颗卫星情况下实现模糊度固定,并且在收星不足情况下,有效地提高了模糊度解算模型观测量的冗余度,可以实现比松组合更高的定位性能。另一方面,由于双模双频紧组合模型直接应用了伪距观测量,而伪距具有噪声较大、精度不足等问题,影响紧组合模糊度固定。为了进一步提高双模双频紧组合定位性能,有必要引入MW宽巷组合模型。首先,MW宽巷组合模型在收星受限情况下实现模糊度固定的能力与双模双频紧组合相当,优于松组合模型。同时,由于GPS和北斗的MW组合波长相差并不大,和L2频段的波长差占比相同,对于控制单差模糊度偏差以及统计宽巷ISB有利,可以帮助固定窄巷紧组合模糊度。其次,可以获得更准确的载波相位ISB估计值,进而提高定位性能。最后,用伪距和载波相位组合构成的MW宽巷观测量,替代双频伪距观测量,有利于降低直接引用伪距带来的误差风险,提高紧组合定位性能。因此,本文研究一种跨系统MW紧组合方法,进一步提高紧组合定位性能。
发明内容
本发明的目的是为了提高定位性能而提供一种跨系统MW紧组合实时动态精密导航定位方法。
本发明的目的是这样实现的:一种跨系统MW紧组合实时动态精密导航定位方法,具体包括以下步骤:
步骤1、获取基站与移动站接收机输出的原始伪距与载波相位观测量;
步骤2、利用步骤1中的原始观测量构造MW紧组合,得到MW组合下观测量LMW;
步骤3、通过站间差分的方式得到卫星q和卫星s宽巷单差观测量;
步骤4、判断卫星q和卫星s是否处于同一系统,若是,则设置参考星为s,构造宽巷双差观测量,转化为同系统内标准双差的形式,然后进行步骤7;否则进行步骤5;
步骤5、若卫星q和卫星s不是同一系统的卫星,引入新约束条件进一步转化,得到跨系统双差观测量;
步骤6、引入步骤5中得到的跨系统双差观测量,计算MW宽巷ISBMW;
步骤7、将步骤4、5、6中得到的MW宽巷组合模型代替双频伪距观测量,得到跨系统固定双频模糊度方程。
本发明还包括这样一些结构特征:
1、步骤1所述原始伪距与载波相位观测量分别为:
式中,上标s和Q分别表示卫星的PRN号和卫星系统,Q表示为G、C分别对应GPS、BDS;下标r和j(j=1,2)分别表示接收机ID和观测值频段号;与分别表示伪距与载波相位观测值减计算值的残差量;xr表示用户位置、对流层等非弥散项增量;表示对应xr的线性化几何设计矩阵;tr与ts分别表示接收机与卫星端钟差;对应频点B1的一阶电离层延迟误差,表示电离层比例系数;与分别表示与频点相关的接收机与卫星端的伪距偏差;表示载波相位波长;表示整周模糊度;与分别表示与频点相关的接收机与卫星端载波相位偏差,包含载波相位硬件偏差与初始相位等;与分别表示伪距与载波相位观测噪声。
2、所述步骤2具体包括以下步骤:
步骤2.1、将L1频段和L2频段的载波相位改为:
式中,L1和L2分别是L1频段和L2频段的载波相位,f1和f2分别是L1频段和L2频段的频率值,将伪距电离层延迟改为:
步骤2.2、构造MW宽巷组合模型,获得新的观测量LMW:
将上述四个式子都带入下式,可化简为将几何距离和电离层消除的形式,
步骤2.3、考虑接收机对双频载波相位和伪距的系统延时偏差ISB,将上式改为:
LMW=λMW·NMW+c·ιMW+εMW。
3、所述步骤5具体包括以下步骤:
步骤5.1、卫星q和卫星s不是同一系统的卫星,设s为GPS,q为北斗,进一步转化,
步骤5.2、采用北斗1号卫星和GPS参考星之间的模糊度固定为0的形式引入新约束条件,将上式转为:
4、所述步骤6具体包括以下步骤:
步骤6.3、对N颗北斗卫星计算得到的ISB平滑得到ISBMW;
5、所述步骤7具体包括将跨系统MW宽巷组合模型代替双频伪距观测量,得到跨系统固定双频模糊度方程:
式中,和ISBMW分别为L1频段的载波相位ISB、L2频段的载波相位ISB以及MW宽巷的ISB,和分别为GPS与GPS参考星的双差模糊度矩阵,和分别为北斗与GPS参考星的双差模糊度矩阵;宽巷观测量由载波相位参与构建,通过宽巷替代伪距观测量进行MW紧组合可提高模糊度固定成功率。
与现有技术相比,本发明的有益效果是:
本发明公开了一种跨系统MW紧组合实时动态精密导航定位方法,该方法在双模双频紧组合的基础上引入MW宽巷模型。保留了双模双频紧组合在收星受限情况下实现模糊度固定的优势,并且使用MW宽巷观测量替代双频伪距观测量,降低直接引用伪距带来的误差风险,提高紧组合定位性能。因此,这种新的方法具有更好的模糊度解算的性能是可以被预见的。
附图说明
图1是跨系统MW紧组合实时动态精密导航定位方法流程图。
具体实施方式
下面结合附图与具体实施方式对本发明作进一步详细描述。
本发明公开了一种跨系统MW紧组合实时动态精密导航定位方法,属于卫星精密导航与定位领域。在双模双频紧组合的基础上引入宽巷模型的形式,保留了双模双频紧组合在收星受限情况下实现模糊度固定的优势,并且可以获得更准确的载波相位ISB估计值,进而提高定位性能。同时,跨系统MW紧组合利用伪距和载波相位组合构成的MW宽巷观测量替代双频伪距观测量,有利于降低直接引用伪距带来的误差风险,提高紧组合定位性能。因此,本文研究的跨系统MW紧组合可以在收星受限时实现模糊度固定,并提高紧组合定位性能。
如附图1所示,是跨系统MW紧组合实时动态精密导航定位方法流程图。下面以北斗系统及GPS系统的双频信号即:B1=1561.098MHZ,B2=1207.14MHZ,L1=1575.42MHZ,L2=1227.60MHZ作为实施例,并结合附对本发明的技术方案作进一步具体的说明。
步骤1、获取基站与移动站接收机输出的原始伪距与载波相位观测量;
式中,上标s和Q分别表示卫星的PRN号和卫星系统,Q可表示为G、C分别对应GPS、BDS;下标r和j(j=1,2)分别表示接收机ID和观测值频段号;与分别表示伪距与载波相位观测值减计算值的残差量;xr表示用户位置、对流层等非弥散项增量;表示对应xr的线性化几何设计矩阵;tr与ts分别表示接收机与卫星端钟差;对应频点B1的一阶电离层延迟误差,表示电离层比例系数;与分别表示与频点相关的接收机与卫星端的伪距偏差;表示载波相位波长;表示整周模糊度;与分别表示与频点相关的接收机与卫星端载波相位偏差,包含载波相位硬件偏差与初始相位等;与分别表示伪距与载波相位观测噪声。
步骤2、利用步骤1中的原始观测量构造MW紧组合,得到MW组合下观测量LMW;
对于北斗或GPS,将L1频段和L2频段的载波相位改为:
式中,L1和L2分别是L1频段和L2频段的载波相位,f1和f2分别是L1频段和L2频段的频率值。类似地,将伪距电离层延迟改为:
构造MW宽巷组合模型,获得新的观测量LMW:
将式(2)、式(3)、式(4)和式(5)都带入式(6),化简后可以将几何距离和电离层消除的形式,
必须考虑接收机对双频载波相位和伪距的系统延时偏差ISB,将式(7)改为:
LMW=λMW·NMW+c·ιMW+εMW (8)
步骤3、利通过站间差分的方式得到卫星q和卫星s宽巷单差观测量;
步骤4、设置参考星为s,构造宽巷双差观测量,若卫星q和卫星s是同一系统的卫星,转化为同系统内标准双差的形式,否则进行步骤5;
进一步,设参考星为s,构造宽巷双差观测量:
步骤5、若卫星q和卫星s不是同一系统的卫星,引入新约束条件进一步转化,得到跨系统双差观测量;
若卫星q和卫星s不是同一系统的卫星,设s为GPS,q为北斗,进一步转化,
式中,是北斗宽巷波长,是GPS宽巷波长,是GPS的宽巷单差模糊度,为GPS和北斗MW组合的波长差。由于系统的不同,不可以直接计算出宽巷混合双差模糊度,采用北斗1号卫星和GPS参考星之间的模糊度固定为0的形式来引入新约束条件,将上式转为:
步骤6、引入步骤5中得到的跨系统双差观测量计算ISBMW;
根据式(13),在单历元中,若北斗卫星j和北斗1号卫星可以成功固定模糊度,则直接引入跨系统宽巷观测量计算ISBMW:
步骤7、将步骤4、5、6中得到的MW宽巷组合模型代替双频伪距观测量,得到跨系统固定双频模糊度方程;
根据式(13),将跨系统MW宽巷组合模型代替双频伪距观测量,得到跨系统固定双频模糊度方程:
本发明的步骤2中构造MW紧组合时,必须考虑接收机对双频载波相位和伪距的系统延时偏差ISB;步骤5利用采用北斗1号卫星和GPS参考星之间的模糊度固定为0的形式来引入新约束条件,可以将MW宽巷组合跨系统模糊度转为北斗和北斗1号卫星的形式,跨系统产生的系统偏差被ISBMW吸收。步骤6中利用北斗卫星双差模糊度辅助,在静态或动态中都可以直接准确、快速地计算ISBMW。步骤7利用步骤4、5、6中得到的MW宽巷组合模型代替双频伪距观测量,得到跨系统固定双频模糊度方程,宽巷观测量是由载波相位参与构建的,通过宽巷替代伪距观测量进行MW紧组合可以进一步提高模糊度固定成功率。
Claims (5)
1.一种跨系统MW紧组合实时动态精密导航定位方法,其特征在于,具体包括以下步骤:
步骤1、获取基站与移动站接收机输出的原始伪距与载波相位观测量;
步骤2、利用步骤1中的原始观测量构造MW宽巷紧组合,得到MW宽巷组合下观测量LMW;
步骤3、通过站间差分的方式得到卫星q和卫星s的MW宽巷单差观测量;
步骤4、判断卫星q和卫星s是否处于同一系统,若是,则设置参考星为s,构造同系统MW宽巷双差观测量,然后进行步骤7;否则进行步骤5;
步骤5、若卫星q和卫星s不是同一系统的卫星,构造跨系统MW宽巷双差观测量;
步骤5.1、卫星q和卫星s不是同一系统的卫星,设s为GPS,q为北斗,进一步转化,
步骤5.2、采用北斗1号卫星和GPS参考星之间的模糊度固定为0的形式引入新约束条件,将上式转为:
步骤6、引入步骤5中得到的跨系统MW宽巷双差观测量,计算MW宽巷ISBMW;
步骤7、将前述步骤中得到的MW宽巷双差观测量代替双频伪距观测量,得到跨系统MW紧组合方程。
2.根据权利要求1所述一种跨系统MW紧组合实时动态精密导航定位方法,其特征在于,步骤1所述原始伪距与载波相位观测量分别为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910387854.2A CN110095796B (zh) | 2019-05-10 | 2019-05-10 | 一种跨系统mw紧组合实时动态精密导航定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910387854.2A CN110095796B (zh) | 2019-05-10 | 2019-05-10 | 一种跨系统mw紧组合实时动态精密导航定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110095796A CN110095796A (zh) | 2019-08-06 |
CN110095796B true CN110095796B (zh) | 2023-03-31 |
Family
ID=67447655
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910387854.2A Active CN110095796B (zh) | 2019-05-10 | 2019-05-10 | 一种跨系统mw紧组合实时动态精密导航定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110095796B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111308528B (zh) * | 2019-12-10 | 2023-02-14 | 哈尔滨工程大学 | 一种北斗/gps紧组合虚拟参考站定位方法 |
CN111239786B (zh) * | 2020-02-28 | 2021-04-13 | 同济大学 | 一种无人驾驶定位测姿的载波相位整周模糊度测定方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9733359B2 (en) * | 2008-01-14 | 2017-08-15 | Trimble Inc. | GNSS signal processing with regional augmentation message |
CN101403792B (zh) * | 2008-11-05 | 2011-06-29 | 武汉大学 | 动态模糊度确定方法 |
US9170335B2 (en) * | 2011-02-14 | 2015-10-27 | Trimble Navigation Limited | GNSS signal processing with ionosphere model for synthetic reference data |
CN102621569A (zh) * | 2012-03-22 | 2012-08-01 | 哈尔滨工程大学 | 一种分布滤波式全球定位及捷联惯导组合导航方法 |
CN103323867B (zh) * | 2013-07-05 | 2015-04-29 | 武汉大学 | 一种北斗二代geo卫星精密定轨系统偏差改进方法 |
CN105204048B (zh) * | 2014-06-24 | 2017-06-23 | 东南大学 | 一种北斗兼容gps/glonass的网络rtk参考站间模糊度快速解算方法 |
CN105158782B (zh) * | 2015-05-29 | 2017-07-28 | 东南大学 | 一种bds和gps观测信息融合的宽巷模糊度解算方法 |
CN105445771A (zh) * | 2015-11-13 | 2016-03-30 | 上海华测导航技术股份有限公司 | 一种单频rtk融合测试分析方法 |
-
2019
- 2019-05-10 CN CN201910387854.2A patent/CN110095796B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN110095796A (zh) | 2019-08-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109764879B (zh) | 一种卫星定轨方法、装置及电子设备 | |
CN109581452B (zh) | 一种gnss参考站载波相位整周模糊度解算方法 | |
CN108363084B (zh) | 利用卫星定位的方法和装置、卫星导航接收机、存储介质 | |
CN111190203B (zh) | 基于非差观测模型和消秩亏的ppp-rtk定位方法 | |
US10802160B2 (en) | Rapid determination of precise position by aiding data | |
CN104102822B (zh) | 一种多频gnss观测值随机特性建模方法 | |
EP2156208B1 (en) | Positioning using a reference station | |
CN108519614A (zh) | 一种gps/bds紧组合载波差分定位方法 | |
CN109001781B (zh) | 一种顾及电离层约束的bds三频模糊度解算方法 | |
Su et al. | Triple-frequency carrier phase precise time and frequency transfer models for BDS-3 | |
CN109884679B (zh) | 一种单模gnss系统跨频点混合双差rtk解算方法 | |
CN109581455B (zh) | 一种bds和gps融合的三频宽巷紧组合定位方法 | |
Wang et al. | FCB estimation with three different PPP models: Equivalence analysis and experiment tests | |
CN116299615B (zh) | 一种实现单北斗实时ppp模糊固定的相位偏差估计方法 | |
CN115933356B (zh) | 一种虚拟原子钟的高精度时间同步系统和方法 | |
CN104898144B (zh) | 一种北斗系统参考站整周模糊度单历元确定方法 | |
CN110095796B (zh) | 一种跨系统mw紧组合实时动态精密导航定位方法 | |
CN116148909B (zh) | 多频多模非组合精密单点定位瞬时窄巷模糊度固定方法 | |
US20240019585A1 (en) | Leo-augmentation-based convergence time shortening method of wide-area uduc ppp-rtk positioning | |
CN109683182A (zh) | 一种引入系统间偏差的gnss多模组合差分定位方法 | |
CN116359968B (zh) | 一种联合北斗二号和北斗三号的三频差分定位方法 | |
CN116540280B (zh) | 多频卫星导航数据状态域改正信息综合处理方法和系统 | |
CN117665871A (zh) | 一种载波距支持下的gnss非差固定解精密钟差确定方法 | |
CN116679329A (zh) | 基于状态域空间的北斗时空基准精确表征方法和系统 | |
Li et al. | Assessment and analysis of the four-satellite QZSS precise point positioning and the integrated data processing with GPS |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |