一种顾及Doppler的北斗修正伪距单点定位方法、系统与设备

技术领域

本发明属于北斗卫星导航系统数据预处理技术领域，具体涉及一种顾及Doppler的北斗修正伪距单点定位方法、系统与设备。

背景技术

目前，由于全球导航卫星系统(GNSS，Global Navigation Satellite System)伪距单点定位模块具有全天时、全天候、体积小、便于携带等优点，伪距单点定位技术已广泛应用于各个领域。

北斗卫星导航系统(BDS,BeiDou Navigation Satellite System)于2020年7月31日正式开通，标志着北斗迈入全球服务新时代，也意味着北斗伪距单点定位技术将具有更加广阔的应用前景

然而，由于伪距观测值受多种误差源的影响，该定位精度难以保证，个别定位误差大于百米，成为制约北斗伪距定位技术推广应用的瓶颈。

发明内容

为了解决现有技术中存在的问题，本发明提供一种顾及Doppler的北斗修正伪距单点定位方法、系统与设备，基于改善北斗伪距观测值的精度，计算北斗伪距单点定位，提高北斗伪距单点定位的精度与可靠性。

为了实现上述目的，本发明采用的技术方案：一种顾及Doppler的北斗修正伪距单点定位方法，具体包括以下步骤：

获取多普勒观测值和原始伪距观测值；

构建顾及Doppler的平滑伪距观测值表达式：

其中，ω为平滑因子；λ为载波波长；P为原始伪距观测值，

向所述平滑伪距观测值表达式中引入权重因子η调整原始伪距观测值与平滑伪距观测值的权重，则修正伪距

其中，η为权重因子；

基于所述修正伪距观测值进行单点定位计算，得到地面接收机的空间三维坐标值，即地面接收机的单点定位值。

将平滑伪距观测值表达式中的平滑因子设定为平滑窗口的倒数ω＝1/k，则平滑伪距观测值的表达式为：

其中T＝t

将卫星高度角的影响引入平滑伪距观测值表达式中，平滑伪距观测值表达式为：

其中，θ为卫星高度角。

平滑伪距的方差为：

其中，

依据误差传播定律，计算修正伪距

F

其中，τ、μ均为替换字符，无实义，将所述权重因子计算值引入所述修正伪距

本发明还提供一种顾及Doppler的北斗修正伪距单点定位系统，包括数据获取模块，数据处理模块以及定位模块，其中数据获取模块用于获取多普勒观测值和原始伪距观测值，数据处理模块基于所述多普勒观测值和原始伪距观测值对原始伪距观测值进行修正，具体如下：

获取多普勒观测值和原始伪距观测值；

构建顾及Doppler的平滑伪距观测值表达式：

其中，ω为平滑因子；λ为载波波长；P为原始伪距观测值，

向所述平滑伪距观测值表达式中引入权重因子η调整原始伪距观测值与平滑伪距观测值的权重，则修正伪距

其中，η为权重因子；

定位模块基于所述修正伪距观测值进行单点定位计算，得到地面接收机的空间三维坐标值，即地面接收机的单点定位值。

为解决上述技术问题，本发明还提出一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现本发明所述顾及Doppler的北斗修正伪距单点定位方法的步骤。

本发明还提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现本发明所述顾及Doppler的北斗修正伪距单点定位方法的步骤。

与现有技术相比，本发明至少具有以下有益效果：充分发挥独立观测值Doppler的作用，削弱北斗原始伪距观测值的误差，修正伪距观测值的精度，权重因子的取值取决于卫星高度角、平滑窗口k和采样间隔T；当卫星高度角与平滑窗口k确定后，权重因子的大小主要由采样间隔T决定，即采样间隔T越大，τ值越大，权重因子η越小，平滑伪距

附图说明

图1a是原始伪距噪声水平随采样间隔变化图，图中，不同颜色(蓝、绿、红)分别对应代表北斗B1、B2和B3频点观测值。

图1b是平滑伪距噪声水平随采样间隔变化图，图中，不同颜色(蓝、绿、红)分别对应代表北斗B1、B2和B3频点观测值。

图1c是修正伪距噪声水平随采样间隔变化图，图中，不同颜色(蓝、绿、红)分别对应代表北斗B1、B2和B3频点观测值。

具体实施方式

在以下的叙述中，为了使读者更好地理解本申请而提出了许多技术细节。但是，本领域的普通技术人员可以理解，即使没有这些技术细节和基于以下各实施方式的种种变化和修改，也可以实现本申请所要求保护的技术方案。

为使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请的实施方式作进一步地详细描述，获取多普勒观测值和原始伪距观测值

步骤1：构建顾及Doppler的平滑伪距观测值表达式：

需指出，式(1)采用递归方式获取平滑伪距

进一步的，平滑因子可设定为平滑窗口的倒数(ω＝1/k)，则平滑伪距观测值表达式为：

需指出，随着平滑窗口的增大，平滑因子变小，原始伪距观测值的比重降低，Doppler观测值的比重变大，其在平滑伪距中发挥的作用增大。

卫星高度角用于表征观测数据质量，将卫星高度角的影响引入平滑伪距观测值表达式中，平滑伪距可表示为

其中，θ为卫星高度角，由于Doppler观测值与时间的相关性可忽略，采用误差传播定律计算平滑伪距的方差值：

需指出，在采样间隔固定，原始伪距与Doppler观测值噪声水平一定的情况下，平滑伪距的方差值主要取决于平滑窗口的大小。由于k既为分子项，又为分母项，故随着平滑窗口的增大，平滑伪距方差值降低，但二者并非线性负相关。

更进一步的，为确定平滑窗口大小，F

需指出，式(5)是关于k的三次方程，即k有3个根。考虑到该函数理论上与X轴仅有一个交点，故该方程的根包括1个实数解与2个共轭复数解。由于平滑窗口必为正整数，需对实数解进行四舍五入取整，则有

需指出，为保证平滑作用起效果，

步骤2：引入权重因子构建修正伪距观测值表达式：

式(7)旨在通过引入权重因子η调整原始伪距观测值P

进一步的，依据误差传播定律，修正伪距的方差值可表示为：

更进一步的，式(8)对k求导并赋值为

式(9)中，权重因子的取值主要取决于卫星高度角、平滑窗口k和采样间隔T。当卫星高度角与平滑窗口k固定后，权重因子的大小主要由采样间隔T决定，即采样间隔T越大，权重因子η越小，平滑伪距

参考图1a、图1b和图1c，其中图示从左至右三个柱体分别表示北斗B1、B2和B3的观测值，原始伪距噪声水平随采样间隔变化图，图1b，平滑伪距噪声水平随采样间隔变化图，图1c是修正伪距噪声水平随采样间隔变化图，由图可知，随着采样间隔的增加，伪距噪声水平均增大，且B1频点伪距噪声水平大于B2与B3频点。当采样间隔不大于2秒时，由于Doppler观测值积分误差较小，平滑伪距噪声水平优于原始伪距观测值。随着采样间隔持续增加，Doppler观测值积分误差显著增大，平滑伪距噪声水平大于原始伪距观测值。由于权重因子的调整，修正伪距噪声水平优于原始伪距与平滑伪距观测值，且在采样间隔较小时，修正效果更佳。

根据已有文献资料可知，Doppler积分累计误差一般不超过2cm/s。当数据采样间隔不超过1s时，Doppler积分累计误差水平优于原始伪距观测值噪声水平，故平滑伪距观测值精度优于原始伪距观测值，其权重将增加。随着采样间隔的增加，Doppler累计积分误差增大，平滑伪距观测值权重将降低。由于采样间隔一般不超过30秒，故权重因子η不会无限接近于0，即平滑伪距

本发明提供一种顾及Doppler的北斗修正伪距单点定位系统，包括数据获取模块，数据处理模块以及定位模块，其中数据获取模块用于获取多普勒观测值和原始伪距观测值，数据处理模块基于所述多普勒观测值和原始伪距观测值对原始伪距观测值进行修正，具体如下：

获取多普勒观测值和原始伪距观测值；

构建顾及Doppler的平滑伪距观测值表达式：

其中，ω为平滑因子；λ为载波波长；P为原始伪距观测值，

向所述平滑伪距观测值表达式中引入权重因子η调整原始伪距观测值与平滑伪距观测值的权重，则修正伪距

其中，η为权重因子；

定位模块基于所述修正伪距观测值进行单点定位计算，得到地面接收机的空间三维坐标值，即地面接收机的单点定位值。

本发明顾及Doppler的北斗伪距观测值修正方法可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。本发明顾及Doppler的北斗伪距观测值修正方法如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。

基于这样的理解，在示例性实施例中，还提供了一种计算机可读存储介质，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于该计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。计算机可读存储介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括电载波信号和电信信号。其中，所述计算机存储介质可以是计算机能够存取的任何可用介质或数据存储设备，包括但不限于磁性存储器(例如软盘、硬盘、磁带、磁光盘(MO)等)、光学存储器(例如CD、DVD、BD、HVD等)、以及半导体存储器(例如ROM、EPROM、EEPROM、非易失性存储器(NANDFLASH)、固态硬盘(SSD))等。

在示例性实施例中，还提供计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述顾及Doppler的北斗伪距观测值修正方法的步骤。处理器可能是中央处理单元(CentralProcessing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital SignalProcessor、DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。

以上所述仅为本发明的优选实施例，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的权利要求范围之内。