一种面向低轨大型星群的多站局部感知联合定轨方法技术

本发明专利技术属于航天测量与控制技术领域，包括主要的数据流和算法设计。本发明专利技术提出的一种面向低轨大型星群的多站局部感知联合定轨方法，利用单个测站局部感知数据进行最小二乘滤波处理，生成法方程系统融入数据中心联合定轨。该方法将各种类型的数据处理过程分散在各测站，数据中心仅通过联合测站法方程系统，就可以获得全局最优轨道状态估计。与传统测定轨方案比较，通过在测控站和数据中心的联合处理，避免了大量测量数据的传输，也显著减小数据传输和中心计算的开销，增强数据中心处理的灵活性和稳健性，满足大型星群的测控管理需求。因此，本方案适用于多个测站以及星间链路支持的大型星群定轨问题。大型星群定轨问题。大型星群定轨问题。

【技术实现步骤摘要】
一种面向低轨大型星群的多站局部感知联合定轨方法


[0001]本专利技术属于航天测量与控制
，涉及一种面向低轨大型星群的基于多站局部感知信息处理的联合定轨方法。

技术介绍

[0002]星群是指一组由多颗卫星组成的卫星系统，以协同工作方式构成一个大的“虚拟航天器”，来完成对地观测、互联网通信、战场态势感知等一系列覆盖范围广、连续性要求高的空间任务。近年来，随着小卫星技术的发展，在近地轨道上(Low Earth Orbit，LEO)多颗卫星联合组网构成编队、星座或星群的运行模式备受青睐，各种大型星群发展迅猛。“依”星、星链等互联网大型星座快速部署，仅星链就规划了42000颗卫星。与传统单星工作模式相比，星群在成本、性能、以及灵活性和可靠性等方面都有巨大优势；但它也给目前以地基主导的测运控模式带来了巨大挑战。
[0003]目前，传统航天测控系统仍然依赖地基S频段测控网，主要采用单颗卫星跟踪定轨方式，即多个或单个测控站同时或者交替跟踪卫星，采集的观测数据直接发送至数据中心，数据中心采用单星定轨方式计算卫星轨道并预报未来跟踪计划，发送至测控站引导设备跟踪。传统的单星定轨方式对于单个卫星或数量很少的小规模星群仍然适用；但对于成百上千颗卫星的大型星群，面向单星定轨的数据中心集中处理模式在数据传输和计算资源等方面开销巨大，不能适应未来对星群多目标测控的任务需求。因此，需要针对大型星群应用研究全新的星群联合定轨方案。

技术实现思路

[0004]为了克服现有技术中存在的问题，本专利技术提出一种面向低轨大型星群的多站局部感知联合定轨方法，利用单个测站局部感知数据进行最小二乘滤波处理，生成法方程系统送入数据中心联合定轨，本方案适用于多个测站以及星间链路支持的大型星群的定轨问题。本专利技术采用广义的星群概念，即卫星编队和星座都是星群的特殊结构。
[0005]本专利技术的技术方案为：
[0006]一种面向低轨大型星群的多站局部感知联合定轨方法，包括步骤如下：
[0007]步骤1，在定轨弧段，测站处理局部感知的所有原始观测数据，构建多源数据融合的线性化观测方程组；
[0008]步骤2，各测站应用局部最小二乘滤波器，构建轨道状态估计的法方程系统，所述的法方程系统包括法矩阵和相应的协方差；其中，对于仅地面观测的卫星，构建独立的法方程系统；对于星间链路测量相关的卫星，构建统一的多星法方程系统；
[0009]步骤3，将各测站法方程系统发送至数据中心，数据中心采用递推最小二乘合并法方程；
[0010]步骤4，求解合并后的线性方程组，得到轨道状态的最优估计。
[0011]进一步地，所述多源数据包括测距、瞬时测速、方位角、俯仰角、星间测距的测量数
据。
[0012]上述测距、瞬时测速、方位角、俯仰角、星间测距的观测矩阵如下：
[0013]测距ρ：
[0014]瞬时测速瞬时测速
[0015][0016]方位角A：
[0017]仰角E：
[0018]星间测距R：
[0019]式中，h
r
(*)表示各类观测矩阵，即观测值对位置r(x,y,z)的偏导数矩阵；M表示地固坐标系到测站坐标系的转换矩阵，HG为惯性坐标系到地固坐标系的转换矩阵；ρ(ρ
x
,ρ
y
,ρ
z
)表示测站坐标系中站心至卫星距离矢量；卫星a、b之间的星间测距表示为
[0020]进一步地，所述对于仅测站观测的卫星a，在测站i中，根据测量关系，构建独立法方程系统如下：
[0021]对卫星a，其待估计状态构建法方程系统
[0022]或
[0023]式中，H为m
×
n维观测矩阵，m、n分别是测站i对卫星a的观测值个数、未知轨道参数的个数；满秩rank(H)＝n，为设计矩阵，包含各类观测数据类型的矩阵；P为m
×
m维观测权重矩阵，由相应的观测数据类型决定，每类观测数据的权重相同；独立法方程系统的法矩阵
[0024]独立法方程系统的法矩阵的协方差
[0025]式中，法矩阵N下标表示测站，上标表示卫星，如表示测站i单站数据生成的卫星a法方程，为后验单位权方差。
[0026]进一步地，所述对于星间链路和测站联合观测的卫星，在测站i中，测量星间链路相关的卫星，构建多星法方程系统包括
[0027]多星法方程系统的法矩阵
[0028]多星法方程系统的协方差
[0029]式中，观测矩阵H包含各类观测数据类型的矩阵，权重矩阵P中每类观测数据的权重相同；表示测站i跟踪测量相应卫星产生的测量数据集和测站i接收到相应卫星的星间链路收端数据的集合。
[0030]进一步地，所述步骤3，具体为：
[0031](1)对于仅测站观测的卫星，采用递推最小二乘估计单星轨道状态，合并多测站法方程，对于卫星a合并法方程
[0032][0033]式中，l为涉及相关测站的数量；
[0034](2)对于星间链路和测站联合观测的卫星，采用递推最小二乘估计多星轨道状态，合并多测站法方程，对于卫星b、...、d合并法方程
[0035][0036]式中，h为涉及相关测站的数量。
[0037]本专利技术的有益效果是：
[0038](1)在测控站处理局部感知观测数据，分散了数据中心的计算量，并且避免了大量的数据传送；
[0039](2)数据中心能够不损失精度，获得全局最优轨道状态估计，只需要联合处理测控站的法方程系统，显著减少了计算量；
[0040](3)数据中心能够灵活处理各测站法方程系统，可根据数据质量剔除测站。因此，本方法通过在测控站和数据中心的联合处理，在确保定轨精度的前提下，将显著减小数据传输和中心计算的开销，满足大型星群的测控管理需求，适用于多个测站以及星间链路支持的大型星群的定轨问题。
附图说明
[0041]图1是传统单星测定轨和星群测定轨示意图。
[0042]图2是传统星群定轨方案。
[0043]图3面向大型星群的局部感知联合定轨方案。
[0044]图4传统定轨方案的定轨精度(与仿真“真”轨道比较)。
[0045]图5本专利技术方案的定轨结果(与仿真“真”轨道比较)。
具体实施方式
[0046]为了更好地说明本专利技术的优点，下面结合附图对本专利技术的具体实施方式作进一步
详细的说明。
[0047]实施例1：
[0048]如图3所示，多站和星间链路联合支持星群测量定轨过程如下(为简化问题，以i、j、k三个测站和5颗卫星a、b、c、d、e为例)：
[0049]地面测站：在定轨弧段，测站i跟踪测量卫星a、b、c，产生测量数据集测站j跟踪测量全部5颗卫星，产生测量数据集测站j跟踪测量全部5颗卫星，产生测量数据集测站k跟踪测量卫星c、d、e，产生测量数据集地基测量元素包括测距ρ、测速方位角A和俯仰角E。
[0050]星间链本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种面向低轨大型星群的多站局部感知联合定轨方法，其特征在于，包括步骤如下：步骤1，在定轨弧段，测站处理局部感知的所有原始观测数据，构建多源数据融合的线性化观测方程组；步骤2，各测站应用局部最小二乘滤波器，构建轨道状态估计的法方程系统，所述的法方程系统包括法矩阵和相应的协方差；其中，对于仅地面观测的卫星，构建独立的法方程系统；对于星间链路测量相关的卫星，构建统一的多星法方程系统；步骤3，将各测站法方程系统发送至数据中心，数据中心采用递推最小二乘合并法方程；步骤4，求解合并后的线性方程组，得到轨道状态的最优估计。2.根据权利要求1所述的面向低轨大型星群的多站局部感知联合定轨方法，其特征在于，所述多源数据包括测距、瞬时测速、方位角、俯仰角、星间测距的测量数据。3.根据权利要求2所述的面向低轨大型星群的多站局部感知联合定轨方法，其特征在于，所述测距、瞬时测速、方位角、俯仰角、星间测距的观测矩阵如下：测距ρ：瞬时测速瞬时测速瞬时测速方位角A：仰角E：星间测距R：式中，h
r
(*)表示各类观测矩阵，即观测值对位置r(x,y,z)的偏导数矩阵；M表示地固坐标系到测站坐标系的转换矩阵，HG为惯性坐标系到地固坐标系的转换矩阵；ρ(ρ
x
,ρ
y
,ρ
z
)表示测站坐标系中站心至卫星距离矢量；卫星a、b之间的星间测距表示为4.根据权利要求1所述的面向低轨大型星群的多站局部感知联合定轨方法，其特征在于，所述对于仅测...

【专利技术属性】
技术研发人员：朱俊，王冲，李恒年，李杰，姜宇，龚轲杰，房亚男，李永华，王婧妍，王蓉晖，
申请(专利权)人：中国西安卫星测控中心，
类型：发明
国别省市：