摘要
为了显著提升航天器编队协同控制的响应速度,满足严格时间窗口的编队任务要求,设计一种新型预定时间保性能控制框架。该框架旨在同时应对编队系统中普遍存在的 3 个关键问题:航天器感知范围受限、执行机构饱和约束及编队成员间的碰撞规避。首先,结合误差转换技术与滑模控制方法,构建了控制李雅普诺夫函数条件,不仅能保障系统在预定时间内实现编队目标,同时保证了预设的瞬态响应性能和稳态精度。其次,通过引入高阶控制障碍函数,实现了对相邻航天器间相对距离的精确管控,在维持编队通信拓扑连通性的同时,有效规避了碰撞风险。进一步,采用二次规划方法求解最优控制量, 在考虑执行机构饱和约束的前提下,实现了预定时间收敛、拓扑连通保持和碰撞规避的多目标协同优化。最后,通过数值仿真对所提控制框架的有效性进行了系统验证,结果表明,该方案在满足多种约束条件下具有良好的控制性能和可靠性。
Abstract
To significantly enhance the response speed of spacecraft formation coordination and meet strict timewindow requirements for formation missions, this paper develops a novel prescribed-time performance-guaranteed control framework. This framework effectively addresses three critical challenges prevalent in formation systems: limited perception range, actuator saturation constraints, and inter-agent collision avoidance requirements. First, by integrating error transformation technique with sliding mode control, a control Lyapunov function condition is constructed. This design not only ensures the system meets strict timing requirements for formation tasks, but also guarantees prescribed transient response characteristics and steady-state performance metrics. Second, through the establishment of high-order control barrier functions, precise regulation of relative distances between adjacent spacecraft is achieved, maintaining formation communication topology connectivity while effectively preventing collision risks. Furthermore, this study employs quadratic programming to solve for optimal control inputs, realizing multi-objective coordinated optimization of prescribed-time convergence, topology connectivity maintenance, and collision avoidance control under actuator saturation constraints. To validate the effectiveness of the proposed control framework, systematic performance verification is conducted through numerical simulations. The simulation results fully demonstrate the reliability and superiority of the proposed control scheme in satisfying all specified constraints.
在过去几十年间,航天器编队飞行( spacecraft formation flying,SFF)系统凭借其显著的成本效益和任务执行灵活性,已成为空间科学和太空探索领域的关键技术范式,在学术界和工程界引发了持续的研究热潮。该系统通过保持特定几何构型来完成既定任务,因而协同控制技术对编队飞行任务具有决定性作用。近年来,研究人员已采用多种控制策略(如自适应控制、滑模控制、鲁棒控制和最优控制等)对航天器协同控制问题展开了深入探索[1]。
在深空探测、对地观测等关键任务中,多航天器编队系统常面临严格的时间窗口约束,对控制系统的收敛时间性能提出了极高要求。传统的渐近稳定控制方法虽然能保证系统最终收敛,却难以胜任具有精确时序要求的任务场景。为此,学者们将有限时间控制理论引入多航天器编队控制领域[2-3]。这类方法借助终端吸引子等机制,显著提升了系统的收敛速度和鲁棒性,但其收敛时间仍受初始状态的影响,在突发工况下可能出现响应延迟。为了提升系统的时间确定性,学者们开展了基于固定时间的多航天器编队飞行控制的研究[4-7],通过构造特殊形式的李雅普诺夫函数,获得了与初始状态无关的收敛时间上界。然而,该方法确定的收敛时间通常依赖控制器参数进行估计,存在较大保守性,且无法实现用户自主设定。针对这一局限,预定时间控制方法应运而生,并取得重要进展。如文献[8-9]基于预定时间理论设计了滑模控制策略,有效解决了考虑输入饱和约束的多航天器编队相对运动控制和基于旋转矩阵的姿态协同控制问题; 文献[10-12] 则通过设计新型预定时间性能函数,分别实现了航天器姿态控制、姿态协同与姿轨协同系统的高精度时间约束,为满足复杂空间任务中的严格时间窗口要求提供了新的技术途径。
然而,上述研究大多忽略了 SFF 系统在机动过程中的避碰问题。确保整个任务期间的安全间距是 SFF 系统的核心需求,特别是在密集编队场景下,亟需开发具备实时避碰能力的控制方法。人工势函数(artificial potential function,APF)法因算法简便、计算效率高等优点,常被视为最具工程应用前景的解决方案之一。文献[13-14] 分别针对欠驱动和全驱动 SFF 系统,提出了基于势能函数的避碰协同控制策略。文献[15] 针对存在空间摄动和推进器故障的多航天器系统,开发了融合扰动观测器的协同控制方案,通过引入指数型人工势函数有效解决了航天器间碰撞规避问题。文献[16] 结合滑模技术攻克了多航天器的固定时间位置协同控制的技术难题。此外,模型预测控制( model predictive control,MPC)作为基于滚动优化的先进控制方法,在 SFF避碰控制中展现出独特优势。文献[17-18] 分别针对航天器的姿态动力学和线性化轨道动力学,设计了基于 MPC 的协同避碰控制策略,该方法虽然在处理多种约束方面具有显著优势,但由于需要在线求解优化问题,对星载计算机的实时计算能力有较高的要求。需要特别指出的是,受星载传感设备的物理限制,航天器的环境感知能力存在固有局限,为此,学术界在连通性保持方面进行了研究。文献 [19]结合 APF 和滑模技术,提出了两种分布式协同控制策略,实现了多障碍物环境下的编队系统的连通性保持和避碰。文献[20-21] 通过引入预设性能函数,分别开发了具备连通性保持和避碰功能的姿轨一体化控制策略和轨道协同控制算法,在确保编队安全性的同时,显著提升了系统的瞬态和稳态性能。然而,上述研究存在关键性局限:理论分析仅能保证系统渐近稳定,收敛时间无法预先确定,难以满足深空探测等任务严格的时间窗口要求。
基于上述研究背景,本文重点研究同时考虑通信连通性保持与避碰的航天器预定时间编队跟踪控制问题,构建控制李雅普诺夫函数条件,保证编队协同误差满足预设的收敛速度、瞬态性能及稳态精度。
1 问题描述
1.1 坐标系与动力学模型
考虑由 N 个跟随航天器构成的多航天器编队系统,定义两类核心坐标系:地心惯性坐标系(earthcentered inertial,ECI),表示为局部垂直局部水平坐标系(local vertical local horizontal,LVLH),表示为该坐标系的原点设置在领航航天器的质心上,坐标系的轴沿着径向指向,与地心至虚拟航天器质心的连线方向重合,轴沿着轨道平面的正法向,轴位于轨道平面内,与轴构成右手正交系,通常指向虚拟航天器运动方向或轨道切向。航天器编队系统的示意图见图1。其中虚拟航天器为跟随航天器提供参考位置和速度信息。为地心到虚拟领航航天器的距离向量,其中分别为真近点角、轨道半长轴和轨道偏心率。为第 i 个跟随航天器相对于虚拟领航者的位置向量,为相应的相对速度向量,则其相对运动动力学模型为
(1)
式中:为第 i 个航天器的质量; fi 为实际的控制输入; 为反对称矩阵; 为重力向量,ri =为第 i 个航天器与地心之间的距离,μ 为地球的重力常量; Di 为非线性项:
式中:真近点角 θ 相对于时间的一阶导数和二阶导数分别为其中为虚拟领航航天器的平均轨道角速度。
图1航天器编队飞行示意图
Fig.1Schematic diagram of the SFF system
1.2 控制目标
针对式(1)描述的 SFF 系统的非线性相对运动学模型,设计领航-跟随编队控制策略。该策略需实现以下 3 个核心目标:首先,系统能够在预设的时间内快速收敛并稳定于期望编队构型,同时跟踪参考轨迹; 其次,在整个瞬态过程中,须严格避免成员间的碰撞; 最后,还需确保编队成员间的通信网络始终保持连通。
1)预设定时间编队跟踪:其中为用户指定的收敛时间;
2)碰撞避免任务:
3)连通性保持任务:。
注释 1 航天器编队协同控制中,碰撞避免与障碍物规避属于不同层级的约束问题,二者在任务需求与理论方法上存在本质差异:碰撞避免是针对编队成员间的安全距离; 障碍物规避是针对航天器与环境障碍物的距离,该约束维度与障碍物运动特性密切相关。本研究聚焦多航天器间碰撞避免的分布式协同控制,通过高阶控制障碍函数与预定时间性能函数联合设计,确保编队系统的安全性。障碍物规避问题由于涉及动态环境感知与障碍物运动预测等复杂因素,将作为独立研究方向在后续工作中深入探讨。
1.3 图论
本文利用无向通信图来描述航天器间的信息交互关系,其中,和分别为边集合和顶点集合。边为第 i 个航天器与第 j 个航天器之间存在双向信息交互。邻接矩阵是一个非负矩阵,若则否则。定义 Laplacian 矩阵:对角线元素非对角线元素。通常假设不存在自循环,也即假设是 N 个跟随者和虚拟领航者之间的通信拓扑,其中领航者仅能发出信息,并不能从其他航天器那里获得信息。若第 i 个航天器能够获得领航者的信息,则 bi 为一个正实数,否则为零。那么相应的 Laplacian 矩阵定义为,其中,。
引理 1 若通信图是连通的,即从领航者到每个跟随航天器都存在信息传输路径,则矩阵是正定的。
2 控制器设计
2.1 预定时间稳定控制李雅普诺夫函数约束
为了确保系统的稳定性,基于反步法构建控制李雅普诺夫函数。记为第 i 个航天器相对于编队中心的期望位置向量。那么,为第 i 个航天器与第 j 个航天器间的期望位置偏差。在这种情况下,第 i 个航天器的跟踪误差为。全局编队跟踪误差为
(2)
定理 1 当编队误差 ei 收敛到原点附近的小邻域内时,系统即可实现期望的编队构型。
证明: 令,则式(2)等价于
(3)
式中:为维的全 1 向量。根据引理 1,矩阵是正定的,所以当编队误差收敛到原点附近的小邻域内时,系统即可实现期望的编队构型。
为了实现编队误差的预定时间收敛,针对第 k 个误差分量构建如下预定时间性能函数:

(4)
式中:k = 1,2,3; 为用户设计的系统收敛时间; 时间变量 t 是性能函数的自变量,用以度量系统启动后所经历的时间。常数决定了用户设定的稳态值。参数的设计需满足:
(5)
式中:是性能函数的初值。若所设计的控制策略能够保证跟踪误差的每个分量都严格地在下列边界内演化,则系统能够实现预设定的瞬态和稳态性能:
(6)
式中:表示编队误差的稳态误差的上边界; 和分别为编队误差的正超调量上边界、负超调量下边界。
为了便于控制器设计,引入误差变换将受约束的误差动力学系统等价转换为无约束系统:
(7)
式中:。记。那么,变量的时间导数可以表示为
(8)
式中:
(9)
结合动力学方程(1),变量ξi 的时间导数为
(10)
式中:。
设计虚拟期望速度为
(11)
那么 ξi 的动力学模型为
(12)
式中:。为了实现控制目标,引入滑模变量:
(13)
式中:为对角正定矩阵。那么转换后的系统动力学模型为
(14)
式中:
结合式(12)和式(13),可得:
(15)
显然,这是一个关于变量的稳定一阶微分方程,其输入项为。因此,若初始条件有界,则的有界性将保证和的有界性。当 t 趋于无穷大时,若趋于零,则向量和同样趋近于零。最重要的是,的有界性意味着跟踪误差的预设性能指标得以实现。
为了驱使滑模变量收敛到零,考虑 Lyapunov 函数:
(16)
对式(16)求时间的导数,并代入运动动力学模型(1),可得:
(17)
为了实现渐近稳定性,设计镇定控制器:
(18)
在控制输入的作用下,Vi 的时间导数为
(19)
根据李雅普诺夫稳定性定理,滑模变量 si 将在时间趋于无穷时收敛到零。因此,Vi 可作为多航天器编队系统的控制李雅普诺夫函数,保证系统预设时间稳定。
2.2 控制障碍函数约束
通过高阶控制障碍函数将避碰和连通性保持约束转换为控制输入约束,确保多航天器系统始终保持安全状态,即成员航天器间的距离被严格地限制在预定义的安全区间内。首先给出关于障碍控制函数的定义和引理。
定义 1 针对如下仿射控制系统:
(20)
式中:是局部 Lipschitz 的,U 是一个控制约束集合,定义集合
(21)
式中:是一个连续可微的函数。如果存在一个 K 类函数 α(·),使得对所有的都有下式成立:
(22)
式中:分别为函数沿着向量场 f 和 g 的李导数。那么函数被称为集合 C 上的控制障碍函数。
引理 2 对于定义 1 中给定的集合若连续可微函数是定义域 D 上的控制障碍函数,则满足的 Lipschitz 连续控制器将使得集合 C 对于模型(1)是前向不变的,其中:
对于任意相邻的航天器 i 和航天器 j,定义二者间的距离为为满足避碰和连通性保持约束,需满足:
(23)
式中:是避碰的最小容许距离,且满足
为了保证系统在编队飞行过程中满足避碰和连通性约束,需要设计零化障碍函数(ZBF)。
碰撞约束 ZBF:
(24)
若则成员航天器间将不会发生碰撞。
连通性保持约束 ZBF:
(25)
若则相邻航天器间的通信连接将不会发生割裂。
接下来,需要设计合适的控制策略,使得由碰撞约束和连通性约束共同定义的安全集合对于系统是前向不变的。由于安全约束相对于航天器动力学系统是二阶可微的,需引入高阶障碍函数来保证上述集合的前向不变性。
对于避碰约束,引入如下函数:
(26)
相应的集合记为:
(27)
基于文献[18] 的分析,若存在常数和使得所有的都有则集合是前向不变的。结合公式(1),不等式等价于:
(28)
其中,
对于连通性保持约束,引入如下函数:
(29)
相应的集合为:
(30)
遵循相同的步骤,是保证集合前向不变的必要条件,即连通性保持条件为:
(31)
其中,变量的表达式同避碰约束。
2.3 基于二次规划的安全控制器设计
基于二次规划框架,将控制障碍函数(CBF)条件表示的安全性约束和控制李雅普诺夫函数(CLF)条件表示的稳定性约束结合起来,通过求解二次规划问题,可以得到满足多目标控制要求的控制策略:

(32)
式中:松弛变量
。
定理 2 若集合非空,则对所有满足的参数,二次规划问题的解存在且唯一。
证明 选择任意且定义:
(33)
(34)
(35)
式中,满足式(32)中的约束,因此若通过调整参数和使得则二次规划问题的解存在。由于式(30)是一个凸的二次规划问题,且可行集合是非空的,因而解是唯一的。
注释 2 如定理 2 所述,二次规划问题的可解性,取决于执行器输出上界能否满足由系统动态所决定的控制需求。当在控制输入边界范围内无法找到满足以下不等式的控制变量时,二次规划问题的可行域可能退化为空集,导致控制器失效。尽管定理 2 从理论上证明了可行解存在的充分条件,但受限于上述关系的复杂性,尚未建立与预定时间、避碰距离及控制参数之间的显式定量关系。从定性上分析,碰撞规避是硬约束,若在控制输入边界内无法找到可行解,可以通过减小的值、减小松弛因子的权重或增加预设定的收敛时间来放松对控制稳定性的约束,以便得到可行解。
(36)
定理 3 若将二次规划问题的解记作,且对于所有其中均有那么控制律可确保编队跟踪误差在预设定时间内从任意初始状态收敛到平衡点。
证明根据李雅普诺夫稳定性理论,若变量则系统(14)具有一致毕竟有界性。这一性质表明,系统变量的有界性可直接推导出变量 ξi(t)的有界性。结合 2.1 小节的理论分析,当变量保持有界时,编队误差将严格遵循预设性能边界演化,并能在预先设定的时间内收敛至稳态域。另外,根据 2.2 小节的理论分析,二次规划问题中控制障碍函数约束将保证集合 C(t)是前向不变的,因而在最优控制律的作用下,系统状态将始终在安全集合 C(t)内演化。综上所述,求解二次规划问题得到的控制策略可确保航天器编队系统在预定时间内实现期望编队构型,同时满足避碰、连通性保持以及执行机构饱和约束。
注释 3 虽然文中定义的编队跟踪误差 ei 是全局的,即对于任意初始跟踪误差 ei(0),只要执行器饱和约束可行(即二次规划问题有解),预定时间性能函数和误差转换技术可保证 ei(0)在时进入预设稳态区域,且不受初始误差幅值的影响。但是对于编队安全问题来讲,只有编队初始状态满足通过高阶控制障碍函数约束的前向不变性,系统才能始终保持安全。综上所述,整个编队系统是半全局稳定的。
3 仿真分析
对提出的领航-跟随编队跟踪控制策略的有效性进行仿真验证。假设编队系统由一个领航者和 4 个跟随者构成,通信拓扑图如图2所示,则其对应的邻接矩阵为
图2通信拓扑
Fig.2Communication topology
假设虚拟领航者在近地椭圆轨道上运行,轨道参数为:半长轴 ac = 7 178 km,偏心率 ec = 0. 01,初始真近点角 θ(0)= 0 rad,地球的引力参数 μ = 3.986 × 1014m3 / s2; 跟随航天器参数:质量 mi =100 kg,安全距离 hsafe = 20 m,最大的通信范围hcomm = 80 m; 4 个跟随航天器的初始位置为:
初始速度υi(0)= [0,0,0 ] Tm / s;
跟随航天器相对于领航者的位置为:
期望的跟踪轨迹预设性能函数的参数选取为:
控制器参数为:
图3为编队协同跟踪误差 ei 的3 个分量随时间变化的动态特性曲线。由图可知,所有的误差分量均在预设定时间 T = 60 s 内快速收敛到指定的稳态集合范围内; 误差演化轨迹严格处于预设性能包络边界内,无超调现象,满足瞬态指标; 收敛后误差稳定维持在 3 cm 以内,验证了控制系统的稳态精度。图4为 4 个跟随航天器的控制输入曲线,结果表明,所有控制输入均严格保持在执行机构的最大输出范围内,未出现饱和现象; 各航天器控制信号连续平滑,无高频抖振。图5为编队系统演化过程中成员航天器间的相对距离变化曲线。由图可知,在整个队形形成和跟踪过程中,任意两个航天器间的最小距离始终满足且距离曲线平滑变化,无突变现象,表明避碰控制策略具有良好的稳定性。另外,具有通信需求的相邻航天器间的距离始终维持在范围内,通信拓扑结构在整个任务过程中保持稳定连接,未出现通信中断情况。综上所述,所构建的协同控制策略能够保证编队系统在同时满足连通性保持、避碰和执行机构饱和限制的多重约束下,仍能够保证编队整体的跟踪性能。
图3编队协同跟踪误差的响应曲线
Fig.3Response curve of formation cooperative tracking error
图44 个跟随航天器的控制输入
Fig.4Control inputs of the four follower spacecraft
为了凸显所提控制策略的优势,选取文献[23] 的研究成果进行仿真对比。该文献结合传统性能函数、滑模控制技术及人工势场函数完成多航天器无碰撞协同跟踪控制任务,所选取的性能函数为
(37)
图5成员航天器间的距离
Fig.5Distance between the follower spacecraft
为了保证对比的公正性,将该性能函数中的参数 l 的值确定为 l = 0. 05,以保证编队系统在两种控制策略作用下的收敛时间大致相同。图6为文献[23]策略下的编队跟踪误差曲线。
图7为对应的控制输入曲线。由图6可知,虽然编队跟踪误差在 t = 60 s 左右收敛到平衡点附近,但是需要多次调整收敛速度。另外,编队跟踪误差的第二个分量 ei2在初始阶段超出了性能包络,会导致控制输入奇异(见图7( b))。由图7可知,其所需的控制输入幅值远大于本文所提出的控制策略,这是由于文献[23]在编队控制策略设计过程中未考虑控制输入饱和问题。综上,本文策略在多约束下的多航天器编队协同控制中展现出卓越的性能。
图74 个跟随航天器的控制输入
Fig.7Control inputs of the four follower spacecraft
4 结论
1)针对复杂约束下(碰撞避免、感知范围有限及执行机构饱和)的多航天器编队跟踪问题,本文提出了一种新型的预定时间保性能控制框架。该框架能在给定时间内确保系统稳定收敛,并可以满足各项性能约束。
2)构建的控制李雅普诺夫函数条件,保证了编队协同跟踪误差的预定时间收敛性及预设精度。同时,约束了瞬态响应特性。
3)通过构建一个二次规划框架,在执行机构饱和约束条件下实现了预定时间收敛、碰撞避免及连通性保持的多目标优化。
4)仿真结果验证了该方法的有效性:在 60 s 内实现误差收敛(稳态误差 <3 cm),控制输入无饱和,航天器间无碰撞且通信连通。
后续研究将进一步考虑空间碎片等动态障碍物的规避问题,提升编队系统在复杂空间环境下的适应性。

