CN111505685A - 一种基于改正系统间偏差的多系统组合rtk模型的定位方法 - Google Patents

一种基于改正系统间偏差的多系统组合rtk模型的定位方法 Download PDF

Info

Publication number
CN111505685A
CN111505685A CN202010296342.8A CN202010296342A CN111505685A CN 111505685 A CN111505685 A CN 111505685A CN 202010296342 A CN202010296342 A CN 202010296342A CN 111505685 A CN111505685 A CN 111505685A
Authority
CN
China
Prior art keywords
difference
intersystem
phase
gnss
satellite
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
CN202010296342.8A
Other languages
English (en)
Other versions
CN111505685B (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.)
National Time Service Center of CAS
Original Assignee
National Time Service Center of CAS
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 National Time Service Center of CAS filed Critical National Time Service Center of CAS
Priority to CN202010296342.8A priority Critical patent/CN111505685B/zh
Publication of CN111505685A publication Critical patent/CN111505685A/zh
Application granted granted Critical
Publication of CN111505685B publication Critical patent/CN111505685B/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
    • 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/40—Correcting position, velocity or attitude
    • G01S19/41—Differential correction, e.g. DGPS [differential GPS]

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

本发明公开了一种基于改正系统间偏差的多系统组合RTK模型的定位方法，首先将基准站和流动站原始的伪距和相位观测值进行系统内和系统间差分，组成系统内和系统间双差观测方程并进行参数估计得到伪距和相位差分系统间偏差(DISB)；然后对DISB进行先验校正来提高系统间差分定位模型的强度，从而提升多GNSS RTK的准确性、统一性和可用性，通过进行参数估计，完成基线解算和模糊度固定，并基于基准站的准确坐标，获得流动站的定位结果。本发明方法通过系统间差分原理，有效地削弱了部分共性误差，增加了组合观测量，提高了参数求解强度和冗余度；同时有效地消除了接收机钟差参数，可以得到高精度的测站位置。

Description

一种基于改正系统间偏差的多系统组合RTK模型的定位方法
技术领域
本发明涉及卫星导航技术领域，具体是涉及一种基于改正系统间偏差的多系统组合RTK模型的定位方法。
背景技术
系统间偏差(Inter-System Biases,ISB)与多频多模全球导航卫星系统(GlobalNavigation Satellite System,GNSS)的数据处理密切相关；ISB除系统时空基准差异外，主要由硬件延迟组成，由GNSS 设备中的不同信号路径产生，这取决于接收机内部的相关性。
如果可以合理处理GNSS差分观测的系统间偏差(Differential Inter-SystemBiases,DISB)，那么使用系统间差分模型可以获得最大化冗余，这对于在严苛观测环境中的差分精密定位至关重要，例如城市地区信号很容易被高层建筑或树木遮挡。
由此可见，对DISB的改正处理有利于提高系统间差分定位模型的强度，从而提升多GNSS RTK的准确性、统一性和可用性，对于多GNSS兼容与互操作具有重要价值，因此，提供一种能够改正系统间偏差的多系统组合模型及方法对于定位测站等具有重要意义。
发明内容
为解决上述技术问题，本发明提供了一种基于改正系统间偏差的多系统组合RTK模型的定位方法，通过系统间差分原理，有效地削弱了部分共性误差，增加了组合观测量，提高了参数求解强度和冗余度；同时有效地消除了接收机钟差参数，可以得到高精度的测站位置。
本发明的技术方案是：一种基于改正系统间偏差的多系统组合RTK模型的定位方法，首先将基准站和流动站原始的伪距和相位观测值进行系统内和系统间差分，组成系统内和系统间双差观测方程并进行参数估计得到伪距和相位差分系统间偏差；然后对DISB进行先验校正来提高系统间差分定位模型的强度，从而提升多GNSS RTK的准确性、统一性和可用性，具体包括以下步骤：
步骤一，数据获取：获取基准站和流动站GNSS伪距和相位观测值、GNSS卫星星历以及基准站准确坐标，基于GNSS卫星星历，对获得的GNSS伪距和相位观测值进行数据质量检查、粗差剔除，删除无卫星星历或观测值不完整的数据，得到干净的数据；
步骤二，观测值差分处理：对于两个独立的GNSS系统，选定某一个系统的某颗卫星作为参考卫星，与该系统内的其他卫星进行差分形成系统内差分观测值，与另外一个系统的每颗卫星进行差分形成系统间差分观测值，并组成GNSS系统内和系统间双差观测方程；
步骤三，差分系统间偏差(DISB)估计模型建立：根据系统内和系统间差分观测值，建立 DISB实时估计模型，采用高度角和观测值噪声确定随机模型，采用最小二乘平差方法进行估计，估计出伪距和相位双差系统间偏差；
步骤四，DISB改正的RTK定位：在系统间差分的RTK定位模型中，将估计出来的伪距和相位 DISB作为先验信息，对系统间差分观测值进行改正，然后再进行参数估计，完成基线解算和模糊度固定，并基于基准站的准确坐标，获得流动站的定位结果。
进一步地，所述步骤二中两个独立的GNSS系统分别记作A系统和B系统，所述GNSS系统内和系统间双差观测方程分别如下所示：
1)GNSS系统内双差观测方程：
GNSS系统A系统内双差伪距和相位观测方程分别表示为：
Figure BDA0002452328360000021
Figure BDA0002452328360000022
其中，P和Φ分别表示以米为单位的伪距和相位观测值；下标b和r分别表示基准站和用户站；下标
Figure 1
表示观测值频率；A表示GNSS系统；1A表示在A系统中选择的参考卫星；
Figure BDA0002452328360000024
表示站星双差运算符，
Figure BDA0002452328360000025
ρ表示卫星和接收机天线间的站星几何距离；λ表示载波频率的波长；N表示以周为单位的整数模糊度；ε和e分别表示伪距和相位中的测量噪声、多径误差以及未建模误差之和；
2)GNSS系统间双差观测方程：
GNSS系统A系统和B系统的系统间双差伪距和相位观测方程分别表示为：
Figure BDA0002452328360000026
Figure BDA0002452328360000027
其中，
Figure BDA0002452328360000031
表示系统间差分伪距偏差(code differential inter-systembiases,DISCB)，
Figure BDA0002452328360000032
表示系统间差分相位偏差(phase differential inter-systembiases,DISPB)。
Figure BDA0002452328360000033
Figure BDA0002452328360000034
其中，IFCB表示伪距频间偏差；δ和
Figure BDA0002452328360000035
表示相位中的硬件延迟和初始相位偏差；1B表示选取的 B系统参考卫星。
进一步地，所述步骤三中采用高度角和观测噪声确定随机模型，其具体形式为：
Figure BDA0002452328360000036
其中σ(θ)表示非差观测值的标准差；θ表示卫星高度角；σ0表示天顶方向的非差观测值的标准差, 各个频点天顶方向的伪距和相位观测值的标准差分别取0.3m和3mm；GPS、BDS、Galileo同类观测值的标准差与GLONASS观测值的标准差的比值设置为1:1.5。
更进一步地，所述步骤三中建立DISB实时估计模型的具体形式为：
根据GNSS系统的A系统和B系统的双频伪距和相位观测值，建立DISB实时估计模型：
Figure BDA0002452328360000037
其中，P和Φ分别表示以米为单位的伪距和相位观测值；下标b和r分别表示基准站和用户站； 1A表示A系统的参考卫星；sA表示A系统的非参考卫星；1B表示B系统的参考卫星；sB表示B系统的非参考卫星；
Figure BDA0002452328360000038
表示A系统观测值频率(i＝1,2)；
Figure BDA0002452328360000039
表示B系统观测值频率(i＝1,2)；
Figure BDA00024523283600000310
表示系统A和B的系统间差分伪距偏差；
Figure BDA00024523283600000311
表示系统A和B的系统间差分相位偏差；ρ表示卫星和接收机天线间的站星几何距离；λ表示载波频率的波长；N表示以周为单位的整数模糊度；ε和e 分别表示伪距和相位中的测量噪声、多径误差以及未建模误差之和。
对应的误差方程可以表示为：
Figure BDA0002452328360000041
其中，V代表残差向量；A是设计矩阵；I是单位矩阵；0是零矩阵；X是基线分量参数(ΔXr，ΔYr，ΔZr)；N是双差模糊度参数；L是常数项。
进一步地，所述步骤四中所述系统间差分的RTK定位模型的具体形式为：
Figure BDA0002452328360000042
其中，
Figure BDA0002452328360000043
表示系统A和B的系统间差分伪距偏差改正数；
Figure BDA0002452328360000044
表示系统A和B的系统间差分相位偏差改正数。
对应的误差方程可以表示为：
Figure BDA0002452328360000045
其中，V代表残差向量；A是设计矩阵；I是单位矩阵；0是零矩阵；X是基线分量参数(ΔXr，ΔYr，ΔZr)；λ表示载波频率的波长；N是双差模糊度参数；L是常数项；
Figure BDA0002452328360000046
表示系统A和B的系统间差分伪距偏差改正数；
Figure BDA0002452328360000047
表示系统A和B的系统间差分相位偏差改正数。
本发明的有益效果是:
(1)本发明采用基于改正系统间偏差的多系统组合RTK模型的定位方法，可以有效减弱共性误差，增加观测量，提高参数求解的强度和性能；通过系统差分原理可以消除一些诸如多路径、接收机钟差、模型偏差等共性误差，同时通过观测值差分组合，增加了观测量，提高了参数求解的模型强度和参数估计性能。
(2)本发明采用基于改正系统间偏差的多系统组合RTK模型的定位方法，可以在严苛的观测条件下，亦可以实现高精度的相对定位；相对于传统系统内差分模型，某一系统必须观测到4颗及以上卫星，才可以组成3组双差观测方程进行参数求解，本发明在极端苛刻的条件下(单系统观测的卫星数量均少于4颗，但多系统卫星总数量大于等于4颗，也可以组成3组双差观测方程)，通过构建系统间差分观测方程，仍然可以进行定位解算，在实际应用中具有重要意义。
附图说明
图1是本发明的多系统组合RTK流程图。
图2是本发明实施例的基线CUAA-CUT3中GPS参考卫星(G08)与GPS卫星(G01、G18)，BDS卫星(C02、C08)所形成的双差伪距和相位观测值。
图3是本发明实施例的基线CUAA-CUT3 GPS-BDS伪距和相位L1-B1和L2-B2 DISB时间序列。
图4是本发明实施例的基线CUAA-CUT3 DISB实时估计与改正模型的定位结果偏差。
具体实施方式
为了便于本领域普通技术人员理解和实施本发明，下面结合附图及实施例对本发明作进一步的详细描述，应当理解，此处所描述的实施例仅用于说明和解释本发明，并不用于限定本发明。
本发明将基准站和流动站原始的伪距和相位观测值直接进行系统间差分，构造系统间差分观测值，可以消除部分共性误差的影响、增加观测量、提高模型强度，对差分系统间偏差参数进行估计，然后将估计出来的差分系统间偏差参数作为先验信息，对系统间差分观测值进行改正，可以得到高精度的GNSS RTK定位结果。
其主要包括以下三大模型，具体如下：
(一)差分系统间偏差实时估计模型
在短基线情况下，忽略电离层等大气延迟误差的影响，顾及DISB的时域稳定性，在连续观测时段内可采用常数模型对其进行多历元实时估计，DISB实时估计模型可以表示为：
Figure BDA0002452328360000061
其中，P和Φ分别表示以米为单位的伪距和相位观测值；下标b和r分别表示基准站和用户站；1A表示A系统的参考卫星；sA表示A系统的非参考卫星；1B表示B系统的参考卫星；sB表示B系统的非参考卫星；
Figure BDA0002452328360000062
表示A系统观测值频率(i＝1,2)；
Figure BDA0002452328360000063
表示B系统观测值频率(i＝1,2)；
Figure BDA0002452328360000064
表示系统A和B 的系统间差分伪距偏差；
Figure BDA0002452328360000065
表示系统A和B的系统间差分相位偏差；ρ表示卫星和接收机天线间的站星几何距离；λ表示载波频率的波长；N表示以周为单位的整数模糊度；ε和e分别表示伪距和相位中的测量噪声、多径误差以及未建模误差之和。
对应的误差方程可以表示为：
Figure BDA0002452328360000066
其中，V代表残差向量；A是设计矩阵；I是单位矩阵；0是零矩阵；X是基线分量参数(ΔXr，ΔYr，ΔZr)；N是双差模糊度参数；L是常数项。
(二)改正DISB的系统间差分定位模型
对于A系统与B系统间差分模型来说，将估计出来的伪距和相位DISB作为先验信息，对系统间差分观测值进行改正，改正DISB的系统间差分定位模型表示为：
Figure BDA0002452328360000071
其中，
Figure BDA0002452328360000072
表示系统A和B的系统间差分伪距偏差改正数；
Figure BDA0002452328360000073
表示系统A和B的系统间差分相位偏差改正数。
对应的误差方程可以表示为：
Figure BDA0002452328360000074
其中，V代表残差向量；A是设计矩阵；I是单位矩阵；0是零矩阵；X是基线分量参数(ΔXr，ΔYr，ΔZr)；λ表示载波频率的波长；N是双差模糊度参数；L是常数项；
Figure BDA0002452328360000075
表示系统A和B的系统间差分伪距偏差改正数；
Figure BDA0002452328360000076
表示系统A和B的系统间差分相位偏差改正数。
(三)系统间差分模型中的随机模型
随机模型描述的是观测量的统计特性，可以通过一个先验的方差-协方差矩阵来表达。随机模型的合理确定是获得正确平差结果的前提，不合理的随机模型可能导致平差结果存在系统性偏差以及各项精度评估指标不可靠。
非差观测值的随机模型采用高度角和观测值噪声确定，
Figure BDA0002452328360000077
其中σ(θ)表示非差观测值的标准差；θ表示卫星高度角；σ0表示天顶方向的非差观测值的标准差，各个频点天顶方向的伪距和相位观测值的标准差分别取0.3m和3mm。
设站间单差伪距和相位观测值方差-协方差矩阵R，进行星间差分的映射矩阵为D，则双差伪距和相位观测值的方差-协方差矩阵可以表示为：
Q＝DRDT (17)
对于系统内差分模型，星间差分的映射矩阵D为：
Figure BDA0002452328360000081
而对于系统间差分模型，映射矩阵D需要考虑两个系统间的星间差分，可表示为：
Figure BDA0002452328360000082
式(18)和(19)中
Figure BDA0002452328360000083
其中
Figure BDA0002452328360000084
表示克罗内克积算子；em表示元素均为1的m维列向量；Im表示m×m的单位矩阵； 0m×n表示元素均为0的m×n维矩阵；mA和mB分别表示A系统和B系统在某一历元所观测到的卫星数。
实施例
本实例采用如下数据：
1、观测文件：从澳大利亚科廷大学GNSS研究中心下载的短基线CUAA-CUT3年积日168-172, 2018，数据网址：http://saegnss2.curtin.edu/ldc/；
2、广播星历文件：当天的brdm全星历文件，下载网址：ftp://cddis.gsfc.nasa.gov/。
首先将基准站和流动站原始的伪距和相位观测值进行系统内和系统间差分，组成系统内和系统间双差观测方程并进行参数估计得到伪距和相位差分系统间偏差；然后对DISB进行先验校正来提高系统间差分定位模型的强度，从而提升多GNSS RTK的准确性、统一性和可用性，如图1所示，其包括数据获取、观测值差分处理、DISB估计模型建立和DISB改正的RTK定位，具体如下：
步骤一，数据获取：获取基准站和流动站站上采集的多模GNSS观测文件和广播星历；
步骤二，观测值差分处理：对于两个独立的GNSS系统，选某一个系统的某颗卫星作为参考卫星，与该系统内的其它卫星进行差分，形成系统内差分观测值；与另外一个系统的每颗卫星进行差分，形成系统间差分观测值；图2为基线CUAA-CUT3中GPS参考卫星(G08)与GPS卫星(G01、 G18)形成的系统内双差伪距和相位观测值，以及与BDS卫星(C02、C08)所形成的系统间双差伪距和相位观测值；
步骤三，DISB估计模型建立：根据系统内和系统间差分观测值，建立DISB实时估计模型，采用高度角和观测值噪声确定随机模型，采用最小二乘平差方法进行DISB的实时估计；图3为基线 CUAA-CUT3 GPS-BDS伪距和相位L1-B1和L2-B2 DISB时间序列，伪距和相位DISB在连续时间内保持稳定；
步骤四，DISB改正的RTK定位：将估计出来的伪距和相位DISB作为先验信息，对系统间差分观测值进行改正，并进行多系统组合RTK定位，实时估计测站坐标；图4为基线CUAA-CUT3 DISB实时估计(左图)与改正模型(右图)的定位结果偏差，结果表明：DISB实时估计模型和改正模型在短基线条件下均可以得到平面优于1厘米、高程优于2厘米的定位精度。

Claims (6)

1.一种基于改正系统间偏差的多系统组合RTK模型的定位方法，其特征在于，包括以下步骤：
步骤一，数据获取：获取基准站和流动站GNSS伪距和相位观测值、GNSS卫星星历以及基准站准确坐标，基于GNSS卫星星历对GNSS伪距和相位观测值进行处理；
步骤二，观测值差分处理：对于两个独立的GNSS系统，选定某一个系统的某颗卫星作为参考卫星，与该系统内的其他卫星进行差分形成系统内差分观测值，与另外一个系统的每颗卫星进行差分形成系统间差分观测值，并组成GNSS系统内和系统间双差观测方程；
步骤三，差分系统间偏差(DISB)估计模型建立：根据系统内和系统间差分观测值，建立DISB实时估计模型，采用高度角和观测值噪声确定随机模型，采用最小二乘平差方法进行估计，估计出伪距和相位双差系统间偏差；
步骤四，DISB改正的RTK定位：在系统间差分的RTK定位模型中，将估计出来的伪距和相位DISB作为先验信息，对系统间差分观测值进行改正，然后再进行参数估计，完成基线解算和模糊度固定，并基于基准站的准确坐标，获得流动站的定位结果。
2.如权利要求1所述的一种基于改正系统间偏差的多系统组合RTK模型的定位方法，其特征在于，所述步骤二中两个独立的GNSS系统分别记作A系统和B系统，所述GNSS系统内和系统间双差观测方程分别如下所示：
1)GNSS系统内双差观测方程：
GNSS系统A系统内双差伪距和相位观测方程分别表示为：
Figure FDA0002452328350000011
Figure FDA0002452328350000012
其中，P和Φ分别表示以米为单位的伪距和相位观测值；下标b和r分别表示基准站和用户站；下标
Figure FDA0002452328350000013
表示观测值频率；A表示GNSS系统；1A表示在A系统中选择的参考卫星；
Figure FDA0002452328350000014
表示站星双差运算符，
Figure FDA0002452328350000015
ρ表示卫星和接收机天线间的站星几何距离；λ表示载波频率的波长；N表示以周为单位的整数模糊度；ε和e分别表示伪距和相位中的测量噪声、多径误差以及未建模误差之和；
2)GNSS系统间双差观测方程：
GNSS系统A系统和B系统的系统间双差伪距和相位观测方程分别表示为：
Figure FDA0002452328350000021
Figure FDA0002452328350000022
其中，
Figure FDA0002452328350000023
表示系统间差分伪距偏差(code differential inter-systembiases,DISCB)，
Figure FDA0002452328350000024
表示系统间差分相位偏差(phase differential inter-systembiases,DISPB)；
Figure FDA0002452328350000025
Figure FDA0002452328350000026
其中，IFCB表示伪距频间偏差；δ和
Figure FDA0002452328350000027
表示相位中的硬件延迟和初始相位偏差；1B表示选取的B系统参考卫星。
3.如权利要求1所述的一种基于改正系统间偏差的多系统组合RTK模型的定位方法，其特征在于，所述步骤三中采用高度角和观测噪声确定随机模型，其具体形式为：
Figure FDA0002452328350000028
其中σ(θ)表示非差观测值的标准差；θ表示卫星高度角；σ0表示天顶方向的非差观测值的标准差,各个频点天顶方向的伪距和相位观测值的标准差分别取0.3m和3mm；GPS、BDS、Galileo同类观测值的标准差与GLONASS观测值的标准差的比值设置为1:1.5。
4.如权利要求1-3任意一项所述的一种基于改正系统间偏差的多系统组合RTK模型的定位方法，其特征在于，所述步骤三中建立DISB实时估计模型的具体形式为：
根据GNSS系统的A系统和B系统的双频伪距和相位观测值，建立DISB实时估计模型：
Figure FDA0002452328350000029
其中，P和Φ分别表示以米为单位的伪距和相位观测值；下标b和r分别表示基准站和用户站；1A表示A系统的参考卫星；sA表示A系统的非参考卫星；1B表示B系统的参考卫星；sB表示B系统的非参考卫星；
Figure FDA00024523283500000210
表示A系统观测值频率(i＝1,2)；
Figure FDA00024523283500000211
表示B系统观测值频率(i＝1,2)；
Figure FDA00024523283500000212
表示系统A和B的系统间差分伪距偏差；
Figure FDA00024523283500000213
表示系统A和B的系统间差分相位偏差；ρ表示卫星和接收机天线间的站星几何距离；λ表示载波频率的波长；N表示以周为单位的整数模糊度；ε和e分别表示伪距和相位中的测量噪声、多径误差以及未建模误差之和。
对应的误差方程可以表示为：
Figure FDA0002452328350000031
其中，V代表残差向量；A是设计矩阵；I是单位矩阵；0是零矩阵；X是基线分量参数(ΔXr，ΔYr，ΔZr)；N是双差模糊度参数；L是常数项。
5.如权利要求1所述的一种基于改正系统间偏差的多系统组合RTK模型的定位方法，其特征在于，所述观测值差分处理的方法为：对于两个独立的GNSS系统，选定某一个系统的某颗卫星作为参考卫星，与该系统内的其他卫星进行差分，形成系统内差分观测值；与另外一个系统的每颗卫星进行差分，形成系统间差分观测值。
6.如权利要求1所述的一种基于改正系统间偏差的多系统组合RTK模型的定位方法，其特征在于，所述步骤四中所述系统间差分的RTK定位模型的具体形式为：
Figure FDA0002452328350000032
其中，
Figure FDA0002452328350000033
表示系统A和B的系统间差分伪距偏差改正数；
Figure FDA0002452328350000034
表示系统A和B的系统间差分相位偏差改正数。
对应的误差方程可以表示为：
Figure FDA0002452328350000041
其中，V代表残差向量；A是设计矩阵；I是单位矩阵；0是零矩阵；X是基线分量参数(ΔXr，ΔYr，ΔZr)；λ表示载波频率的波长；N是双差模糊度参数；L是常数项；
Figure FDA0002452328350000042
表示系统A和B的系统间差分伪距偏差改正数；
Figure FDA0002452328350000043
表示系统A和B的系统间差分相位偏差改正数。
CN202010296342.8A 2020-04-15 2020-04-15 一种基于改正系统间偏差的多系统组合rtk模型的定位方法 Active CN111505685B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010296342.8A CN111505685B (zh) 2020-04-15 2020-04-15 一种基于改正系统间偏差的多系统组合rtk模型的定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010296342.8A CN111505685B (zh) 2020-04-15 2020-04-15 一种基于改正系统间偏差的多系统组合rtk模型的定位方法

Publications (2)

Publication Number Publication Date
CN111505685A true CN111505685A (zh) 2020-08-07
CN111505685B CN111505685B (zh) 2022-03-15

Family

ID=71869206

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010296342.8A Active CN111505685B (zh) 2020-04-15 2020-04-15 一种基于改正系统间偏差的多系统组合rtk模型的定位方法

Country Status (1)

Country Link
CN (1) CN111505685B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112462399A (zh) * 2020-10-22 2021-03-09 中国科学院空天信息创新研究院 一种基于接收机的伪距偏差计算及改正的方法及装置
CN112485814A (zh) * 2020-11-17 2021-03-12 中国人民解放军战略支援部队航天工程大学 Glonass测距码实时校准的rtk定位方法及系统
CN113064189A (zh) * 2021-03-26 2021-07-02 同济大学 一种新模式差分rtk定位方法
CN113433576A (zh) * 2021-06-28 2021-09-24 中国科学院国家授时中心 一种gnss与v-slam融合定位方法及系统
CN113640846A (zh) * 2021-09-17 2021-11-12 深圳市电咖测控科技有限公司 一种用于确定gnss接收器位置的方法及系统
CN113865592A (zh) * 2021-09-09 2021-12-31 河海大学 适用于多频gnss精密导航定位中的多路径参数化方法及存储介质
CN114553982A (zh) * 2022-02-23 2022-05-27 国汽大有时空科技(安庆)有限公司 一种应用于gnss的n-rtk差分增强服务的数据编码方法及数据传输方法
CN114859389A (zh) * 2022-04-18 2022-08-05 华力智芯(成都)集成电路有限公司 一种gnss多系统抗差自适应融合rtk解算方法
CN115639583A (zh) * 2022-12-21 2023-01-24 深圳华大北斗科技股份有限公司 基于Android系统的高精度定位装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101843004B1 (ko) * 2017-11-15 2018-03-29 한국 천문 연구원 다중 위성 항법 시스템의 시스템 바이어스를 이용한 통합 측위 장치 및 그 방법
CN108519614A (zh) * 2018-03-16 2018-09-11 东南大学 一种gps/bds紧组合载波差分定位方法
CN109597105A (zh) * 2018-12-13 2019-04-09 东南大学 一种顾及载波系统间偏差的gps/glonass紧组合定位方法
CN109683182A (zh) * 2018-12-24 2019-04-26 中国电子科技集团公司第二十研究所 一种引入系统间偏差的gnss多模组合差分定位方法
CN110018507A (zh) * 2019-05-08 2019-07-16 中国科学院国家授时中心 一种基于星座间作差的组合精密单点定位方法及系统
CN110646823A (zh) * 2019-10-12 2020-01-03 上海交通大学 一种基于Helmet后验定权法的GPS\BDS紧组合精密单点定位方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101843004B1 (ko) * 2017-11-15 2018-03-29 한국 천문 연구원 다중 위성 항법 시스템의 시스템 바이어스를 이용한 통합 측위 장치 및 그 방법
CN108519614A (zh) * 2018-03-16 2018-09-11 东南大学 一种gps/bds紧组合载波差分定位方法
CN109597105A (zh) * 2018-12-13 2019-04-09 东南大学 一种顾及载波系统间偏差的gps/glonass紧组合定位方法
CN109683182A (zh) * 2018-12-24 2019-04-26 中国电子科技集团公司第二十研究所 一种引入系统间偏差的gnss多模组合差分定位方法
CN110018507A (zh) * 2019-05-08 2019-07-16 中国科学院国家授时中心 一种基于星座间作差的组合精密单点定位方法及系统
CN110646823A (zh) * 2019-10-12 2020-01-03 上海交通大学 一种基于Helmet后验定权法的GPS\BDS紧组合精密单点定位方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
ABBAS ABOLFATHI MOMTAZ ET AL.: "Multi-GNSS Constellation Fusion Based on Statistical Features of Positioning Error", 《ELECTRICAL ENGINEERING (ICEE), IRANIAN CONFERENCE ON》 *
CIRO GIOIA ET AL.: "Inter-system Bias: Stability and Impact on Multi-constellation Positioning", 《2015 IEEE METROLOGY FOR AEROSPACE (METROAEROSPACE)》 *
YIJUN TIAN ET AL.: "Estimating inter-system biases for tightly combined Galileo/BDS/GPS RTK", 《ADVANCES IN SPACE RESEARCH》 *
周英东等: "BDS_Galileo紧组合系统间偏差估计与模糊度固定效果分析", 《测绘科学技术学报》 *
张小红等: "BeiDou B2/Galileo E5b短基线紧组合相对定位模型及性能评估", 《测绘学报》 *
楼益栋等: "GPS/BDS混合双差分RTK定位方法及结果分析", 《大地测量与地球动力学》 *
王彦恒等: "顾及系统间偏差的BDS_GPS单钟差定位法", 《测绘工程》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112462399A (zh) * 2020-10-22 2021-03-09 中国科学院空天信息创新研究院 一种基于接收机的伪距偏差计算及改正的方法及装置
CN112485814B (zh) * 2020-11-17 2023-06-30 中国人民解放军战略支援部队航天工程大学 Glonass测距码实时校准的rtk定位方法及系统
CN112485814A (zh) * 2020-11-17 2021-03-12 中国人民解放军战略支援部队航天工程大学 Glonass测距码实时校准的rtk定位方法及系统
CN113064189A (zh) * 2021-03-26 2021-07-02 同济大学 一种新模式差分rtk定位方法
CN113433576A (zh) * 2021-06-28 2021-09-24 中国科学院国家授时中心 一种gnss与v-slam融合定位方法及系统
CN113433576B (zh) * 2021-06-28 2023-09-01 中国科学院国家授时中心 一种gnss与v-slam融合定位方法及系统
CN113865592A (zh) * 2021-09-09 2021-12-31 河海大学 适用于多频gnss精密导航定位中的多路径参数化方法及存储介质
CN113865592B (zh) * 2021-09-09 2024-05-10 河海大学 适用于多频gnss精密导航定位中的多路径参数化方法及存储介质
CN113640846A (zh) * 2021-09-17 2021-11-12 深圳市电咖测控科技有限公司 一种用于确定gnss接收器位置的方法及系统
CN114553982A (zh) * 2022-02-23 2022-05-27 国汽大有时空科技(安庆)有限公司 一种应用于gnss的n-rtk差分增强服务的数据编码方法及数据传输方法
CN114553982B (zh) * 2022-02-23 2024-02-09 国汽大有时空科技(安庆)有限公司 一种应用于gnss的n-rtk差分增强服务的数据编码方法及数据传输方法
CN114859389A (zh) * 2022-04-18 2022-08-05 华力智芯(成都)集成电路有限公司 一种gnss多系统抗差自适应融合rtk解算方法
CN115639583A (zh) * 2022-12-21 2023-01-24 深圳华大北斗科技股份有限公司 基于Android系统的高精度定位装置

Also Published As

Publication number Publication date
CN111505685B (zh) 2022-03-15

Similar Documents

Publication Publication Date Title
CN111505685B (zh) 一种基于改正系统间偏差的多系统组合rtk模型的定位方法
CN111751853B (zh) 一种gnss双频载波相位整周模糊度解算方法
CN106990424B (zh) 一种双天线gps测姿方法
KR101151782B1 (ko) 단독 위치측정장치 및 단독 위치측정방법
CN106842268B (zh) 双gnss接收机载波相位双差整周模糊度浮点解向量估计方法
CN108415049A (zh) 提高网络rtk双差宽巷模糊度固定正确率的方法
CN110389364B (zh) 定位方法、装置、计算机设备和存储介质
CN109196381A (zh) 通过辅助数据对精确位置的快速确定
US8638257B2 (en) Ultra short baseline GNSS receiver
CN107703526B (zh) 基线测向方法、装置和系统
CN111965673A (zh) 基于多gnss的单频精密单点定位算法的时间频率传递方法
CN111308528A (zh) 一种北斗/gps紧组合虚拟参考站定位方法
CN109669196B (zh) 一种顾及基线形变的多天线gnss载波相位精密测姿方法
CN111538039B (zh) 未知点位基准站精确坐标确定方法
CN112146557A (zh) 一种基于gnss的实时桥梁变形监测系统及方法
CN112394379B (zh) 双天线联合卫星导航定位方法和装置
CN116840876B (zh) 适用于相位失真的卫星导航天线的双天线定向方法及装置
CN114935770B (zh) 一种多历元加快精密单点定位收敛速度的方法及装置
CN110568464B (zh) 基于bds/gnss多模芯片的精密定位方法及装置
CN115079236A (zh) 通过低轨增强缩短广域非差非组合ppp-rtk定位收敛时间的方法
CN114879222A (zh) 一种基于自适应随机模型的全球电离层建模方法
CN109143286B (zh) 一种顾及非模型化误差的卫星导航定位方法
CN116359968B (zh) 一种联合北斗二号和北斗三号的三频差分定位方法
CN115308781B (zh) 基于bdgim辅助的相位平滑伪距高精度时间传递方法
CN115932920A (zh) 对流层延迟的插值方法

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