异构集群系统分布式自适应输出时变编队跟踪控制

田 磊 董希旺 ,2 赵启伦 李清东 吕金虎 ,2 任 章 ,2

集群系统中个体之间通过局部的感知和简单的行为交互,可以在宏观上涌现出更为复杂与高级的群体行为,即集群智能.生物界中有很多集群智能的例子:鸟群利用上洗气流编队飞行,可以降低飞行阻力,实现远距离大范围迁徙;狼群通过协同围捕,利用时空优势,能够以弱胜强,捕食大型猎物;蚁群通过协同搬运,可以充分弥补个体力量薄弱的劣势,以小博大,提高觅食效率.军事领域也有很多通过集群协同方式提升作战效能的例子:多个微纳卫星通过协同可以构建一个虚拟卫星,实现深空探测及高精度三维成像等复杂功能,完成单个卫星无法完成的任务;多架低成本无人机通过协同方式可以实施蜂群作战模,使敌方高价值防空体系的火力通道陷入饱和,进而大幅度提高突防能力、毁伤能力和作战效费比.正是因为集群智能可以实现上述1+1≫2的效果,近年来备受学者关注.编队控制作为集群智能的重要基础性研究分支,也涌现出了较为丰富的研究成果.早期的编队控制方法主要是集中式,即存在中心节点,且对通信拓扑要求高.随着一致性理论的发展,分布式编队控制方法以其不依赖中心节点、具备良好的可扩展性等优点,取得了长足的发展[1].研究的方向从固定编队[2]扩展到时变编队[3−4],研究的内容从普通编队控制[5−6]发展到编队跟踪控制[7−8]、编队合围控制[9−10]等,研究的模型也逐渐丰富,包括了线性系统[11]、非线性系统[12−13]、同构系统[14]、异构系统[15−16]等,研究的通信拓扑也涵盖无向拓扑、有向拓扑、连续性拓扑、非连续性拓扑、拓扑固定、拓扑切换等[17−20].不仅如此,研究的编队稳定性问题也从渐近稳定问题逐渐发展到有限时间稳定问题甚至固定时间稳定问题[21−22].随着研究的深入,如何尽可能减少编队控制所需的全局信息成为分布式控制的研究重点,如文献[3]、[8]、[10]、[17]、[19]中设计控制协议参数时需要计算与通信拓扑相关的拉普拉斯矩阵的特征值,文献[16]、[17]、[19]、[21]中构造控制协议时需要用到领导者输入的上界值.此外,文献[2−22]中时变编队函数是需要设定的全局信息,也并非分布式实时可调控的.

基于上述分析可知,编队控制方法在理论上尚有许多可以研究和改进的地方.同时,由于理论研究通常受实际应用驱动,根据不同的任务场景,实际应用中期望集群系统的智能体模型既可以是同构的(无人机蜂群作战任务场景),也可以是异构的(空地/空海协同搜索任务场景),而考虑到强对抗博弈和强电磁干扰带来的影响,还期望通信拓扑是有向切换的,并且实际应用中阵风、紊流等未知外部扰动也是不可忽略的因素.需要指出的是,在研究编队跟踪问题时,跟随者跟踪的领导者通常分为合作目标和非合作目标.所谓合作目标是指领导者与跟随者之间存在协同合作关系,通常情况下,领导者为己方可控目标,其控制输入大多是已知可获取的.所谓非合作目标是指领导者与跟随者之间不存在协同合作关系,通常情况下,领导者为非己方不可控目标,其控制输入大多是未知不易获取的.经过文献调研发现,目前尚无研究成果能够解决高阶异构集群系统在领导者存在已知或未知控制输入、领导者和跟随者均存在未知扰动、有向通信拓扑存在切换等多种因素并存的情况下的输出时变编队跟踪控制问题.

受上述理论研究和工程应用价值的激励,本文研究了高阶异构集群系统的时变编队跟踪控制问题.主要创新有以下两点:1)集群系统中所有智能体均采用异构模型,同时研究了领导者存在已知和未知控制输入两种情况.需要说明的是,考虑到异构集群系统中智能体的动力学模型可能完全不同,当领导者的控制输入是已知时,可以在领导者层和跟随者层构造对自身状态的观测器,基于该观测器构造的控制协议不需要利用领导者的全部状态,只需要利用领导者的输出信息.而在领导者的控制输入是未知时,无法通过领导者的输出信息得到其全部状态信息,进而需要某个跟随者能够获取领导者的全部状态信息.2) 在上述领导者存在已知或未知控制输入的情况下,进一步考虑领导者和跟随者存在未知扰动以及有向通信拓扑存在切换等多种因素并存对集群系统的影响,结合观测器理论、自适应理论和滑模控制理论,构造了完全分布式的控制协议,所有跟随者均不需要获取领导者输入的上界值,也不需要利用与通信拓扑相关的拉普拉斯矩阵的特征值,此外,协议中还针对时变编队函数设计了分布式观测器,使其不再是全局信息,成为分布式实时可调控的.由于本文研究的集群系统均采用有向通信拓扑,因此极大降低了集群系统对通信带宽需求.综上所述,本文研究的内容兼具理论研究与工程应用价值.考虑到有人−无人机协同飞行(领导者为合作目标,见图1)和空地协同围捕打击(领导者为非合作目标,见图2)等任务场景下,集群系统中个体间动力学模型差异较大即异构模型,同时强电磁干扰时的通信限制和强对抗博弈时外界环境的不确定干扰为集群系统的协同控制带来了一定的困难,因此本文研究的方法能够为实战场景下的分布式 协同控制提供有力的理论支撑和技术保障.

图1 有人−无人飞行器协同飞行示意图Fig.1 The illustration of cooperative flight strategy for manned and unmanned aerial vehicles

图2 空地协同围捕目标示意图Fig.2 The illustration of hunting target by air-ground collaboration

1 预备知识和问题描述

1.1 图论知识

集群系统的通信拓扑关系对应的数学基础即为图论知识中关于拉普拉斯矩阵的相关知识.假设集群系统中含有q个智能体,智能体g和智能体k之间存在着从k到g的有向通信,则wgk=1,否则wgk=0.定义邻接矩阵为W=[wgk]∈Rq×q,定义度矩阵为D=,g=1,2,···,q},定义拉普拉斯矩阵为L=D−W.如果在集群系统的通信拓扑中,某个智能体r到其他所有智能体都存在至少一条有向连接,则称集群系统的通信拓扑中存在一个以 智能体r为根节点的生成树.

1.2 问题描述

假设集群系统中含有 1 个领导者和N个跟随者,不失一般性,所有智能体采用如下异构模型:

通过上述异构模型可知,集群系统中所有智能体的动力学模型可以完全不同.为实现集群系统的输出编队,假定所有智能体的输出具有相同的维度,即pi=p(i=1,2,···,1+N).同时为方便后续表达,领导者标号为1,跟随者标号为2,3,···,1+N.

假设 1[16].集群系统的通信拓扑包含一个以领导者为根节点的生成树,且领导者与跟随者之间的通信方式是单向的,方向是从领导者到跟随者,不同跟随者之间的通信方式可以是单向的,也可以是双向的.

定义集群系统的通信拓扑在固定时间间隔内保持不变,该时间间隔被称为驻留时间,记为α,其最小值被称为最小驻留时间,记为αmin.定义t=t0=0为初始时刻,t=tβ(β=0,1,2,···) 为通信拓扑发生变化的时刻,称为切换时刻.如果假设1 成立,则系统在 [tβ,tβ+1)内的拉普拉斯矩阵Lβ有如下形式:

为了前后表达方式的连贯性,后续带有上标β的变量都表示其位于由区间 [tβ,tβ+1) 定义的时间间隔内.

假设 2.矩阵方程组(3)存在对解矩阵 (Rg,Sg),同时矩阵方程(4)存在解矩阵Tg.

其中,g=2,3,···,1+N.

注1.假设2 要求方程(3)和(4)有解,其可解性是异构集群系统实现输出编队跟踪的必要条件,其物理意义是领导者的已知或未知输入能够通过跟随者的控制输入ug(t) 进行补偿.该假设与文献[16]中的假设3 和假设4 类似,由参考文献[16]中的注释3 可知,式(3)的可解性是异构集群系统能够实现输出调节的必要条件[23−27],式(4)的可解性类似采用鲁棒控制方法时匹配条件[28−29],因此假设2 是合理的.

假设 3.所有智能体的控制输入和未知扰动是有界的,即满足‖ui‖1≤χ,‖di‖1≤,χ和为未知正实数,同时未知扰动满足limt→∞di(t)存在.

注2.实际物理系统中智能体的控制输入和未知扰动是有界的,但是其界限值通常不易获取,假设3 要求扰动的极限存在是由于下文中设计控制协议时借鉴了文献[30]中扰动观测器的设计方法,此假设为文献[30]中扰动观测器能够收敛的前提条件.需要说明的是,针对不同于假设2 的扰动形式,下文中的注释4 给出了相应改进方案,因此假设3是合理的.

定义 1.如果通信拓扑如假设1 所述,且对于任意有界初始状态,等式(5)成立,则称异构集群系统能够实现有向拓扑切换条件下的输出时变编队跟踪.

其中,g=2,3,···,1+N,zg(h(t)) 为由h(t) 定义的时变向量,可以表征跟随者之间输出编队的队形以及整个编队队形与领导者之间的相对关系.

由定义1 可知,集群系统的输出时变编队跟踪队形刻画为以h(t)为自变量zg(h(t)) 为因变量的复合向量函数,与现有研究时变编队的文献[2−22]不同,h(t) 不是固化的全局信息,而是由领导者或与领导者存在直接拓扑连接的跟随者调控的时变向量.后续将构造分布式观测器对其进行估计,估计值记为.显然,如果,则.假设输出是1 维的,根据式(5) 所描述的等式关系,图3 给出了二维XOY平面内集群系统输出时变编队跟踪示意图,其中.在图1中,δg(h(t))表示跟随者g与领导者1 在二维XOY平面内相对距离,分别表示跟随者g与领导者1 在OX轴和OY轴方向上的相对距离,h(t) 为调整上述相对距离的时变信息,例如当跟随者绕领导者做旋转运动时,h(t) 可表示旋转半径和旋转角速度信息(见数值仿真部分).

图3 集群系统输出时变编队跟踪示意图Fig.3 The illustration of time-varying output formation tracking for swarm systems

引理 1[17].对于任意正定实对称矩阵J1,J2∈Rn×n和任意向量a∈Rn,以下不等式关系成立:

其中,λmax(J)表示矩阵J的最大特征值,下同.

引理 2[31].如果假设1 成立,式(2)中分块矩阵的所有特征值均具有正实部.同时,存在正定对角矩阵使Eβ >0,记表示的最小特征值.

2 控制协议的设计

本节将给出控制协议的设计方法,考虑到集群系统中领导者存在已知或未知控制输入,领导者和跟随者存在未知扰动,有向通信拓扑存在切换等多种因素影响,且领导者的控制输入在已知和未知两种情况下控制协议的设计并不相同,因此对这两种情况进行分开讨论.

2.1 领导者存在已知输入时控制协议的设计

当领导者的控制输入是已知可设计时,通常对应领导者为合作目标的情况,则对领导者设计如下控制协议:

其中,式(7a)和(7b)为具有未知扰动时领导者对自身状态的观测器,式(7c)为基于该观测器设计的控制器,表示领导者对自身状态x1(t) 的估计值,表示领导者对其未知扰动d1(t) 的估计值,表示根据任务需求和领导者模型设计的控制输入,以实现对领导者输出轨迹的实时调控,H1、由下文中的算法1 确定.对跟随者g ∈{2,3,···,1+N}设计如下控制协议:

算法 1.通过下列步骤确定控制协议(7)和(8)中的待定参数:

步骤1.设计Rg、Sg和Ug使等式(3)和(4)成立.

步骤2.选取矩阵Hi和H¯i使得矩阵A¯i是Hurwitz 的,其中i=1,2,···,1+N.

步骤3.选取θ >0,求下列Riccati 方程的正定实对称解阵P.

其中,In1为n1阶单位阵.由于 (A1,In1) 是可控的,所以存在唯一的解使得等式(10)成立[16].

步骤4.选取ϑg >0,ρg >0,ϕg >0,µ>0,σg >0,选取使得Ag+是Hurwitz 的,计算

步骤5.求下列Lyapunov 方程的正定实对称解阵Qg.

其中,Ing为ng阶单位阵.

步骤6.选取和满足,且为非奇异矩阵.

证明.首先证明观测器的收敛性.有关自身状态观测器(7a)和(7b)以及式(8a)和(8b)的收敛性证明过程详见文献[30],本文不再赘述.根据文献[30]中的定理,定义dg(t),可知.现证明分布式观测器(8c)～(8i)以及式(8j)和(8k)的收敛性.

构造Lyapunov 函数如下:

对V(t) 求导可得

令ζ(t)=V2(t)+V3(t)+V4(t),根据算法1 中的步骤3,可将式(25)进一步化为

更进一步由式(26)可得

因为系统通信拓扑在t=tβ发生变化,如果β ≥1,则由引理1 可知

结合式(27)和(29)可得

同理根据式(26)可以推出

进一步根据递归方法可得

其中

讨论以下两种情况:

由t−t0=t−tβ+tβ−tβ−1+···+t1−t0≥t−tβ+βαmin且 lnς1−θαmin<0可得limt→∞V(t)≤ζ(t),矛盾,所以上述假设不成立.

进一步可得

接下来证明在控制器(8m)～(8q)的作用下异构集群系统能够实现输出时变编队跟踪.将控制器(8m)～(8q)代入式(1)中可得

由算法1 中的步骤6 和等式(12)可知

进而式(44)可化简为

考虑如下Lyapunov 函数:

进而可得

2.2 领导者存在未知输入时控制协议的设计

当领导者的控制输入是未知不可设计时,通常对应领导者是非合作目标的情况,则只对跟随者g ∈{2,3,···,1+N}设计如下控制协议:

其中,各符号变量表示的含义保持不变.与控制协议(8)对比分析可知,式(54d)中利用的是领导者的全部状态信息,而式(8d)中利用的是领导者根据输出信息对其自身状态信息的估计值,造成上述协议不同的原因是当领导者的控制输入是未知不可设计时,无法通过领导者的输出信息估计领导者的全部状态信息,因此需要某个跟随者能够获取领导者的全部状态信息.相应的,式(54c)相比于式(8c),去掉对领导者自身状态观测器估计误差的滑模补偿项.需要说明的是,此时h(t) 为由某个能够获取领导者全部状态信息的跟随者调控的时变向量,控制协议(54)中的其他待定参数仍由算法1 确定.值得注意的是,当领导者的控制输入是未知不可设计时,根据式(1)可知,此时领导者的未知输入和未知扰动对跟随者来说可整体看作领导者的未知输入,因此可以得到下列定理.

定理2.对于任意给定的有界初始状态,如果集群系统中的跟随者采用控制协议(54),通信拓扑的最小驻留时间满足 lnς1−θαmin<0 和lnς2−2µαmin<0,同时等式(12)成立,则异构集群系统能够在有向拓扑切换条件下实现对存在未知输入领导者的输出时变编队跟踪.

由于控制协议(8)和(54)的主要不同体现在跟随者对领导者状态信息的观测方法上,因此定理2 的证明过程类似定理1,此处略.由符号函数sgn(·)可知,控制协议(8)和(54)采用滑模控制方法来消除领导者的控制输入、领导者的未知扰动以及跟随者自身的未知扰动对集群系统的影响,因此抖振现象必然存在.现对控制协议(8)和(54)进行改进,参考文献[16]提供的方法,分别用下列连续向量函数、cg(t) 代替控制协议(8) 和(54)中的、fg(t),进而使控制输入产生较为平滑的连续性变化,能够很好地抑制滑模的抖振现象.

推论1.如果在控制协议(8) 中采用代替控制、fg(t),则异构集群系统对存在已知输入领导者的输出时变编队跟踪误差(定义为‖yg(t)−zg(h(t))−y1(t)‖1)是一致最终有界的.

推论2.如果在控制协议(54) 中采用代替控制,则异构集群系统对存在未知输入领导者的输出时变编队跟踪误差(定义为‖yg(t)−zg(h(t))−y1(t)‖1)是一致最终有界的.

考虑到篇幅原因,本文略去推论1 和推论2 的证明过程,感兴趣的读者可结合本文定理1 以及文献[16]中定理2 进行推导证明.由推论1 和推论2可知 limt→∞‖yg(t)−zg(h(t))−y1(t)‖1≤Θg,Θg与、∂g的选取有关.通过选取恰当的∂g可使‖yg(t)−zg(h(t))−y1(t)‖1收敛到0 附近较小的有界区间内,以满足实际应用需求.

由于控制协议(8)和(54)中均采用滑模控制方法,分析自适应律(8g)～(8i)、(8p)、(54f)、(54g)、(54m)可知,自适应增益φg(t) 、ηg(t) 、υg(t) 和κg(t)会一直缓慢地变大,因此可对自适应律(8g)～(8i)、(8p)、(54f)、(54g)、(54m)增加如下约束:

注4.由式(1)可知,本文研究的集群系统模型考虑比较简单的匹配扰动,同时要求未知扰动满足假设3,因此设计的自身状态观测器(7a)、(7b)、(8a)、(8b)以及(54a)、(54b)借鉴文献[30]的研究成果,其结构相对简单.如果模型考虑不匹配扰动或者观测器匹配条件不满足等更加复杂的条件时,对于未知扰动的假设条件也会发生相应的变化,可根据扰动形式的不同参考文献[32−33]设计更加复杂的全维或降维观测器.需要指出的是,当集群系统模型中扰动形式更加复杂时,除了对观测器(7a)、(7b)、(8a)、(8b)以及(54a)、(54b)进行改进外,分布式观测器(8c)～(8i)、(54c)～(54g)以及控制器(8m)～(8q)、(54j)～(54n)也要进行一定的改进,例如模型考虑不匹配扰动,即式(1)写为

则分布式观测器(8c)～(8i)、(54c)～(54g)以及控制器(8m)～(8q)、(54j)～(54n)中需要增加针对不匹配扰动Didi(t) 的滑模补偿项,此处不作详细论述,后续将进一步开展相关研究.

注5.本文提供的控制方法具有较强的普适性.首先本文研究的系统模型为异构模型,如果系统模型为同构模型,所有智能体具有相同的动力学模型,满足Ai=A ∈Rn×n、Bi=B ∈Rn×m、Ci=C ∈Rp×n、可以验证假设2 能够成立,此时Rg=In、Ug=Im、Sg=0m×n,因此文献[2−4],[6],[8],[11−14],[17],[20−21]研究的同构系统模型可以看作异构系统模型的特例.其次,本文研究的集群系统通信拓扑考虑有向拓扑切换的情况,而文献[10−12]中考虑的无向拓扑固定、文献[17]中考虑的无向拓扑切换以及文献[4−9],[13],[15−16],[18−19],[21−22]中采用有向拓扑固定等情况均可以看作有向拓扑切换的特例.再次,本文研究的集群系统考虑考虑领导者存在已知或未知控制输入,领导者和跟随者均存在未知扰动的情况,而文献[16],[19]中研究的领导者不存在输入或领导者和跟随者不存在扰动等情况都可以看作本文研究问题的特例.最后,本文对时变编队函数设计了分布式观测器,进而使时变编队函数不再是全局信息,文献[2−22]中的时变编队函数均可看作是本文中时变编队复合函数zg(h(t))的特例.

3 数值仿真

由于定理1 和定理2 的主要区别在于领导者的控制输入是已知可设计的还是未知不可设计的,对应控制协议(8)和控制协议(54)的观测器设计方法存在些许差异,在仿真层面上主要体现在领导者的运动轨迹是否可控,对于跟随者是否能够实现对领导者的输出时变编队跟踪没有明显影响,考虑到篇幅原因,在验证定理1 和定理2 时,对集群系统模型以及控制协议(8)和(54)的相同参数部分采用同样的选取方法,并对控制协议(8)和(54)中自适应律(8g)～(8i)、(8p)、(54f)～(54g)、(54m)增加式(58)中的约束.假定集群系统含有7 个智能体,其中包括1 个领导者(对合作目标和非合作目标两种情况分开讨论)和6 个跟随者,领导者标号为1,跟随者标号为2～7,所有智能体采用如式(1)所示的异构模型,设定如下:

仿真中仅考虑XOY平面运动,集群系统在OX轴和OY轴方向上的控制可视为两个独立方向的解耦控制,所有智能体在OX轴和OY轴方向上的初始状态在 (−8,8)中随机选取,未知扰动在(0,1)中随机选取,当领导者是合作目标时,领导者采用控制协议(7),其中=−0.1 cos(t),跟随者采用控制协议(8),当领导者是非合作目标时,领导者的未知控制输入假定为−2 cos(t),跟随者采用控制协议(54),根据算法1以及式(58),控制协议(7)、(8)和(54)中的其他参数确定如下:

表征输出时变编队跟踪队形的复合函数及相关参数设计如下:

通过分析上述复合函数可知,仿真中期望跟随者在以领导者为圆心、为半径的圆弧上绕领导者运动,同时跟随者之间保持固定的相对距离,领导者或与领导者直接通信的跟随者能够对圆弧半径以及跟随者绕领导者的转动角速度的实时调控,当t∈[0,22) 时,圆弧半径为8,转动角速度为1 rad/s,当t∈[22,+∞) 时,圆弧半径增大为15,转动角速度增加为2 rad/s.图4 给出了集群系统的通信拓扑关系,设定通信拓扑在驻留时间α=4s保持不变,到达切换时刻后在G1～G4之间切换.接下来对领导者是合作目标和非合作目标两种情况的图例进行分别说明.

图4 集群系统的通信拓扑关系Fig.4 The communication topologies of swarm systems

1)领导者是合作目标

图5 给出了控制协议(7)和(8)中的观测器的估计误差(以所有智能体对自身状态和领导者调控的时变向量的估计误差为例),定义所有智能体对自身状态估计误差的二范数其中i=1,2,···,7,定义领导者调控的时变向量估计误差的二范数其中g=2,3,···,7.从图5中可以看出,控制协议(7)和(8)中的观测器是收敛的.需要说明的是,图5中在t=22 s时出现非连续的变化,原因是由于在t=22 s时编队半径瞬间增大且编队转动角速度瞬间增加.图6 给出了跟随者控制协议(8)中自适应增益(以和为例) 的变化曲线,其中g=2,3,···,7,由于对上述自适应增益增加了(58)中的约束,同时在t=22 s 时编队半径瞬间增大且编队转动角速度瞬间增加,因此图6 中的自适应增益曲线存在部分不可微点.图7 给出了跟随者采用控制协议(8) 时在OX轴和OY轴方向上控制输入的变化曲线,需要说明的是,控制输入在t=22 s 时发生非连续变化,也是由于该时刻编队半径瞬间增大且编队转动角速度瞬间增加造成的.图8 给出了集群系统在35s 内的输出轨迹,并给出了部分时刻领导者和跟随者的输出分布.图9给出了集群系统XOY平面内的输出时变编队跟踪误差,定义如下:

图5 控制协议(7)和(8)中观测器的估计误差Fig.5 The estimation errors of observers in control protocol (7) and control protocol (8)

图6 跟随者控制协议(8)中的自适应增益变化曲线Fig.6 The varying curve of adaptive gains in followers′ control protocol (8)

图7 跟随者控制协议(8)中的控制输入变化曲线Fig.7 The varying curve of control input in followers′ control protocol (8)

图8 当领导者是合作目标时集群系统的输出轨迹以及部分时刻的输出信息Fig.8 The trajectory of output and the output information at partial moments associated with swarm systems when leader is cooperative target

其中,g=2,3,···,7.从图9 中可以看出输出时变编队跟踪误差在编队函数连续变化时快速收敛到0,在非连续变化瞬间会短暂增大(对应t=22 s 时编队半径瞬间增大且编队转动角速度增加),随即快速收敛到0.

图9 当领导者是合作目标时集群系统输出时变编队跟踪误差Fig.9 The error of time-varying output formation tracking related to swarm systems when leader is cooperative target

2)领导者是非合作目标

考虑到领导者是合作目标和非合作目标两种情况下,仿真中集群系统参数选取的差异主要体现在领导者控制输入的不同,为更直观地进行仿真对比,当领导者是非合作目标时,仿真结果仅展示集群系统的输出分布和编队跟踪误差.图10 给出了集群系统在35s 内的输出轨迹,并给出了部分时刻领导者和跟随者的输出分布.与图8 中领导者相对规律的轨迹相比,图10 模拟了非合作领导者进行恶意机动的场景.图11 给出了集群系统XOY平面内的输出时变编队跟踪误差 Υg.与图9 类似,从图11中可以看出输出时变编队跟踪误差在编队函数连续变化时快速收敛到0,在非连续变化瞬间会短暂增大(对应t=22 s 时编队半径瞬间增大且编队转动角速度增加),随即快速收敛到0.因此定理1 和定理 2 通过数值仿真得到了验证.

图10 当领导者是非合作目标时集群系统的输出轨迹以及部分时刻的输出信息Fig.10 The trajectory of output and the output information at partial moments associated with swarm systems when leader is noncooperative target

图11 当领导者是非合作目标时集群系统输出时变编队跟踪误差Fig.11 The error of time-varying output formation tracking related to swarm systems when leader is noncooperative target

4 结论

本文研究了高阶异构集群系统输出时变编队跟踪控制问题.基于观测器理论、自适应控制理论、滑模控制理论等构造了分布式的控制协议,该控制协议分为观测器和控制器两个部分,观测器包括对自身状态的观测器和跟随者对领导者状态的分布式观测器,还包括针对时变编队函数的分布式观测器,基于上述观测器构造了控制器,摆脱了对领导者的控制输入上界值、与通信拓扑相关的拉普拉斯矩阵特征值以及时变编队函数等全局信息的依赖,使本文的控制方法成为完全分布式的控制方法,因此本文提供的控制方法具有较高的理论研究和实际应用意义.未来会进一步将相关理论成果在无人机和无人车等实际物理系统上进行试验验证.