CN115390109A - 北斗卫星集中式星座自主导航方法 - Google Patents
北斗卫星集中式星座自主导航方法 Download PDFInfo
- Publication number
- CN115390109A CN115390109A CN202210985343.2A CN202210985343A CN115390109A CN 115390109 A CN115390109 A CN 115390109A CN 202210985343 A CN202210985343 A CN 202210985343A CN 115390109 A CN115390109 A CN 115390109A
- Authority
- CN
- China
- Prior art keywords
- satellite
- algorithm
- centralized
- constellation
- orbit
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 83
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 181
- 238000001914 filtration Methods 0.000 claims abstract description 41
- 238000005259 measurement Methods 0.000 claims description 48
- 238000012795 verification Methods 0.000 claims description 23
- 238000012545 processing Methods 0.000 claims description 19
- 238000004088 simulation Methods 0.000 claims description 19
- 230000002457 bidirectional effect Effects 0.000 claims description 16
- 238000004891 communication Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000012360 testing method Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 description 14
- 238000013461 design Methods 0.000 description 8
- 230000007704 transition Effects 0.000 description 6
- 238000012937 correction Methods 0.000 description 5
- 230000007547 defect Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000001133 acceleration Effects 0.000 description 3
- 238000004873 anchoring Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
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/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/393—Trajectory determination or predictive tracking, e.g. Kalman filtering
-
- 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/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/09—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing processing capability normally carried out by 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/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/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- 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/396—Determining accuracy or reliability of position or pseudorange measurements
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供了一种北斗卫星集中式星座自主导航方法，所述北斗卫星集中式星座自主导航方法包括：采用推广卡尔曼滤波算法估计，代替地面批处理算法的最小二乘批处理估计，避免计算导航星座各星精密轨道时的实时性差；结合算法工程实现流程，使集中式导航算法具备星上使用条件。
Description
技术领域
本发明涉及卫星导航技术领域，特别涉及一种北斗卫星集中式星座自主导航方法。
背景技术
新一代北斗三号卫星将于2020年全面建成，为全球提供导航服务。相比于上一代北斗二号卫星，新一代北斗三号卫星均增加有星间链路载荷，以实现北斗导航星座中双星间的测距与通信。与GPS卫星采用的UHF链路不同，北斗卫星采用的为Ka星间链路。相比于UHF链路，Ka链路建链速度更快，抗干扰能力更强，是一种性能更为优越的星间链路设计。基于Ka星间链路的星间测距及通信能力，新一代北斗卫星均设计有卫星自主导航功能。通过实时采集星间、星地双向测距信息，可实现对卫星动力学预报轨道的修正，进而实现卫星自主导航。但与GPS卫星一样，北斗卫星也采用分布式算法以实现卫星星间测距自主导航。在分布式算法中，算法只获取与本星相关的星间测距值，仅根据星间测距值确定自身导航状态。由于分布式算法是以假设解耦为前提的，其精度将不可避免地遭受损失。
与分布式算法估计方法不同，集中式算法是采集整网卫星星间测距信息，并基于测距信息确定所有星座卫星的导航信息，可得到全局最优解。因而，从统计学角度分析，集中式算法精度远优于分布式算法。地面即是通过集中式批处理算法实现卫星轨道的精密确定。近年来，本领域技术人员均探讨了集中式导航算法的在轨实现方法。例如，为降低集中式自主定轨算法的计算量，提出了一种以长期预报星历为基础，利用观测信息实时修正广播星历变量的集中式自主定轨方法。再比如，提出了一种快速、稳定的运动学集中式定轨算法。该算法不基于轨道动力学模型，利用高阶多项式对真实轨道与长期预报星历的差异进行描述。上述虽分别设计了星上可用的集中式导航算法。但均以星上处理能力受限为前提，将算法进行了相应的简化。而随着北斗三号卫星在星载计算机和星间通信技术的突破，Ka星间链路可以实现卫星与卫星之间100k的通信速率，而星上使用的龙芯1E300处理器的峰值频率为200MHz，内存为512M。
发明内容
本发明的目的在于提供一种北斗卫星集中式星座自主导航方法，以解决现有的北斗卫星星座自主集中式导航算法因为星上处理能力受限被简化的问题。
为解决上述技术问题，本发明提供一种北斗卫星集中式星座自主导航方法，所述北斗卫星集中式星座自主导航方法包括：
采用推广卡尔曼滤波算法估计，代替地面批处理算法的最小二乘批处理估计，避免计算导航星座各星精密轨道时的实时性差；
结合算法工程实现流程，使集中式导航算法具备星上使用条件。
可选的，在所述的北斗卫星集中式星座自主导航方法中，还包括：
评估推广卡尔曼滤波算法是否受星上处理能力制约，将推广卡尔曼滤波算法注入星上使用的龙芯1E300处理器中进行仿真测试，在评估推广卡尔曼滤波算法精度的同时验证算法的工程可行性；
将推广卡尔曼滤波算法注入于龙芯1E300处理器中，以星上使用的龙芯1E300处理器为硬件基础，开展集中式导航算法的仿真验证；
通过在龙芯1E300处理器中验证集中式导航算法，对推广卡尔曼滤波算法的精度、稳定性及在轨应用可行性均进行评估。
可选的，在所述的北斗卫星集中式星座自主导航方法中，还包括：
北斗三号卫星能够星载计算机通信和星间通信，Ka星间链路实现卫星与卫星之间100k的通信速率，星上使用的龙芯1E300处理器的峰值频率为200MHz，内存为512M；
根据北斗卫星星间测距能力，3s完成双星间的相互测距，30s内完成本星与最多10颗卫星的相互测距；
北斗卫星Ka星间链路最高实现卫星与卫星之间100kb的通信速率，将整网测距信息传入主卫星中。
可选的，在所述的北斗卫星集中式星座自主导航方法中，还包括：
针对北斗卫星整网集中式定轨算法在轨计算，通过基于推广卡尔曼滤波的整网集中式导航算法方法、基于北斗卫星Ka星间链路载荷、以及龙芯1E300处理器进行北斗卫星整网集中式定轨算法的星上计算流程；
将北斗卫星整网集中式定轨算法注入于龙芯1E300处理器中，并开展集中式导航算法的仿真验证，以获得利用整网集中式导航算法；
且周期内，龙芯1E300处理器用时满足算法处理星上使用要求需求。
可选的，在所述的北斗卫星集中式星座自主导航方法中，
分布于各个子卫星上的信息采集模块根据星座建链表完成整个星座的卫星间的相互测距，获取星间双向测距信息；
所述子卫星根据星座路由表，将所述星间双向测距信息传入所述主卫星中；
承载于主卫星上的集中式导航算法及可行性验证模块建立星间测距观测模型生成模块、轨道预报算法模块及滤波估计模块，其中：
所述轨道预报算法模块依据卫星动力学模型得到星座卫星预报轨道信息；
所述星间测距观测模型生成模块对星间双向测距信息进行处理，得到算法所需的观测量并构建观测方程，形成整网星间观测信息；
所述滤波估计模块利用扩展卡尔曼滤波算法将所述整网星间测距信息作为模型，修正所述星座卫星预报轨道信息，更新整个星座的定轨信息；
所述主卫星将所述整个星座的定轨信息实时发送至各个子卫星，执行各个卫星轨道信息的更新，完成集中式自主导航算法的一次定轨；
所述集中式导航算法及可行性验证模块对所述集中式自主导航算法的精度及工程可行性进行评估，并对整网卫星建链测距用时及主卫星的星上处理器运算处理时长进行评估。
可选的，在所述的北斗卫星集中式星座自主导航方法中，所述轨道预报算法模块建立所述卫星动力学模型，并根据所述卫星动力学模型生成卫星外推轨道，所述卫星动力学模型模块将所述卫星外推轨道发送至所述星间测距观测模型生成模块；
所述星间测距观测模型生成模块将所述外推轨道作为一步预报；
所述卫星动力学模型模块还用于建立轨道状态转移矩阵预报模型，并将所述轨道状态转移矩阵预报模型发送至所述星间测距观测模型生成模块。
可选的，在所述的北斗卫星集中式星座自主导航方法中，所述卫星动力学模型由公式(1)表示：
其中，
为卫星在惯性坐标系中的位置矢量；
为卫星在惯性坐标系中的速度矢量；
为卫星受各种作用力而产生的加速度。
可选的，在所述的北斗卫星集中式星座自主导航方法中，所述星间双向测距信息由公式(2)表示：
其中，ρAB为卫星A发B收的星间距离测量值，ρBA为卫星B发A收的星间距离测量值，
与
为星间测距理论值，δtA为卫星A钟差，δtB为卫星B钟差；dErrAB与dErrBA为星间测量误差，所述星间测量误差包括收发时延误差、相位中心偏差、相对论效应误差，以及电离层延迟误差；
由公式(2)计算所述观测量，所述观测量由公式(3)表示：
根据所述观测量与两个卫星的位置关系，以及两个卫星的速度信息关系，得到公式(4)：
其中，
为卫星A的位置向量，
为卫星B的位置向量；
将公式(4)右式进行一阶泰勒公式展开，得到所述观测方程，所述观测方程由公式(5)表示：
其中，
为卫星A的预报位置矢量，
为卫星B的预报位置矢量，
为卫星A的一阶位置信息改正量，
为卫星B的一阶位置信息改正量，disAB为卫星A与卫星B之间的轨道预报距离，ε为观测噪声，计算方法如公式(6)所示；
可选的，在所述的北斗卫星集中式星座自主导航方法中，所述滤波估计模块利用扩展卡尔曼滤波算法将所述整网星间测距信息作为模型，修正所述星座卫星预报轨道信息包括：
设定：
所述扩展卡尔曼滤波过程由式(7)～(14)表示：
初始化：
状态评估：
测量更新：
其中，Xk，Xk+1表示第k，k+1时刻卫星轨道滤波更新结果；
为基于Xk得到的第k+1时刻的轨道预报信息，Φk,k+1为k时刻到k+1时刻的轨道状态转移矩阵，H为观测方程中观测矩阵，P为状态协方差矩阵，Z为算法观测量，Q为过程噪声协方差矩阵，R为观测噪声协方差矩阵。
可选的，在所述的北斗卫星集中式星座自主导航方法中，所述集中式导航算法及可行性验证模块承载于所述主卫星的星上处理器中；
所述主卫星的星上处理器的峰值频率为200MHz，内存为512M。
可选的，在所述的北斗卫星集中式星座自主导航方法中，所述集中式导航算法及可行性验证模块根据所述北斗星座的卫星运行特点，构建24颗北斗MEO星座卫星仿真轨道；
所述集中式导航算法及可行性验证模块设置仿真周期为30天，星间测距周期为5分钟，星间测距精度为0.1米；
所述集中式导航算法及可行性验证模块设置整网卫星单星平均链路数为8条；
在一个所述星间测距周期内，所述地面锚固站与4颗在轨卫星建链；
所述集中式自主导航算法注入于所述主卫星的星上处理器上运行，以评估算法周期内所需时间。
在本发明提供的北斗卫星集中式星座自主导航方法中，通过将推广卡尔曼滤波算法代替地面传统的批处理算法，并合理设计算法工程实现流程，使集中式导航算法具备了星上使用条件，实现了集中式导航算法的在轨使用；本发明不但建立了算法设计方法，还对算法的在轨实现流程进行了设计，充分评估了集中式导航算法的应用可行性，克服了现有技术中星载处理器数据处理能力限制制约星上使用全星座自主导航算法的缺陷。
本发明选择扩展卡尔曼滤波算法估计方法，实现了集中式导航算法的星上使用，避免一般情况下在地面采用集中式导航算法计算导航星座各星精密轨道时采用最小二乘批处理估计方法，克服了最小二乘批处理方法计算导航星座各星精密轨道时存在着实时性差，计算量大的缺陷。
本发明为评估算法是否受星上处理能力制约，特将算法注入星上使用的龙芯1E300处理器中进行仿真测试，在评估算法精度的同时也验证算法的工程可行性。本发明将算法注入于龙芯1E300处理器中，以星上使用的龙芯1E300处理器为硬件基础，开展集中式导航算法的仿真验证。龙芯1E300处理器是由中科龙芯公司最新研发的宇航级处理器。处理器已在北斗卫星CA34，CA35中使用。龙芯1E300处理器的。通过在龙芯1E300处理器中验证集中式导航算法，可对算法的精度，稳定性及在轨应用可行性均进行评估。
进一步的，随着北斗三号卫星在星载计算机和星间通信技术的突破，Ka星间链路可以实现卫星与卫星之间100k的通信速率，而星上使用的龙芯1E300处理器的峰值频率为200MHz，内存为512M。本发明根据北斗卫星目前星间测距能力，3s即可完成双星间的相互测距，因而30s内即可完成本星与最多10颗卫星的相互测距。且北斗卫星Ka星间链路最高可以实现卫星与卫星之间100kb的通信速率，也为将整网测距信息传入主卫星中提供了可能。
附图说明
图1是本发明一实施例的北斗卫星星座集中式自主导航方法流程示意图；
图2是本发明一实施例的北斗卫星集中式星座自主导航方法的整网卫星轨道预报平均URE误差示意图；
图3是本发明一实施例的北斗卫星集中式星座自主导航方法的卫星整网集中式分布式导航算法轨道平均URE误差示意图。
具体实施方式
以下结合附图和具体实施例对本发明提出的北斗卫星集中式星座自主导航方法作进一步详细说明。根据下面说明和权利要求书，本发明的优点和特征将更清楚。需说明的是，附图均采用非常简化的形式且均使用非精准的比例，仅用以方便、明晰地辅助说明本发明实施例的目的。
本发明的核心思想在于提供一种星上可用的北斗卫星集中式星座自主导航方法，以解决现有的北斗卫星星座自主集中式导航算法因为星上处理能力受限被简化的问题。
为实现上述思想，本发明提供了一种星上可用的北斗卫星集中式星座自主导航方法，所述星上可用的北斗卫星集中式星座自主导航方法包括分布于各个子卫星上的信息采集模块、承载于主卫星上的集中式导航算法及可行性验证模块，其中：所述信息采集模块根据星座建链表完成整个星座的卫星间的相互测距，获取星间双向测距信息；所述子卫星根据星座路由表，将所述星间双向测距信息传入所述主卫星中；所述集中式导航算法及可行性验证模块包括星间测距观测模型生成模块、轨道预报算法模块及滤波估计模块，其中：所述轨道预报算法模块依据卫星动力学模型得到星座卫星预报轨道信息；所述星间测距观测模型生成模块对星间双向测距信息进行处理，得到算法所需的观测量并构建观测方程，形成整网星间观测信息；所述滤波估计模块利用扩展卡尔曼滤波算法将所述整网星间测距信息作为模型，修正所述星座卫星预报轨道信息，更新整个星座的定轨信息；所述主卫星将所述整个星座的定轨信息实时发送至各个子卫星，执行各个卫星轨道信息的更新，完成集中式自主导航算法的一次定轨；所述集中式导航算法及可行性验证模块对所述集中式自主导航算法的精度及工程可行性进行评估，并对整网卫星建链测距用时及主卫星的星上处理器运算处理时长进行评估。
<实施例一>
本实施例提供一种北斗卫星集中式星座自主导航方法，如图1所示，所述北斗卫星集中式星座自主导航方法包括分布于各个子卫星上的信息采集模块、承载于主卫星上的集中式导航算法及可行性验证模块，其中：所述信息采集模块根据星座建链表完成整个星座的卫星间的相互测距，获取星间双向测距信息；所述子卫星根据星座路由表，将所述星间双向测距信息传入所述主卫星中；所述集中式导航算法及可行性验证模块包括星间测距观测模型生成模块、轨道预报算法模块及滤波估计模块，其中：所述轨道预报算法模块依据卫星动力学模型得到星座卫星预报轨道信息；所述星间测距观测模型生成模块对星间双向测距信息进行处理，得到算法所需的观测量并构建观测方程，形成整网星间观测信息；所述滤波估计模块利用扩展卡尔曼滤波算法将所述整网星间测距信息作为模型，修正所述星座卫星预报轨道信息，更新整个星座的定轨信息；所述主卫星将所述整个星座的定轨信息实时发送至各个子卫星，执行各个卫星轨道信息的更新，完成集中式自主导航算法的一次定轨；所述集中式导航算法及可行性验证模块对所述集中式自主导航算法的精度及工程可行性进行评估，并对整网卫星建链测距用时及主卫星的星上处理器运算处理时长进行评估。
具体的，要使集中式导航算法在轨稳定运行，除建立算法设计方法外，还需对算法的在轨实现流程加以设计。具体而言，若要在北斗卫星上实现集中式定轨算法，需首先基于星座建链表完成整网卫星间的相互测距。而后，依据星座路由表，将测距值传入星座卫星中所设的主节点卫星中。接着，主卫星运行集中式导航算法，利用整网星间测距信息得到整网星座卫星轨道信息。最后，主卫星将所更新的导航信息传入各个子卫星中，实现各卫星轨道信息的更新。
根据北斗卫星目前星间测距能力，3s即可完成双星间的相互测距，因而30s内即可完成本星与最多10颗卫星的相互测距。且北斗卫星Ka星间链路最高可以实现卫星与卫星之间100kb的通信速率，也为将整网测距信息传入主卫星中提供了可能。因而，若规定算法周期为5min，可设计如图1在轨实现流程。
进一步的，本实施例设计的基于推广卡尔曼滤波的集中式导航算法是以卫星动力学模型产生的外推轨道作为一步预报，星间双向测距值作为观测量，推广卡尔曼滤波法作为估计方法的导航算法。
具体的，在所述的卫星星座自主导航方法中，所述卫星动力学模型由公式(1)表示：
其中，
为卫星在惯性坐标系中的位置矢量；
为卫星在惯性坐标系中的速度矢量；
为卫星受各种作用力而产生的加速度。
由于北斗卫星为中高轨卫星，卫星主要受地球二体引力，非球形摄动力，日月三体引力，光压摄动力作用。因而，对于卫星所受加速度
可利用表1中轨道信息预报模型构建。同时，由于算法采用扩展卡尔曼滤波算法，卫星动力学模型中还需对轨道状态转移矩阵预报模型进行设计。如表1所示，轨道状态转移矩阵仅考虑二体中心引力与J2非球形摄动模型。
表1卫星动力学建模方法
具体的，在所述的卫星星座自主导航方法中，所述星间双向测距信息由公式(2)表示：
其中，ρAB为卫星A发B收的星间距离测量值，ρBA为卫星B发A收的星间距离测量值，
与
为星间测距理论值，δtA为卫星A钟差，δtB为卫星B钟差；dErrAB与dErrBA为星间测量误差，所述星间测量误差包括收发时延误差、相位中心偏差、相对论效应误差，以及电离层延迟误差；
由公式(2)计算所述观测量，所述观测量由公式(3)表示：
根据所述观测量与两个卫星的位置关系，以及两个卫星的速度信息关系，得到公式(4)：
其中，
为卫星A的位置向量，
为卫星B的位置向量；
将公式(4)右式进行一阶泰勒公式展开，得到所述观测方程，所述观测方程由公式(5)表示：
其中，
为卫星A的预报位置矢量，
为卫星B的预报位置矢量，
为卫星A的一阶位置信息改正量，
为卫星B的一阶位置信息改正量，disAB为卫星A与卫星B之间的轨道预报距离，计算方法如公式(6)所示；
具体的，在地面采用集中式导航算法计算导航星座各星精密轨道时，一般采用的是最小二乘批处理估计方法。最小二乘批处理方法存在着实时性差，计算量大等缺点。因而，为实现集中式导航算法星上使用，本实施例选择推广卡尔曼滤波算法估计方法。包括：
设定：
所述扩展卡尔曼滤波过程由式(7)～(14)表示：
初始化：
状态评估：
测量更新：
其中，Xk，Xk+1表示第k，k+1时刻卫星轨道滤波更新结果；
为基于Xk得到的第k+1时刻的轨道预报信息，Φk,k+1为k时刻到k+1时刻的轨道状态转移矩阵，H为观测方程中观测矩阵，P为状态协方差矩阵，Z为算法观测量，Q为过程噪声协方差矩阵，R为观测噪声协方差矩阵。
由于现有技术均指出星载处理器数据处理能力的限制是制约星上使用全星座自主导航算法的主要因素，因而本发明将算法注入于龙芯1E300处理器中，以星上使用的龙芯1E300处理器为硬件基础，开展集中式导航算法的仿真验证。龙芯1E300处理器是由中科龙芯公司最新研发的宇航级处理器。处理器已在北斗卫星CA34，CA35中使用。龙芯1E300处理器的峰值频率可达200MHz，可用内存达512M。通过在龙芯1E300处理器中验证集中式导航算法，可对算法的精度，稳定性及在轨应用可行性均进行评估。
仿真条件依据北斗星座卫星运行特点，构建24颗北斗MEO星座卫星仿真轨道。设定仿真周期为30天，星间测距周期为5min，星间测距精度为0.1m。整网卫星单星平均链路数为8条。有1个地面锚固站参与定轨，周期内锚固站与4颗在轨卫星建链。算法注入于龙芯1E300处理器上运行，以评估算法周期内所需时间。
首先，由于集中式自主导航算法是以星上动力学一步轨道预报信息为基础的导航算法，为更好评估算法定轨精度，对预报轨道精度作以分析。图2为采用所述轨道动力学的一步预报方法，30天整网卫星轨道预报URE平均误差。如图2所示，由于预报模型与真实轨道存在偏差，导致卫星轨道预报误差缓慢发散，到30天时，误差可达330m。
具体的，依据上述仿真条件，对算法定轨精度进行验证。为充分评估集中式算法精度，本发明也在相同的仿真条件下，对采用分布式算法的估计精度进行评估。图3给出了集中式与分布式整网卫星平均URE误差仿真结果。如图3所示，集中式定轨算法整网卫星平均URE误差仅为0.18m，远优于分布式导航算法2.8m的平均URE轨道误差。最后，由于算法在星上使用的龙芯1E300处理器中运行，算法的工程可行性也相应得到了验证。经仿真评估，周期内，龙芯1E300处理器57s可生成整网卫星定轨结果，可满足算法的工程实现方法的要求。
<实施例二>
本实施例还提供了一种北斗卫星集中式星座自主导航方法，如图1所示，所述北斗卫星集中式星座自主导航方法包括：分布于各个子卫星上的信息采集模块根据星座建链表完成整个星座的卫星间的相互测距，获取星间双向测距信息；所述子卫星根据星座路由表，将所述星间双向测距信息传入所述主卫星中；承载于主卫星上的集中式导航算法及可行性验证模块建立星间测距观测模型生成模块、轨道预报算法模块及滤波估计模块，其中：所述轨道预报算法模块依据卫星动力学模型得到星座卫星预报轨道信息；所述星间测距观测模型生成模块对星间双向测距信息进行处理，得到算法所需的观测量并构建观测方程，形成整网星间观测信息；所述滤波估计模块利用扩展卡尔曼滤波算法将所述整网星间测距信息作为模型，修正所述星座卫星预报轨道信息，更新整个星座的定轨信息；所述主卫星将所述整个星座的定轨信息实时发送至各个子卫星，执行各个卫星轨道信息的更新，完成集中式自主导航算法的一次定轨；所述集中式导航算法及可行性验证模块对所述集中式自主导航算法的精度及工程可行性进行评估，并对整网卫星建链测距用时及主卫星的星上处理器运算处理时长进行评估。
综上所述，在本发明提供的北斗卫星星座自主导航方法中，不但建立了算法设计方法，还对算法的在轨实现流程进行了设计，实现了集中式导航算法在轨稳定运行。本发明将扩展卡尔曼滤波算法代替地面传统的批处理算法，并合理设计算法工程实现流程，实现了集中式导航算法的在轨使用。
另外，本发明选择扩展卡尔曼滤波算法估计方法，实现了集中式导航算法的星上使用，避免一般情况下在地面采用集中式导航算法计算导航星座各星精密轨道时采用最小二乘批处理估计方法，克服了最小二乘批处理方法计算导航星座各星精密轨道时存在着实时性差，计算量大的缺陷。
最后，本发明对北斗卫星整网集中式定轨算法在轨实现方法开展研究，设计了基于推广卡尔曼滤波的整网集中式导航算法方法并基于北斗卫星Ka星间链路载荷及龙芯1E300处理器设计了算法的星上实现流程。将算法注入于龙芯1E300处理器中，并开展集中式导航算法的仿真验证表明，利用本发明设计的整网集中式导航算法，整网卫星URE精度为0.18m,远优于分布式导航算法整网URE精度2.8m。且周期内，龙芯1E300处理器用时57s钟即可完成算法处理，可充分满足星上使用要求。
综上，上述实施例对北斗卫星星座自主导航方法的不同构型进行了详细说明，当然，本发明包括但不局限于上述实施中所列举的构型，任何在上述实施例提供的构型基础上进行变换的内容，均属于本发明所保护的范围。本领域技术人员可以根据上述实施例的内容举一反三。
上述描述仅是对本发明较佳实施例的描述，并非对本发明范围的任何限定，本发明领域的普通技术人员根据上述揭示内容做的任何变更、修饰，均属于权利要求书的保护范围。
Claims (5)
1.一种北斗卫星集中式星座自主导航方法，其特征在于，所述北斗卫星集中式星座自主导航方法包括：
采用推广卡尔曼滤波算法估计，代替地面批处理算法的最小二乘批处理估计，避免计算导航星座各星精密轨道时的实时性差；
结合算法工程实现流程，使集中式导航算法具备星上使用条件。
2.如权利要求1所述的北斗卫星集中式星座自主导航方法，其特征在于，还包括：
评估推广卡尔曼滤波算法是否受星上处理能力制约，将推广卡尔曼滤波算法注入星上使用的龙芯1E300处理器中进行仿真测试，在评估推广卡尔曼滤波算法精度的同时验证算法的工程可行性；
将推广卡尔曼滤波算法注入于龙芯1E300处理器中，以星上使用的龙芯1E300处理器为硬件基础，开展集中式导航算法的仿真验证；
通过在龙芯1E300处理器中验证集中式导航算法，对推广卡尔曼滤波算法的精度、稳定性及在轨应用可行性均进行评估。
3.如权利要求2所述的北斗卫星集中式星座自主导航方法，其特征在于，还包括：
北斗三号卫星能够星载计算机通信和星间通信，Ka星间链路实现卫星与卫星之间100k的通信速率，星上使用的龙芯1E300处理器的峰值频率为200MHz，内存为512M；
根据北斗卫星星间测距能力，3s完成双星间的相互测距，30s内完成本星与最多10颗卫星的相互测距；
北斗卫星Ka星间链路最高实现卫星与卫星之间100kb的通信速率，将整网测距信息传入主卫星中。
4.如权利要求1所述的北斗卫星集中式星座自主导航方法，其特征在于，还包括：
针对北斗卫星整网集中式定轨算法在轨计算，通过基于推广卡尔曼滤波的整网集中式导航算法方法、基于北斗卫星Ka星间链路载荷、以及龙芯1E300处理器进行北斗卫星整网集中式定轨算法的星上计算流程；
将北斗卫星整网集中式定轨算法注入于龙芯1E300处理器中，并开展集中式导航算法的仿真验证，以获得利用整网集中式导航算法；
且周期内，龙芯1E300处理器用时满足算法处理星上使用要求需求。
5.如权利要求1所述的北斗卫星集中式星座自主导航方法，其特征在于，还包括：
分布于各个子卫星上的信息采集模块根据星座建链表完成整个星座的卫星间的相互测距，获取星间双向测距信息；
所述子卫星根据星座路由表，将所述星间双向测距信息传入所述主卫星中；
承载于主卫星上的集中式导航算法及可行性验证模块建立星间测距观测模型生成模块、轨道预报算法模块及滤波估计模块，其中：
所述轨道预报算法模块依据卫星动力学模型得到星座卫星预报轨道信息；
所述星间测距观测模型生成模块对星间双向测距信息进行处理，得到算法所需的观测量并构建观测方程，形成整网星间观测信息；
所述滤波估计模块利用扩展卡尔曼滤波算法将所述整网星间测距信息作为模型，修正所述星座卫星预报轨道信息，更新整个星座的定轨信息；
所述主卫星将所述整个星座的定轨信息实时发送至各个子卫星，执行各个卫星轨道信息的更新，完成集中式自主导航算法的一次定轨；
所述集中式导航算法及可行性验证模块对所述集中式自主导航算法的精度及工程可行性进行评估，并对整网卫星建链测距用时及主卫星的星上处理器运算处理时长进行评估。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210985343.2A CN115390109A (zh) | 2020-04-30 | 2020-04-30 | 北斗卫星集中式星座自主导航方法 |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210985343.2A CN115390109A (zh) | 2020-04-30 | 2020-04-30 | 北斗卫星集中式星座自主导航方法 |
| CN202010360123.1A CN111522036B (zh) | 2020-04-30 | 2020-04-30 | 星上可用的北斗卫星集中式星座自主导航系统及导航方法 |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010360123.1A Division CN111522036B (zh) | 2020-04-30 | 2020-04-30 | 星上可用的北斗卫星集中式星座自主导航系统及导航方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115390109A true CN115390109A (zh) | 2022-11-25 |
Family
ID=71911334
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010360123.1A Active CN111522036B (zh) | 2020-04-30 | 2020-04-30 | 星上可用的北斗卫星集中式星座自主导航系统及导航方法 |
| CN202210985343.2A Pending CN115390109A (zh) | 2020-04-30 | 2020-04-30 | 北斗卫星集中式星座自主导航方法 |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010360123.1A Active CN111522036B (zh) | 2020-04-30 | 2020-04-30 | 星上可用的北斗卫星集中式星座自主导航系统及导航方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (2) | CN111522036B (zh) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115342818A (zh) * | 2020-08-12 | 2022-11-15 | 中国科学院微小卫星创新研究院 | 星间链路导航方法 |
| CN112014869B (zh) * | 2020-08-12 | 2023-11-28 | 中国科学院微小卫星创新研究院 | 基于天文导航的星间链路自主导航方法及系统 |
| CN112213747B (zh) * | 2020-09-29 | 2023-12-29 | 中国科学院微小卫星创新研究院 | 粗轨航天器与北斗三号卫星星间链路信号双向捕获方法 |
| CN112394381B (zh) * | 2020-09-30 | 2021-07-23 | 中国人民解放军军事科学院国防科技创新研究院 | 基于球卫星的全自主月面导航和数据通信方法 |
| CN112764057B (zh) * | 2020-12-30 | 2021-10-22 | 中国人民解放军32039部队 | 一种同步卫星外测数据的质量监控方法、装置和电子设备 |
| CN112987044B (zh) * | 2021-02-25 | 2021-12-10 | 中国科学院微小卫星创新研究院 | 一种星座星历在轨更新方法 |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6760664B1 (en) * | 2001-06-25 | 2004-07-06 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Autonomous navigation system based on GPS and magnetometer data |
| US7372400B2 (en) * | 2005-11-07 | 2008-05-13 | The Boeing Company | Methods and apparatus for a navigation system with reduced susceptibility to interference and jamming |
| US20110238308A1 (en) * | 2010-03-26 | 2011-09-29 | Isaac Thomas Miller | Pedal navigation using leo signals and body-mounted sensors |
| CN104407321B (zh) * | 2014-11-25 | 2017-04-26 | 西安电子科技大学 | 基于量子测距的编队卫星相对轨道估计方法及其设备 |
| CN106403960B (zh) * | 2016-09-23 | 2019-07-09 | 南京航空航天大学 | 一种适用于自组织星座的星座导航方法 |
| CN107421550B (zh) * | 2017-07-25 | 2020-08-28 | 北京航空航天大学 | 一种基于星间测距的地球-Lagrange联合星座自主定轨方法 |
| CN109459017B (zh) * | 2018-12-11 | 2022-07-01 | 上海航天控制技术研究所 | 一种利用外部基准辅助的星座自主导航方法 |
| CN110208831A (zh) * | 2019-07-09 | 2019-09-06 | 中国人民解放军61540部队 | 一种实现北斗三号卫星定轨和时间同步的方法 |
| CN110426720B (zh) * | 2019-08-13 | 2023-03-28 | 中国人民解放军61540部队 | 星间测量实现geo卫星机动后快速恢复方法 |
| CN110955969A (zh) * | 2019-11-27 | 2020-04-03 | 中国科学院微小卫星创新研究院 | 基于Windows平台的导航卫星星间信息处理单元模拟器 |
-
2020
- 2020-04-30 CN CN202010360123.1A patent/CN111522036B/zh active Active
- 2020-04-30 CN CN202210985343.2A patent/CN115390109A/zh active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| CN111522036A (zh) | 2020-08-11 |
| CN111522036B (zh) | 2022-07-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111522036B (zh) | 星上可用的北斗卫星集中式星座自主导航系统及导航方法 | |
| CN111522037B (zh) | 星座同轨道面卫星自主导航方法及导航系统 | |
| CN106338753B (zh) | 一种基于地面站/星间链路/gnss联合测量的地球同步轨道星座定轨方法 | |
| CN102230969B (zh) | 一种卫星星座星间链路的长时间自主维持方法 | |
| CN110426720B (zh) | 星间测量实现geo卫星机动后快速恢复方法 | |
| Choi et al. | Onboard orbit determination using GPS observations based on the unscented Kalman filter | |
| CN102679985A (zh) | 一种应用星间跟踪的航天器星座分散化自主导航方法 | |
| Lightsey et al. | Real-time navigation for Mars missions using the Mars network | |
| CN103033188A (zh) | 基于综合孔径观测的导航卫星自主时间同步方法 | |
| CN113608427A (zh) | 一种集中式天基时间基准建立方法 | |
| CN113253312A (zh) | 联合卫星导航方法、系统、电子设备及存储介质 | |
| CN112014869A (zh) | 基于天文导航的星间链路自主导航方法及系统 | |
| CN111487660B (zh) | 一种高精度实时微纳卫星集群导航方法 | |
| CN104035110A (zh) | 应用于多模卫星导航系统中的快速卡尔曼滤波定位方法 | |
| El-Mowafy et al. | The impact of orbital and clock errors on positioning from LEO constellations and proposed orbital solutions | |
| CN112731504B (zh) | 对月球探测器自主定轨的方法及装置 | |
| Iiyama et al. | Autonomous and decentralized orbit determination and clock offset estimation of lunar navigation satellites using GPS signals and inter-satellite ranging | |
| CN117890936B (zh) | 一种低轨卫星在轨实时星间时间传递方法和系统 | |
| Zhou et al. | Emerging errors in the orientation of the constellation in GNSS autonomous orbit determination based on inter-satellite link measurements | |
| CN114740541B (zh) | 基于主从星测速模式的小行星重力场反演方法及系统 | |
| CN117111112A (zh) | 一种gnss和lnss数据融合处理方法及系统 | |
| Zhao et al. | PPP augmentation and real-time precise orbit determination for LEO satellites | |
| CN115639582A (zh) | 一种基于共视法授时的geo卫星轨道机动期间轨道测定方法 | |
| CN114280650A (zh) | 基于leo获取gnss卫星频间偏差的方法和系统 | |
| CN112578418B (zh) | 一种面向导航星座测控管理的天地联合轨道计算方法 |
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 |