具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

如图1所示，本发明提供了一种基于GNSS信号非连续跟踪的导航方法，在同一采集时间过顶的卫星超过四颗时，通过组合算法确定每一次采集对应的GNSS卫星组合，且不同采集时间点对应的GNSS卫星组合不同，从而实现多卫星导航并避免其中一个卫星出错，而导致导航出错的问题，另外，根据采集时间点与卫星过顶时间点的对比，可筛选出最靠近目标时区的卫星序列，进一步的提高导航的精确性。

具体包括以下步骤：

步骤100、批量生成GNSS卫星轨迹数据，并基于卫星轨迹数据确定多GNSS卫星在目标时区的过顶预报时间，其中，过顶预报时间为根据GNSS卫星的时间-位置数据确定的经过目标时区头顶的时间。

在步骤100中，基于卫星轨迹数据确定多GNSS卫星在目标时区的过顶预报时间的实现步骤为：

1、通过两行轨道根数TLE确定所有GNSS卫星运动轨迹的轨道参数，并利用SGP4模型将TLE轨道参与轨道运算模型结合测算卫星实时状态，得到所有GNSS卫星的时间-位置数据。

TLE数据只是提供菜颗卫星的轨道参数，卫星实时状态的测算需要将E轨道参数具体的轨道运算模型或算法结合，才能得到较为精确的结果，SGP4,SDP4(Simplified Deep-space Perturbation Version 4)为简化扰动模型，简化扰动系列模型考虑了地球形状、大气阻力、宇宙辐射，以及来自其他太空物体的干扰，计算地心惯性坐标系下.卫星、空间碎片等人造太空物体的时间-位置和时间-速率等物理量。简化扰动模型通常代指SGP4模型，这是因为SGP4模型是使用最广泛的简化扰动模型，许多学者对它的精度和误差分析也进行了大量的研究，通常认为TLE结合SGP4模型计算卫星位置和速度是比较可取的。

2、将时间-位置数据写成CZML数据文件，利用循环算法获取连续时间点的位置坐标。

Python是一种开源的、面向对象的脚本语言。这种编程语言风格简洁、简单易学、强制格式规范，还有强大的学习社区支持，灵活的可扩展性和功能丰富的库。因此，本文采用Python语言将TlE数据与SGP4模型结合，对导航卫星的时间-位置，时间速率进行计算，并将计算结果存入CZML格式的数据文件中。

因此卫星轨道数据(以及其他飞行物的轨迹数据)的生成分为两个步骤:首先由TLE数据与SGP4模型得到卫星的时间-位置数据；再将时间位置数据写成CZML数据文件。

卫星的轨道是由一系列连续的时间位置样本点数据组成的，将这一系列样本点连接起来形成卫星的椭圆(圆)轨道，因此轨道数据的获取就分成两个部分位置数据和时间数据，卫星的位置坐标可通过上文提到的SGP4模型计算得出，接下来，需要获取连续时间点的位置坐标，并把这些时间、位置数据存储成CZML格式兼容的数据。CZML中“position”"属性提供的“carisian”子属性是以米为单位的三维笛卡尔坐标，恰好与SPG4模型生成的坐标属同以坐标系。因此可直接将上文生成的数据存储到“arisian”属性中。一段时间内的连续时间-位置坐标数据的生成需要通过循环实现，例如生成某颗卫星10分钟内的时间-位置数据：首先需要指定时间的起止点，然后指定单位间隔(例如隔3秒钟)来生成每坐标数据的生成需要通过循环实现，例如生成某颗卫星10分钟内的时间-位置数据，接着通过一个循环实现一段时间内时间-位置(如10分钟内共200次数据采样)坐标的计算，并在每次计算时记录当前时刻:最后将计算所得的时刻-坐标数据存储到"carisian”属性中。

即在上述步骤中，CZML数据文件为GNSS卫星轨道CZML数据，通过确定生成GNSS卫星的时间-位置样本的时间间隔t，同时利用CZML提供的内插属性，将时间-位置样本数据之间通过高阶拉格朗日算法内插填充样本数据间的时间间隔t内，其中，每个GNSS卫星的CZML数据只记录该GNSS卫星在额定秒数内的轨迹状况，每个CZML数据文件都需要Python调用SPG4库，根据额定秒数和时间间隔t计算得到该GNSS卫星一系列的时间-位置数据。

3、根据连续时间点的位置坐标获取每个GNSS卫星经过目标时区的过顶预报时间。

卫星的过面时间是卫星经过观测站或观测者头顶的时间，过顶预报是供据架测者所处的位置(经度、纬度、高程)、气压、温度和卫星星历信息预报卫星未来的过顶时间，卫量的过顶预报一般是先计算某时刻卫星的位置-速度，然后算出该时刻卫星相对于目标时区的仰角是否满足所设定的观测约束条件，如果满足就认为该时刻属于卫星的过顶预报时间。

再由从时间-位置数据内筛选GNSS卫星经过目标时区的过顶预报时间，并将过顶预报时间和对应的GNSS卫星组成卫星序列。

本实施方式利用作为卫星轨道数据源的两行轨道数据，常用的SGP4轨道计算模型，并使用Python脚本语言将这两种资源结合生成CZML格式的卫星轨道数据，进一步的，还可以通过目标物CZML轨道数据可视化操作，来建立多卫星环绕动态图，具体的实现过程可参考现有技术中的轨迹数据的流式传输与可视化操作文献。

步骤200、按照过顶预报时间的时间顺序预测依次在目标时区内过顶的卫星序列，且根据GNSS卫星的导航定位需求确定单次跟踪的卫星接收数量。

由于GNSS卫星比较多，且环绕地球转动，因此每个卫星的过顶时间并不是完全相同的，因此卫星序列为多行多列式，卫星序列的同一行为同一过顶时间点在目标时区内过顶的GNSS卫星，卫星序列的不同行从上到下表示按照时间顺序依次经过目标时区的GNSS卫星排列顺序。

GNSS卫星的导航定位需求确定单次跟踪的卫星接收数量为四个，同一时间点在目标时区内过顶的GNSS卫星大于等于卫星接收数量。

也就是说，当卫星序列的一行表示同一过顶时间点在目标时区内过顶的GNSS卫星，在卫星序列一行包含的GNSS卫星数量大于卫星接收数量时，那在非连续的采集跟踪时，则需要对GNSS卫星进行选择，来保证跟踪采集的稳定工作，特别在多次非连续的采集跟踪在卫星序列中对应的GNSS卫星相同时，则更需要确定每次采集跟踪的卫星，本实施方式通过组合算法确定卫星序列某一行的GNSS卫星组合，且一般来说，每次采集跟踪的GNSS卫星组合均不相同，可以避免某一卫星信息不准确时，导致多次采集跟踪出现错误，且提高导航精确性，具体的实现过程如下述步骤300所得。

步骤300、确定非连续采集的采集时间列表，每次采集时间点基于卫星接收数量通过组合筛选方式从卫星序列中确定用于跟踪的GNSS卫星组合，组合筛选方式用于在至少两次采集时间点对应的卫星序列相同且卫星序列的数量大于卫星接收数量时，通过组合算法确定每次采集时间点对应的GNSS卫星组合。

在步骤300中，先根据采集时间点与卫星序列中的两行数据的过顶时间点计算采集时间点的跟踪偏向，再根据跟踪偏向选择一行GNSS卫星作为目标筛选卫星。

每次采集时间点基于卫星接收数量通过组合筛选方式从目标筛选卫星组中确定用于跟踪的GNSS卫星组合，其中，在多次采集时间点对应的目标筛选卫星相同时，利用组合筛选方式使得每次的采集时间点对应的GNSS卫星组合不同。

在同一采集时间过顶的卫星超过四颗时，通过组合算法确定每一次采集对应的GNSS卫星组合，且不同采集时间点对应的GNSS卫星组合不同，从而实现多卫星导航并避免其中一个卫星出错，而导致导航出错的问题。

进一步的，计算采集时间点的跟踪偏向的实现步骤为：

计算卫星序列中两行GNSS卫星的过顶时间点确定两个过顶时间点的中心时间点，中心时间点将两行GNSS卫星分为上一行GNSS卫星的过顶时间点与中心时间点之间的上阈值，以及下一行GNSS卫星的过顶时间点与中心时间点之间的下阈值。

判断采集时间点处于上阈值或下阈值，来选择其中一行GNSS卫星作为目标筛选卫星。

确定目标筛选卫星的过程即筛选靠近目标时区的GNSS卫星，越靠近目标时区的GNSS卫星的导航信息比较精确，因此本实施方式通过对比采集时间点靠近上一行GNSS卫星的过顶时间点，还是靠近下一行GNSS卫星的过顶时间点，来表示哪一行的GNSS卫星越来越靠近目标时区，从而一方面提高导航的精确性，另一方面，提高GNSS天线接收数据的条理性和稳定性。

步骤400、通过GNSS天线按照采集时间列表的采集时间点依次非连续性采集和跟踪筛选出的GNSS卫星组合信号，将捕获得到的粗略的码相位和载波多普勒频移参数输入到跟踪环路以便对GNSS卫星组合信号进行跟踪。

步骤500、对GNSS卫星组合信号实现稳定跟踪并得到观测量，从观测量提取得到伪距和伪距率信息并进行导航定位解算。

GNSS天线接收来自GNSS卫星组合的GNSS信号；基于GNSS卫星组合的至少一个的轨道或位置，以及GNSS卫星组合的至少一个的时钟偏移获得精确卫星信息；在获得的GNSS信号之后，识别周跳所可能影响的至少一个GNSS信号的子集，并估算对确定GNSS天线的位置或GNSS天线位置的变化有用的参数，以使用不属于周跳影响子集的至少一些获得的GNSS信号和精确卫星信息。

也即是，非连续地采集和跟踪GNSS卫星信号是指：每隔一预设周期对GNSS信号进行采集与跟踪，每次采集的GNSS信号中频数据段的长度；

其中，Lcode为GNSS信号一个伪码周期的码片数，f_c ^code为GNSS信号伪码码率，n_i为跟踪第i颗GNSS卫星信号时的相关积分时长对应的GNSS信号伪码周期数，s_s ^IF为GNSS信号中频频率，为INS辅助估算得到的第i颗GNSS卫星信号在采集的GNSS信号中频数据段起始(第i个)中频样点处的伪码相位值；

为第i颗GNSS卫星信号的伪码多普勒频率估算值，f_car为GNSS信号载波频率，为INS辅助估算得到的第i颗GNSS.卫星信号的载波多普勒频率。

以上实施例仅为本申请的示例性实施例，不用于限制本申请，本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内，对本申请做出各种修改或等同替换，这种修改或等同替换也应视为落在本申请的保护范围内。