CN112305918A - 超螺旋观测器下的多智能体系统滑模容错一致性控制算法 - Google Patents

超螺旋观测器下的多智能体系统滑模容错一致性控制算法 Download PDF

Info

Publication number
CN112305918A
CN112305918A CN202011257346.1A CN202011257346A CN112305918A CN 112305918 A CN112305918 A CN 112305918A CN 202011257346 A CN202011257346 A CN 202011257346A CN 112305918 A CN112305918 A CN 112305918A
Authority
CN
China
Prior art keywords
fault
agent
observer
actuator
consistency
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011257346.1A
Other languages
English (en)
Inventor
杨蒲
胡旭凯
张芷晴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN202011257346.1A priority Critical patent/CN112305918A/zh
Publication of CN112305918A publication Critical patent/CN112305918A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种超螺旋观测器下的多智能体系统滑模容错一致性控制算法。考虑二阶非线性多智能体系统速度信息不可测且发生执行器加性故障,首先提出了一种超螺旋二阶滑模故障观测器,可以有效地对智能体速度状态信息和执行器故障信息进行估计,针对初始误差问题和峰值问题,通过改进观测器结构并引入增益系数,使得观测器估计更加快速准确。然后基于积分滑模控制,提出了一种容错控制方法,根据智能体之间的相对状态信息定义了一致性误差变量,基于一致性误差变量设计滑模面。最后结合观测器估计值以及智能体之间的通信结构信息,设计了一种滑模容错一致性控制律。该控制律可以使得多智能体系统在受到执行器加性故障和固有非线性动态影响下,仍能在有限时间内实现一致性目标,从而有效的避免了多智能体系统因执行器故障而导致一致性任务无法完成的问题。本发明可用于一类含有执行器加性故障和固有非线性动态的多智能体系统的一致性控制实现问题。

Description

超螺旋观测器下的多智能体系统滑模容错一致性控制算法
技术领域
本发明涉及一种超螺旋观测器下的多智能体系统滑模容错一致性控制算法,属于多智能体系统的容错控制技术领域。
背景技术
近年来,随着人工智能技术的快速发展,研究人员在智能算法领域取得了许多创新进展和成果。然而,由于长时间的频繁使用,智能设备会不可避免地出现故障。如果没有预先设定补偿故障的措施,一旦故障发生就可能影响到整个任务的执行进度,严重时甚至会导致财产损失或人员伤亡。因此,容错控制算法应运而生,并成为了当代控制领域研究的焦点和热点。一般来说,根据控制系统利用冗余处理故障的方式不同,容错控制可以划分为被动容错控制和主动容错控制。被动容错控制利用预先设置的系统冗余来达到系统对于故障的不敏感作用,它借鉴了鲁棒控制的设计思想,不改变控制器参数或结构,也不需要故障诊断模块的帮助,只对特定故障进行抑制;主动容错控制一般在故障发生后通过在线调节控制器参数或结构来达到容错目的,它往往需要借助故障诊断所得到的故障信息,然后进行控制律的调整或者重构。
自多智能体系统的概念诞生以来,多智能体系统的分布式协同控制方法便因其与传统集中控制相比更好的鲁棒性、实用性和灵活度而越发受到关注,诸如多机器人系统和多飞行器系统等在现实中已经得到了广泛应用。而多智能体系统的一致性作为其中最基本也是最重要的问题之一,也吸引了大量研究人员投入到对其的研究之中。所谓一致性,是指多智能体系统在适当的一致性控制策略作用下,各个智能体能够通过彼此之间的信息交互使得其最终状态都趋向于同一个期望,它在多智能体系统的编队控制、集群控制等问题上都有所体现,可见多智能体系统一致性控制的研究具有重要的价值。
然而,对于多智能体系统这样一个复杂而庞大的研究对象,故障的发生就显得更为频繁而危险。以由四旋翼无人飞行器组成的多智能体系统为例,受到老化、磨损或突发状况影响,四旋翼飞行器可能会发生执行器故障、传感器故障、通信故障等多种类型的故障,这些故障不仅会作用在故障飞行器本身,还会对系统内的其它飞行器造成影响,进而妨碍到整个系统的一致性控制。因此,在对多智能体系统一致性控制问题的研究中引入合适的容错控制器以补偿故障带来的影响就显得十分重要。
值得注意的是,目前多数工作在研究中都假定多智能体的位置和速度等信息是可测的,而在实际工程系统中,考虑到检测工具的成本以及智能体的负载问题,多智能体系统中未必会配置有完备的传感器设备,从而导致部分状态不可测的情况发生。针对这一现实问题,现有的研究成果中往往会采用结合观测器的方法来对不可测状态进行观测和预估。多智能体系统研究领域中利用状态观测器的合理构建来设计效果良好的控制方案已经取得了一定进展,结合观测器设计容错控制算法以解决多智能体系统在状态信息不可测和执行器故障情况下的一致性问题,不仅在理论和实际上都具有可行性且容易实现。
发明内容
发明目的:针对上述研究背景,提出一种新型的超螺旋观测器下的多智能体系统滑模容错一致性控制算法。设计了一种能够准确估计出执行器加性故障和不可测速度信息的超螺旋二阶滑模观测器;为了更好地确保观测的快速性,设计了能够很好处理峰值问题的增益系数,并且能够使系统在初始误差过大的情况下实现快速收敛;结合智能体之间的相对状态信息和观测器估计值,设计了滑模容错一致性控制律,通过设计积分滑模面取代线性滑模面,避免了滑模趋近过程的失稳,使得多智能体系统在容错的基础上仍能在有限时间内实现一致性控制任务。
技术方案:一种新型的超螺旋观测器下的多智能体系统滑模容错一致性控制算法。针对系统中存在的执行器加性故障和不可测速度信息,设计了带有增益系数的超螺旋二阶滑模观测器,使得观测器和系统在初始值相差较大的情况下仍能够快速收敛并准确获得估计值,并且有效避免了峰值问题带来的影响;根据智能体之间的相对状态信息定义了一致性误差变量,以此为基础设计积分滑模面以提高控制系统的鲁棒性;在滑模容错一致性控制律中引入双曲正切项,从而遏制了抖振对系统的影响,使设计的控制器能够发挥更好的控制效果。
一种超螺旋观测器下的多智能体系统滑模容错一致性控制算法,包括如下具体步骤:
步骤1)确定本发明所适用的多智能体系统模型及其参数,包括如下步骤:
步骤1.1)确定领导智能体的状态空间模型如式(1)所示:
Figure BSA0000224257420000021
其中,x0(t)、v0(t)∈Rm分别表示领导者的位置和速度状态;u0(t)∈Rm为控制输入。
步骤1.2)确定带有执行器加性故障和固有非线性动态的跟随智能体的状态空间模型如式(2)所示:
Figure BSA0000224257420000031
其中,xi(t)、vi(t)∈Rm分别表示第i个智能体的位置和速度状态;ui(t)∈Rm为控制输入;fi(xi(t),vi(t),t)∈Rm表示第i个智能体的固有非线性动态;βi(t)∈Rm代表第i个智能体时变时产生的执行器加性故障,βi(t)=0表示无故障发生,βi(t)≠0表示第i个智能体的执行器发生了加性故障;
对于固有非线性动态fi(xi(t),vi(t),t)和执行器加性故障βi(t)以及其一阶微分,满足下面不等式(3)和(4):
Figure BSA0000224257420000032
其中
Figure BSA0000224257420000033
为非负常数;
Figure BSA0000224257420000034
其中
Figure BSA0000224257420000035
为非负常数;
步骤2)确定本发明所适用的多智能体系统的通信拓扑结构:
考虑由一个领导者和n个跟随者组成的多智能体系统的通信拓扑结构,领导者标记为0,跟随者标记为i=1,2,...,n;该结构可利用图
Figure BSA0000224257420000036
来对包括领导者和跟随者在内所有节点之间的信息交互进行描述,其中,
Figure BSA0000224257420000037
表示图中各节点的集合,
Figure BSA0000224257420000038
表示边集合;定义子图G={V,E}以表示跟随者之间的通信拓扑网络,其中,V={1,2,...,n}表示子图中各节点的集合,
Figure BSA0000224257420000039
表示边集合;定义Ni={j|(i,j)∈E,i≠j}表示与节点i相邻的各节点集合,A=(aij)为子图G的邻接矩阵,如果(i,j)∈E,i≠j,即表示节点i和j之间有直接的路径,那么aij>0;否则,aij=0;定义节点i的入度为
Figure BSA0000224257420000041
则子图G的入度矩阵可表示为D=diag{d1,...,dn}。由此可得到子图的拉普拉斯矩阵表示为L=[lij]n×n=D-A;
对于领导者和跟随者之间的邻接矩阵,可定义为B=diag(b1,b2,...,bn);当跟随者i能够直接获得来自领导者的信息时,bi>0;否则,bi=0;
步骤3)构造超螺旋二阶滑模故障观测器,包括如下步骤:
步骤3.1)首先为每个跟随智能体构造如式(5)所示的超螺旋观测器:
Figure BSA0000224257420000042
其中,
Figure BSA0000224257420000043
Figure BSA0000224257420000044
分别表示第i个智能体的速度估计和执行器加性故障估计,α为辅助变量,λ1和λ2为待定参数;设计增益系数ηi(t),可以保证系统在初始误差较大时也能快速收敛至终态,其表达式如式(6)所示:
Figure BSA0000224257420000045
其中,
Figure BSA0000224257420000046
和σi为大于0的常数;由于增益系数随时间而快速衰减,因此也避免了系统初始值与估计值不同而导致的峰值问题;
步骤3.2)定义速度和故障状态的估计误差如式(7)所示:
Figure BSA0000224257420000047
对式(7)求导,可以得到系统的误差动态如式(8)所示:
Figure BSA0000224257420000048
步骤3.3)在式(4)成立的条件下,对于状态空间模型(2)和对应设计的观测器(5),如果参数λ1和λ2满足
Figure BSA0000224257420000049
则估计误差(7)最终一致有界且系统的误差动态(8)是有限时间稳定的,即:
Figure BSA0000224257420000051
Figure BSA0000224257420000052
其中,T为有限收敛时间;V1(0)为选取的李雅普诺夫函数在t=0时的值,且有
Figure BSA0000224257420000053
Figure BSA0000224257420000054
且P和Q均为正定矩阵;λmin{·}和λmax{·}分别表示矩阵的最小和最大特征根;通过为λ1和λ2选取合适的值,可以将估计误差调整到极小的值;
步骤4)设计容错控制算法,包括如下步骤:
步骤4.1)根据第i个智能体获取的邻居信息,定义一致性跟踪位置误差变量exi(t)和速度误差变量evi(t),其形式如下:
Figure BSA0000224257420000055
定义
Figure BSA0000224257420000056
Figure BSA0000224257420000057
式(13)可改写为如式(14)所示的向量形式:
Figure BSA0000224257420000058
其中,Im为m维的单位矩阵;定义
Figure BSA0000224257420000059
Figure BSA0000224257420000061
可得一致性跟踪误差系统如式(15)所示:
Figure BSA0000224257420000062
步骤4.2)设计积分滑模面函数如式(16)所示
Figure BSA0000224257420000063
其中
Figure BSA0000224257420000064
si(t),i=1,...,n为滑模变量;k1和k2为正常数;
步骤4.3)为跟随者智能体设计如式(17)所示的容错一致性控制算法:
Figure BSA0000224257420000065
其中γ和ε为正常数;通过理论推导可以得到(1)和(2)组成的多智能体系统在容错控制律(17)下能够实现执行器故障情况下的有限时间一致性。
有益效果:针对存在执行器加性故障和固有非线性动态,且速度不可测的多智能体系统,设计了一种新型的超螺旋观测器下的多智能体系统滑模容错一致性控制算法。设计了超螺旋二阶滑模观测器对故障和速度信息进行快速准确的预估。根据估计值,结合智能体之间的交互信息以及设计的积分滑模面,设计了滑模容错一致性控制算法,可以实现多智能体系统在执行器加性故障和非线性动态影响下的有限时间一致性。总的来说,具有如下具体优点:
①通过设计超螺旋二阶滑模观测器,同时对系统中的执行器加性故障和不可测速度进行有限时间内的准确估计,通过改进观测器结构,进一步提高预估的快速性,从而能够在初始预估值和实际值误差较大的场合既实现故障和状态的精准估计又可以节省设备资源,提高硬件的冗余度。
②所设计的滑模容错一致性控制律引入了积分滑模面和双曲正切项,避免了趋近过程中的失稳并削弱了抖振问题,有利于保证系统的鲁棒性。
③利用本发明所设计容错控制算法,能够解决一类速度状态不可测的多智能体系统在受到执行器故障和固有非线性动态影响下的有限时间一致性实现问题。这表明该算法的设计不仅节省了状态检测设备的开销,同时还增强了整个系统的鲁棒性和容错能力。
本发明所提出的超螺旋观测器下的多智能体系统滑模容错一致性控制算法,具有一定的应用意义,易于实现,实时性好,准确性高,能够有效提高控制系统安全性且可操作性强,节省时间,效率更高,可广泛应用于一类速度状态不可测的多智能体系统在执行器故障和固有非线性动态下的有限时间一致性实现问题。
附图说明
图1是本发明方法的流程图;
图2是Quanser公司研制的实验装置Qball-X4四旋翼直升机以及它的结构简图;
图3是多Qball-X4四旋翼直升机系统的通信拓扑网络;
图4是第2架Qball-X4四旋翼直升机的故障以及它的上限平方估计值曲线;
图5是第3架Qball-X4四旋翼直升机的故障以及它的上限平方估计值曲线;
图6是第4架Qball-X4四旋翼直升机的故障以及它的上限平方估计值曲线;
图7是多Qball-X4四旋翼直升机系统的位置跟踪误差。
图8是多Qball-X4四旋翼直升机系统的速度跟踪误差。
具体实施方式
下面结合附图对本发明做更进一步的解释。
如图1所示,针对系统中存在的执行器加性故障和不可测速度信息,设计了带有增益系数的超螺旋二阶滑模观测器,使得观测器和系统在初始值相差较大的情况下仍能够快速收敛并准确获得估计值,并且有效避免了峰值问题带来的影响;根据智能体之间的相对状态信息定义了一致性误差变量,以此为基础设计积分滑模面以提高控制系统的鲁棒性;在滑模容错一致性控制律中引入双曲正切项,从而遏制了抖振对系统的影响,使设计的控制器能够发挥更好的控制效果。
一种超螺旋观测器下的多智能体系统滑模容错一致性控制算法,包括如下具体步骤:
步骤1)确定本发明所适用的多智能体系统模型及其参数,包括如下步骤:
步骤1.1)确定领导智能体的状态空间模型如式(1)所示:
Figure BSA0000224257420000071
其中,x0(t)、v0(t)∈Rm分别表示领导者的位置和速度状态;u0(t)∈Rm为控制输入。
步骤1.2)确定带有执行器加性故障和固有非线性动态的跟随智能体的状态空间模型如式(2)所示:
Figure BSA0000224257420000081
其中,xi(t)、vi(t)∈Rm分别表示第i个智能体的位置和速度状态;ui(t)∈Rm为控制输入;fi(xi(t),vi(t),t)∈Rm表示第i个智能体的固有非线性动态;βi(t)∈Rm代表第i个智能体时变时产生的执行器加性故障,βi(t)=0表示无故障发生,βi(t)≠0表示第i个智能体的执行器发生了加性故障;
对于固有非线性动态fi(xi(t),vi(t),t)和执行器加性故障βi(t)以及其一阶微分,满足下面不等式(3)和(4):
Figure BSA0000224257420000082
其中
Figure BSA0000224257420000083
为非负常数;
Figure BSA0000224257420000084
其中
Figure BSA0000224257420000085
为非负常数;
步骤2)确定本发明所适用的多智能体系统的通信拓扑结构:
考虑由一个领导者和n个跟随者组成的多智能体系统的通信拓扑结构,领导者标记为0,跟随者标记为i=1,2,...,n;该结构可利用图
Figure BSA0000224257420000086
来对包括领导者和跟随者在内所有节点之间的信息交互进行描述,其中,
Figure BSA0000224257420000087
表示图中各节点的集合,
Figure BSA0000224257420000088
表示边集合;定义子图G={V,E}以表示跟随者之间的通信拓扑网络,其中,V={1,2,...,n}表示子图中各节点的集合,
Figure BSA0000224257420000089
表示边集合;定义Ni={j|(i,j)∈E,i≠j}表示与节点i相邻的各节点集合,A=(aij)为子图G的邻接矩阵,如果(i,j)∈E,i≠j,即表示节点i和j之间有直接的路径,那么aij>0;否则,aij=0;定义节点i的入度为
Figure BSA0000224257420000091
则子图G的入度矩阵可表示为D=diag{d1,...,dn}。由此可得到子图的拉普拉斯矩阵表示为L=[lij]n×n=D-A;
对于领导者和跟随者之间的邻接矩阵,可定义为B=diag(b1,b2,...,bn);当跟随者i能够直接获得来自领导者的信息时,bi>0;否则,bi=0;
步骤3)构造超螺旋二阶滑模故障观测器,包括如下步骤:
步骤3.1)首先为每个跟随智能体构造如式(5)所示的超螺旋观测器:
Figure BSA0000224257420000092
其中,
Figure BSA0000224257420000093
Figure BSA0000224257420000094
分别表示第i个智能体的速度估计和执行器加性故障估计,α为辅助变量,λ1和λ2为待定参数;设计增益系数ηi(t),可以保证系统在初始误差较大时也能快速收敛至终态,其表达式如式(6)所示:
Figure BSA0000224257420000095
其中,
Figure BSA0000224257420000096
和σi为大于0的常数;由于增益系数随时间而快速衰减,因此也避免了系统初始值与估计值不同而导致的峰值问题;
步骤3.2)定义速度和故障状态的估计误差如式(7)所示:
Figure BSA0000224257420000097
对式(7)求导,可以得到系统的误差动态如式(8)所示:
Figure BSA0000224257420000098
步骤3.3)在式(4)成立的条件下,对于状态空间模型(2)和对应设计的观测器(5),如果参数λ1和λ2满足
Figure BSA0000224257420000099
则估计误差(7)最终一致有界且系统的误差动态(8)是有限时间稳定的,即:
Figure BSA0000224257420000101
Figure BSA0000224257420000102
其中,T为有限收敛时间;V1(0)为选取的李雅普诺夫函数在t=0时的值,且有
Figure BSA0000224257420000103
Figure BSA0000224257420000104
且P和Q均为正定矩阵;λmin{·}和λmax{·}分别表示矩阵的最小和最大特征根;通过为λ1和λ2选取合适的值,可以将估计误差调整到极小的值;
步骤4)设计容错控制算法,包括如下步骤:
步骤4.1)根据第i个智能体获取的邻居信息,定义一致性跟踪位置误差变量exi(t)和速度误差变量evi(t),其形式如下:
Figure BSA0000224257420000105
定义
Figure BSA0000224257420000106
Figure BSA0000224257420000107
式(13)可改写为如式(14)所示的向量形式:
Figure BSA0000224257420000108
其中,Im为m维的单位矩阵;定义
Figure BSA0000224257420000109
Figure BSA0000224257420000111
可得一致性跟踪误差系统如式(15)所示:
Figure BSA0000224257420000112
步骤4.2)设计积分滑模面函数如式(16)所示
Figure BSA0000224257420000113
其中
Figure BSA0000224257420000114
si(t),i=1,...,n为滑模变量;k1和k2为正常数;
步骤4.3)为跟随者智能体设计如式(17)所示的容错一致性控制算法:
Figure BSA0000224257420000115
其中γ和ε为正常数;通过理论推导可以得到(1)和(2)组成的多智能体系统在容错控制律(17)下能够实现执行器故障情况下的有限时间一致性。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
下面以实际案例仿真说明实施方案的有效性。
采用由加拿大Quanser公司研制的Qball-X4四旋翼直升机飞行控制系统作为应用研究对象。由包含标记为0的领导者和标记为i,i=1,...,5的跟随者的6架Qball-X4四旋翼直升机组成的多智能体系统,其通信拓扑结构如图2所示。设定通信拓扑图中每条边的权值均为1,则根据通信拓扑结构,可计算得出拉普拉斯矩阵L和邻接矩阵B的具体表达式为:
Figure BSA0000224257420000116
对于该多智能体系统的模型,领导智能体的状态空间模型描述如下:
Figure BSA0000224257420000121
其中,领导智能体的初始状态值为x0(0)=0,v0(0)=1,控制输入为u0=sin(t/15);带有执行器加性故障和固有非线性动态的跟随智能体i,i=1,...,5的状态空间模型描述如下:
Figure BSA0000224257420000122
其中固有非线性动态fi(xi(t),vi(t),t)=sin(xi(t))+sin(vi(t)),且满足||fi(xi(t),vi(t),t)||≤2;跟随智能体的初始状态值分别为x1(0)=1,x2(0)=1.2,x3(0)=0.5,x4(0)=-1,x5(0)=-1.5,v1(0)=1.5,v2(0)=1.2,v3(0)=2,v4(0)=0.5,v5(0)=0。
为了说明本发明设计的容错控制算法在解决带有执行器故障的多智能体系统的-致性控制问题时的有效性和优越性,假设跟随智能体2和4发生执行器加性故障,其余飞行器无故障。其中,飞行器2的执行器故障描述为:
β2(t)=0.2+0.1sin(0.5πt)
飞行器4的执行器故障描述为:
β4(t)=0.1(1-e-0.2t)+0.2sin(0.5πt)cos(0.3t)
显然,这些执行器故障函数取值存在上限,可以取
Figure BSA0000224257420000123
因此,根据步骤3.3),可以将超螺旋二阶滑模观测器中的待定参数设置为λ1=5,λ2=13;时变增益系数ηi(t)中的参数设置为
Figure BSA0000224257420000124
σ=5。为了更好的体现本发明设计的观测器在故障估计时的快速性和准确性,我们选取相同的故障情况,分别得出本观测器和传统观测器在进行故障估计时的对比曲线(见图3和图4)。通过图像可知,对比传统观测器,本发明设计的观测器的故障观测曲线表现更为平滑且峰值更低,能够有效弥补初始误差过大带来的影响;与此同时,在故障估计的耗时方面,对于两种不同的执行器故障函数,本发明设计的观测器的估计时间均小于传统观测器,从而表明所设计的观测器在处理执行器加性故障时能够实现更快速的故障估计。
选取滑模面系数为k1=2,k2=1;容错控制器中的参数定为γ=0.5,ε=0.15。图5和图6为多智能体系统在控制器作用下各个智能体的位置和速度轨迹曲线。可以看到,尽管部分跟随智能体受到不同程度的执行器加性故障的影响,但在固定拓扑结构的基础上,通过及时的信息交互,仍能在有限时间内快速收敛误差,追踪上领导者。图7为各个智能体的控制律曲线,可以看出在本发明所设计的滑模容错一致性控制律的作用下,控制信号在收敛过程中曲线平滑。图8为智能体2在一般的容错控制算法和本发明设计的容错控制算法下的控制律曲线的对比,可以看出,由于本发明设计的容错控制算法能够更好地抑制抖振,因此控制律曲线收敛速度更快,波动更为平稳,即多智能体系统在本发明设计地控制器地控制下更为稳定,效果更为优秀。综上,对于发生执行器加性故障的多智能体系统,本案例仿真的控制方法是行之有效的。

Claims (1)

1.一种超螺旋观测器下的多智能体系统滑模容错一致性控制算法,其特征在于,包括如下具体步骤:
步骤1)确定本发明所适用的多智能体系统模型及其参数,包括如下步骤:
步骤1.1)确定领导智能体的状态空间模型如式(1)所示:
Figure FSA0000224257410000011
其中,x0(t)、v0(t)∈Rm分别表示领导者的位置和速度状态;u0(t)∈Rm为控制输入。
步骤1.2)确定带有执行器加性故障和固有非线性动态的跟随智能体的状态空间模型如式(2)所示:
Figure FSA0000224257410000012
其中,xi(t)、vi(t)∈Rm分别表示第i个智能体的位置和速度状态;ui(t)∈Rm为控制输入;fi(xi(t),vi(t),t)∈Rm表示第i个智能体的固有非线性动态;βi(t)∈Rm代表第i个智能体时变时产生的执行器加性故障,βi(t)=0表示无故障发生,βi(t)≠0表示第i个智能体的执行器发生了加性故障;
对于固有非线性动态fi(xi(t),vi(t),t)和执行器加性故障βi(t)以及其一阶微分,满足下面不等式(3)和(4):
Figure FSA0000224257410000013
其中
Figure FSA0000224257410000014
为非负常数;
Figure FSA0000224257410000015
其中
Figure FSA0000224257410000016
为非负常数;
步骤2)确定本发明所适用的多智能体系统的通信拓扑结构:
考虑由一个领导者和n个跟随者组成的多智能体系统的通信拓扑结构,领导者标记为0,跟随者标记为i=1,2,...,n;该结构可利用图
Figure FSA0000224257410000017
来对包括领导者和跟随者在内所有节点之间的信息交互进行描述,其中,
Figure FSA0000224257410000018
表示图中各节点的集合,
Figure FSA0000224257410000021
表示边集合;定义子图G={V,E}以表示跟随者之间的通信拓扑网络,其中,V={1,2,...,n}表示子图中各节点的集合,
Figure FSA0000224257410000022
表示边集合;定义Ni={j|(i,j)∈E,i≠j}表示与节点i相邻的各节点集合,A=(aij)为子图G的邻接矩阵,如果(i,j)∈E,i≠j,即表示节点i和j之间有直接的路径,那么aij>0;否则,aij=0;定义节点i的入度为
Figure FSA0000224257410000023
则子图G的入度矩阵可表示为D=diag{d1,...,dn}。由此可得到子图的拉普拉斯矩阵表示为L=[lij]n×n=D-A;
对于领导者和跟随者之间的邻接矩阵,可定义为B=diag(b1,b2,...,bn);当跟随者i能够直接获得来自领导者的信息时,bi>0;否则,bi=0;
步骤3)构造超螺旋二阶滑模故障观测器,包括如下步骤:
步骤3.1)首先为每个跟随智能体构造如式(5)所示的超螺旋观测器:
Figure FSA0000224257410000024
其中,
Figure FSA0000224257410000025
Figure FSA0000224257410000026
分别表示第i个智能体的速度估计和执行器加性故障估计,α为辅助变量,λ1和λ2为待定参数;设计增益系数ηi(t),可以保证系统在初始误差较大时也能快速收敛至终态,其表达式如式(6)所示:
Figure FSA0000224257410000027
其中,
Figure FSA0000224257410000028
和σi为大于0的常数;由于增益系数随时间而快速衰减,因此也避免了系统初始值与估计值不同而导致的峰值问题;
步骤3.2)定义速度和故障状态的估计误差如式(7)所示:
Figure FSA0000224257410000029
对式(7)求导,可以得到系统的误差动态如式(8)所示:
Figure FSA0000224257410000031
步骤3.3)在式(4)成立的条件下,对于状态空间模型(2)和对应设计的观测器(5),如果参数λ1和λ2满足
Figure FSA0000224257410000032
则估计误差(7)最终一致有界且系统的误差动态(8)是有限时间稳定的,即:
Figure FSA0000224257410000033
Figure FSA0000224257410000034
其中,T为有限收敛时间;V1(0)为选取的李雅普诺夫函数在t=0时的值,且有
Figure FSA0000224257410000035
Figure FSA0000224257410000036
且P和Q均为正定矩阵;λmin{·}和λmax{·}分别表示矩阵的最小和最大特征根;通过为λ1和λ2选取合适的值,可以将估计误差调整到极小的值;
步骤4)设计容错控制算法,包括如下步骤:
步骤4.1)根据第i个智能体获取的邻居信息,定义一致性跟踪位置误差变量exi(t)和速度误差变量evi(t),其形式如下:
Figure FSA0000224257410000037
定义
Figure FSA0000224257410000041
Figure FSA0000224257410000042
式(13)可改写为如式(14)所示的向量形式:
Figure FSA0000224257410000043
其中,Im为m维的单位矩阵;定义
Figure FSA0000224257410000044
Figure FSA0000224257410000045
可得一致性跟踪误差系统如式(15)所示:
Figure FSA0000224257410000046
步骤4.2)设计积分滑模面函数如式(16)所示
Figure FSA0000224257410000047
其中
Figure FSA0000224257410000048
si(t),i=1,...,n为滑模变量;k1和k2为正常数;
步骤4.3)为跟随者智能体设计如式(17)所示的容错一致性控制算法:
Figure FSA0000224257410000049
其中γ和ε为正常数;通过理论推导可以得到(1)和(2)组成的多智能体系统在容错控制律(17)下能够实现执行器故障情况下的有限时间一致性。
CN202011257346.1A 2020-11-11 2020-11-11 超螺旋观测器下的多智能体系统滑模容错一致性控制算法 Pending CN112305918A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011257346.1A CN112305918A (zh) 2020-11-11 2020-11-11 超螺旋观测器下的多智能体系统滑模容错一致性控制算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011257346.1A CN112305918A (zh) 2020-11-11 2020-11-11 超螺旋观测器下的多智能体系统滑模容错一致性控制算法

Publications (1)

Publication Number Publication Date
CN112305918A true CN112305918A (zh) 2021-02-02

Family

ID=74326674

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011257346.1A Pending CN112305918A (zh) 2020-11-11 2020-11-11 超螺旋观测器下的多智能体系统滑模容错一致性控制算法

Country Status (1)

Country Link
CN (1) CN112305918A (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113093554A (zh) * 2021-04-14 2021-07-09 西北工业大学 一种基于Time Base Generator的固定时间动态区域跟踪控制方法
CN113241794A (zh) * 2021-05-28 2021-08-10 合肥工业大学 一种基于多智能体的孤岛微电网自适应控制方法
CN113391550A (zh) * 2021-05-19 2021-09-14 山东大学 多智能体储能电池一致性控制方法及系统
CN113472242A (zh) * 2021-07-05 2021-10-01 江南大学 基于多智能体的抗干扰自适应模糊滑膜协同控制方法
CN113741192A (zh) * 2021-09-06 2021-12-03 杭州电子科技大学 基于可切换拓扑的时滞多智能体系统约束容错控制方法
CN113885499A (zh) * 2021-10-08 2022-01-04 四川大学 腔体内检测用机器人轨迹容错控制方法
CN114326664A (zh) * 2021-12-22 2022-04-12 同济大学 一种非线性多智能体的容错控制器设计方法及存储介质
CN115840426A (zh) * 2022-11-15 2023-03-24 华东交通大学 一种基于事件的多智能体系统主动抗扰一致性控制方法

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113093554A (zh) * 2021-04-14 2021-07-09 西北工业大学 一种基于Time Base Generator的固定时间动态区域跟踪控制方法
CN113391550A (zh) * 2021-05-19 2021-09-14 山东大学 多智能体储能电池一致性控制方法及系统
CN113241794A (zh) * 2021-05-28 2021-08-10 合肥工业大学 一种基于多智能体的孤岛微电网自适应控制方法
CN113241794B (zh) * 2021-05-28 2022-08-26 合肥工业大学 一种基于多智能体的孤岛微电网自适应控制方法
CN113472242A (zh) * 2021-07-05 2021-10-01 江南大学 基于多智能体的抗干扰自适应模糊滑膜协同控制方法
CN113741192A (zh) * 2021-09-06 2021-12-03 杭州电子科技大学 基于可切换拓扑的时滞多智能体系统约束容错控制方法
CN113741192B (zh) * 2021-09-06 2024-05-07 杭州电子科技大学 基于可切换拓扑的时滞多智能体系统约束容错控制方法
CN113885499A (zh) * 2021-10-08 2022-01-04 四川大学 腔体内检测用机器人轨迹容错控制方法
CN114326664A (zh) * 2021-12-22 2022-04-12 同济大学 一种非线性多智能体的容错控制器设计方法及存储介质
CN114326664B (zh) * 2021-12-22 2023-08-29 同济大学 一种非线性多智能体的容错控制器设计方法及存储介质
CN115840426A (zh) * 2022-11-15 2023-03-24 华东交通大学 一种基于事件的多智能体系统主动抗扰一致性控制方法
CN115840426B (zh) * 2022-11-15 2024-05-14 华东交通大学 一种基于事件的多智能体系统主动抗扰一致性控制方法

Similar Documents

Publication Publication Date Title
CN112305918A (zh) 超螺旋观测器下的多智能体系统滑模容错一致性控制算法
CN109557818B (zh) 具有多故障的多智能体跟踪系统的滑模容错控制方法
CN109116736B (zh) 基于滑模的线性多智能体系统执行器故障的容错控制方法
CN109445447B (zh) 一种多智能体编队跟踪控制方法及系统
CN108803316B (zh) 针对多智能体系统执行器故障的主动容错控制方法
Yu et al. Decentralized fractional-order backstepping fault-tolerant control of multi-UAVs against actuator faults and wind effects
CN109407520B (zh) 基于滑模控制的二阶多智能体系统的容错一致性控制算法
CN110058519B (zh) 一种基于快速自适应技术的主动编队容错控制方法
Yang et al. Distributed formation control of nonholonomic autonomous vehicle via RBF neural network
CN108427401B (zh) 一种具有联合连通拓扑的飞行控制系统协同故障诊断方法
CN108333949A (zh) 针对多智能体系统执行器部分失效故障的滑模容错控制方法
CN109986562B (zh) 一种平面机械臂系统及其分散式容错控制方法
CN112527016B (zh) 一种智能集群一体化容错时变编队控制方法及系统
CN109884902B (zh) 一种基于区间观测器无人机编队系统故障检测方法
CN110497415B (zh) 一种基于干扰观测器的多机械臂系统的一致控制方法
CN112947086B (zh) 一种无人机和无人车组成的异构多智能体系统编队控制中执行器故障的自适应补偿方法
CN112650290B (zh) 带有扰动补偿的高层消防无人机编队优化控制方法及系统
Liu et al. Adaptive distributed finite-time formation control for multi-UAVs under input saturation without collisions
CN114237041A (zh) 一种基于预设性能的空-地协同固定时间容错控制方法
CN116466588A (zh) 多智能体系统的有限时间时变编队跟踪控制方法及系统
CN117055605A (zh) 多无人机姿态控制方法及系统
Li et al. Learning-observer-based adaptive tracking control of multiagent systems using compensation mechanism
CN114326664A (zh) 一种非线性多智能体的容错控制器设计方法及存储介质
CN116610106B (zh) 一种多无人机容错跟踪控制方法
CN111806728B (zh) 具有执行器饱和的航天器多智能体姿态同步容错控制方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination