CN113687392A - 一种基于gnss信号非连续跟踪的导航方法 - Google Patents

一种基于gnss信号非连续跟踪的导航方法 Download PDF

Info

Publication number
CN113687392A
CN113687392A CN202110965457.6A CN202110965457A CN113687392A CN 113687392 A CN113687392 A CN 113687392A CN 202110965457 A CN202110965457 A CN 202110965457A CN 113687392 A CN113687392 A CN 113687392A
Authority
CN
China
Prior art keywords
satellite
gnss
time
satellites
tracking
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
CN202110965457.6A
Other languages
English (en)
Other versions
CN113687392B (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.)
Shenzhen Lexin Information Technology Co ltd
Original Assignee
Shenzhen Dianka Measurement And Control 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 Shenzhen Dianka Measurement And Control Technology Co ltd filed Critical Shenzhen Dianka Measurement And Control Technology Co ltd
Priority to CN202110965457.6A priority Critical patent/CN113687392B/zh
Publication of CN113687392A publication Critical patent/CN113687392A/zh
Application granted granted Critical
Publication of CN113687392B publication Critical patent/CN113687392B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
    • G01S19/256—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to timing, e.g. time of week, code phase, timing offset
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
    • G01S19/258—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to the satellite constellation, e.g. almanac, ephemeris data, lists of satellites in view
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13—Receivers
    • G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/27—Acquisition or tracking or demodulation of signals transmitted by the system creating, predicting or correcting ephemeris or almanac data within the receiver
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining 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/42—Determining position

Abstract

本发明公开了一种基于GNSS信号非连续跟踪的导航方法，批量生成GNSS卫星轨迹数据，并基于卫星轨迹数据确定多GNSS卫星在目标时区的过顶预报时间；按照过顶预报时间的时间顺序预测依次在目标时区内过顶的卫星序列，且根据GNSS卫星的导航定位需求确定单次跟踪的卫星接收数量；确定非连续采集的采集时间列表，每次采集时间点基于卫星接收数量通过组合筛选方式从卫星序列中确定用于跟踪的GNSS卫星组合；采集GNSS卫星组合信号，将捕获得到的粗略的码相位和载波多普勒频移参数输入到跟踪环路；对GNSS卫星组合信号实现稳定跟踪并得到观测量，从观测量提取得到伪距和伪距率信息并进行导航定位解算；本发明实现多卫星非连续导航，且通过组合筛选卫星提高导航的精确性。

Description

一种基于GNSS信号非连续跟踪的导航方法
技术领域
本发明涉及GNSS导航技术领域，具体涉及一种基于GNSS信号非连续跟踪的导航方法。
背景技术
全球导航卫星系统(GNSS)包括全球定位系统GPS(美国)，GLONASS(俄罗斯)，伽俐略(欧洲)和COMPASS(中国)，GNSS典型地使用多个绕地球飞行的卫星。这多个卫星形成了卫星的星座。GNSS接收器检测到在卫星所广播的电磁信号上调制的伪随机噪声(PRN)代码。也将这个代码称为测距码。代码检测包括将在所述广播信号上调制的位序列与接收器所生成的即将检测到的代码的复制品比较。基于用于一系列卫星中每一个的代码到达时间的检测，GNSS接收器估算其位置，定位包括地理定位，即地球表面的定位。
全球导航卫星系统是能在地球表面或近地空间的任何地点为用户提供全天候的3维坐标和速度以及时间信息的空基无线电导航定位系统。因此，通俗一点说，如果你除了要知道经纬度还想知道高度的话，那么，必须对收到4颗卫星才能准确定位。
然而当卫星数量超过4颗，或者卫星经过目标时区的过顶时间差比较相近，那么在使用GNSS信号完成非连续跟踪导航时，则无法依据现有的GNSS接收方法稳定接收GNSS信号，造成接收操作紊乱，甚至导航出错的问题。
发明内容
本发明的目的在于提供一种基于GNSS信号非连续跟踪的导航方法，以解决现有技术中当卫星数量超过4颗，或者卫星经过目标时区的过顶时间差比较相近时，无法依据现有的GNSS接收方法稳定接收GNSS信号，造成接收操作紊乱，甚至导航出错的问题。
为解决上述技术问题，本发明具体提供下述技术方案：
一种基于GNSS信号非连续跟踪的导航方法，包括以下步骤：
步骤100、批量生成GNSS卫星轨迹数据，并基于卫星轨迹数据确定多GNSS卫星在目标时区的过顶预报时间，其中，所述过顶预报时间为根据所述GNSS卫星的时间-位置数据确定的经过所述目标时区头顶的时间；
步骤200、按照过顶预报时间的时间顺序预测依次在所述目标时区内过顶的卫星序列，且根据GNSS卫星的导航定位需求确定单次跟踪的卫星接收数量；
步骤300、确定非连续采集的采集时间列表，每次采集时间点基于所述卫星接收数量通过组合筛选方式从卫星序列中确定用于跟踪的GNSS卫星组合，所述组合筛选方式用于在至少两次所述采集时间点对应的卫星序列相同且卫星序列的数量大于所述卫星接收数量时，通过组合算法确定每次采集时间点对应的GNSS卫星组合；
步骤400、通过GNSS天线按照采集时间列表的采集时间点依次非连续性采集和跟踪筛选出的GNSS卫星组合信号，将捕获得到的粗略的码相位和载波多普勒频移参数输入到跟踪环路以便对所述GNSS卫星组合信号进行跟踪；
步骤500、对所述GNSS卫星组合信号实现稳定跟踪并得到观测量，从所述观测量提取得到伪距和伪距率信息并进行导航定位解算。
作为本发明的一种优选方案，在步骤100中，基于卫星轨迹数据确定多GNSS卫星在目标时区的过顶预报时间的实现步骤为：
步骤101、通过两行轨道根数TLE确定所有GNSS卫星运动轨迹的轨道参数，并利用SGP4模型将TLE轨道参与轨道运算模型结合测算卫星实时状态，得到所有GNSS卫星的时间-位置数据；
步骤102、将时间-位置数据写成CZML数据文件，利用循环算法获取连续时间点的位置坐标；
步骤103、根据连续时间点的位置坐标获取每个所述GNSS卫星经过目标时区的过顶预报时间。
作为本发明的一种优选方案，在步骤102中，所述CZML数据文件为GNSS卫星轨道CZML数据，通过确定生成GNSS卫星的时间-位置样本的时间间隔t，同时利用CZML提供的内插属性，将时间-位置样本数据之间通过高阶拉格朗日算法内插填充样本数据间的时间间隔t内，其中，每个GNSS卫星的CZML数据只记录该GNSS卫星在额定秒数内的轨迹状况，每个CZML数据文件都需要Python调用SPG4库，根据额定秒数和时间间隔t计算得到该GNSS卫星一系列的时间-位置数据。
作为本发明的一种优选方案，从所述时间-位置数据内筛选GNSS卫星经过所述目标时区的过顶预报时间，并将过顶预报时间和对应的GNSS卫星组成卫星序列。
作为本发明的一种优选方案，在步骤200中，所述卫星序列为多行多列式，所述卫星序列的同一行为在同一过顶时间点在所述目标时区内过顶的所述GNSS卫星，所述卫星序列的不同行从上到下表示按照时间顺序依次经过所述目标时区的GNSS卫星排列顺序；
所述GNSS卫星的导航定位需求确定单次跟踪的卫星接收数量为四个，同一时间点在所述目标时区内过顶的所述GNSS卫星大于等于所述卫星接收数量。
作为本发明的一种优选方案，在步骤300中，先根据所述采集时间点与所述卫星序列中的两行数据的过顶时间点计算所述采集时间点的跟踪偏向，再根据跟踪偏向选择一行所述GNSS卫星作为目标筛选卫星。
作为本发明的一种优选方案，所述每次采集时间点基于所述卫星接收数量通过组合筛选方式从目标筛选卫星组中确定用于跟踪的GNSS卫星组合，其中，在多次所述采集时间点对应的所述目标筛选卫星相同时，利用所述组合筛选方式使得每次的采集时间点对应的所述GNSS卫星组合不同。
作为本发明的一种优选方案，在所述步骤300中，计算所述采集时间点的跟踪偏向的实现步骤为：
计算所述卫星序列中两行所述GNSS卫星的过顶时间点确定两个所述过顶时间点的中心时间点，所述中心时间点将两行所述GNSS卫星分为上一行所述GNSS卫星的过顶时间点与中心时间点之间的上阈值，以及下一行所述GNSS卫星的过顶时间点与中心时间点之间的下阈值；
判断所述采集时间点处于所述上阈值或下阈值，来选择其中一行所述GNSS卫星作为目标筛选卫星。
作为本发明的一种优选方案，在步骤400中，GNSS天线接收来自所述GNSS卫星组合的GNSS信号；基于所述GNSS卫星组合的至少一个的轨道或位置，以及所述GNSS卫星组合的至少一个的时钟偏移获得精确卫星信息；在获得的GNSS信号之后，识别周跳所可能影响的至少一个GNSS信号的子集，并估算对确定GNSS天线的位置或所述GNSS天线位置的变化有用的参数，以使用不属于所述周跳影响子集的至少一些获得的GNSS信号和精确卫星信息。
本发明与现有技术相比较具有如下有益效果：
本发明在同一采集时间过顶的卫星超过四颗时，通过组合算法确定每一次采集对应的GNSS卫星组合，且不同采集时间点对应的GNSS卫星组合不同，从而实现多卫星导航并避免其中一个卫星出错，而导致导航出错的问题，另外，根据采集时间点与卫星过顶时间点的对比，可筛选出最靠近目标时区的卫星序列，进一步的提高导航的精确性。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案，下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是示例性的，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图引伸获得其它的实施附图。
图1为本发明实施例提供的非连续跟踪导航方法的流程示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图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信号中频数据段的长度；
Figure BDA0003223661170000091
其中，Lcode为GNSS信号一个伪码周期的码片数，fc code为GNSS信号伪码码率，ni为跟踪第i颗GNSS卫星信号时的相关积分时长对应的GNSS信号伪码周期数，ss IF为GNSS信号中频频率，
Figure BDA0003223661170000092
为INS辅助估算得到的第i颗GNSS卫星信号在采集的GNSS信号中频数据段起始(第i个)中频样点处的伪码相位值；
Figure BDA0003223661170000093
为第i颗GNSS卫星信号的伪码多普勒频率估算值，fcar为GNSS信号载波频率，
Figure BDA0003223661170000094
为INS辅助估算得到的第i颗GNSS.卫星信号的载波多普勒频率。
以上实施例仅为本申请的示例性实施例，不用于限制本申请，本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内，对本申请做出各种修改或等同替换，这种修改或等同替换也应视为落在本申请的保护范围内。

Claims (9)

1.一种基于GNSS信号非连续跟踪的导航方法，其特征在于，包括以下步骤：
步骤100、批量生成GNSS卫星轨迹数据，并基于卫星轨迹数据确定多GNSS卫星在目标时区的过顶预报时间，其中，所述过顶预报时间为根据所述GNSS卫星的时间-位置数据确定的经过所述目标时区头顶的时间；
步骤200、按照过顶预报时间的时间顺序预测依次在所述目标时区内过顶的卫星序列，且根据GNSS卫星的导航定位需求确定单次跟踪的卫星接收数量；
步骤300、确定非连续采集的采集时间列表，每次采集时间点基于所述卫星接收数量通过组合筛选方式从卫星序列中确定用于跟踪的GNSS卫星组合，所述组合筛选方式用于在至少两次所述采集时间点对应的卫星序列相同且卫星序列的数量大于所述卫星接收数量时，通过组合算法确定每次采集时间点对应的GNSS卫星组合；
步骤400、通过GNSS天线按照采集时间列表的采集时间点依次非连续性采集和跟踪筛选出的GNSS卫星组合信号，将捕获得到的粗略的码相位和载波多普勒频移参数输入到跟踪环路以便对所述GNSS卫星组合信号进行跟踪；
步骤500、对所述GNSS卫星组合信号实现稳定跟踪并得到观测量，从所述观测量提取得到伪距和伪距率信息并进行导航定位解算。
2.根据权利要求1所述的一种基于GNSS信号非连续跟踪的导航方法，其特征在于：在步骤100中，基于卫星轨迹数据确定多GNSS卫星在目标时区的过顶预报时间的实现步骤为：
步骤101、通过两行轨道根数TLE确定所有GNSS卫星运动轨迹的轨道参数，并利用SGP4模型将TLE轨道参与轨道运算模型结合测算卫星实时状态，得到所有GNSS卫星的时间-位置数据；
步骤102、将时间-位置数据写成CZML数据文件，利用循环算法获取连续时间点的位置坐标；
步骤103、根据连续时间点的位置坐标获取每个所述GNSS卫星经过目标时区的过顶预报时间。
3.根据权利要求2所述的一种基于GNSS信号非连续跟踪的导航方法，其特征在于：在步骤102中，所述CZML数据文件为GNSS卫星轨道CZML数据，通过确定生成GNSS卫星的时间-位置样本的时间间隔t，同时利用CZML提供的内插属性，将时间-位置样本数据之间通过高阶拉格朗日算法内插填充样本数据间的时间间隔t内，其中，每个GNSS卫星的CZML数据只记录该GNSS卫星在额定秒数内的轨迹状况，每个CZML数据文件都需要Python调用SPG4库，根据额定秒数和时间间隔t计算得到该GNSS卫星一系列的时间-位置数据。
4.根据权利要求3所述的一种基于GNSS信号非连续跟踪的导航方法，其特征在于，从所述时间-位置数据内筛选GNSS卫星经过所述目标时区的过顶预报时间，并将过顶预报时间和对应的GNSS卫星组成卫星序列。
5.根据权利要求2所述的一种基于GNSS信号非连续跟踪的导航方法，其特征在于，在步骤200中，所述卫星序列为多行多列式，所述卫星序列的同一行为在同一过顶时间点在所述目标时区内过顶的所述GNSS卫星，所述卫星序列的不同行从上到下表示按照时间顺序依次经过所述目标时区的GNSS卫星排列顺序；
所述GNSS卫星的导航定位需求确定单次跟踪的卫星接收数量为四个，同一时间点在所述目标时区内过顶的所述GNSS卫星大于等于所述卫星接收数量。
6.根据权利要求5所述的一种基于GNSS信号非连续跟踪的导航方法，其特征在于：在步骤300中，先根据所述采集时间点与所述卫星序列中的两行数据的过顶时间点计算所述采集时间点的跟踪偏向，再根据跟踪偏向选择一行所述GNSS卫星作为目标筛选卫星。
7.根据权利要求6所述的一种基于GNSS信号非连续跟踪的导航方法，其特征在于：所述每次采集时间点基于所述卫星接收数量通过组合筛选方式从目标筛选卫星组中确定用于跟踪的GNSS卫星组合，其中，在多次所述采集时间点对应的所述目标筛选卫星相同时，利用所述组合筛选方式使得每次的采集时间点对应的所述GNSS卫星组合不同。
8.根据权利要求6所述的一种基于GNSS信号非连续跟踪的导航方法，其特征在于，在所述步骤300中，计算所述采集时间点的跟踪偏向的实现步骤为：
计算所述卫星序列中两行所述GNSS卫星的过顶时间点确定两个所述过顶时间点的中心时间点，所述中心时间点将两行所述GNSS卫星分为上一行所述GNSS卫星的过顶时间点与中心时间点之间的上阈值，以及下一行所述GNSS卫星的过顶时间点与中心时间点之间的下阈值；
判断所述采集时间点处于所述上阈值或下阈值，来选择其中一行所述GNSS卫星作为目标筛选卫星。
9.根据权利要求1所述的一种基于GNSS信号非连续跟踪的导航方法，其特征在于，在步骤400中，GNSS天线接收来自所述GNSS卫星组合的GNSS信号；基于所述GNSS卫星组合的至少一个的轨道或位置，以及所述GNSS卫星组合的至少一个的时钟偏移获得精确卫星信息；在获得的GNSS信号之后，识别周跳所可能影响的至少一个GNSS信号的子集，并估算对确定GNSS天线的位置或所述GNSS天线位置的变化有用的参数，以使用不属于所述周跳影响子集的至少一些获得的GNSS信号和精确卫星信息。
CN202110965457.6A 2021-08-23 2021-08-23 一种基于gnss信号非连续跟踪的导航方法 Active CN113687392B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110965457.6A CN113687392B (zh) 2021-08-23 2021-08-23 一种基于gnss信号非连续跟踪的导航方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110965457.6A CN113687392B (zh) 2021-08-23 2021-08-23 一种基于gnss信号非连续跟踪的导航方法

Publications (2)

Publication Number Publication Date
CN113687392A true CN113687392A (zh) 2021-11-23
CN113687392B CN113687392B (zh) 2024-04-12

Family

ID=78581201

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110965457.6A Active CN113687392B (zh) 2021-08-23 2021-08-23 一种基于gnss信号非连续跟踪的导航方法

Country Status (1)

Country Link
CN (1) CN113687392B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115032671A (zh) * 2022-08-11 2022-09-09 成都国星宇航科技股份有限公司 一种低轨卫星跟踪预报时间段计算方法和装置

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101424731A (zh) * 2007-10-31 2009-05-06 中国科学院微电子研究所 全球定位系统接收机信号缺失情况下快速重捕和定位的方法
CN102498414A (zh) * 2009-09-19 2012-06-13 天宝导航有限公司 用以估计轨道的gnss信号处理
US8228230B1 (en) * 2010-04-07 2012-07-24 Stmicroelectronics S.R.L. Satellite orbit prediction method, receiving apparatus and navigation satellite system employing the method
US20120306689A1 (en) * 2011-06-01 2012-12-06 Andrew Thomas Yule Hybrid satellite positioning with prediction
CN105044745A (zh) * 2015-07-15 2015-11-11 中国人民解放军理工大学 一种圆轨道低轨卫星过顶剩余可见时长预测方法
US20160259061A1 (en) * 2015-03-06 2016-09-08 Gatekeeper Systems, Inc. Low-energy consumption location of movable objects
CN105954776A (zh) * 2016-05-14 2016-09-21 四川中卫北斗科技有限公司 一种导航信号接收方法和接收机
CN109459776A (zh) * 2018-10-08 2019-03-12 上海交通大学 基于gnss信号非连续跟踪的gnss/ins深组合导航方法
CN111751789A (zh) * 2020-06-30 2020-10-09 北京无线电测量研究所 人造卫星经过雷达探测范围的预报方法、系统、介质及设备

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101424731A (zh) * 2007-10-31 2009-05-06 中国科学院微电子研究所 全球定位系统接收机信号缺失情况下快速重捕和定位的方法
CN102498414A (zh) * 2009-09-19 2012-06-13 天宝导航有限公司 用以估计轨道的gnss信号处理
US8228230B1 (en) * 2010-04-07 2012-07-24 Stmicroelectronics S.R.L. Satellite orbit prediction method, receiving apparatus and navigation satellite system employing the method
US20120306689A1 (en) * 2011-06-01 2012-12-06 Andrew Thomas Yule Hybrid satellite positioning with prediction
US20160259061A1 (en) * 2015-03-06 2016-09-08 Gatekeeper Systems, Inc. Low-energy consumption location of movable objects
CN105044745A (zh) * 2015-07-15 2015-11-11 中国人民解放军理工大学 一种圆轨道低轨卫星过顶剩余可见时长预测方法
CN105954776A (zh) * 2016-05-14 2016-09-21 四川中卫北斗科技有限公司 一种导航信号接收方法和接收机
CN109459776A (zh) * 2018-10-08 2019-03-12 上海交通大学 基于gnss信号非连续跟踪的gnss/ins深组合导航方法
CN111751789A (zh) * 2020-06-30 2020-10-09 北京无线电测量研究所 人造卫星经过雷达探测范围的预报方法、系统、介质及设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ROBERT ODOLINSKI ET AL.: "Low-Cost, High-Precision, Single-Frequency GPS-BDS RTK Postioning", 《GPS SOLUTIONS》, vol. 21, no. 3, pages 1315 - 1330, XP036261692, DOI: 10.1007/s10291-017-0613-x *
刘婧: "GNSS接收机基带处理电路低功耗设计综述", 《数字技术与应用》, no. 4, pages 166 - 167 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115032671A (zh) * 2022-08-11 2022-09-09 成都国星宇航科技股份有限公司 一种低轨卫星跟踪预报时间段计算方法和装置

Also Published As

Publication number Publication date
CN113687392B (zh) 2024-04-12

Similar Documents

Publication Publication Date Title
EP2486418B1 (en) Improvements in or relating to radio navigation
US8525727B2 (en) Position and velocity uncertainty metrics in GNSS receivers
US8825375B2 (en) Snap-to-road using wireless access point data
CN102033236B (zh) 一种卫星导航位置速度联合估计方法
RU2432584C2 (ru) Способ определения координат мобильного приемника спутниковой радионавигационной системы (срнс)
KR100712237B1 (ko) 지역적 전리층 오차 모델링 방법 및 그를 이용한 오차 보정방법
CN107607971B (zh) 基于gnss共视时间比对算法的时间频率传递方法及接收机
CN113848579B (zh) Ins辅助gnss定位的粗差剔除方法及系统
CN109633718B (zh) 一种归一化加权最小二乘导航定位方法
CN108181633A (zh) 一种gnss时间频率传递接收机及接收方法
CN113687392B (zh) 一种基于gnss信号非连续跟踪的导航方法
US20170097422A1 (en) Method and system for positioning and timing of a radionavigation receiver
Bisnath et al. Innovation: Examining precise point positioning now and in the future
Beuchert et al. Snappergps: Algorithms for energy-efficient low-cost location estimation using gnss signal snapshots
US20230126547A1 (en) Machine learning in gnss receivers for improved velocity outputs
Akopian et al. A network aided iterated LS method for GPS positioning and time recovery without navigation message decoding
Schubert et al. Integrity of navigation for land users: Study concept and simulator architecture
Innerkofler et al. GNSS radio occultation excess-phase processing for climate applications including uncertainty estimation
Le et al. Use of global and regional ionosphere maps for single-frequency precise point positioning
RU2614039C2 (ru) Способ определения показателя надежности, связанного с траекторией движения подвижного объекта
Vani et al. Evaluation of GPS data for navigational solution and error reduction using kalman filter
Spangenberg Safe navigation for vehicles
KR20190006560A (ko) 위성 신호 처리를 위한 방법 및 시스템
ABDELGAFAR Atmospheric water vapor maps generations from stochastic interpolation of GNSS Zenith Tropospheric Delays
Wang et al. Multi Constellation Fast Selection Method Based on Equivalent Distance Error

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20240307

Address after: 518000, Building A, Xinglongfu Maker Park, Yousong Road, Yousong Community, Longhua Street, Longhua District, Shenzhen City, Guangdong Province 607W

Applicant after: Shenzhen Lexin Information Technology Co.,Ltd.

Country or region after: China

Address before: 518000 317, mingzesheng complex building, No. 12, Industrial Park Road, Tongsheng community, Dalang street, Longhua District, Shenzhen, Guangdong

Applicant before: Shenzhen Dianka measurement and Control Technology Co.,Ltd.

Country or region before: China

GR01 Patent grant
GR01 Patent grant