一种航天器编队轨迹跟踪分布式控制方法及相关设备技术

本发明专利技术提供了一种航天器编队轨迹跟踪分布式控制方法及相关设备，包括：构建目标航天器编队每个航天器的相对运动PH方程；分别将每个航天器的相对运动PH方程进行转换，得到目标航天器的相对运动轨迹跟踪误差PH方程作为目标航天器编队的跟踪误差系统；基于跟踪误差系统，构造将目标航天器编队中各航天器的相对位置误差耦合到跟踪误差系统的期望哈密顿函数；根据期望哈密顿函数设计目标航天器编队的控制律并计算得到目标航天器编队的分布式协同控制律的值；根据分布式协同控制律的值对目标航天器编队进行轨迹跟踪分布式控制；航天器编队能够准确、快速地达到期望的运动轨迹并保持构型，实现了航天器集群网络化的分布式协同控制。制。制。

【技术实现步骤摘要】
一种航天器编队轨迹跟踪分布式控制方法及相关设备


[0001]本专利技术涉及航天器编队轨迹跟踪控制
，特别涉及一种航天器编队轨迹跟踪分布式控制方法及相关设备。

技术介绍

[0002]航天器编队是指由多个航天器根据指定任务需求按一定拓扑链接组成的航天器集群，通过编队内各成员航天器相互合作，从而实现一颗大“虚拟航天器”的功能。编队航天器在轨飞行时，为了实现不同的任务，需要对编队运行轨迹进行跟踪控制，从而达到预定构型。目前已有多种控制方法应用到航天器编队控制，如非线性自适应控制、滑模方法、线性二次调节、IDA
‑
PBC(Interconnection and Damping Assignment Passivity
‑
Based Control，基于图论和互联阻尼分配的无源控制)、脉冲控制、模型预测控制和快速搜索随机树等方法，编队协同方法主要有领导跟随法，虚拟结构法和蜂拥控制等方法。
[0003]近年来，采用能量思想的Port
‑
Hamiltonian(PH，端口哈密顿)系统理论逐步兴起并快速发展，为复杂非线性系统建模与控制提供了有潜力的求解方案。另外PH理论从能量的角度出发对系统实现控制，解决了李亚普洛夫函数的选取问题，用其建立的模型更加符合工程实际应用。在航天器编队控制领域，端口哈密顿系统方法得到了初步研究。Chang Liu等将圆形限制性三体问题为为端口哈密顿系统形式，基于能量整形和阻尼注入设计了控制策略，保证了渐近稳定性，并且能够任意设置平衡点。Ewoud Vos等建立了端口哈密顿系统为的赤道平面二维二体轨道动力学模型，利用广义正则变换方法获得相对目标轨道的误差系统PH模型，然后设计了基于虚拟弹簧阻尼的内部控制律，将每个卫星引导到目标轨道，以及基于虚拟弹簧阻尼交互的分布式相对运动控制律，使多个卫星在轨道上均匀分布。该文是多卫星端口哈密顿系统协同控制的首次尝试，针对的是特殊星座位置保持控制问题；但由于航天器编队轨迹跟踪是时变的，而时变端口哈密顿系统一般不满足无源性。
[0004]IDA
‑
PBC是一种针对非线性系统的控制律设计方法，自从Ortega.R等在2001年提出该方法以来，已经应用到非常多的领域，在航天器控制方面也初步得到了应用。如Najmeh Javanmardi等应用IDA
‑
PBC控制方法以及压缩分析设计了跟踪控制律。但是该文的方法仅是基于领导/追随者的一对一跟踪，跟踪控制中缺乏相互协调和反馈机制，如何实现集群网络化的分布式协同控制，该研究并未提供解决思路。

技术实现思路

[0005]本专利技术提供了一种航天器编队轨迹跟踪分布式控制方法及相关设备，其目的是为了实现航天器集群网络化的分布式协同控制，并减小各航天器间的相对距离误差。
[0006]为了达到上述目的，本专利技术提供了一种航天器编队轨迹跟踪分布式控制方法，包括：
[0007]步骤1，构建目标航天器编队中每个航天器的相对运动PH方程；相对运动PH方程用于描述每个航天器的相对运动轨迹；
[0008]步骤2，分别将每个航天器的相对运动PH方程进行转换，得到目标航天器编队的相对运动轨迹跟踪误差PH方程，相对运动轨迹跟踪误差PH方程作为所述目标航天器编队的跟踪误差系统；相对运动轨迹跟踪误差PH方程用于描述相对运动轨迹的跟踪误差；
[0009]步骤3，基于跟踪误差系统，构造将目标航天器编队中各航天器的相对位置误差耦合到跟踪误差系统的期望哈密顿函数；
[0010]步骤4，根据期望哈密顿函数设计目标航天器编队的控制律，并将目标航天器编队的实际状态和期望状态代入控制律进行计算，得到目标航天器编队的分布式协同控制律的值；
[0011]步骤5，根据分布式协同控制律的值对目标航天器编队进行轨迹跟踪分布式控制。
[0012]进一步来说，每个航天器的相对运动PH方程为：
[0013][0014][0015]其中，是x的导数，为航天器速度和加速度，J为航天器的互联矩阵，H(x,t)为能量函数，x为航天器的状态变量，t为时间，G无具体物理含义，u为航天器的输入，I为单位矩阵，y为航天器的输出。
[0016]进一步来说，目标航天器编队的相对运动轨迹跟踪误差PH方程为：
[0017][0018][0019]其中，为航天器速度误差和加速度误差，I
n
为n阶单位矩阵，n为航天器数量，为矩阵Kronecker积，为航天器互联矩阵，为由多个航天器跟踪误差系统组成的航天器编队跟踪误差系统的哈密顿函数，为编队跟踪误差系统状态，为系数矩阵，无具体物理含义，为编队误差系统输出，为系统控制律。
[0020]进一步来说，分布式拓扑相对位置误差的期望哈密顿函数为
[0021][0022]其中，n为航天器数量，a
ij
为邻接矩阵元素，a
ij
≠0，k
p
为互联系数，为航天器i与航天器j之间相对距离保持的误差，x
i
为第i个航天器的误差，为第j个航天器的误差。
[0023]进一步来说，步骤4包括：
[0024]根据由多个航天器跟踪误差系统组成的编队跟踪误差系统的哈密顿函数可得：
[0025][0026]其中，为目标航天器编队跟踪误差系统状态；
[0027]分布式拓扑相对位置误差的期望哈密顿函数为：
[0028][0029]根据若系统中期望哈密顿函数H
d
(x)、对称矩阵R
d
(x)≥0和反对称矩阵J
d
(x)满足偏微分方程的定理，求解匹配方程为：
[0030][0031]其中，L为拉普拉斯矩阵，A为邻接矩阵，J
12
为矩阵J的元素，
[0032]根据和的任意性，满足该方程的条件为J
11
＝R
11
，J
12
＝I3，可设
[0033][0034][0035]根据对于系统，若期望哈密顿函数H
d
(x)、对称矩阵R
d
(x)≥0和反对称矩阵J
d
(x)满足偏微分方程的定理，可知系统的控制律为：
[0036][0037]将
代入中，求得目标航天器编队的分布式协同控制律为：
[0038][0039]其中，k
d
为阻尼系数。
[0040]进一步来说，根据目标航天器编队的分布式协同控制律，可求得目标航天器编队的跟踪误差系统中第i个航天器的控制律为：
[0041][0042]根据跟踪误差系统中第i个航天器的控制律，可得目标航天器编队中第i个航天器的控制律为：
[0043][0044]其中，
[0045]本专利技术还提供了一种航天器编队轨迹跟踪分布式控制装置，包括：
[0046]构建模块，用于构建目标航天器编队中每个航天器的相对运动PH方程；相对本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种航天器编队轨迹跟踪分布式控制方法，其特征在于，包括：步骤1，构建目标航天器编队中每个航天器的相对运动PH方程；所述相对运动PH方程用于描述每个航天器的相对运动轨迹；步骤2，分别将每个航天器的相对运动PH方程进行转换，得到所述目标航天器编队的相对运动轨迹跟踪误差PH方程，所述相对运动轨迹跟踪误差PH方程作为所述目标航天器编队的跟踪误差系统；所述相对运动轨迹跟踪误差PH方程用于描述相对运动轨迹的跟踪误差；步骤3，基于所述跟踪误差系统，构造将所述目标航天器编队中各航天器的相对位置误差耦合到所述跟踪误差系统的期望哈密顿函数；步骤4，根据所述期望哈密顿函数设计所述目标航天器编队的控制律，并将目标航天器编队的实际状态和期望状态代入所述控制律进行计算，得到所述目标航天器编队的分布式协同控制律的值；步骤5，根据所述分布式协同控制律的值对所述目标航天器编队进行轨迹跟踪分布式控制。2.根据权利要求1所述的航天器编队轨迹跟踪分布式控制方法，其特征在于，每个所述航天器的相对运动PH方程为：航天器的相对运动PH方程为：其中，是x的导数，为航天器速度和加速度，J为航天器的互联矩阵，H(x,t)为能量函数，x为航天器的状态变量，t为时间，G无具体物理含义，u为航天器的输入，I为单位矩阵，y为航天器的输出。3.根据权利要求2所述的航天器编队轨迹跟踪分布式控制方法，其特征在于，所述目标航天器编队的相对运动轨迹跟踪误差PH方程为：航天器编队的相对运动轨迹跟踪误差PH方程为：其中，为航天器速度误差和加速度误差，I
n
为n阶单位矩阵，n为航天器数量，为矩阵Kronecker积，为航天器互联矩阵，为由多个航天器跟踪误差系统组成的航天器编队跟踪误差系统的哈密顿函数，为编队跟踪误差系统状态，为系数矩阵，无具体物理含义，
为编队误差系统输出，为系统控制律。4.根据权利要求3所述的航天器编队轨迹跟踪分布式控制方法，其特征在于，所述分布式拓扑相对位置误差的期望哈密顿函数为式拓扑相对位置误差的期望哈密顿函数为其中，n为航天器数量，a
ij
为邻接矩阵元素，a
ij
≠0，k
p
为互联系数，为航天器i与航天器j之间相对距离保持的误差，x
i
为第i个航天器的误差，为第j个航天器的误差。5.根据权利要求4所述的航天器编队轨迹跟踪分布式控制方法，其特征在于，所述步骤4包括：根据由多个航天器跟踪误差系统组成的航天器编队跟踪误差系统的哈密顿函数可得：其中，为编队误差系统状态；所述分布式拓扑相对位置误差的期望哈密顿函数为：根据若系统中期望哈密顿...

【专利技术属性】
技术研发人员：陈琪锋，刘俊，廖宇新，戴明哲，魏才盛，金华，
申请(专利权)人：中南大学，
类型：发明
国别省市：