一种卫星定位信号解算方法及装置与流程

本发明属于卫星定位技术领域，具体涉及一种卫星定位信号解算方法及装置。

背景技术：

近年来基于全球导航卫星系统gnss(theglobalnavigationsatellitesystem)的定位方法由于其高精度的定位结果广泛应用于各领域，要获得高精度的定位结果需要使用gnss载波相位测量值，而在实际应用场景中，由于卫星信号的中断、信噪比过低、卫星高度角偏低、电离层变化比较剧烈和多径误差较大等原因，gnss载波相位测量值的整周计数很容易发生跳变或中断(即周跳)。载波相位测量值只有在正确探测修复周跳的条件下才具有可用性，才能得到高精度的定位结果。因此探测周跳是gnss定位中非常重要的一个环节。

目前，gnss周跳探测已有多种方法，其中，多项式拟合法、高次差法、三差法等难以用于动态情况；mw组合法(melbourne和wubbena于1985年提的组合观测值算法)、电离层残差法及三频组合法无法探测单频周跳，而伪距相位组合法探测不到小周跳，而利用惯导位置和速度进行周跳探测的方法会受惯性导航系统(ins，简称惯导)的累积误差的影响。

在使用惯导位置和速度进行周跳探测时，需要使用加速度信息。惯性测量单元imu的加速度计输出移动站的加速度信息，现有技术会对其积分得到速度信息，再根据相邻时刻的速度平均值与时间间隔的乘积来计算出位置增量。这种方法由于使用积分存在较大的累积误差。

技术实现要素：

本发明主要针对现有技术中周跳探测不准的问题，发明了一种卫星定位信号解算方法及装置，避免了imu累积误差对周跳的影响，提高了周跳的准确度，从而提高了定位的精度。

本发明的上述技术问题是通过以下技术方案得以实施的：一种卫星定位信号解算装置，其特征在于，包括：

惯性测量模块，用于获取加速度计测量值；

计算模块，用于使用所述加速度计测量值计算姿态信息，还用于使用姿态信息和加速度计测量值计算移动站的加速度，使用所述移动站的加速度计算周跳检验量；

比较模块，用于比较阈值s与周跳检验量；

定位解算模块，当周跳检验量小于s时，定位解算模块进行定位解算；当周跳检验量大于s时，重置定位解算模块的整周模糊度，进行定位解算。

作为优选，还包括位置检测模块，用于获取移动站的位置信息，所述位置信息包括当地水平坐标系下的重力向量。

一种卫星定位信号解算方法，采用上述的卫星定位信号解算装置，其特征在于，解算方法如下：

使用所述惯性测量模块获取加速度计测量值；

使用所述计算模块根据所述加速度计测量值计算姿态信息；

使用所述计算模块根据姿态信息及惯性测量单元的加速度计测量值计算装置的加速度，根据装置加速度计算周跳检验量；

使用所述比较模块比较阈值s与周跳检验量，当周跳检验量小于s时，定位解算模块进行定位解算；

使用所述比较模块比较阈值s与周跳检验量，当周跳检验量大于s时，重置整周模糊度，进行定位解算。

作为优选，使用所述位置检测模块获取移动站的位置信息，所述位置信息包括当地水平坐标系下的重力向量；

作为优选，所述计算模块根据基于位置信息以及惯性测量单元的加速度计测量值作为观测向量建立观测方程，确定姿态信息。

作为优选，所述阈值s在0到1之间。

作为优选，所述计算模块使用计算周跳检验量，其中，p、q为不同卫星，u为移动站，r为基站，e为方向余弦，和分别为卫星p、q和移动站的加速度，λ为载波相位波长，为载波相位双差的加速变化率。

作为优选，所述观测方程为以imu的加速度计测量值作为测量向量，姿态信息作为状态向量滤波得到姿态信息。

作为优选，采用rtk定位技术进行卫星定位。

综上所述，本发明与现有技术相比具有如下优点：

本发明使用imu加速度计测量值建立观测向量，通过滤波得到了较为准确的姿态信息，避免使用imu直接输出姿态产生累积误差。

本发明设置了周跳检验量和周跳检验阈值，当周跳检验量超过所设阈值时，通过重置模糊度操作改善了定位结果。

本发明通过避免imu累积误差的影响得到跟准确的姿态信息。

本发明通过避免imu累积误差的影响准确探测出周跳。

附图说明

图1是根据本申请的一些实施例的卫星定位信号解算流程图。

图2是本申请的一些实施例的静态条件载波相位双差值变化图。

图3是本申请的一些实施例的静态条件周跳检验量变化图。

图4是本申请的一些实施例的动态条件载波相位双差值变化图。

图5是本申请的一些实施例的动态条件周跳检验量变化图。

具体实施方式

在下面的详细描述中，通过示例阐述了本申请的许多具体细节，以便提供对相关披露的透彻理解。然而，对于本领域的普通技术人员来讲，本申请显而易见的可以在没有这些细节的情况下实施。应当理解的是，本申请中使用“系统”、“装置”、“单元”和/或“模块”术语，是用于区分在顺序排列中不同级别的不同部件、元件、部分或组件的一种方法。然而，如果其他表达式可以实现相同的目的，这些术语可以被其他表达式替换。

应当理解的是，当设备、单元或模块被称为“在……上”、“连接到”或“耦合到”另一设备、单元或模块时，其可以直接在另一设备、单元或模块上，连接或耦合到或与其他设备、单元或模块通信，或者可以存在中间设备、单元或模块，除非上下文明确提示例外情形。例如，本申请所使用的术语“和/或”包括一个或多个相关所列条目的任何一个和所有组合。

参看下面的说明以及附图，本申请的这些或其他特征和特点、操作方法、结构的相关元素的功能、部分的结合以及制造的经济性可以被更好地理解，其中说明和附图形成了说明书的一部分。然而，可以清楚地理解，附图仅用作说明和描述的目的，并不意在限定本申请的保护范围。可以理解的是，附图并非按比例绘制。

本申请中使用了多种结构图用来说明根据本申请的实施例的各种变形。应当理解的是，前面或下面的结构并不是用来限定本申请。本申请的保护范围以权利要求为准。

下面结合附图和实施例对本发明进一步说明。

实施例1：

图1所示，一种卫星定位信号解算方法及装置，其中，卫星定位信号解算装置包括惯性测量模块、计算模块、比较模块、定位解算模块、位置检测模块。

所述惯性测量模块用于获取加速度计测量值；所述计算模块用于使用所述加速度计测量值计算姿态信息；所述计算模块还用于使用姿态信息和加速度计测量值计算移动站的加速度，根据所述移动站加速度计算周跳检验量。

所述比较模块用于比较阈值s与周跳检验量，当周跳检验量小于s时，定位解算模块进行定位解算；当周跳检验量大于s时，重置定位解算模块的整周模糊度，进行定位解算。

所述位置检测模块根据基于位置的状态空间模型，以惯性测量单元的加速度计测量值作为观测向量建立观测方程，确定姿态信息。

所述位置检测模块还用于获取移动站所处位置的重力向量。

本实施例，使用rtk(real-timekinematic，实时动态)载波相位差分技术进行卫星定位，具体步骤为：

1.通过当地水平坐标系里导航方程的状态空间模型得到方程

式中，g^l为当地水平坐标系下的重力向量。为b系到l系的坐标变换矩阵，和为角速度和相应的反对称矩阵，f^b为加速度计的比例测量值，和v^l为rtk输出的速度和加速度。具体地：

其中，仅有姿态信息为未知量，因此将f^b作为观测向量，姿态角(俯仰角p、横滚角r和航向角y)为状态向量的观测方程，有：其中，

上式中，p,r,y分别为俯仰角、横滚角和航向角。

2.基于观测方程建立滤波模型其中b代表imu，u代表移动站，e为方向余弦，为imu加速度计测量值，为移动站的加速度。以imu的加速度计测量值作为测量向量，姿态角为状态向量滤波；最终得到姿态信息

3.基于步骤2得到的较为准确姿态信息结合imu的加速度计测量值获得准确的移动站加速度，进一步基于双差加速变化率的计算值和测量值构造周跳检验量：

其中，p、q为不同卫星，u代表移动站，r代表基站，e为方向余弦，和分别为卫星p、q和移动站的加速度，λ为载波相位波长，为载波相位双差的加速变化率。

4.进一步设置周跳检验阈值为0.4，当周跳检验量没有超过所设阈值时，直接进行定位解算，当周跳检验量超过所设阈值时，则重置模糊度，进一步进行rtk定位解算。

本发明的优势在于，使用imu加速度计测量值建立观测向量，通过滤波得到了较为准确的姿态信息，避免使用imu直接输出姿态产生累积误差。

本发明的优势还在于，设置了周跳检验量和周跳检验阈值，当周跳检验量超过所设阈值时，通过重置模糊度操作改善了定位结果。

图2，图3，图4，图5为本实施例的实验数据图。

在静态条件下，载波相位测量值加入3mm测量误差时，在第50，150历元处给第二颗卫星分别加入半周、一周跳变，在第80历元处给第三颗卫星加入2周跳变时，载波相位双差值如图2所示，周跳检验量的变化情况如图3所示。

由图2，图3可知，在静态条件下，本发明可以准确的探测出半周及其以上周跳发生的位置。

在动态条件下，载波相位测量值加入3mm测量误差时，在第50，150历元处给第二颗卫星分别加入半周、一周跳变，在第80历元处给第三颗卫星加入2周跳变时，载波相位双差值如图4所示，周跳检验量的变化情况如图5所示。

由图4，图5可知，在动态条件下，本发明可以准确的探测出半周及其以上周跳发生的位置。

本发明相比于现有技术，具有如下有益效果：

一、通过避免imu累积误差的影响得到跟准确的姿态信息。

二、通过避免imu累积误差的影响准确探测出周跳。

应当理解的是，本申请的上述具体实施方式仅仅用于示例性说明或解释本申请的原理，而不构成对本申请的限制。因此，在不偏离本申请的精神和范围的情况下所做的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。此外，本申请所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。