一种用于多航天器追逃博弈任务的动态匹配方法组成比例

本发明专利技术公开的一种用于多航天器追逃博弈任务的动态匹配方法，属于航空航天技术领域。本发明专利技术基于航天器动力学模型，将多对多问题转化为若干个一对一子问题，即将多对多追逃博弈问题拆解为顶层优化和底层优化分别解算，在顶层优化中，将追逃双方的状态信息转化为对应的权重矩阵，进而解决逃逸器和追踪器间的最优匹配问题；在底层优化中，求解追逃双方最优控制策略。本发明专利技术通过评估追逃双方的机动能力、位置和速度就能生成有效的任务匹配策略。根据本发明专利技术求解确定的追逃双方的最优控制策略，实现追逃双方动态最优控制，使追踪方拦截的逃逸器数目最多，且终端距离较小，实现追逃双方的最优追击。本发明专利技术具有控制方法简单、计算量小、实时性强等特点。时性强等特点。时性强等特点。

【技术实现步骤摘要】
一种用于多航天器追逃博弈任务的动态匹配方法


[0001]本专利技术涉及一种多航天器追逃博弈的任务匹配方法，尤其涉及一种考虑航天器状态变化的动态任务匹配方法，适用于躲避或追踪多颗卫星或空间碎片的任务，属于航空航天


技术介绍

[0002]近年来，全球航天技术迅猛发展，航天发射频率持续提高，商业低轨卫星巨型星座开启井喷式发展浪潮，对航天器多对多追逃博弈的研究提出了迫切需求。目前，对于航天器追逃博弈问题的研究大多局限于一追一逃和多追一逃的范畴，鲜少涉及多追多逃问题。
[0003]现有的一追一逃问题多基于二人零和博弈模型来求解双方策略，有学者研究了椭圆轨道上侦察卫星和非合作目标之间的微分博弈问题和近地轨道的双航天器追逃问题，提出的主要算法有间接式启发算法、半直接控制参数化方法、基于多重打靶的混合算法和灵敏度法等。对于不完全信息下的轨道追逃博弈问题，针对博弈中的规避者，有学者提出了一种基于估计的最优策略和一种改进的在线参数估计算法。随着航天技术和追逃理论的不断发展，对航天器追逃博弈问题的研究逐渐向多追一逃和追逃防三方博弈拓展深化，应用的主要算法有间接启发算法、水平集法、分布式在线任务规划算法和集中式任务规划算法等。
[0004]随着航天技术不断发展，未来的航天器追逃博弈将呈现集群化、智能化及对抗化的特点，现有的轨道博弈形式难以满足实际的任务需求，有必要开展多对多航天器追逃博弈研究，这也是追逃理论深入发展的必然走向。此外，空间中遍布各个轨道的巨型星座、废弃航天器及其碎片严重挤占了地球轨道资源，严重威胁国家安全和整个人类的太空资源安全。为缓解太空资源紧张，保障航天器安全运行，有必要开展空间碎片清除任务，还需针对巨型星座可能实施的监视和拦截行为做出应对措施。对于空间碎片清除，通过一次发射任务释放多颗小卫星去清理多颗碎片的性价比较高；对于巨型星座实施的监视和拦截等行为，该场景往往涉及多个航天器。这是典型的军民两用任务，它们均涉及了多目标最优匹配与轨道优化问题，值得深入研究。为了使研究更具普适性，本专利技术将碎片和巨型星座均视为具有机动能力的非合作目标，转而研究多对多航天器追逃博弈问题，提出了一种有效解决多航天器任务匹配的方法。

技术实现思路

[0005]为了解决多航天器追逃博弈的任务分配问题，针对航天器状态会随着时间发生变化进而影响追逃策略，本专利技术的主要目的是提供一种用于多航天器追逃博弈任务的动态分配方法，基于航天器动力学模型，将多对多问题转化为若干个一对一子问题，进而求解追逃双方的最优控制策略，即本专利技术仅通过评估追逃双方航天器的机动能力、位置和速度信息，给定任务时间后，即能够给出追逃双方的动态匹配方案和最优控制策略，具有控制方法简单、计算量小、实时性强等特点。本专利技术能够保证追方航天器获得整体最优的追踪效果。
[0006]本专利技术的目的是通过下述技术方案实现的：
[0007]本专利技术公开的一种用于多航天器追逃博弈任务的动态匹配方法，包括如下步骤：
[0008]步骤一：建立简化追逃博弈模型涉及的相关坐标系，选择与追逃航天器距离较近的一点作为参考点，参考点沿Kepler轨道围绕地球运行。以参考点为原点建立参考坐标系Sxyz，Sx轴由地心指向参考点，Sy轴位于轨道平面内指向参考点运动方向，与Sx轴垂直，Sz轴满足右手定则。
[0009]从地球到参考点的距离远大于追逃双方的相对距离，采用下述方程来描述追逃航天器相对参考点的运动，形式如下：
[0010][0011][0012][0013][0014][0015]其中，下标k取E或P，分别表示逃逸器和追踪器；X为航天器的状态量，包含位置和速度信息，x,y,z表示航天器相对参考点的位置分量；u为航天器的控制量输入，u
x
,u
y
,u
z
表示相应的加速度分量；ω表示参考点的轨道角速度；A和B分别为状态量和控制量的系数矩阵。
[0016]步骤二：将追逃双方的机动能力、相对距离和相对速度量化为权重矩阵。
[0017](1)构造判断矩阵：追逃双方的机动能力、相对距离和相对速度是追踪器选择追击目标时需要考虑的因素。将所述三个因素两两比较，根据表1所示的1
‑
9标度方法来构造判断矩阵C＝(c
ij
)3×3。其中，c
ij
表示因素i相对j的判断，c
ji
＝1/c
ij
表示因素j相对i的判断；
[0018]表1 1
‑
9标度方法
[0019][0020](2)层次单排序：计算判断矩阵最大特征根λ的归一化特征向量w，其各元素表示各因素对于其他因素的相对重要性权值。
[0021](3)一致性检测：定义一致性比率CR＝CI/RI和一致性指标CI＝(λ
‑
n)/(n
‑
1)。RI为
随机一致性指标，其取值与因素个数n有关。当CR小于0.1时，判定判断矩阵的不一致程度在容许范围之内，检验通过，否则需重新构造判断矩阵，即重复步骤二(1)。
[0022](4)建立因素权重：将追逃双方的机动能力、相对距离和相对速度转化为对应的权重值。
[0023][0024][0025][0026]其中，取值下标i,j为逃逸器和追踪器编号；μ
a
,μ
d
,μ
v
分别为追逃双方机动能力、距离和相对速度的权重值；v为航天器速度；a为加速度幅值；r为追逃双方相对距离；D1和D2分别表示追逃双方的最近距离和最远距离，由双方位置信息得出；k
d
为距离计算系数，根据经验取值；k
v
为速度方向系数，当逃逸器数目少于追踪器时取0.25，其他情形取1。在公式(6)中，机动能力权重在0到1之间，若逃逸器的机动能力相对追踪器越强，则追踪难度越高，该权重值越接近0，公式(8)与此类似。公式(7)符合追踪器就近选择逃逸目标的原则，双方距离越近，距离权值越接近1。
[0027](5)计算权重矩阵：根据逃逸器和追踪器的机动能力、相对距离和相对速度建立隶属度函数，权重矩阵的各元素T
ij
计算方式如下：
[0028][0029]其中，μ
ij
为步骤(4)计算所得的因素权重矩阵。
[0030]由于在限定时间内，追逃双方的机动能力、相对距离和相对速度大小均能够影响最终的追逃效果，需要比较三个因素的重要程度。由于航天器均采用小推力方式进行机动，在给定时间内速度的改变较小，判定机动能力的重要程度最弱。至此，至于剩下两个因素中哪一个更重要，给定追逃航天器的数目和状态信息，按照相对速度、相对距离、机动能力和相对距离、相对速度、机动能力两种重要性序列分别生成权重矩阵和匹配方案，对比追逃结果的优劣。结果显示，三个因素的重要性由高到低排序为：相对速度、相对距离、机动能力，由此可产生6种不同的符合条件的判断矩阵，本专利技术根据其在实例中的应用结果来选取最优的判断本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种用于多航天器追逃博弈任务的动态匹配方法，其特征在于，包括如下步骤：步骤一：建立简化追逃博弈模型涉及的相关坐标系，选择与追逃航天器距离较近的一点作为参考点，参考点沿Kepler轨道围绕地球运行；以参考点为原点建立参考坐标系Sxyz，Sx轴由地心指向参考点，Sy轴位于轨道平面内指向参考点运动方向，与Sx轴垂直，Sz轴满足右手定则；从地球到参考点的距离远大于追逃双方的相对距离，采用下述方程来描述追逃航天器相对参考点的运动，形式如下：相对参考点的运动，形式如下：相对参考点的运动，形式如下：相对参考点的运动，形式如下：相对参考点的运动，形式如下：其中，下标k取E或P，分别表示逃逸器和追踪器；X为航天器的状态量，包含位置和速度信息，x,y,z表示航天器相对参考点的位置分量；u为航天器的控制量输入，u
x
,u
y
,u
z
表示相应的加速度分量；ω表示参考点的轨道角速度；A和B分别为状态量和控制量的系数矩阵；步骤二：将追逃双方的机动能力、相对距离和相对速度量化为权重矩阵；(1)构造判断矩阵：追逃双方的机动能力、相对距离和相对速度是追踪器选择追击目标时需要考虑的因素；将所述三个因素两两比较，根据表1所示的1
‑
9标度方法来构造判断矩阵C＝(c
ij
)3×3；其中，c
ij
表示因素i相对j的判断，c
ji
＝1/c
ij
表示因素j相对i的判断；表1 1
‑
9标度方法(2)层次单排序：计算判断矩阵最大特征根λ的归一化特征向量w，其各元素表示各因素对于其他因素的相对重要性权值；(3)一致性检测：定义一致性比率CR＝CI/RI和一致性指标CI＝(λ
‑
n)/(n
‑
1)；RI为随机一致性指标，其取值与因素个数n有关；当CR小于0.1时，判定判断矩阵的不一致程度在容许
范围之内，检验通过，否则需重新构造判断矩阵，即重复步骤二(1)；(4)建立因素权重：将追逃双方的机动能力、相对距离和相对速度转化为对应的权重值；值；值；其中，取值下标i,j为逃逸器和追踪器编号；μ
a
,μ
d
,μ
v
分别为追逃双方机动能力、相对距离和相对速度的权重值；v为航天器速度；a为加速度幅值；r为追逃双方相对距离；D1和D2分别表示追逃双方的最近距离和最远距离，由双方位置信息得出；k
d
为距离计算系数，根据经验取值；k
v
为速度方向系数；在公式(6)中，机动能力权重在0到1之间，若逃逸器的机动能力相对追踪器越强，则追踪难度越高，该权重值越接近0，公式(8)与此类似；公式(7)符合追踪器就近选择逃逸目标的原则，双方距离越近，距离权值越接近1；(5)计算权重矩阵：根据逃逸器和追踪器的机动能力、相对距离和相对速度建立隶属度函数，权重矩阵的各元素T
ij
计算方式如下：其中，μ
ij
为步骤(4)计算所得的因素权重矩阵；由于在限定时间内，追逃双方的机动能力、相对距离和相对速度大小均能够影响最终的追逃效果，需要比较三个因素的重要程度；由于航天器均采用小推力方式进行机动，在给定时间内速度的改变较小，判定机动能力的重要程度最弱；至此，至于剩下两个因素中哪一个更重要，给定追逃航天器的数目和状态信息，按照相对速度、相对距离、机动能力和相对距离、相对速度、机动能力两种重要性序列分别生成权重矩阵和匹配方案，对比追逃结果的优劣；三个因素的重要性由高到低排序为：相对速度、相对距离、机动能力，由此产生6种不同的符合条件的判断矩阵，根据应用结果选取最优的判断矩阵；步骤三：基于步骤二得到的权重矩阵，将多航天器追逃问题转换为带权二分图的最优匹配问题，为各追踪器匹配逃逸器；多对多航天器追逃博弈问题的分为等数量追逃、少追多逃和多追少逃问题三类情形进行匹配，直到找到相等子图的完全匹配为止；对于二分图的基本概念：若图中的顶点分为两个不相交的子集，且图中的每条边所关联的两个顶点分别属于不同的顶点集，则此图为一个二分图，在二分图中，一个匹配是一个边的集合，其中任意两条边都没有公共顶点，定义匹配点、匹配边、未匹配点、未匹配边等概念；二分图的最大匹配是指，在一个图的所有匹配中，所含匹配边数最多的匹配；可以通过寻找增广路径来增加匹配中的匹配边和匹配点，从而求解出最大匹配；在二分图中，从一个未匹配点出发，依次经过未匹配边、匹配边、未匹配边
……
的交替顺序形成的路径称为交替
路径；若交替路径最终到达了另一个未匹配...

【专利技术属性】
技术研发人员：张甜甜，祁瑞，
申请(专利权)人：北京理工大学，
类型：发明
国别省市：