CN114296119A - 一种精密单点定位方法、装置、电子设备和存储介质 - Google Patents

一种精密单点定位方法、装置、电子设备和存储介质 Download PDF

Info

Publication number
CN114296119A
CN114296119A CN202210050395.0A CN202210050395A CN114296119A CN 114296119 A CN114296119 A CN 114296119A CN 202210050395 A CN202210050395 A CN 202210050395A CN 114296119 A CN114296119 A CN 114296119A
Authority
CN
China
Prior art keywords
frequency
satellite
epoch
data
jth
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
CN202210050395.0A
Other languages
English (en)
Other versions
CN114296119B (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.)
Guangzhou Asensing Technology Co Ltd
Original Assignee
Guangzhou Asensing 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 Guangzhou Asensing Technology Co Ltd filed Critical Guangzhou Asensing Technology Co Ltd
Priority to CN202210050395.0A priority Critical patent/CN114296119B/zh
Publication of CN114296119A publication Critical patent/CN114296119A/zh
Application granted granted Critical
Publication of CN114296119B publication Critical patent/CN114296119B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本申请实施例提供一种精密单点定位方法、装置、电子设备和存储介质，其中，方法包括：获取卫星的观测数据和卫星的差分数据，所述差分数据包括电离层延迟数据和湿对流层延迟数据；组建非组合PPP观测方程，所述非组合PPP观测方程包括电离层参数和湿对流层参数；根据所述电离层延迟数据和所述湿对流层延迟数据计算所述电离层参数的电离层延迟值和所述湿对流层参数的湿对流层延迟值；利用所述电离层延迟值和所述湿对流层延迟值对所述非组合PPP观测方程的所述电离层参数和所述湿对流层参数进行约束；求解所述非组合PPP观测方程，得到定位坐标参数信息。实施上述实施例，可以减少求解非组合PPP观测方程时的收敛速度，快速获取坐标参数信息。

Description

一种精密单点定位方法、装置、电子设备和存储介质
技术领域
本申请涉及定位技术领域，具体而言，涉及一种精密单点定位方法、装置、电子设备和计算机可读存储介质。
背景技术
目前，全球导航卫星系统(Global Navigation Satellite System，GNSS)定位技术的主要方案有实时差分(real-time kinematic，RTK)定位技术和精密单点定位(precisepoint positioning，PPP)技术。常采用的RTK技术依赖于附近的参考站，距离限制一般不能大于15公里，且需要通过网络接收差分数据。PPP技术不依赖附近参考站，且可以通过卫星信号实时接收卫星精密轨道钟差等信息，但其收敛速度慢，通常需要20-30分钟的连续观测后才能获得厘米级的定位信息，数十分钟的收敛时间严重影响该技术的应用。
发明内容
本申请实施例的目的在于提供一种精密单点定位方法、装置、电子设备和计算机可读存储介质，缩短PPP技术的收敛时间，提高收敛速度，快速得到高精度位置坐标信息。
第一方面，本申请实施例提供了一种精密单点定位方法，包括：
获取卫星的观测数据和卫星的差分数据，所述差分数据包括电离层延迟数据和湿对流层延迟数据；
组建非组合PPP观测方程，所述非组合PPP观测方程包括电离层参数和湿对流层参数；
根据所述电离层延迟数据和所述湿对流层延迟数据计算所述电离层参数的电离层延迟值和所述湿对流层参数的湿对流层延迟值；
利用所述电离层延迟值和所述湿对流层延迟值对所述非组合PPP观测方程的所述电离层参数和所述湿对流层参数进行约束；
求解所述非组合PPP观测方程，得到定位坐标参数信息。
在上述实现过程中，与现有技术不同的是，获取的差分数据包括电离层延迟数据和湿对流层延迟数据，非组合PPP观测方程包括电离层参数和湿对流层参数；根据电离层延迟数据和湿对流层数据计算电离层参数的电离层延迟值和湿对流层参数的湿对流层延迟值，利用电离层延迟值和湿对流层延迟值对非组合PPP观测方程的电离层参数和湿对流层参数进行约束。基于上述实施方式，能够提高求解非组合PPP观测方程中的收敛速度，从而快速地获取定位坐标参数信息。
进一步地，所述非组合PPP观测方程为：
Figure BDA0003473784780000021
其中，Ps,j为卫星S的第j(j＝1,2)频率的伪距观测数据；
Ls,j为所述卫星S的第j(j＝1,2)频率的载波相位观测数据；
Figure BDA0003473784780000022
为接收机r的天线相位中心和所述卫星S相位中心间的几何距离；
c为真空中的光速；
dtr,j为所述接收机r的第j(j＝1,2)频率对应的钟差；
dts为所述卫星S的精密轨道钟差；
T为第j(j＝1,2)频率的所述湿对流层延迟参数；
Figure BDA0003473784780000023
为第j(j＝1,2)频率的所述电离层延迟参数；
γ为第j(j＝1,2)频率的平方比
Figure BDA0003473784780000024
fj为第j(j＝1,2)频率的频率值；
br,j为接收机r的伪距硬件；
Figure BDA0003473784780000031
为所述卫星S的延迟数据；
Br,j为所述接收机r的相位硬件延迟数据；
Figure BDA0003473784780000032
为所述卫星S的相位硬件延迟数据；
Figure BDA0003473784780000033
为所述卫星S的第j(j＝1,2)频率的波长；
Figure BDA0003473784780000034
为所述卫星S的第j(j＝1,2)频率的相位模糊度；
εp,j为所述卫星S的第j(j＝1,2)频率的伪距观测数据的观测噪声；
εL,j为所述卫星S的第j(j＝1,2)频率的载波相位观测数据的观测噪声。
进一步地，所述电离层延迟值和所述湿对流层延迟值的计算公式为：
Figure BDA0003473784780000035
其中，
Figure BDA0003473784780000036
为所述卫星S的第j(j＝1,2)频率的电离层延迟数据；
Figure BDA0003473784780000037
为所述卫星S的第j(j＝1,2)频率的湿对流层延迟数据；
εI,j为j(j＝1,2)频率的所述电离层延迟数据的精度误差；
εT,j为j(j＝1,2)频率的所述湿对流层延迟数据的精度误差；
Figure BDA0003473784780000038
为所述电离层延迟值；
T为所述湿对流层延迟值。
进一步地，所述求解所述非组合PPP观测方程，得到定位坐标参数信息的步骤，包括：
对移动星和参考星的星间单差宽巷模糊度固定成单差宽巷整数模糊度；
采用lambda搜索方法固定移动星和参考星的星间单差窄巷模糊度，得到单差窄巷整数模糊度；
利用所述单差宽巷整数模糊度和单差窄巷整数模糊度更新所述非组合PPP观测方程中的模糊度，得到模糊度固定的所述坐标参数信息。
进一步地，在所述组建非组合PPP观测方程的步骤之前，还包括：
对所述观测数据进行预处理。
进一步地，所述对所述观测数据进行预处理的步骤，包括：
根据每个所述卫星的伪距观测数据获取每个所述卫星的伪距频间偏差；
根据所有所述卫星的伪距频间偏差剔除所有所述卫星的伪距观测数据中的异常伪距观测数据。
进一步地，所述对所述观测数据进行预处理的步骤，包括：
根据每个所述卫星的伪距观测数据和多普勒观测数据获取每个所述卫星的多普勒探测伪距检验值；
根据每个所述卫星的多普勒探测伪距检验值剔除所有所述卫星的伪距观测数据中的异常伪距观测数据。
进一步地，所述根据每个所述卫星的伪距观测数据获取每个所述卫星的伪距频间偏差的步骤，包括：
根据每个所述卫星的伪距观测数据获取每个所述卫星的第一频率对应的伪距观测数据和第二频率对应的伪距观测数据；
获取每个所述卫星的第一频率对应的伪距观测数据和第二频率对应的伪距观测数据的差值，得到每个所述卫星的伪距频间偏差。
进一步地，所述根据所有所述卫星的伪距频间偏差剔除所有所述卫星的伪距观测数据中的异常伪距观测数据的步骤，包括：
获取所有所述卫星的伪距频间偏差的中位数；
计算每个所述卫星的伪距频间偏差和所有所述卫星的伪距频间偏差的中位数的差值，得到每个所述卫星的第一检验值；
剔除所述第一检验值大于第一预设阈值的所述卫星对应的伪距观测数据。
进一步地，通过以下公式获取所述多普勒探测伪距检验值：
ΔPi,j,k＝Pi,j,k-Pi,j,k-1-0.5*λi,j(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；
Di,j,k为k历元卫星i的j(j＝1,2)频率的多普勒观测数据；
Pi,j,k-1为k-1历元卫星i的j(j＝1,2)频率的伪距观测数据；
Di,j,k-1为k-1历元卫星i的j(j＝1,2)频率的多普勒观测数据；
λi,j为卫星i的第j(j＝1,2)频率的波长值；
tk为k历元的时间；
tk-1为历元的时间；
ΔPi,j,k为k历元卫星i的j(j＝1,2)频率的多普勒探测伪距检验值。
进一步地，所述对所述观测数据进行预处理的步骤，还包括：
根据以下公式对每个所述卫星的所述伪距观测数据进行平滑处理，得到多普勒平滑伪距；
PSi,j,k＝ωkPi,j,k+(1-ωk)(PSi,j,k-1+0.5*λi,j(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，其中，PSi,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒平滑伪距；
ωk为k历元对应的平滑因子；
Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；
PSi,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒平滑伪距；
λi,j为卫星i的第j(j＝1,2)频率的波长值；
Di,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒观测数据；
Di,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒观测数据；
tk为k历元的时间；
tk-1为k-1历元的时间。
进一步地，所述对所述卫星的观测数据进行预处理的步骤，包括：
判断每个载波相位观测数据的整周模糊度是否发生跳变；
若是，对整周模糊度跳变的所述载波相位观测数据进行标记。
进一步地，所述判断每个载波相位观测数据的整周模糊度是否发生跳变的步骤包括：
采用历元间差分算法对所述载波相位观测数据进行单频周跳探测，得到跳变结果。
进一步地，在所述采用历元间差分算法对所述载波相位观测数据进行单频周跳探测，得到跳变结果的步骤之前，还包括：利用以下公式获取所述载波相位观测数据的周跳探测值：
ΔLi,j,k＝Li,j,k-Li,j,k-1-0.5*(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的周跳探测值；
Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Di,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒观测数据；
Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Di,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒观测数据；
tk为k历元的时间；
tk-1为k-1历元的时间；
若|ΔLi,j,k|大于第二预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；
若|ΔLi,j,k|小于或等于所述第二预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
进一步地，所述采用历元间差分算法对所述载波相位观测数据进行单频周跳探测，得到跳变结果的步骤，包括：
利用以下公式对所述载波相位观测数据进行单频周跳探测：
ΔLi,j,k＝Li,j,k-Li,j,k-1-(Pi,j,k-Pi,j,k-1)/λi,j；
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的单频周跳探测值；
Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；
Pi,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的伪距观测数据；
λi,j为卫星i的第j(j＝1,2)频率的波长值；
若|ΔLi,j,k|大于第三预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；
若|ΔLi,j,k|小于或等于所述第三预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
进一步地，所述采用历元间差分算法对所述载波相位观测数据进行单频周跳探测，得到跳变结果的步骤，包括：
利用以下公式对所述载波相位观测数据进行单频周跳探测：
ΔLi,j,k＝Li,j,k-3*Li,j,k-1+3*Li,j,k-2-Li,j,k-3：
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的单频周跳探测值；
Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Li,j,k-2为k-2历元所述卫星i的第j(j＝1,2)频率的载波相位观测数据；
Li,j,k-3为k-3历元所述卫星i的第j(j＝1,2)频率的载波相位观测数据；
λi,j为卫星i的第j(j＝1,2)频率的波长值；
若|ΔLi,j,k|大于第四预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；
若|ΔLi,j,k|小于或等于所述第四预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
进一步地，所述观测数据通过导航卫星发送到接收机，所述差分数据通过通信卫星发送到所述接收机。
第二方面，本申请实施例还提供一种精密单点定位装置，包括：
获取模块，用于获取卫星的观测数据和卫星的差分数据，所述差分数据包括电离层延迟数据和湿对流层延迟数据；
方程组建模块，用于组建非组合PPP观测方程，所述非组合PPP观测方程包括电离层参数和湿对流层参数；
延迟值获取模块，用于根据所述电离层延迟数据和所述湿对流层延迟数据计算所述电离层参数的电离层延迟值和所述湿对流层参数的湿对流层延迟值；
约束模块，用于利用所述电离层延迟值和所述湿对流层延迟值对所述非组合PPP观测方程的所述电离层参数和所述湿对流层参数进行约束；
求解模块，用于求解所述非组合PPP观测方程，得到定位坐标参数信息。
第三方面，本申请实施例提供的一种电子设备，包括：存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如第一方面任一项所述的方法的步骤。
第四方面，本申请实施例提供的一种计算机可读存储介质，所述计算机可读存储介质上存储有指令，当所述指令在计算机上运行时，使得所述计算机执行如第一方面任一项所述的方法。
本申请公开的其他特征和优点将在随后的说明书中阐述，或者，部分特征和优点可以从说明书推知或毫无疑义地确定，或者通过实施本申请公开的上述技术即可得知。
为使本申请的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本申请实施例的技术方案，下面将对本申请实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本申请的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
图1为本申请实施例提供的精密单点定位方法的流程示意图；
图2为本申请实施例提供的求解非组合PPP观测方程的流程示意图；
图3为本申请实施例提供的对观测数据进行预处理的流程示意图；
图4为本申请实施例提供的获取每个卫星的多普勒探测伪距检验值的流程示意图；
图5为本申请实施例提供的剔除所有卫星的伪距观测数据中的异常伪距观测数据的流程示意图；
图6为本申请实施例提供的对观测数据进行预处理的另一流程示意图；
图7为本申请实施例提供的精密单点定位装置的结构组成结构图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行描述。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本申请的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。
实施例1
参见图1，本申请实施例提供一种精密单点定位方法，包括：
S1：获取卫星的观测数据和卫星的差分数据，差分数据包括电离层延迟数据和湿对流层延迟数据；
S2：组建非组合PPP观测方程，该非组合PPP观测方程包括电离层参数和湿对流层参数；
S3：根据电离层延迟数据和湿对流层延迟数据计算电离层参数的电离层延迟值和湿对流层参数的湿对流层延迟值；
S4：利用电离层延迟值和湿对流层延迟值对非组合PPP观测方程的电离层参数和湿对流层参数进行约束；
S5：求解非组合PPP观测方程，得到定位坐标参数信息。
与现有技术不同的是，获取的差分数据包括电离层延迟数据和湿对流层延迟数据，非组合PPP观测方程包括电离层参数和湿对流层参数；根据电离层延迟数据和湿对流层数据计算电离层参数的电离层延迟值和湿对流层参数的湿对流层延迟值，利用电离层延迟值和湿对流层延迟值对非组合PPP观测方程的电离层参数和湿对流层参数进行约束。基于上述实施方式，能够提高求解非组合PPP观测方程中的收敛速度，从而快速地获取定位坐标参数信息。
在一种可能的实施方式中，非组合PPP观测方程为：
Figure BDA0003473784780000101
其中，Ps,j为卫星S的第j(j＝1,2)频率的伪距观测数据；Ls,j为卫星S的第j(j＝1,2)频率的载波相位观测数据；
Figure BDA0003473784780000111
为接收机r的天线相位中心和卫星S相位中心间的几何距离；c为真空中的光速；dtr,j为接收机r的第j(j＝1,2)频率对应的钟差；dts为卫星S的精密轨道钟差；T为第j(j＝1,2)频率的湿对流层延迟参数；
Figure BDA0003473784780000112
为第j(j＝1,2)频率的电离层延迟参数；γ为第j(j＝1,2)频率的平方比
Figure BDA0003473784780000113
fj为第j(j＝1,2)频率的频率值；br,j为接收机r的伪距硬件；
Figure BDA0003473784780000114
为卫星S的延迟数据；Br,j为接收机r的相位硬件延迟数据；
Figure BDA0003473784780000115
为卫星S的相位硬件延迟数据；
Figure BDA0003473784780000116
为卫星S的第j(j＝1,2)频率的波长；
Figure BDA0003473784780000117
为卫星S的第j(j＝1,2)频率的相位模糊度；εp,j为卫星S的第j(j＝1,2)频率的伪距观测数据的观测噪声；εL,j为卫星S的第j(j＝1,2)频率的载波相位观测数据的观测噪声。
上述实施例中，将卫星的观测数据和差分数据进行组建，得到非组合PPP观测方程，由于非组合PPP观测方程中包括多种观测数据和差分数据，因此，基于上述公式获得的定位信息相比现有技术获得的定位信息更加准确。
在一种可能的实施方式中，电离层延迟值和湿对流层延迟值的计算公式为：
Figure BDA0003473784780000118
其中，
Figure BDA0003473784780000119
为卫星S的第j(j＝1,2)频率的电离层延迟数据；
Figure BDA00034737847800001110
为卫星S的第j(j＝1,2)频率的湿对流层延迟数据；εI,j为j(j＝1,2)频率的电离层延迟数据的精度误差；εT,j为j(j＝1,2)频率的湿对流层延迟数据的精度误差；
Figure BDA00034737847800001111
为电离层延迟值；T为湿对流层延迟值。
基于上述实施方式，基于电离层延迟数据和湿对流层延迟数据，分别获得电离层延迟值和湿对流层延迟值，基于上述两个延迟值，可以实现对非组合PPP观测方程的约束，能够减少非组合PPP观测方程的收敛时间。
参见图2，本申请实施例提供一种求解非组合PPP观测方程的实现方式，过程包括：
S51：对移动星和参考星的星间单差宽巷模糊度固定成单差宽巷整数模糊度；
S52：采用lambda搜索方法固定移动星和参考星的星间单差窄巷模糊度，得到单差窄巷整数模糊度；
S53：利用单差宽巷整数模糊度和单差窄巷整数模糊度更新非组合PPP观测方程中的模糊度，得到模糊度固定的坐标参数信息。
上述实施例中，参考星为所有数据对应的卫星中高度角最高的卫星，其他卫星即是移动星，移动星的模糊度减去参考星的模糊度即是单差模糊度。
优选地，通过以下公式计算移动星和参考星的星间单差宽巷模糊度：WL＝N0,1-N0,2-(Nr,1-Nr,2)，N**分别是移动星和参考星的第1和第2频率的相位模糊度。再通过以下公式WL0＝N0,1-N0,2-(Nr,1-Nr,2)作为模糊度虚拟观测方程约束更新非组合PPP观测方程中的模糊度参数。
基于上述实施方式，可以获取固定模糊度固定的坐标参数信息。
在一种可能的实施方式中，在S5之前，还包括：
对观测数据进行预处理。
本申请实施例提供多种数据预处理方法，需要说明的是，可以使用以下的一种或多种方法对卫星的观测数据和卫星的差分数据进行预处理，本申请对此不做限制。
在一种可能的实施方式中，对观测数据进行预处理的步骤，包括：根据每个卫星的伪距观测数据获取每个卫星的伪距频间偏差；根据所有卫星的伪距频间偏差剔除所有卫星的伪距观测数据中的异常伪距观测数据。
参见图3，在一种可能的实施方式中，对观测数据进行预处理的步骤，包括：
S61：根据每个卫星的伪距观测数据和多普勒观测数据获取每个卫星的多普勒探测伪距检验值；
S62：根据每个卫星的多普勒探测伪距检验值剔除所有卫星的伪距观测数据中的异常伪距观测数据。
粗差数据对定位具有致命性的影响，由于卫星信号的传输距离较长，信号在传输的过程中会发生变化，所获得的数据中会存在粗差数据，因此，获取每个卫星的伪距频间偏差，根据所有导航卫星的伪距频间偏差剔除所有卫星的伪距观测数据中的异常伪距观测数据能够消除终端本身的系统偏差。
参见图4，在一种可能的实施方式中，S61包括以下子步骤：
S611：根据每个卫星的伪距观测数据获取每个卫星的第一频率对应的伪距观测数据和第二频率对应的伪距观测数据；
S612：获取每个卫星的第一频率对应的伪距观测数据和第二频率对应的伪距观测数据的差值，得到每个卫星的伪距频间偏差。
上述实现过程中，对每个卫星的第一频率的数据和第二频率的数据进行采样，根据第一频率对应的伪距观测数据和第二频率对应的伪距观测数据的差值作为每个导航卫星的伪距频间偏差。
参见图5，在一种可能的实施方式中，S62包括以下子步骤：
S621：获取所有卫星的伪距频间偏差的中位数；
S622：计算每个卫星的伪距频间偏差和所有卫星的伪距频间偏差的中位数的差值，得到每个卫星的第一检验值；
S623：剔除第一检验值大于第一预设阈值的卫星对应的伪距观测数据。
采样信号频率间可能会存在系统偏差，需要扣除系统偏差的影响，中位数能够衡量采集到的所有卫星数据的平均水平，计算每个卫星的伪距频间偏差和所有卫星的伪距频间偏差的中位数的差值，得到第一检验值，将该第一检验值和第一预设阈值进行比较，根据比较结果能够获取到异常的伪距观测数据，进一步将该卫星数据删除。
进一步地，本申请实施例提供一种公式，能够快速地对该伪距观测数据是否是异常卫星数据进行判断，该公式如下：
ΔPi＝dPi-dP0；
其中，dPi为颗卫星i的伪距频间偏差，dP0为所有卫星的伪距频间偏差的中位数，ΔPi为卫星i和所有卫星的伪距频间偏差的中位数的差值，δ1为第一预设阈值；若|ΔPi|≤δ1；则判定该卫星的伪距观测数据发生正常，若|ΔPi|＞δ1；则判定该卫星的伪距观测数据发生异常，将该卫星的伪距观测数据删除，δ1为第一预设阈值。示例性地，δ1可以取为50米。
本申请实施例提供第二种预处理的方式，该方法是针对异常伪距观测数据。方法如下：
通过以下公式获取多普勒探测伪距检验值：
ΔPi,j,k＝Pi,j,k-Pi,j,k-1-0.5*λi,j(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；Di,j,k为k历元卫星i的j(j＝1,2)频率的多普勒观测数据；Pi,j,k-1为k-1历元卫星i的j(j＝1,2)频率的伪距观测数据；Di,j,k-1为k-1历元卫星i的j(j＝1,2)频率的多普勒观测数据；λi,j为卫星i的第j(j＝1,2)频率的波长值；tk为k历元的时间；tk-1为历元的时间；ΔPi,j,k为k历元卫星i的j(j＝1,2)频率的多普勒探测伪距检验值。
示例性地，若|ΔPi,j,k|≤δ2，则认为k历元中的卫星i第j频率的伪距观测数据和多普勒观测数据正常，若|ΔPi,j,k|＞δ2，则认为k历元中的卫星i第j(j＝1,2)频率的伪距观测数据和多普勒观测数据异常，将k历元中的卫星i第j频率的伪距观测数据和多普勒观测数据剔除。示例性地，δ2＝50米。
本申请实施例提供第三种预处理的方法，包括：
根据以下公式对每个卫星的伪距观测数据进行平滑处理，得到多普勒平滑伪距；
PSi,j,k＝ωkPi,j,k+(1-ωk)(PSi,j,k-1+0.5*λi,j(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，PSi,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒平滑伪距；ωk为k历元对应的平滑因子；Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；PSi,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒平滑伪距；λi,j为卫星i的第j(j＝1,2)频率的波长值；Di,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒观测数据；Di,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒观测数据；tk为k历元的时间；tk-1为k-1历元的时间。示例性地，当k大于60时，ωk＝1/60.0，否则，
Figure BDA0003473784780000151
通过上述公式，可以利用多普勒观测数据，对伪距观测数据进行处理，将多普勒平滑伪距代替原来的多普勒伪距离，从而减少多普勒观测数据中的误差，提高采样精度。
需要说明的是，如果是上述两种实施方式进行结合，那么可以先对多普勒观测数据进行剔除，再计算多普勒平滑伪距。
参见图6，本申请实施例提供第四种预处理的方法，包括以下步骤：
S63：判断每个载波相位观测数据的整周模糊度是否发生跳变；若是，执行S64；
S64：对整周模糊度跳变的载波相位观测数据进行标记。
上述实施例中，标记的目的是为了后续更新模糊度参数，具体参见S53。
优选地，可以采用历元间差分算法对载波相位观测数据进行单频周跳探测，得到跳变结果。
在一种可能的实施方式中，在采用历元间差分算法对载波相位观测数据进行单频周跳探测，得到跳变结果的步骤之前，还包括：利用以下公式获取载波相位观测数据的周跳探测值：
ΔLi,j,k＝Li,j,k-Li,j,k-1-0.5*(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的周跳探测值；Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Di,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒观测数据；Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Di,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒观测数据；tk为k历元的时间；tk-1为k-1历元的时间；若|ΔLi,j,k|大于第二预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；若|ΔLi,j,k|小于或等于第二预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
若|ΔLi,j,k|≤δ3，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测值未发生跳变；|ΔLi,j,k|＞δ3，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测值发生跳变。其中，δ3为第二预设阈值，优选地，第二预设阈值为0.5周。
在一种可能的实施方式中，采用历元间差分算法对载波相位观测数据进行单频周跳探测，得到跳变结果的步骤，包括：
利用以下公式对载波相位观测数据进行单频周跳探测：
ΔLi,j,k＝Li,j,k-Li,j,k-1-(Pi,j,k-Pi,j,k-1)/λi,j；
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的单频周跳探测值；Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；Pi,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的伪距观测数据；λi,j为卫星i的第j(j＝1,2)频率的波长值；若|ΔLi,j,k|大于第三预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；若|ΔLi,j,k|小于或等于第三预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
上述过程可以利用如下公式表示：
若|ΔLi,j,k|≤δ4，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测值未发生跳变；若|ΔLi,j,k|＞δ4，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测值发生跳变，其中，δ4为第三预设阈值，优选地，第三预设阈值取为100周。
在一种可能的实施方式中，还可以利用下述方式对载波相位观测值进行单频周跳探测：
ΔLi,j,k＝Li,j,k-3*Li,j,k-1+3*Li,j,k-2-Li,j,k-3：
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的单频周跳探测值；Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Li,j,k-2为k-2历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Li,j,k-3为k-3历元卫星i的第j(j＝1,2)频率的载波相位观测数据；λi,j为卫星i的第j频率的波长值；若|ΔLi,j,k|大于第四预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；若|ΔLi,j,k|小于或等于第四预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
上述过程可以用公式表示如下：若|ΔLi,j,k|≤δ5，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测值未发生跳变；若|ΔLi,j,k|＞δ5，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测值发生跳变，其中，δ5为第四预设阈值，优选地，第四预设阈值取为3周。
可以理解的是，如果载波相位观测值发生周跳，则将该载波相位观测值进行标记，在解非组合PPP观测方程中重置模糊度参数。
在一种可能的实施方式中，观测数据通过导航卫星发送到接收机，差分数据通过通信卫星发送到接收机。
现有技术中，常采用的RTK技术依赖于附近的参考站，距离限制一般不能大于15公里，且需要通过网络接收差分数据。通过上述实施方式，可以不依赖参考站和网络通讯，在1分钟内提供厘米级的定位信息。
实施例2
参见图7，本申请实施例提供一种精密单点定位装置，包括：
获取模块1，用于获取卫星的观测数据和卫星的差分数据，差分数据包括电离层延迟数据和湿对流层延迟数据；
方程组建模块2，用于组建非组合PPP观测方程，非组合PPP观测方程包括电离层参数和湿对流层参数；
延迟值获取模块3，用于根据电离层延迟数据和湿对流层延迟数据计算电离层参数的电离层延迟值和湿对流层参数的湿对流层延迟值；
约束模块4，用于利用电离层延迟值和湿对流层延迟值对非组合PPP观测方程的电离层参数和湿对流层参数进行约束；
求解模块5，用于求解非组合PPP观测方程，得到定位坐标参数信息。
在一种可能的实施方式中，方程组建模块2还用于组建以下非组合PPP观测方程：
Figure BDA0003473784780000191
其中，Ps,j为卫星S的第j(j＝1,2)频率的伪距观测数据；Ls,j为卫星S的第j(j＝1,2)频率的载波相位观测数据；
Figure BDA0003473784780000192
为接收机r的天线相位中心和卫星S相位中心间的几何距离；c为真空中的光速；dtr,j为接收机r的第j(j＝1,2)频率对应的钟差；dts为卫星S的精密轨道钟差；T为第j(j＝1,2)频率的湿对流层延迟参数；
Figure BDA0003473784780000193
为第j(j＝1,2)频率的电离层延迟参数；γ为第j(j＝1,2)频率的平方比
Figure BDA0003473784780000194
fj为第j(j＝1,2)频率的频率值；br,j为接收机r的伪距硬件；
Figure BDA0003473784780000195
为卫星S的延迟数据；Br,j为接收机r的相位硬件延迟数据；
Figure BDA0003473784780000196
为卫星S的相位硬件延迟数据；
Figure BDA0003473784780000197
为卫星S的第j(j＝1,2)频率的波长；
Figure BDA0003473784780000198
为卫星S的第j(j＝1,2)频率的相位模糊度；εp,j为卫星S的第j(j＝1,2)频率的伪距观测数据的观测噪声；εL,j为卫星S的第j(j＝1,2)频率的载波相位观测数据的观测噪声。
在一种可能的实施方式中，延迟值获取模块3还用于通过以下公式获取电离层延迟值和湿对流层延迟值：
Figure BDA0003473784780000199
其中，
Figure BDA00034737847800001910
为卫星S的第j(j＝1,2)频率的电离层延迟数据；
Figure BDA00034737847800001911
为卫星S的第j(j＝1,2)频率的湿对流层延迟数据；εI,j为j(j＝1,2)频率的电离层延迟数据的精度误差；εT,j为j(j＝1,2)频率的湿对流层延迟数据的精度误差；
Figure BDA00034737847800001912
为电离层延迟值；T为湿对流层延迟值。
在一种可能的实施方式中，求解模块5还用于对移动星和参考星的星间单差宽巷模糊度固定成单差宽巷整数模糊度；采用lambda搜索方法固定移动星和参考星的星间单差窄巷模糊度，得到单差窄巷整数模糊度；利用单差宽巷整数模糊度和单差窄巷整数模糊度更新非组合PPP观测方程中的模糊度，得到模糊度固定的坐标参数信息。
在一种可能的实施方式中，装置还包括预处理模块，用于对观测数据进行预处理。
在一种可能的实施方式中，预处理模块还用于根据每个卫星的伪距观测数据获取每个卫星的伪距频间偏差；根据所有卫星的伪距频间偏差剔除所有卫星的伪距观测数据中的异常伪距观测数据。
在一种可能的实施方式中，预处理模块还用于根据每个卫星的伪距观测数据和多普勒观测数据获取每个卫星的多普勒探测伪距检验值；根据每个卫星的多普勒探测伪距检验值剔除所有卫星的伪距观测数据中的异常伪距观测数据。
在一种可能的实施方式中，预处理模块还用于根据每个卫星的伪距观测数据获取每个卫星的第一频率对应的伪距观测数据和第二频率对应的伪距观测数据；获取每个卫星的第一频率对应的伪距观测数据和第二频率对应的伪距观测数据的差值，得到每个卫星的伪距频间偏差。
在一种可能的实施方式中，预处理模块还用于获取所有卫星的伪距频间偏差的中位数；计算每个卫星的伪距频间偏差和所有卫星的伪距频间偏差的中位数的差值，得到每个卫星的第一检验值；剔除第一检验值大于第一预设阈值的卫星对应的伪距观测数据。
在一种可能的实施方式中，预处理模块还用于通过以下公式获取多普勒探测伪距检验值：ΔPi,j,k＝Pi,j,k-Pi,j,k-1-0.5*λi,j(-Di,j,k-Di,j,k-1)/(tk-tk-1)；其中，Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；Di,j,k为k历元卫星i的j(j＝1,2)频率的多普勒观测数据；Pi,j,k-1为k-1历元卫星i的j(j＝1,2)频率的伪距观测数据；Di,j,k-1为k-1历元卫星i的j(j＝1,2)频率的多普勒观测数据；λi,j为卫星i的第j(j＝1,2)频率的波长值；tk为k历元的时间；tk-1为历元的时间；ΔPi,j,k为k历元卫星i的j(j＝1,2)频率的多普勒探测伪距检验值。
在一种可能的实施方式中，预处理模块还用于根据以下公式对每个卫星的伪距观测数据进行平滑处理，得到多普勒平滑伪距：
PSi,j,k＝ωkPi,j,k+(1-ωk)(PSi,j,k-1+0.5*λi,j(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，PSi,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒平滑伪距；ωk为k历元对应的平滑因子；Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；PSi,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒平滑伪距；λi,j为卫星i的第j(j＝1,2)频率的波长值；Di,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒观测数据；Di,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒观测数据；tk为k历元的时间；tk-1k-1历元的时间。
在一种可能的实施方式中，预处理模块还用于判断每个载波相位观测数据的整周模糊度是否发生跳变；若是，对整周模糊度跳变的载波相位观测数据进行标记。
在一种可能的实施方式中，预处理模块还用于采用历元间差分算法对载波相位观测数据进行单频周跳探测，得到跳变结果。
在一种可能的实施方式中，预处理模块还用于利用以下公式获取载波相位观测数据的周跳探测值：
ΔLi,j,k＝Li,j,k-Li,j,k-1-0.5*(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的周跳探测值；Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Di,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒观测数据；Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Di,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒观测数据；tk为k历元的时间；tk-1为k-1历元的时间；若ΔLi,j,k大于第二预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；若ΔLi,j,k小于或等于第二预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
在一种可能的实施方式中，预处理模块还用于利用以下公式对载波相位观测数据进行单频周跳探测：
ΔLi,j,k＝Li,j,k-Li,j,k-1-(Pi,j,k-Pi,j,k-1)/λi,j；
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的单频周跳探测值；Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；Pi,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的伪距观测数据；λi,j为卫星i的第j(j＝1,2)频率的波长值；若ΔLi,j,k大于第三预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；若ΔLi,j,k小于或等于第三预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
在一种可能的实施方式中，预处理模块还用于利用以下公式对载波相位观测数据进行单频周跳探测：
ΔLi,j,k＝Li,j,k-3*Li,j,k-1+3*Li,j,k-2-Li,j,k-3；
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的单频周跳探测值；Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Li,j,k-2为k-2历元卫星i的第j(j＝1,2)频率的载波相位观测数据；Li,j,k-3为k-3历元卫星i的第j(j＝1,2)频率的载波相位观测数据；λi,j为卫星i的第j频率的波长值；若ΔLi,j,k大于第四预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；若ΔLi,j,k小于或等于第四预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
在一种可能的实施方式中，观测数据通过导航卫星发送到接收机，差分数据通过通信卫星发送到接收机。
本申请实施例还提供一种电子设备，电子设备可以包括处理器、通信接口、存储器和至少一个通信总线。其中，通信总线用于实现这些组建直接的连接通信。其中，本申请实施例中电子设备的通信接口用于与其他节点设备进行信令或数据的通信。处理器可以是一种集成电路芯片，具有信号的处理能力。
上述的处理器可以是通用处理器，包括中央处理器(Central Processing Unit，CPU)、网络处理器(Network Processor，NP)等；还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组建。可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
存储器可以是，但不限于，随机存取存储器(Random Access Memory，RAM)，只读存储器(Read Only Memory，ROM)，可编程只读存储器(Programmable Read-Only Memory，PROM)，可擦除只读存储器(Erasable Programmable Read-Only Memory，EPROM)，电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory，EEPROM)等。存储器中存储有计算机可读取指令，当计算机可读取指令由处理器执行时，电子设备可以执行上述方法实施例涉及的各个步骤。
可选地，电子设备还可以包括存储控制器、输入输出单元。
存储器、存储控制器、处理器、外设接口、输入输出单元各元件相互之间直接或间接地电性连接，以实现数据的传输或交互。例如，这些元件相互之间可通过一条或多条通信总线实现电性连接。处理器用于执行存储器中存储的可执行模块，例如电子设备包括的软件功能模块或计算机程序。
输入输出单元用于提供给用户创建任务以及为该任务创建启动可选时段或预设执行时间以实现用户与服务器的交互。输入输出单元可以是，但不限于，鼠标和键盘等。
本申请实施例还提供一种计算机可读存储介质，计算机可读存储介质上存储有指令，当指令在计算机上运行时，计算机程序被处理器执行时实现方法实施例的方法，为避免重复，此处不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本申请的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。
另外，在本申请各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。
功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-Only Memory，ROM)、随机存取存储器(Random Access Memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上仅为本申请的实施例而已，并不用于限制本申请的保护范围，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
以上，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以权利要求的保护范围为准。
需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、物品或者设备中还存在另外的相同要素。

Claims (20)

1.一种精密单点定位方法，其特征在于，包括：
获取卫星的观测数据和卫星的差分数据，所述差分数据包括电离层延迟数据和湿对流层延迟数据；
组建非组合PPP观测方程，所述非组合PPP观测方程包括所述观测数据、电离层参数和湿对流层参数；
根据所述电离层延迟数据和所述湿对流层延迟数据计算所述电离层参数的电离层延迟值和所述湿对流层参数的湿对流层延迟值；
利用所述电离层延迟值和所述湿对流层延迟值对所述非组合PPP观测方程的所述电离层参数和所述湿对流层参数进行约束；
求解所述非组合PPP观测方程，得到定位坐标参数信息。
2.根据权利要求1所述的精密单点定位方法，其特征在于，所述非组合PPP观测方程为：
Figure FDA0003473784770000011
其中，Ps,j为卫星S的第j(j＝1,2)频率的伪距观测数据；
Ls,j为所述卫星S的第j(j＝1,2)频率的载波相位观测数据；
Figure FDA0003473784770000012
为接收机r的天线相位中心和所述卫星S相位中心间的几何距离；
c为真空中的光速；
dtr,j为所述接收机r的第j(j＝1,2)频率对应的钟差；
dts为所述卫星S的精密轨道钟差；
T为第j(j＝1,2)频率的所述湿对流层延迟参数；
Figure FDA0003473784770000013
为第j(j＝1,2)频率的所述电离层延迟参数；
γ为第j(j＝1,2)频率的平方比
Figure FDA0003473784770000014
fj为第j(j＝1,2)频率的频率值；
br,j为接收机r的伪距硬件；
Figure FDA0003473784770000021
为所述卫星S的延迟数据；
Br,j为所述接收机r的相位硬件延迟数据；
Figure FDA0003473784770000022
为所述卫星S的相位硬件延迟数据；
Figure FDA0003473784770000023
为所述卫星S的第j(j＝1,2)频率的波长；
Figure FDA0003473784770000024
为所述卫星S的第j(j＝1,2)频率的相位模糊度；
εp,j为所述卫星S的第j(j＝1,2)频率的伪距观测数据的观测噪声；
εL,j为所述卫星S的第j(j＝1,2)频率的载波相位观测数据的观测噪声。
3.根据权利要求1所述的精密单点定位方法，其特征在于，所述电离层延迟值和所述湿对流层延迟值的计算公式为：
Figure FDA0003473784770000025
其中，
Figure FDA0003473784770000026
为所述卫星S的第j(j＝1,2)频率的电离层延迟数据；
Figure FDA0003473784770000027
为所述卫星S的第j(j＝1,2)频率的湿对流层延迟数据；
εI,j为j(j＝1,2)频率的所述电离层延迟数据的精度误差；
εT,j为j(j＝1,2)频率的所述湿对流层延迟数据的精度误差；
Figure FDA0003473784770000028
为所述电离层延迟值；
T为所述湿对流层延迟值。
4.根据权利要求1至3任一项所述的精密单点定位方法，其特征在于，所述求解所述非组合PPP观测方程，得到定位坐标参数信息的步骤，包括：
对移动星和参考星的星间单差宽巷模糊度固定成单差宽巷整数模糊度；
采用lambda搜索方法固定移动星和参考星的星间单差窄巷模糊度，得到单差窄巷整数模糊度；
利用所述单差宽巷整数模糊度和单差窄巷整数模糊度更新所述非组合PPP观测方程中的模糊度，得到模糊度固定的所述坐标参数信息。
5.根据权利要求1所述的精密单点定位方法，其特征在于，在所述组建非组合PPP观测方程的步骤之前，还包括：
对所述观测数据进行预处理。
6.根据权利要求5所述的精密单点定位方法，其特征在于，所述对所述观测数据进行预处理的步骤，包括：
根据每个所述卫星的伪距观测数据获取每个所述卫星的伪距频间偏差；
根据所有所述卫星的伪距频间偏差剔除所有所述卫星的伪距观测数据中的异常伪距观测数据。
7.根据权利要求5所述的精密单点定位方法，其特征在于，所述对所述观测数据进行预处理的步骤，包括：
根据每个所述卫星的伪距观测数据和多普勒观测数据获取每个所述卫星的多普勒探测伪距检验值；
根据每个所述卫星的多普勒探测伪距检验值剔除所有所述卫星的伪距观测数据中的异常伪距观测数据。
8.根据权利要求6所述的精密单点定位方法，其特征在于，所述根据每个所述卫星的伪距观测数据获取每个所述卫星的伪距频间偏差的步骤，包括：
根据每个所述卫星的伪距观测数据获取每个所述卫星的第一频率对应的伪距观测数据和第二频率对应的伪距观测数据；
获取每个所述卫星的第一频率对应的伪距观测数据和第二频率对应的伪距观测数据的差值，得到每个所述卫星的伪距频间偏差。
9.根据权利要求8所述的精密单点定位方法，其特征在于，所述根据所有所述卫星的伪距频间偏差剔除所有所述卫星的伪距观测数据中的异常伪距观测数据的步骤，包括：
获取所有所述卫星的伪距频间偏差的中位数；
计算每个所述卫星的伪距频间偏差和所有所述卫星的伪距频间偏差的中位数的差值，得到每个所述卫星的第一检验值；
剔除所述第一检验值大于第一预设阈值的所述卫星对应的伪距观测数据。
10.根据权利要求7所述的精密单点定位方法，其特征在于，通过以下公式获取所述多普勒探测伪距检验值：
ΔPi,j,k＝Pi,j,k-Pi,j,k-1-0.5*λi,j(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；
Di,j,k为k历元卫星i的j(j＝1,2)频率的多普勒观测数据；
Pi,j,k-1为k-1历元卫星i的j(j＝1,2)频率的伪距观测数据；
Di,j,k-1为k-1历元卫星i的j(j＝1,2)频率的多普勒观测数据；
λi,j为卫星i的第j(j＝1,2)频率的波长值；
tk为k历元的时间；
tk-1为历元的时间；
ΔPi,j,k为k历元卫星i的j(j＝1,2)频率的多普勒探测伪距检验值。
11.根据权利要求5至10任一项所述的精密单点定位方法，其特征在于，所述对所述观测数据进行预处理的步骤，还包括：
根据以下公式对每个所述卫星的伪距观测数据进行平滑处理，得到多普勒平滑伪距；
PSi,j,k＝ωkPi,j,k+(1-ωk)(PSi,j,k-1+0.5*λi,j(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，PSi,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒平滑伪距；
ωk为k历元对应的平滑因子；
Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；
PSi,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒平滑伪距；
λi,j为卫星i的第j(j＝1,2)频率的波长值；
Di,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒观测数据；
Di,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒观测数据；
tk为k历元的时间；
tk-1为k-1历元的时间。
12.根据权利要求5至11任一项所述的精密单点定位方法，其特征在于，所述对所述卫星的观测数据进行预处理的步骤，包括：
判断每个载波相位观测数据的整周模糊度是否发生跳变；
若是，对整周模糊度跳变的所述载波相位观测数据进行标记。
13.根据权利要求12所述的精密单点定位方法，其特征在于，所述判断每个载波相位观测数据的整周模糊度是否发生跳变的步骤包括：
采用历元间差分算法对所述载波相位观测数据进行单频周跳探测，得到跳变结果。
14.根据权利要求13所述的精密单点定位方法，其特征在于，在所述采用历元间差分算法对所述载波相位观测数据进行单频周跳探测，得到跳变结果的步骤之前，还包括：利用以下公式获取所述载波相位观测数据的周跳探测值：
ΔLi,j,k＝Li,j,k-Li,j,k-1-0.5*(-Di,j,k-Di,j,k-1)/(tk-tk-1)；
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的周跳探测值；
Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Di,j,k为k历元卫星i的第j(j＝1,2)频率的多普勒观测数据；
Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Di,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的多普勒观测数据；
tk为k历元的时间；
tk-1为k-1历元的时间；
若|ΔLi,j,k|大于第二预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；
若|ΔLi,j,k|小于或等于所述第二预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
15.根据权利要求13所述的精密单点定位方法，其特征在于，所述采用历元间差分算法对所述载波相位观测数据进行单频周跳探测，得到跳变结果的步骤，包括：
利用以下公式对所述载波相位观测数据进行单频周跳探测：
ΔLi,j,k＝Li,j,k-Li,j,k-1-(Pi,j,k-Pi,j,k-1)/λi,j；
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的单频周跳探测值；
Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Pi,j,k为k历元卫星i的第j(j＝1,2)频率的伪距观测数据；
Pi,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的伪距观测数据；
λi,j为卫星i的第j(j＝1,2)频率的波长值；
若|ΔLi,j,k|大于第三预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；
若|ΔLi,j,k|小于或等于所述第三预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
16.根据权利要求13所述的精密单点定位方法，其特征在于，所述采用历元间差分算法对所述载波相位观测数据进行单频周跳探测，得到跳变结果的步骤，包括：
利用以下公式对所述载波相位观测数据进行单频周跳探测：
ΔLi,j,k＝Li,j,k-3*Li,j,k-1+3*Li,j,k-2-Li,j,k-3；
其中，ΔLi,j,k为k历元卫星i的第j(j＝1,2)频率的单频周跳探测值；
Li,j,k为k历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Li,j,k-1为k-1历元卫星i的第j(j＝1,2)频率的载波相位观测数据；
Li,j,k-2为k-2历元所述卫星i的第j(j＝1,2)频率的载波相位观测数据；
Li,j,k-3为k-3历元所述卫星i的第j(j＝1,2)频率的载波相位观测数据；
λi,j为卫星i的第j(j＝1,2)频率的波长值；
若|ΔLi,j,k|大于第四预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据发生跳变；
若|ΔLi,j,k|小于或等于所述第四预设阈值，则判定k历元卫星i的第j(j＝1,2)频率的载波相位观测数据未发生跳变。
17.根据权利要求1所述的精密单点定位方法，其特征在于，所述观测数据通过导航卫星发送到接收机，所述差分数据通过通信卫星发送到所述接收机。
18.一种精密单点定位装置，其特征在于，包括：
获取模块，用于获取卫星的观测数据和卫星的差分数据，所述差分数据包括电离层延迟数据和湿对流层延迟数据；
方程组建模块，用于组建非组合PPP观测方程，所述非组合PPP观测方程包括电离层参数和湿对流层参数；
延迟值获取模块，用于根据所述电离层延迟数据和所述湿对流层延迟数据计算所述电离层参数的电离层延迟值和所述湿对流层参数的湿对流层延迟值；
约束模块，用于利用所述电离层延迟值和所述湿对流层延迟值对所述非组合PPP观测方程的所述电离层参数和所述湿对流层参数进行约束；
求解模块，用于求解所述非组合PPP观测方程，得到定位坐标参数信息。
19.一种电子设备，其特征在于，包括：存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如权利要求1-17任一项所述的精密单点定位方法的步骤。
20.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有指令，当所述指令在计算机上运行时，使得所述计算机执行如权利要求1-17任一项精密单点定位所述的方法。
CN202210050395.0A 2022-01-17 2022-01-17 一种精密单点定位方法、装置、电子设备和存储介质 Active CN114296119B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210050395.0A CN114296119B (zh) 2022-01-17 2022-01-17 一种精密单点定位方法、装置、电子设备和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210050395.0A CN114296119B (zh) 2022-01-17 2022-01-17 一种精密单点定位方法、装置、电子设备和存储介质

Publications (2)

Publication Number Publication Date
CN114296119A true CN114296119A (zh) 2022-04-08
CN114296119B CN114296119B (zh) 2023-10-20

Family

ID=80976701

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210050395.0A Active CN114296119B (zh) 2022-01-17 2022-01-17 一种精密单点定位方法、装置、电子设备和存储介质

Country Status (1)

Country Link
CN (1) CN114296119B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115826018A (zh) * 2023-02-22 2023-03-21 广州导远电子科技有限公司 一种模糊度固定方法、装置、接收机及存储介质
CN115826009A (zh) * 2023-02-22 2023-03-21 广州导远电子科技有限公司 一种周跳探测方法、存储介质和接收机

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160077213A1 (en) * 2014-09-15 2016-03-17 Fugro N.V. Integer ambiguity-fixed precise point positioning method and system
CN108415049A (zh) * 2018-04-19 2018-08-17 千寻位置网络有限公司 提高网络rtk双差宽巷模糊度固定正确率的方法
CN109581452A (zh) * 2018-12-18 2019-04-05 辽宁工程技术大学 一种gnss参考站载波相位整周模糊度解算方法
CN111025346A (zh) * 2019-11-18 2020-04-17 广州南方卫星导航仪器有限公司 一种快速估计gnss精密卫星钟差的方法及存储介质
CN112835082A (zh) * 2021-01-05 2021-05-25 广州星际互联科技有限公司 一种gnss区域增强电离层和对流层大气产品质量指标计算方法、电子设备、存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160077213A1 (en) * 2014-09-15 2016-03-17 Fugro N.V. Integer ambiguity-fixed precise point positioning method and system
CN108415049A (zh) * 2018-04-19 2018-08-17 千寻位置网络有限公司 提高网络rtk双差宽巷模糊度固定正确率的方法
CN109581452A (zh) * 2018-12-18 2019-04-05 辽宁工程技术大学 一种gnss参考站载波相位整周模糊度解算方法
CN111025346A (zh) * 2019-11-18 2020-04-17 广州南方卫星导航仪器有限公司 一种快速估计gnss精密卫星钟差的方法及存储介质
CN112835082A (zh) * 2021-01-05 2021-05-25 广州星际互联科技有限公司 一种gnss区域增强电离层和对流层大气产品质量指标计算方法、电子设备、存储介质

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
伍冠滨 等: ""基于非差非组合PPP-RTK的大气改正模型及其性能验证"" *
周锋: ""多系统GNSS非差非组合精密单点定位相关理论和方法研究"" *
潘宗鹏: ""实时精密单点定位及模糊度固定"" *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115826018A (zh) * 2023-02-22 2023-03-21 广州导远电子科技有限公司 一种模糊度固定方法、装置、接收机及存储介质
CN115826009A (zh) * 2023-02-22 2023-03-21 广州导远电子科技有限公司 一种周跳探测方法、存储介质和接收机
CN115826009B (zh) * 2023-02-22 2023-04-28 广州导远电子科技有限公司 一种周跳探测方法、存储介质和接收机

Also Published As

Publication number Publication date
CN114296119B (zh) 2023-10-20

Similar Documents

Publication Publication Date Title
AU2004281011B2 (en) Method for using three GPS frequencies to resolve carrier-phase integer ambiguities
US7961143B2 (en) Partial search carrier-phase integer ambiguity resolution
US7692578B2 (en) GNSS signal processing with clock difference modeling
US6127968A (en) On-the-fly RTK positioning system with single frequency receiver
JP5421903B2 (ja) 部分探索搬送波位相整数アンビギュイティ決定
EP2067053B1 (en) Method for using three gps frequencies to resolve whole-cycle carrier-phase ambiguities
Li Cycle slip detection and ambiguity resolution algorithms for dual-frequency GPS data processing
CN114296119B (zh) 一种精密单点定位方法、装置、电子设备和存储介质
CN112731496B (zh) 一种面向智能终端的gnss精密单点定位数据质量控制方法
CN112540393A (zh) 海上精密单点定位保护水平质量检核方法及系统
AMAMI Enhancing Stand-Alone GPS Code Positioning Using Stand-Alone Double Differencing Carrier Phase Relative Positioning
AU2012200450B2 (en) Method for using three GPS frequencies to resolve whole-cycle carrier-phase ambiguities
CN113093237A (zh) Ssr轨钟改正数质量因子实时评估方法、装置、设备及介质
CN115877414B (zh) 一种固定解坐标的检核方法、装置、存储介质及电子设备
CN115826018B (zh) 一种模糊度固定方法、装置、接收机及存储介质
CN117331106A (zh) 基于速度辅助的高精度手机定位方法及系统
CN117908060A (zh) 一种定位方法及装置
Nikiforov Asymptotically efficient estimation of a nonlinear model of the heteroscedasticity and the calibration of measurement systems
CN112987048A (zh) 一种适用于Andriod智能终端的高精度定位方法、智能终端以及存储介质
CN115993620A (zh) 模糊度固定方法及其系统
CN115902958A (zh) 接收机间相位频间偏差标定方法、装置及设备
Ya'acob et al. Ionospheric Correction and Ambiguity Resolution in DGPS with Single Frequency

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