CN116359968A - 一种联合北斗二号和北斗三号的三频差分定位方法 - Google Patents

一种联合北斗二号和北斗三号的三频差分定位方法 Download PDF

Info

Publication number
CN116359968A
CN116359968A CN202310517617.XA CN202310517617A CN116359968A CN 116359968 A CN116359968 A CN 116359968A CN 202310517617 A CN202310517617 A CN 202310517617A CN 116359968 A CN116359968 A CN 116359968A
Authority
CN
China
Prior art keywords
beidou
difference
double
wide lane
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.)
Granted
Application number
CN202310517617.XA
Other languages
English (en)
Other versions
CN116359968B (zh
Inventor
邵慧超
郭向欣
叶瑾
涂振
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Leador Spatial Information Technology Co ltd
Original Assignee
Leador Spatial Information Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Leador Spatial Information Technology Co ltd filed Critical Leador Spatial Information Technology Co ltd
Priority to CN202310517617.XA priority Critical patent/CN116359968B/zh
Publication of CN116359968A publication Critical patent/CN116359968A/zh
Application granted granted Critical
Publication of CN116359968B publication Critical patent/CN116359968B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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/421—Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

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

本发明公开了一种联合北斗二号和北斗三号的三频差分定位方法。该方法选择相同中心频点的北斗二号系统（B1I,B2I,B3I）三频信号原始测量值和北斗三号系统（B1I,B2b,B3I）三频信号原始测量值进行联合组合，依次求得三频单差载波相位测量值、三频双差载波相位测量值，利用超宽巷、宽巷、窄巷线性组合以及LAMBDA算法求解整周模糊度，根据最小二乘算法进行三频差分基线解算得到高精度定位结果。相对单独北斗二号或者单独北斗三号的差分定位方法，本发明提出的方法将北斗二号和北斗三号联合组成线性组合，方法简单、计算量较小，同时提高了北斗卫星导航系统观测数据的利用率，也显著增加了差分定位的固定率，进而提升了三频差分定位的精度。

Description

一种联合北斗二号和北斗三号的三频差分定位方法
技术领域
本发明属于卫星导航定位技术领域，具体涉及一种联合北斗二号和北斗三号的三频差分定位方法。
背景技术
目前，北斗二号提供B1I、B2I和B3I三个公开服务信号，北斗三号提供B1I、B1C、B2a、B2b和B3I五个公开服务信号。三频信号组合观测具有长波长、弱电离层延迟、观测噪声小等特点，有利于模糊度的准确快速固定，提升卫星导航差分定位精度。北斗三号五个频点信号为三频信号组合提供更多的组合选择，再结合已有的北斗二号的三频信号，为提高导航定位的精度带来了更优选择的可能性，但也带来了可选组合模式过多、计算量大、计算模型复杂等问题。此外，针对单卫星系统进行差分定位，当同步观测卫星颗数少时，观测模型冗余不足，导致固定率降低。
发明内容
针对现有技术不足，本发明提供一种联合北斗二号和北斗三号的三频差分定位方法，包括以下步骤：
步骤1，获取北斗二号、北斗三号三频信号原始测量值；
步骤2，计算单差载波相位测量值；
步骤3，基于单差载波相位测量值，计算双差载波相位测量值；
步骤4，基于双差载波相位测量值和三频信号线性组合，建立超宽巷线性组合，并进行超宽巷模糊度固定解算；
步骤5，建立宽巷线性组合，并基于超宽巷模糊度固定解后的双差伪距组合测量值进行宽巷模糊度固定解算；
步骤6，建立窄巷线性组合，基于宽巷模糊度固定解后的双差伪距组合测量值，利用LAMBDA算法解算窄巷模糊度；
步骤7，基于最小二乘算法求解最优基线向量校正量，得到高精度定位结果。
而且，所述步骤2中流动站接收机
Figure SMS_1
和基站接收机/>
Figure SMS_2
之间对卫星/>
Figure SMS_3
的单差载波相位测量值/>
Figure SMS_4
为：
Figure SMS_5
（1）
式中，
Figure SMS_9
表示波长，/>
Figure SMS_10
表示流动站/>
Figure SMS_13
和基站/>
Figure SMS_7
之间的几何距离，/>
Figure SMS_12
表示电离层延迟，/>
Figure SMS_14
表示对流层延迟，/>
Figure SMS_16
表示载波频率，/>
Figure SMS_6
表示接收机钟差，/>
Figure SMS_11
表示流动站/>
Figure SMS_15
和基站
Figure SMS_17
的单差整周模糊度，/>
Figure SMS_8
表示单差量测噪声。
而且，所述步骤3中基于公式（1）单差载波相位测量值方程，在流动站接收机
Figure SMS_18
和基站接收机/>
Figure SMS_19
之间对卫星/>
Figure SMS_20
和卫星/>
Figure SMS_21
之间的单差载波相位测量值进行差分，由于
Figure SMS_22
，因此得到双差载波相位测量值/>
Figure SMS_23
：
Figure SMS_24
（2）
式中，
Figure SMS_26
表示波长，/>
Figure SMS_31
、/>
Figure SMS_34
分别为基线向量在基站对卫星/>
Figure SMS_27
和卫星/>
Figure SMS_28
观测的投影方向，/>
Figure SMS_29
为基线向量，/>
Figure SMS_32
表示流动站接收机/>
Figure SMS_25
和基站接收机/>
Figure SMS_30
的双差整周模糊度，/>
Figure SMS_33
表示双差量测噪声。
而且，所述步骤4中基于公式（2）双差载波相位测量值方程，针对北斗二号系统（B1I, B2I, B3I）三频信号和北斗三号系统（B1I, B2b, B3I）三频信号，得到某一时刻的三频双差载波相位测量值
Figure SMS_35
分别为：
Figure SMS_36
（3）
Figure SMS_37
（4）
Figure SMS_38
（5）
式中，下标“1”“2”“3”分别表示B1、B2、B3三个载波频率，
Figure SMS_39
表示波长，/>
Figure SMS_40
表示几何距离，/>
Figure SMS_41
表示卫星星历误差，/>
Figure SMS_42
表示对流层延迟，/>
Figure SMS_43
表示电离层延迟，/>
Figure SMS_44
表示整周模糊度，/>
Figure SMS_45
表示量测噪声。
基于公式（3）-公式（5），对三频双差载波测量值进行三频信号线性组合
Figure SMS_46
的表达式为：
Figure SMS_47
（6）
式中，系数
Figure SMS_48
既可以是整数，也可以是分数，将组合系数记成/>
Figure SMS_49
。
基于公式（6）设置组合系数(0, -1, 1)为超宽巷线性组合，由于超宽巷为长波长观测值，中短基线下的组合方程残余误差小于该波长，因此对超宽巷模糊度取整能够得到超宽巷模糊度固定解，即：
Figure SMS_50
（7）
式中，
Figure SMS_51
表示超宽巷模糊度，下标/>
Figure SMS_52
表示超宽巷组合，/>
Figure SMS_53
是超宽巷的波长，/>
Figure SMS_54
是超宽巷双差载波组合测量值，/>
Figure SMS_55
是超宽巷双差伪距组合测量值。
而且，所述步骤5中基于公式（6）设置组合系数(1, -1, 0)为宽巷线性组合，联合公式（7）超宽巷模糊度固定解后的修正双差伪距组合测量值，对宽巷模糊度取整得到宽巷模糊度固定解，则有：
Figure SMS_56
（8）
式中，
Figure SMS_57
表示宽巷模糊度，下标/>
Figure SMS_58
表示宽巷组合，/>
Figure SMS_59
是宽巷的波长，/>
Figure SMS_60
是宽巷双差载波组合测量值，/>
Figure SMS_61
是超宽巷修正双差伪距组合观测值。
而且，所述步骤6中为了得到高精度厘米级的定位精度，基于公式（6）设置组合系数(1, 1, 0)为窄巷线性组合，联合公式（8）宽巷模糊度固定解后的修正双差伪距组合测量值，则有：
Figure SMS_62
（9）
式中，
Figure SMS_63
表示窄巷模糊度，下标/>
Figure SMS_64
表示窄巷组合，/>
Figure SMS_65
是窄巷的波长，/>
Figure SMS_66
是窄巷双差载波组合测量值，/>
Figure SMS_67
是宽巷修正双差伪距组合观测值。
由于窄巷波长小，组合观测噪声小，直接取整窄巷模糊度容易产生错误，因此采用最小二乘模糊度降相关算法LAMBDA解算窄巷模糊度。
而且，所述步骤7中联合公式（9），双差载波相位的相对定位方程的矩阵形式为：
Figure SMS_68
（10）
式中，
Figure SMS_69
是双差载波相位测量值，/>
Figure SMS_70
是基线向量校正量，/>
Figure SMS_71
和/>
Figure SMS_72
是常系数矩阵，
Figure SMS_73
是窄巷模糊度。
基于最小二乘原理求得最优解
Figure SMS_74
使得测量残差的加权平方和最小，得到基线向量校正量/>
Figure SMS_75
，进而得到最优基线，再结合已知基站坐标，即可得到最终的高精度定位结果。
与现有技术相比，本发明具有如下优点：
相对单独北斗二号或者单独北斗三号的差分定位方法，本发明提出的方法将北斗二号和北斗三号联合组成线性组合，方法简单、计算量较小，同时提高了北斗卫星导航系统观测数据的利用率，也显著增加了差分定位的固定率，可达99.8%，进而提升了三频差分定位的精度。
附图说明
图1为本发明实施例的流程图。
具体实施方式
本发明提供一种联合北斗二号和北斗三号的三频差分定位方法，选择相同中心频点的北斗二号系统（B1I, B2I, B3I）三频信号原始测量值和北斗三号系统（B1I, B2b,B3I）三频信号原始测量值进行联合组合，依次求得三频单差载波相位测量值、三频双差载波相位测量值，利用超宽巷、宽巷、窄巷线性组合以及LAMBDA算法求解整周模糊度，根据最小二乘算法进行三频差分基线解算得到高精度定位结果。
下面结合附图和实施例对本发明的技术方案作进一步说明。
如图1所示，本发明实施例的流程包括以下步骤：
步骤1，获取北斗二号、北斗三号三频信号原始测量值。
北斗二号提供B1I、B2I和B3I三个公开服务信号，其中B1I中心频率为1561.098MHz，B2I中心频率为1207.14 MHz，B3I中心频率为1268.52 MHz。北斗三号提供B1I、B1C、B2a、B2b和B3I五个公开服务信号，其中B1C中心频率为1575.42 MHz，B2a中心频率为1176.45 MHz，B2b中心频率为1207.14 MHz。利用接收机获取每颗观测卫星的北斗二号系统（B1I, B2I, B3I）三频信号原始测量值和北斗三号系统（B1I, B2b, B3I）三频信号原始测量值的载波相位、伪距等原始测量值。
步骤2，计算单差载波相位测量值。
流动站接收机
Figure SMS_76
和基站接收机/>
Figure SMS_77
之间对卫星/>
Figure SMS_78
的单差载波相位测量值/>
Figure SMS_79
为：
Figure SMS_80
（1）
式中，
Figure SMS_83
表示波长，/>
Figure SMS_85
表示流动站/>
Figure SMS_90
和基站/>
Figure SMS_84
之间的几何距离，/>
Figure SMS_87
表示电离层延迟，/>
Figure SMS_88
表示对流层延迟，/>
Figure SMS_92
表示载波频率，/>
Figure SMS_81
表示接收机钟差，/>
Figure SMS_86
表示流动站/>
Figure SMS_89
和基站
Figure SMS_91
的单差整周模糊度，/>
Figure SMS_82
表示单差量测噪声。
步骤3，基于单差载波相位测量值，计算双差载波相位测量值。
基于单差载波相位测量值方程（公式（1）），在流动站接收机
Figure SMS_93
和基站接收机/>
Figure SMS_94
之间对卫星/>
Figure SMS_95
和卫星/>
Figure SMS_96
之间的单差载波相位测量值进行差分，即在站间和星间各自求一次差分，接收机钟差/>
Figure SMS_97
可以被双差进一步消除。由于/>
Figure SMS_98
，则双差载波相位测量值/>
Figure SMS_99
计算公式可以表示为：
Figure SMS_100
（2）
式中，
Figure SMS_103
表示波长，/>
Figure SMS_104
、/>
Figure SMS_109
分别为基线向量在基站对卫星/>
Figure SMS_102
和卫星/>
Figure SMS_106
观测的投影方向，/>
Figure SMS_108
为基线向量，/>
Figure SMS_110
表示流动站接收机/>
Figure SMS_101
和基站接收机/>
Figure SMS_105
的双差整周模糊度，/>
Figure SMS_107
表示双差量测噪声。
步骤4，基于双差载波相位测量值和三频信号线性组合，建立超宽巷线性组合，并进行超宽巷模糊度固定解算。
在相同测量噪声量的条件下，组合测量值的载波波长越长，对其载波相位整周模糊度的求解就越有利。基于双差载波相位测量值方程（公式（2）），针对北斗二号系统（B1I,B2I, B3I）三频信号和北斗三号系统（B1I, B2b, B3I）三频信号，得到某一时刻的三频双差载波相位测量值
Figure SMS_111
分别为：
Figure SMS_112
（3）
Figure SMS_113
（4）
Figure SMS_114
（5）
式中，下标“1”“2”“3”分别表示B1、B2、B3三个载波频率，
Figure SMS_115
表示波长，/>
Figure SMS_116
表示几何距离，/>
Figure SMS_117
表示卫星星历误差，/>
Figure SMS_118
表示对流层延迟，/>
Figure SMS_119
表示电离层延迟，/>
Figure SMS_120
表示整周模糊度，/>
Figure SMS_121
表示量测噪声。
通过对这些多频信号的测量值进行线性组合，尤其是超宽巷、宽巷组合方式，可以创造出具有长波长的组合测量值，从而促进整周模糊度能够被准确又快速地求解。基于公式（3）-公式（5），对三频双差载波测量值进行三频信号线性组合
Figure SMS_122
的表达式可以写成：
Figure SMS_123
（6）
式中，系数
Figure SMS_124
既可以是整数，也可以是分数，将组合系数记成/>
Figure SMS_125
。
基于公式（6）设置组合系数(0, -1, 1)为超宽巷线性组合，由于超宽巷为长波长观测值，中短基线下的组合方程残余误差小于该波长，因此对超宽巷模糊度取整能够得到超宽巷模糊度固定解，即：
Figure SMS_126
（7）
式中，
Figure SMS_127
表示超宽巷模糊度，下标/>
Figure SMS_128
表示超宽巷组合，/>
Figure SMS_129
是超宽巷的波长，/>
Figure SMS_130
是超宽巷双差载波组合测量值，/>
Figure SMS_131
是超宽巷双差伪距组合测量值。
步骤5，建立宽巷线性组合，并基于超宽巷模糊度固定解后的双差伪距组合测量值进行宽巷模糊度固定解算。
基于公式（6）设置组合系数(1, -1, 0)为宽巷线性组合，联合公式（7）超宽巷模糊度固定解后的修正双差伪距组合测量值，对宽巷模糊度取整得到宽巷模糊度固定解，则有：
Figure SMS_132
（8）
式中，
Figure SMS_133
表示宽巷模糊度，下标/>
Figure SMS_134
表示宽巷组合，/>
Figure SMS_135
是宽巷的波长，/>
Figure SMS_136
是宽巷双差载波组合测量值，/>
Figure SMS_137
是超宽巷修正双差伪距组合观测值。
步骤6，建立窄巷线性组合，基于宽巷模糊度固定解后的双差伪距组合测量值，利用LAMBDA算法解算窄巷模糊度。
虽然超宽巷和宽巷组合的长波长便于模糊度固定，但组合观测噪声大，为了得到高精度厘米级的定位精度，基于公式（6）设置组合系数(1, 1, 0)为窄巷线性组合，联合公式（8）宽巷模糊度固定解后的修正双差伪距组合测量值，则有：
Figure SMS_138
（9）
式中，
Figure SMS_139
表示窄巷模糊度，下标/>
Figure SMS_140
表示窄巷组合，/>
Figure SMS_141
是窄巷的波长，/>
Figure SMS_142
是窄巷双差载波组合测量值，/>
Figure SMS_143
是宽巷修正双差伪距组合观测值。
由于窄巷波长小，组合观测噪声小，直接取整窄巷模糊度容易产生错误，因此采用最小二乘模糊度降相关算法（Leastsquare AMBiguity Decorrelation Adjustment，LAMBDA）解算窄巷模糊度。
步骤7，基于最小二乘算法求解最优基线向量校正量，得到高精度定位结果。
联合公式（9），双差载波相位的相对定位方程的矩阵形式为：
Figure SMS_144
（10）
式中，
Figure SMS_145
是双差载波相位测量值，/>
Figure SMS_146
是基线向量校正量，/>
Figure SMS_147
和/>
Figure SMS_148
是常系数矩阵，
Figure SMS_149
是窄巷模糊度。
基于最小二乘原理求得最优解
Figure SMS_150
使得测量残差的加权平方和最小，得到基线向量校正量/>
Figure SMS_151
，进而得到最优基线，再结合已知基站坐标，即可得到最终的高精度定位结果。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施案例，做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (8)

1.一种联合北斗二号和北斗三号的三频差分定位方法，其特征在于，包括以下步骤：
步骤1，获取北斗二号、北斗三号三频信号原始测量值；
步骤2，计算单差载波相位测量值；
步骤3，基于单差载波相位测量值，计算双差载波相位测量值；
步骤4，基于双差载波相位测量值和三频信号线性组合，建立超宽巷线性组合，并进行超宽巷模糊度固定解算；
步骤5，建立宽巷线性组合，并基于超宽巷模糊度固定解后的双差伪距组合测量值进行宽巷模糊度固定解算；
步骤6，建立窄巷线性组合，基于宽巷模糊度固定解后的双差伪距组合测量值，利用LAMBDA算法解算窄巷模糊度；
步骤7，基于最小二乘算法求解最优基线向量校正量，得到高精度定位结果。
2.如权利要求1所述的一种联合北斗二号和北斗三号的三频差分定位方法，其特征在于：步骤2中流动站接收机
Figure QLYQS_1
和基站接收机/>
Figure QLYQS_2
之间对卫星/>
Figure QLYQS_3
的单差载波相位测量值/>
Figure QLYQS_4
为：
Figure QLYQS_5
（1）
式中，
Figure QLYQS_7
表示波长，/>
Figure QLYQS_9
表示流动站/>
Figure QLYQS_11
和基站/>
Figure QLYQS_13
之间的几何距离，/>
Figure QLYQS_15
表示电离层延迟，
Figure QLYQS_16
表示对流层延迟，/>
Figure QLYQS_17
表示载波频率，/>
Figure QLYQS_6
表示接收机钟差，/>
Figure QLYQS_8
表示流动站/>
Figure QLYQS_10
和基站/>
Figure QLYQS_12
的单差整周模糊度，/>
Figure QLYQS_14
表示单差量测噪声。
3.如权利要求1所述的一种联合北斗二号和北斗三号的三频差分定位方法，其特征在于：步骤3中基于公式（1）单差载波相位测量值方程，在流动站接收机
Figure QLYQS_18
和基站接收机/>
Figure QLYQS_19
之间对卫星/>
Figure QLYQS_20
和卫星/>
Figure QLYQS_21
之间的单差载波相位测量值进行差分，由于/>
Figure QLYQS_22
，因此得到双差载波相位测量值/>
Figure QLYQS_23
：
Figure QLYQS_24
（2）
式中，
Figure QLYQS_26
表示波长，/>
Figure QLYQS_28
、/>
Figure QLYQS_30
分别为基线向量在基站对卫星/>
Figure QLYQS_31
和卫星/>
Figure QLYQS_32
观测的投影方向，
Figure QLYQS_33
为基线向量，/>
Figure QLYQS_34
表示流动站接收机/>
Figure QLYQS_25
和基站接收机/>
Figure QLYQS_27
的双差整周模糊度，/>
Figure QLYQS_29
表示双差量测噪声。
4.如权利要求1所述的一种联合北斗二号和北斗三号的三频差分定位方法，其特征在于：步骤4中基于公式（2）双差载波相位测量值方程，针对北斗二号系统B1I、B2I、B3I三频信号和北斗三号系统B1I、B2b、B3I三频信号，得到某一时刻的三频双差载波相位测量值
Figure QLYQS_35
分别为：
Figure QLYQS_36
（3）
Figure QLYQS_37
（4）
Figure QLYQS_38
（5）
式中，下标“1”“2”“3”分别表示B1、B2、B3三个载波频率，
Figure QLYQS_39
表示波长，/>
Figure QLYQS_40
表示几何距离，
Figure QLYQS_41
表示卫星星历误差，/>
Figure QLYQS_42
表示对流层延迟，/>
Figure QLYQS_43
表示电离层延迟，/>
Figure QLYQS_44
表示整周模糊度，/>
Figure QLYQS_45
表示量测噪声；
基于公式（3）-公式（5），对三频双差载波测量值进行三频信号线性组合
Figure QLYQS_46
的表达式为：
Figure QLYQS_47
（6）
式中，系数
Figure QLYQS_48
为整数或者分数，将组合系数记成/>
Figure QLYQS_49
。
5.如权利要求4所述的一种联合北斗二号和北斗三号的三频差分定位方法，其特征在于：步骤4中基于公式（6）设置组合系数(0, -1, 1)为超宽巷线性组合，由于超宽巷为长波长观测值，中短基线下的组合方程残余误差小于该波长，因此对超宽巷模糊度取整能够得到超宽巷模糊度固定解，即：
Figure QLYQS_50
（7）
式中，
Figure QLYQS_51
表示超宽巷模糊度，下标/>
Figure QLYQS_52
表示超宽巷组合，/>
Figure QLYQS_53
是超宽巷的波长，/>
Figure QLYQS_54
是超宽巷双差载波组合测量值，/>
Figure QLYQS_55
是超宽巷双差伪距组合测量值。
6.如权利要求1所述的一种联合北斗二号和北斗三号的三频差分定位方法，其特征在于：步骤5中基于公式（6）设置组合系数(1, -1, 0)为宽巷线性组合，联合公式（7）超宽巷模糊度固定解后的修正双差伪距组合测量值，对宽巷模糊度取整得到宽巷模糊度固定解，则有：
Figure QLYQS_56
（8）
式中，
Figure QLYQS_57
表示宽巷模糊度，下标/>
Figure QLYQS_58
表示宽巷组合，/>
Figure QLYQS_59
是宽巷的波长，/>
Figure QLYQS_60
是宽巷双差载波组合测量值，/>
Figure QLYQS_61
是超宽巷修正双差伪距组合观测值。
7.如权利要求1所述的一种联合北斗二号和北斗三号的三频差分定位方法，其特征在于：步骤6中为了得到高精度厘米级的定位精度，基于公式（6）设置组合系数(1, 1, 0)为窄巷线性组合，联合公式（8）宽巷模糊度固定解后的修正双差伪距组合测量值，则有：
Figure QLYQS_62
（9）
式中，
Figure QLYQS_63
表示窄巷模糊度，下标/>
Figure QLYQS_64
表示窄巷组合，/>
Figure QLYQS_65
是窄巷的波长，
Figure QLYQS_66
是窄巷双差载波组合测量值，/>
Figure QLYQS_67
是宽巷修正双差伪距组合观测值；
由于窄巷波长小，组合观测噪声小，直接取整窄巷模糊度容易产生错误，因此采用最小二乘模糊度降相关算法LAMBDA解算窄巷模糊度。
8.如权利要求1所述的一种联合北斗二号和北斗三号的三频差分定位方法，其特征在于：步骤7中联合公式（9），双差载波相位的相对定位方程的矩阵形式为：
Figure QLYQS_68
（10）
式中，
Figure QLYQS_69
是双差载波相位测量值，/>
Figure QLYQS_70
是基线向量校正量，/>
Figure QLYQS_71
和/>
Figure QLYQS_72
是常系数矩阵，/>
Figure QLYQS_73
是窄巷模糊度；
基于最小二乘原理求得最优解
Figure QLYQS_74
使得测量残差的加权平方和最小，得到基线向量校正量/>
Figure QLYQS_75
，进而得到最优基线，再结合已知基站坐标，即可得到最终的高精度定位结果。
CN202310517617.XA 2023-05-10 2023-05-10 一种联合北斗二号和北斗三号的三频差分定位方法 Active CN116359968B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310517617.XA CN116359968B (zh) 2023-05-10 2023-05-10 一种联合北斗二号和北斗三号的三频差分定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310517617.XA CN116359968B (zh) 2023-05-10 2023-05-10 一种联合北斗二号和北斗三号的三频差分定位方法

Publications (2)

Publication Number Publication Date
CN116359968A true CN116359968A (zh) 2023-06-30
CN116359968B CN116359968B (zh) 2023-09-12

Family

ID=86909843

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310517617.XA Active CN116359968B (zh) 2023-05-10 2023-05-10 一种联合北斗二号和北斗三号的三频差分定位方法

Country Status (1)

Country Link
CN (1) CN116359968B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116736352A (zh) * 2023-08-01 2023-09-12 深圳市中车智联科技有限公司 一种基于北斗差分定位的移动交通设备追踪系统及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006108227A1 (en) * 2005-04-11 2006-10-19 Yanming Feng Improved phase ambiguity resolution using three gnss signals
CN105699999A (zh) * 2016-02-01 2016-06-22 武汉大学 一种固定北斗地基增强系统基准站窄巷模糊度的方法
US20170269216A1 (en) * 2016-03-18 2017-09-21 Deere & Company Rapid determination of precise position by aiding data
US20180252819A1 (en) * 2015-02-26 2018-09-06 Profound Positioning Inc. Method and system for performing precise point positioning (ppp) ambiguity resolution using gnss triple frequency signals
CN109001781A (zh) * 2018-08-01 2018-12-14 太原理工大学 一种顾及电离层约束的bds三频模糊度解算方法
CN112987059A (zh) * 2021-04-16 2021-06-18 湖南联智科技股份有限公司 一种基于三频模糊度解算的整周模糊度计算方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2006108227A1 (en) * 2005-04-11 2006-10-19 Yanming Feng Improved phase ambiguity resolution using three gnss signals
US20180252819A1 (en) * 2015-02-26 2018-09-06 Profound Positioning Inc. Method and system for performing precise point positioning (ppp) ambiguity resolution using gnss triple frequency signals
CN105699999A (zh) * 2016-02-01 2016-06-22 武汉大学 一种固定北斗地基增强系统基准站窄巷模糊度的方法
US20170269216A1 (en) * 2016-03-18 2017-09-21 Deere & Company Rapid determination of precise position by aiding data
CN109001781A (zh) * 2018-08-01 2018-12-14 太原理工大学 一种顾及电离层约束的bds三频模糊度解算方法
CN112987059A (zh) * 2021-04-16 2021-06-18 湖南联智科技股份有限公司 一种基于三频模糊度解算的整周模糊度计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HUICHAO SHAO 等: "Analysis of Triple-Frequency Multiple Combination Differential Positioning of BDS-3", 《CHINA SATELLITE NAVIGATION CONFERENCE (CSNC 2021) PROCEEDINGS》, pages 427 - 438 *
张冠显 等: "改进的北斗三频RTK整周模糊度固定方法", 传感器与微系统, no. 06 *
赵才新 等: "一种改进的北斗三频基线解算方法", 《测绘科学》, vol. 44, no. 2, pages 18 - 25 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116736352A (zh) * 2023-08-01 2023-09-12 深圳市中车智联科技有限公司 一种基于北斗差分定位的移动交通设备追踪系统及方法
CN116736352B (zh) * 2023-08-01 2023-12-26 深圳市中车智联科技有限公司 一种基于北斗差分定位的移动交通设备追踪系统及方法

Also Published As

Publication number Publication date
CN116359968B (zh) 2023-09-12

Similar Documents

Publication Publication Date Title
CN109581452B (zh) 一种gnss参考站载波相位整周模糊度解算方法
CN108415049B (zh) 提高网络rtk双差宽巷模糊度固定正确率的方法
CN110275186B (zh) Leo卫星增强的gnss电离层归一化与融合建模方法
CN109001781B (zh) 一种顾及电离层约束的bds三频模糊度解算方法
CN108802782B (zh) 一种惯导辅助的北斗三频载波相位整周模糊度求解方法
CN107728171B (zh) 基于粒子滤波的gnss相位系统间偏差实时追踪和精密估计方法
CN109581455B (zh) 一种bds和gps融合的三频宽巷紧组合定位方法
CN111505685A (zh) 一种基于改正系统间偏差的多系统组合rtk模型的定位方法
CN116148909B (zh) 多频多模非组合精密单点定位瞬时窄巷模糊度固定方法
CN109884679B (zh) 一种单模gnss系统跨频点混合双差rtk解算方法
CN116299615B (zh) 一种实现单北斗实时ppp模糊固定的相位偏差估计方法
CN116359968B (zh) 一种联合北斗二号和北斗三号的三频差分定位方法
CN112146557A (zh) 一种基于gnss的实时桥梁变形监测系统及方法
CN115079236B (zh) 通过低轨增强缩短广域非差非组合ppp-rtk定位收敛时间的方法
CN111983650B (zh) 一种基于gnss的高精度时间传递方法
CN110568464B (zh) 基于bds/gnss多模芯片的精密定位方法及装置
CN115933356A (zh) 一种虚拟原子钟的高精度时间同步系统和方法
CN110646823A (zh) 一种基于Helmet后验定权法的GPS\BDS紧组合精密单点定位方法
CN110095796B (zh) 一种跨系统mw紧组合实时动态精密导航定位方法
CN113156477B (zh) 一种基于Android智能手机的高精度RTK定位方法
CN112485813B (zh) Glonass测站间非组合测距码频间偏差校正方法及系统
CN115061170A (zh) 短距离大高差环境网络rtk方法
CN110208841B (zh) 一种改进的面向非重叠频率的gnss紧组合方法
CN111323748B (zh) 差分定位方法及系统
CN111830545A (zh) Bds+gps双系统高精度定位双频rtk方法

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