CN105158783B - 一种实时动态差分定位方法及其设备 - Google Patents

一种实时动态差分定位方法及其设备 Download PDF

Info

Publication number
CN105158783B
CN105158783B CN201510520435.3A CN201510520435A CN105158783B CN 105158783 B CN105158783 B CN 105158783B CN 201510520435 A CN201510520435 A CN 201510520435A CN 105158783 B CN105158783 B CN 105158783B
Authority
CN
China
Prior art keywords
solution
carrier wave
double
float
integer ambiguity
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
Application number
CN201510520435.3A
Other languages
English (en)
Other versions
CN105158783A (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.)
SHANGHAI HIGH GAIN INFORMATION TECHNOLOGY Co Ltd
Original Assignee
SHANGHAI HIGH GAIN 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 SHANGHAI HIGH GAIN INFORMATION TECHNOLOGY Co Ltd filed Critical SHANGHAI HIGH GAIN INFORMATION TECHNOLOGY Co Ltd
Priority to CN201510520435.3A priority Critical patent/CN105158783B/zh
Publication of CN105158783A publication Critical patent/CN105158783A/zh
Application granted granted Critical
Publication of CN105158783B publication Critical patent/CN105158783B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/42Determining position
    • G01S19/43Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining 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/40Correcting position, velocity or attitude
    • G01S19/41Differential 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

本发明公开一种实时动态差分定位方法及其设备,包括:获取基准站和流动站共同接收的多个参考卫星的双频伪距观测值和双频载波相位观测值;将所述宽巷组合方程和无电离层组合方程进行等效代换得到关于所述载波L1和所述载波L2的整周模糊度浮点解的载波相位双差方程;根据所述载波相位双差方程和由所述基准站和流动站的双频伪距观测值建立的双差伪距方程,确定载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度的浮点解;整个定位过程只用到了一次模糊度搜素,有效减少了模糊度搜索的次数,在保证定位精度不损失的条件下提高了实时动态差分定位的效率。

Description

一种实时动态差分定位方法及其设备
技术领域
本发明涉及全球导航卫星系统(GNSS)卫星定位方法,尤其涉及一种实时动态差分定位方法及其设备。
背景技术
随着导航系统的发展,高精度定位的应用范围越来越广泛。在卫星导航系统中,定位的基本原理是从导航系统中获得两种基本的观测量,即伪距和载波相位。伪距是通过解码卫星播报的星历从而获得一种近似距离,该距离由于码片的局限性只能够得到米级的精度。而载波相位观测量要比伪距精确的多,以1575MHZ频率下的载波为例,每个波长仅有19cm,因此准确的获得载波信息是高精度定位、定向的基本方法。在接收机准确的捕获卫星信号获得实时的载波相位值后,定位的问题即转变成了求解载波相位中存在着的整周模糊度。整周模糊度是一个未知的固定整数值,它描述了接收机到卫星之间存在的整数个完整电磁波。而整周模糊度一旦固定,载波相位就能作为没有模糊度的精确的伪距测量,因而模糊度的解算和确定是高精度定位、定向中的核心问题之一。
现有技术中模糊度的解算和确定都是先进行宽巷组合观测量固定模糊度后再用无电离层组合观测量进行相对定位解算,这是因为宽巷观测量组合后波长最大因此有利于模糊度搜索,之后无电离层组合方程解算消除了电离层误差,因此会提高相对定位的精度。但是在实际的应用中发现以上模糊度的解算和确定的整个过程中,都不止一次的进行了模糊度搜索的计算,每多一次模糊度的搜索都会明显的增加定位解算的时间,因此,亟需一种可以减少模糊度搜索次数、缩短解算时间的实时动态差分定位方法。
发明内容
本发明实施例提供一种实时动态差分定位方法及其设备,用以解决现有技术中多次模糊度搜索占用时间较长,定位解算效率低的问题。
本发明方法包括:一种实时动态差分定位方法,包括:获取基准站和流动站共同接收的多个参考卫星的双频伪距观测值和双频载波相位观测值;根据所述双频载波相位观测值建立宽巷组合方程和无电离层组合方程;根据宽巷组合模糊度浮点解和无电离层组合模糊度浮点解与所述参考卫星使用的载波L1和载波L2的整周模糊度浮点解之间的关系,将所述宽巷组合方程和无电离层组合方程进行等效代换得到关于所述载波L1和所述载波L2的整周模糊度浮点解的载波相位双差方程;根据所述载波相位双差方程和由所述基准站和流动站的双频伪距观测值建立的双差伪距方程,确定载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度的浮点解;根据所述载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度浮点解进行一次模糊度搜索得到的载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解,并根据得到的所述载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解确定流动站到基准站的相对位置矢量。
基于同样的发明构思,本发明实施例还提供一种实时动态差分定位的设备,包括:获取单元,用于获取基准站和流动站共同接收的多个参考卫星的双频伪距观测值和双频载波相位观测值;方程建立单元,用于根据所述双频载波相位观测值建立宽巷组合方程和无电离层组合方程;载波相位双差方程确定单元,用于根据根据宽巷组合模糊度浮点解和无电离层组合模糊度浮点解与载波的整周模糊度浮点解之间的关系,将所述宽巷组合方程和无电离层组合方程进行等效代换得到关于所述载波L1和所述载波L2的整周模糊度浮点解的载波相位双差方程;第一确定单元,用于根据所述载波相位双差方程和由所述基准站和流动站的双频伪距观测值建立的双差伪距方程,确定载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度的浮点解;第二确定单元,用于根据所述载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度浮点解进行一次模糊度搜索得到的载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解,并根据得到的所述载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解确定流动站到基准站的相对位置矢量。
本发明实施例一方面结合使用双频载波相位观测值建立宽巷组合方程和无电离层组合方程,并根据宽巷组合模糊度浮点解和无电离层组合模糊度浮点解与所述参考卫星使用的载波L1和载波L2的整周模糊度浮点解之间的关系,将所述宽巷组合方程和无电离层组合方程进行等效代换得到关于所述载波L1和所述载波L2的整周模糊度浮点解的载波相位双差方程;另一方面所述载波相位双差方程结合根据所述基准站和流动站的双频伪距观测值建立的双差伪距方程得到的观测方程组,通过求解该观测方程组得到所述确定载波L1、载波L2的双差整周模糊度的浮点解,最终根据所述载波L1、载波L2的双差整周模糊度浮点解进行一次模糊度搜索得到的载波L1、载波L2双差整周模糊度整数解确定流动站到基准站的相对位置矢量。可见,整个定位过程只用到了一次模糊度搜素,有效减少了模糊度搜索的次数,在保证定位精度不损失的条件下提高了实时动态差分定位的效率。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简要介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为现有技术求解宽巷模糊度浮点解的方法流程示意图;
图2为本发明实施例提供一种实时动态差分定位方法流程示意图;
图3为本发明实施例提供的一种基站站和流动站接收四个参考卫星的测量数据示意图;
图4为本发明实施例提供一种实时动态差分定位步骤流程图;
图5为本发明实施例提供一种实时动态差分定位的设备示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述,显然,所描述的实施例仅仅是本发明一部份实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
载波相位和伪距是卫星定位系统的两大基本观测量。由于载波相位的分辨率远高于伪距,故利用载波相位的载波相位差分定位技术被广泛地应用于高精度静态或动态的定位计算中。载波相位观测量实质上是同一时刻接收机所产生的基准信号的相位和接收到的来自卫星的载波信号的相位之差,其一般表达式为:
其中,为相位差的不足整周的小数部分,为起始历元t0时刻至观测历元t时刻的相位变化整周数,为相应起始历元t0时刻的初始时刻整周模糊度,ε(t)为各种误差的影响。
考虑到若卫星j在历元t0被接收机i锁定,则由接收机自动连续记数,故(1)式中为已知观测量。显然,要实现正确定位必须确定出初始相位整周数亦称初始时刻整周模糊度。理论上讲,初始时刻整周模糊度应为整数,但由于相位观测量中含有许多误差,如卫星与接收机钟差、大气延迟、观测噪声等,使得确定并非易事。可是一旦确定了初始时刻整周模糊度,只要同步观测的卫星数不少于4颗,即使只观测1个历元也能获得唯一的定位结果。因此,如果能寻求一种方法,它可以只用较少的观测数据就能很快地解出相位模糊度,则不仅能较大地提高作业效率,而且在动态定位计算应用中也不必担心信号的失真和相位的周跳,因此,如何准确而快速地确定整周模糊度一直是人们竞相研究的热点课题。
如图1所示,现有技术首先利用载波与伪距的线性组合建立方程求解宽巷模糊度浮点解,根据宽巷模糊度的长波特性,通过多历元求解平均值确定宽巷模糊度整数解;接着采用电离层无关组合建立方程分离出L1模糊度,继而确定L2模糊度;最后进行模糊度解算的正确性检验
可见,该方法进行两次模糊度搜索的计算,每进行一次模糊度的搜索都会明显的增加定位解算的时间,因此模糊度搜索在整体解算过程中占用的时间较长。为了解决这一问题,本发明实施例提供一种实时动态差分定位方法流程示意图,如图2所示,该流程包括:
步骤S101,获取基准站和流动站共同接收的多个参考卫星的双频伪距观测值和双频载波相位观测值。
步骤S102,根据所述双频载波相位观测值建立宽巷组合方程和无电离层组合方程。
步骤S103,根据宽巷组合模糊度浮点解和无电离层组合模糊度浮点解与所述参考卫星使用的载波L1和载波L2的整周模糊度浮点解之间的关系,将所述宽巷组合方程和无电离层组合方程进行等效代换得到关于所述载波L1和所述载波L2的整周模糊度浮点解的载波相位双差方程。
步骤S104,根据所述载波相位双差方程和由所述基准站和流动站的双频伪距观测值建立的双差伪距方程,确定载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度的浮点解。
步骤S105,根据所述载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度浮点解进行一次模糊度搜索得到的载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解,并根据得到的所述载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解确定流动站到基准站的相对位置矢量。
在步骤S101执行之前,如图3所示,所述流动站至少从全球定位系统GPS、北斗卫星导航系统BDS和格洛纳斯导航系统GLONASS中的一个系统中选择出N个参考卫星,其中,所述所述流动站与基准站选择的N个参考卫星相同;所述流动站的接收机接收N个参考卫星的卫星信号以获取所述N个参考卫星的双频伪距观测值和双频载波相位观测值,同时接收所述基准站发送的所述N个参考卫星的双频伪距观测值和双频载波相位观测值,可见,所述流动站可以从不同的导航系统中选择出多个参考卫星,因此所述流动站选择的参考卫星具有如下特点,即所述多个参考卫星中每个参考卫星对应的不同定位系统部分相同或全不相同,这样通过结合使用多个系统的参考卫星,避免了某个系统因为信号干扰造成的影响,提高了定位的准确性。
在步骤S103中,所述宽巷组合整周模糊度浮点解和无电离层组合整周模糊度浮点解与载波L1、载波L2的整周模糊度浮点解之间的关系,包括:
其中,宽巷组合的整周模糊度的浮点解,分别是载波L1和载波L2双差载波观测值的整周模糊度浮点解,为无电离层组合的整周模糊度的浮点解,f1、f2分别是载波L1和载波L2频率,c为光速,为载波相位组合双差值,表示双差算子,表示单差算子。
具体地,对于短基线(基线长度小于10km)情况,利用宽巷组合L5与载波L1、L2观测值的关系,可以将载波L2的双差观测方程中的双差模糊度替换,得到公式三:
而对于中长基线(基线长度大于10km)情况,利用宽巷组合L5、无电离层组合L3与L1、L2观测值的关系,可以将无电离层组合中的模糊度参数用宽巷组合模糊度参数和L1模糊度参数替换。具体推导如下:
无电离层组合观测方程为:
其中,无电离层组合的整周模糊度的浮点解,f1、f2分别是L1和L2载波频率,分别是L1和L2的双差载波观测值,c为光速,为载波相位组合双差值。
在步骤S102中,结合上述模糊度浮点解之间的关系,对现有的宽巷组合方程和无电离层组合方程作相应的等效变换,其中:
无电离层组合方程为:
其中,首先做如下定义:
宽巷组合方程为:
其中,首先做如下定义:
BWL表示矩阵系数[0 -λx2 0 λx2],表示第i个参考卫星的基准站的载波相位观测值和流动站的载波相位观测值的双差测量常数向量,可表示为分别是L1和L2的双差载波观测值,AP是基准站和流动站间的相对位置向量系数矩阵,BWL是基准站和流动站间的整周模糊度向量系数矩阵;基准站和流动站间的载波相位双差测量残差;a:基准站和流动站间的相对位置向量;b:基准站和流动站间的整周模糊度向量。
在此基础之上,我们可以得到载波相位双差方程如下:
考虑到上述载波相位双差方程中的未知数a、b的个数多于方程个数,因此还需要联立双差伪距方程得到观测方程组,所述观测方程组如下:
其中,Lρ表示基准站和流动站间的伪距双差测量常数向量,Vρ表示基准站和流动站间的伪距双差测量残差,Pρ表示基准站和流动站间的伪距双差测量权系数矩阵,q表示基准站和流动站间的伪距双差测量权系数因子,其大小依赖于伪距测量精度及误差水平,典型值为10-4~10-6,其余参数代表含义有上述公式相同。
上述公式[8]可表示为如下的矩阵形式:
公式[9]加上系数矩阵,利用最小二乘法可得:
上式可简化为:
则上式可简记为:
所述结果向量b中包括载波L1、载波L2的双差模糊度浮点解,根据所述载波L1、载波L2的双差模糊度浮点解使用LAMBDA搜索得出多组模糊度和一个表征模糊度可靠性的Ratio门限值,并根据多组模糊度和一个表征模糊度可靠性的Ratio门限值确定所述载波L1的双差整周模糊度整数解和所述载波L2的双差整周模糊度整数解;根据所述双差整周模糊度整数解确定流动站到基准站的相对位置矢量。
可见,本发明实施例的实时动态定位过程如图4所示,流动站和基准站根据双频伪距观测值和双频载波相位观测值分别组建三个方程:无电离层组合观测方程、宽巷组合观测方程和双差伪距观测方程,通过三个方程的联合求解求得载波L1、载波L2的双差模糊度浮点解,然后通过LAMBDA搜索确定双差整周模糊度整数解,根据所述双差整周模糊度整数解确定流动站到基准站的相对位置矢量。这样,整个定位过程只用到了一次模糊度搜素,有效减少了模糊度搜索的次数,在保证定位精度不损失的条件下提高了实时动态差分定位的效率。
基于相同的技术构思,本发明实施例还提供一种实时动态差分定位的设备,该设备可执行上述方法实施例,该设备可以安装在流动站上,也可以独立于基站站和流动站而单独安装。本发明实施例提供的设备如图5所示,包括:获取单元401,方程建立单元402,载波相位双差方程确定单元403,第一确定单元404,第二确定单元405,其中:
获取单元401,用于获取基准站和流动站共同接收的多个参考卫星的双频伪距观测值和双频载波相位观测值;
方程建立单元402,用于根据所述双频载波相位观测值建立宽巷组合方程和无电离层组合方程;
载波相位双差方程确定单元403,用于根据根据宽巷组合模糊度浮点解和无电离层组合模糊度浮点解与载波的整周模糊度浮点解之间的关系,将所述宽巷组合方程和无电离层组合方程进行等效代换得到关于所述载波L1和所述载波L2的整周模糊度浮点解的载波相位双差方程;
第一确定单元404,用于根据所述载波相位双差方程和由所述基准站和流动站的双频伪距观测值建立的双差伪距方程,确定载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度的浮点解;
第二确定单元405,用于根据所述载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度浮点解进行一次模糊度搜索得到的载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解,并根据得到的所述载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解确定流动站到基准站的相对位置矢量。
在获取单元401获取到测量值之前,还需要至少从全球定位系统GPS、北斗卫星导航系统BDS和格洛纳斯导航系统GLONASS系统中的一个系统中选择出N个参考卫星,其中,所述基准站和流动站选择的N个参考卫星相同,并且接收N个参考卫星的卫星信号以获取所述N个参考卫星的双频伪距观测值和双频载波相位观测值,同时接收所述基准站发送的所述N个参考卫星的双频伪距观测值和双频载波相位观测值,可见,所述流动站可以从不同的导航系统中选择出多个参考卫星,因此所述流动站选择的参考卫星具有如下特点,即所述多个参考卫星中每个参考卫星对应的不同定位系统部分相同或全不相同,这样通过结合使用多个系统的参考卫星,避免了某个系统因为信号干扰造成的影响,提高了定位的准确性。
所述宽巷组合整周模糊度浮点解和无电离层组合整周模糊度浮点解与载波L1、载波L2的整周模糊度浮点解之间的关系,包括:
其中,宽巷组合的整周模糊度的浮点解,分别是载波L1和载波L2双差载波观测值的整周模糊度浮点解,为无电离层组合的整周模糊度的浮点解,f1、f2分别是载波L1和载波L2频率,c为光速,为载波相位组合双差值。
所述载波相位双差方程为:
所述第二确定单元405具体用于:根据所述载波L1、载波L2的双差模糊度浮点解使用LAMBDA搜索得出多组模糊度和一个表征模糊度可靠性的Ratio门限值,并根据多组模糊度和一个表征模糊度可靠性的Ratio门限值确定所述载波L1的双差整周模糊度整数解和所述载波L2的双差整周模糊度整数解。
这样,整个定位过程只用到了一次模糊度搜素,有效减少了模糊度搜索的次数,在保证定位精度不损失的条件下提高了实时动态差分定位的效率。
尽管已描述了本发明的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例作出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (10)

1.一种实时动态差分定位方法,其特征在于,包括:
获取基准站和流动站共同接收的多个参考卫星的双频伪距观测值和双频载波相位观测值;
根据所述双频载波相位观测值建立宽巷组合方程和无电离层组合方程;
根据宽巷组合模糊度浮点解和无电离层组合模糊度浮点解与所述参考卫星使用的载波L1和载波L2的整周模糊度浮点解之间的关系,将所述宽巷组合方程和无电离层组合方程进行等效代换得到关于所述载波L1和所述载波L2的整周模糊度浮点解的载波相位双差方程;
根据所述载波相位双差方程和由所述基准站和流动站的双频伪距观测值建立的双差伪距方程,确定载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度的浮点解;
根据所述载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度浮点解进行一次模糊度搜索得到载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解,并根据得到的所述载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解确定流动站到基准站的相对位置矢量。
2.如权利要求1所述的方法,其特征在于,所述多个参考卫星中每个参考卫星对应的不同定位系统部分相同或全不相同。
3.如权利要求1所述的方法,其特征在于,所述宽巷组合整周模糊度浮点解和无电离层组合整周模糊度浮点解与载波L1、载波L2的整周模糊度浮点解之间的关系,包括:
公式一:
公式二:
其中,为宽巷组合的整周模糊度的浮点解,分别是载波L1和载波L2双差载波观测值的整周模糊度浮点解,为无电离层组合的整周模糊度的浮点解,f1、f2分别是载波L1和载波L2频率,c为光速,为载波相位组合双差值。
4.如权利要求1所述的方法,其特征在于,所述载波相位双差方程为:
其中,A:基准站和流动站间的无电离层和宽巷的组合的相对位置向量系数矩阵,B:基准站和流动站间的无电离层和宽巷的组合的整周模糊度向量系数矩阵,a:基准站和流动站间的相对位置向量,b:基准站和流动站间的整周模糊度向量,基准站和流动站间的无电离层和宽巷的组合的载波相位双差测量残差,基准站和流动站间的无电离层和宽巷的组合的载波相位双差测量常数向量,基准站和流动站间的无电离层和宽巷的组合的载波相位双差测量权系数矩阵,P:基础权系数矩阵,n为解算星数,lin、col表示矩阵P的行和列,为测量噪声方差,ALC是基准站和流动站间的无电离层的相对位置向量系数矩阵,AWL是基准站和流动站间的宽巷的的相对位置向量系数矩阵,为基准站和流动站间的无电离层的载波相位双差测量常数向量,为基准站和流动站间的宽巷的载波相位双差测量常数向量,BLC是基准站和流动站间的无电离层的整周模糊度向量系数矩阵,BWL为基准站和流动站间的宽巷的整周模糊度向量系数矩阵,是基准站和流动站间的无电离层的载波相位双差测量残差,是基准站和流动站间的宽巷的载波相位双差测量残差,PLC是基准站和流动站间的无电离层的载波相位双差测量权系数矩阵,PWL是基准站和流动站间的宽巷的载波相位双差测量权系数矩阵。
5.如权利要求1所述的方法,其特征在于,根据所述载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度浮点解进行一次模糊度搜索得到的载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解,包括:
根据所述载波L1、载波L2的双差模糊度浮点解使用LAMBDA搜索得出多组模糊度和一个表征模糊度可靠性的Ratio门限值,并根据多组模糊度和一个表征模糊度可靠性的Ratio门限值确定所述载波L1的双差整周模糊度整数解和所述载波L2的双差整周模糊度整数解。
6.一种实时动态差分定位的设备,其特征在于,包括:
获取单元,用于获取基准站和流动站共同接收的多个参考卫星的双频伪距观测值和双频载波相位观测值;
方程建立单元,用于根据所述双频载波相位观测值建立宽巷组合方程和无电离层组合方程;
载波相位双差方程确定单元,用于根据宽巷组合模糊度浮点解和无电离层组合模糊度浮点解与载波的整周模糊度浮点解之间的关系,将所述宽巷组合方程和无电离层组合方程进行等效代换得到关于载波L1和载波L2的整周模糊度浮点解的载波相位双差方程;
第一确定单元,用于根据所述载波相位双差方程和由所述基准站和流动站的双频伪距观测值建立的双差伪距方程,确定所述载波L1的双差整周模糊度的浮点解和所述载波L2的双差整周模糊度的浮点解;
第二确定单元,用于根据所述载波L1的双差整周模糊度的浮点解和载波L2的双差整周模糊度浮点解进行一次模糊度搜索得到载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解,并根据得到的所述载波L1的双差整周模糊度整数解和载波L2的双差整周模糊度整数解确定流动站到基准站的相对位置矢量。
7.如权利要求6所述的设备,其特征在于,所述多个参考卫星中每个参考卫星对应的不同定位系统部分相同或全不相同。
8.如权利要求6所述的设备,其特征在于,根据宽巷组合模糊度浮点解和无电离层组合模糊度浮点解与所述载波L1和所述载波L2的整周模糊度浮点解之间的关系,包括:
公式一:
公式二:
其中,宽巷组合的整周模糊度的浮点解,分别是载波L1和载波L2双差载波观测值的整周模糊度浮点解,为无电离层组合的整周模糊度的浮点解,f1、f2分别是载波L1和载波L2频率,c为光速,为载波相位组合双差值。
9.如权利要求6所述的设备,其特征在于,所述载波相位双差方程为:
其中,A:基准站和流动站间的无电离层和宽巷的组合的相对位置向量系数矩阵,B:基准站和流动站间的无电离层和宽巷的组合的整周模糊度向量系数矩阵,a:基准站和流动站间的相对位置向量,b:基准站和流动站间的整周模糊度向量,基准站和流动站间的无电离层和宽巷的组合的载波相位双差测量残差,基准站和流动站间的无电离层和宽巷的组合的载波相位双差测量常数向量,基准站和流动站间的无电离层和宽巷的组合的载波相位双差测量权系数矩阵,P:基础权系数矩阵,n为解算星数,lin、col表示矩阵P的行和列,为测量噪声方差,ALC是基准站和流动站间的无电离层的相对位置向量系数矩阵,AWL是基准站和流动站间的宽巷的的相对位置向量系数矩阵,为基准站和流动站间的无电离层的载波相位双差测量常数向量,为基准站和流动站间的宽巷的载波相位双差测量常数向量,BLC是基准站和流动站间的无电离层的整周模糊度向量系数矩阵,BWL为基准站和流动站间的宽巷的整周模糊度向量系数矩阵,是基准站和流动站间的无电离层的载波相位双差测量残差,是基准站和流动站间的宽巷的载波相位双差测量残差,PLC是基准站和流动站间的无电离层的载波相位双差测量权系数矩阵,PWL是基准站和流动站间的宽巷的载波相位双差测量权系数矩阵。
10.如权利要求6所述的设备,其特征在于,所述第二确定单元具体用于:
根据所述载波L1、载波L2的双差模糊度浮点解使用LAMBDA搜索得出多组模糊度和一个表征模糊度可靠性的Ratio门限值,并根据多组模糊度和一个表征模糊度可靠性的Ratio门限值确定所述载波L1的双差整周模糊度整数解和所述载波L2的双差整周模糊度整数解。
CN201510520435.3A 2015-08-21 2015-08-21 一种实时动态差分定位方法及其设备 Active CN105158783B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510520435.3A CN105158783B (zh) 2015-08-21 2015-08-21 一种实时动态差分定位方法及其设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510520435.3A CN105158783B (zh) 2015-08-21 2015-08-21 一种实时动态差分定位方法及其设备

Publications (2)

Publication Number Publication Date
CN105158783A CN105158783A (zh) 2015-12-16
CN105158783B true CN105158783B (zh) 2018-06-29

Family

ID=54799699

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510520435.3A Active CN105158783B (zh) 2015-08-21 2015-08-21 一种实时动态差分定位方法及其设备

Country Status (1)

Country Link
CN (1) CN105158783B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109541660A (zh) * 2018-11-12 2019-03-29 广东星舆科技有限公司 一种vrs切换的方法、接收机、可读介质及云服务器

Families Citing this family (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105467415B (zh) * 2016-01-18 2017-09-29 中国民航大学 一种基于差分气压高度约束的小型无人机rtk相对定位方法
CN106840160B (zh) * 2016-12-07 2019-12-20 中国人民解放军63921部队 深空探测器无线电干涉测量差分相位整周模糊度解算方法
CN106814379B (zh) * 2017-01-19 2019-07-05 湖南北云科技有限公司 一种云rtk定位方法及系统
CN107329147B (zh) * 2017-06-08 2020-01-14 广州海格通信集团股份有限公司 卫星导航系统的定位方法和系统及地面接收系统
CN110703193B (zh) * 2018-07-10 2022-05-06 中移物联网有限公司 一种位置监测方法及装置、设备、存储介质
CN109085626B (zh) * 2018-09-05 2021-04-20 上海海积信息科技股份有限公司 一种定位方法及装置
CN109116394B (zh) * 2018-09-10 2020-08-07 中国科学院国家授时中心 一种适用于不同长度基线的实时动态定位方法
CN111142142B (zh) * 2018-11-06 2023-04-28 千寻位置网络有限公司 Rtk参考站切换的处理方法及其装置
CN111381264B (zh) * 2018-12-27 2023-08-15 北京六分科技有限公司 网络rtk中长基线模糊度固定方法和平台
CN110208835B (zh) * 2019-05-21 2023-05-05 哈尔滨工程大学 一种基于消电离层组合的跨系统紧组合差分定位方法
CN110133696A (zh) * 2019-05-27 2019-08-16 青岛杰瑞自动化有限公司 一种基于浮动平台的高精度定位方法及系统
CN110133585A (zh) * 2019-06-27 2019-08-16 江苏芯盛智能科技有限公司 双频双动态定位方法、装置、定位设备及运载工具
CN110646820B (zh) * 2019-09-20 2021-11-30 广州市中海达测绘仪器有限公司 Rtk定位数据的质检方法、装置、设备和存储介质
CN110687558B (zh) * 2019-12-11 2020-05-12 深圳微品致远信息科技有限公司 一种基于北斗动动定位的独立双基线解算验证方法
CN111290005B (zh) * 2020-03-04 2021-03-05 高维时空(北京)网络有限公司 载波相位的差分定位方法、装置、电子设备及存储介质
CN111538056B (zh) * 2020-05-27 2022-09-09 北京自动化控制设备研究所 动态精密单点定位解算方法
CN113687399B (zh) * 2021-10-26 2022-02-18 中移(上海)信息通信科技有限公司 一种定位方法、系统、终端设备及卫星
CN114397684A (zh) * 2022-01-21 2022-04-26 腾讯科技(深圳)有限公司 模糊度的固定方法及相关设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050101248A1 (en) * 2003-10-28 2005-05-12 Trimble Navigation Limited, A California Corporation Ambiguity estimation of GNSS signals for three or more carriers
CN101403792A (zh) * 2008-11-05 2009-04-08 武汉大学 动态模糊度确定方法
CN101770033A (zh) * 2010-02-08 2010-07-07 东南大学 连续运行参考站系统站间整周模糊度网络固定方法
CN103576175A (zh) * 2013-11-06 2014-02-12 西南交通大学 一种双频多星座gnss整周模糊度otf解算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7586440B2 (en) * 2004-08-25 2009-09-08 The Ritsumeikan Trust Independent positioning device and independent positioning method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050101248A1 (en) * 2003-10-28 2005-05-12 Trimble Navigation Limited, A California Corporation Ambiguity estimation of GNSS signals for three or more carriers
CN101403792A (zh) * 2008-11-05 2009-04-08 武汉大学 动态模糊度确定方法
CN101770033A (zh) * 2010-02-08 2010-07-07 东南大学 连续运行参考站系统站间整周模糊度网络固定方法
CN103576175A (zh) * 2013-11-06 2014-02-12 西南交通大学 一种双频多星座gnss整周模糊度otf解算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
GNSS 整周模糊度实时解算方法研究;李娜;《中国优秀硕士学位论文全文数据库 基础科学辑》;20140715(第07期);正文第9页倒数第3段-第12页第1段、第31页第1段-第32页第1段、第33页倒数第2段-第34页 *
一种用于长距离网络 RTK 基准站模糊度固定的非组合方法;张明等;《测绘科学技术学报》;20150228;第32卷(第1期);正文第33页左侧栏-右侧栏第2段、倒数第1段 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109541660A (zh) * 2018-11-12 2019-03-29 广东星舆科技有限公司 一种vrs切换的方法、接收机、可读介质及云服务器

Also Published As

Publication number Publication date
CN105158783A (zh) 2015-12-16

Similar Documents

Publication Publication Date Title
CN105158783B (zh) 一种实时动态差分定位方法及其设备
CN109116394A (zh) 一种适用于不同长度基线的实时动态定位方法
CN102230971B (zh) Gps多天线测姿方法
CN110531392A (zh) 一种基于ppp算法的高精度定位方法和系统
CN104102822B (zh) 一种多频gnss观测值随机特性建模方法
CN105242293B (zh) 一种全球导航卫星系统的高精度厘米级定位方法
CN106526629B (zh) 卫星导航设备及其定向方法及装置
CN108196281A (zh) 一种基于位置域曲线约束的单频动态周跳探测与修复方法
CN105425261B (zh) 基于GPS/Beidou2/INS的组合导航与定位方法
CN107703526B (zh) 基线测向方法、装置和系统
CN105158778B (zh) 多系统联合实施载波相位差分故障卫星剔除方法及其系统
CN109782313A (zh) 地面精密单点定位数据处理方法
CN1967282B (zh) 全站仪与gps单频实时动态组合测量方法及其系统
CN111751853A (zh) 一种gnss双频载波相位整周模糊度解算方法
CN105589087A (zh) 一种高精度实时卫星定位装置及其方法
CN105849589A (zh) 全球导航卫星系统、定位终端、定位方法以及记录介质
US11125886B2 (en) Method and apparatus for multipath mitigation in GNSS
CN103728641A (zh) 基于北斗系统的三天线阵高精度定向方法及系统
CN106569242A (zh) 固定参考星的gnss单差处理方法
CN110058282A (zh) 一种基于双频gnss智能手机的ppp高精度定位方法
RU2012150814A (ru) Система и способ определения курсового угла транспортного средства
CN104459722B (zh) 一种基于多余观测分量的整周模糊度可靠性检验方法
CN111538039B (zh) 未知点位基准站精确坐标确定方法
Quddus et al. Validation of map matching algorithms using high precision positioning with GPS
CN109975849B (zh) 一种基线向量的确定方法、服务器和计算机存储介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A real-time dynamic differential positioning method and device thereof

Effective date of registration: 20220826

Granted publication date: 20180629

Pledgee: China Merchants Bank Co.,Ltd. Shanghai Branch

Pledgor: SHANGHAI HIGH GAIN INFORMATION TECHNOLOGY Co.,Ltd.

Registration number: Y2022310000205

PC01 Cancellation of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20230818

Granted publication date: 20180629

Pledgee: China Merchants Bank Co.,Ltd. Shanghai Branch

Pledgor: SHANGHAI HIGH GAIN INFORMATION TECHNOLOGY Co.,Ltd.

Registration number: Y2022310000205

PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A real-time dynamic differential positioning method and its equipment

Effective date of registration: 20230829

Granted publication date: 20180629

Pledgee: China Merchants Bank Co.,Ltd. Shanghai Branch

Pledgor: SHANGHAI HIGH GAIN INFORMATION TECHNOLOGY Co.,Ltd.

Registration number: Y2023310000499