CN113325719A - 时变时滞多智能体的分布式脉冲一致性控制方法及系统 - Google Patents

时变时滞多智能体的分布式脉冲一致性控制方法及系统 Download PDF

Info

Publication number
CN113325719A
CN113325719A CN202110667285.4A CN202110667285A CN113325719A CN 113325719 A CN113325719 A CN 113325719A CN 202110667285 A CN202110667285 A CN 202110667285A CN 113325719 A CN113325719 A CN 113325719A
Authority
CN
China
Prior art keywords
time
agent
pulse
lag
model
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.)
Granted
Application number
CN202110667285.4A
Other languages
English (en)
Other versions
CN113325719B (zh
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.)
Jiangnan University
Original Assignee
Jiangnan University
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 Jiangnan University filed Critical Jiangnan University
Priority to CN202110667285.4A priority Critical patent/CN113325719B/zh
Publication of CN113325719A publication Critical patent/CN113325719A/zh
Application granted granted Critical
Publication of CN113325719B publication Critical patent/CN113325719B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G05B13/042Adaptive 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 in which a parameter or coefficient is automatically adjusted to optimise the performance

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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种时变时滞多智能体的分布式脉冲一致性控制方法及系统,包括以下步骤:构建含有多个独立智能体的时变时滞非线性多智能模型,基于时变时滞非线性多智能模型获得目标状态模型;基于脉冲牵制控制方法构建分布式脉冲牵制控制器;根据智能体的时变时滞非线性多智能模型和分布式脉冲牵制控制器,定义误差向量并建立对应的误差系统模型;基于误差系统模型构建李雅普诺夫函数并证明其导数值有上界;求得脉冲发生时刻李雅普诺夫函数值的关系;利用脉冲比较原理和参数变分方法进行上界分析,得出带有时变时滞的非线性多智能体系统能够实现一致性。其系统的一致性误差将会以指数形式收敛到零,正确率高。

Description

时变时滞多智能体的分布式脉冲一致性控制方法及系统
技术领域
本发明涉及信息技术领域,尤其是指一种时变时滞多智能体的分布式脉冲一致性控制方法及系统。
背景技术
过去的几十年以来,非线性多智能体一致性问题由于在无人汽车驾驶,无人机编队,通信网络,无线传感网络的控制等方面有着十分广阔的应用前景,吸引了国内外大量学者的目光。在多智能体系统中,每个独立的智能体都拥有一定的独立性和一定的学习能力。但多智能体是一个讲究协调的系统,需要多智能体中的每一个智能体达到共同的状态,从而实现效率的最大化。
随着人们对于非线性多智能体研究的不断深入,脉冲牵制控制开始更多地被应用到非线性多智能体一致性研究领域。基于脉冲微分方程的脉冲牵制控制方法可以使得每个智能体之间的状态信息只在脉冲时刻进行传输,大大减少了智能体间信息传输的数量。由于脉冲输入本质上属于不连续的信号,而智能体之间的信息交换是连续的且智能体一般是分布式的,因此通常采用分布式脉冲牵制控制。很明显的是,分布式脉冲牵制控制可以很好地减少控制成本,并且对于实际生活中大型的多智能体来说,实际应用效率更高。
从系统的理论观点来看,时变时滞的存在使得任何存在于现实世界上的系统的过去状态不可避免地对系统当前的状态产生一定影响,即系统演化趋势不仅仅依赖于当前的状态,也取决于过去某一时刻或者过去若干个时刻的状态。在系统建模中,如果没有很好地描述或者考虑系统中时变时滞的存在,会产生设计缺陷或者得到不正确的分析结论。
因此,对于非线性多智能体的一致性问题,考虑系统中时变时滞带来的影响是至关重要的。
发明内容
为此,本发明所要解决的技术问题在于克服现有技术中如果没有很好地描述或者考虑系统中时变时滞的存在,会产生设计缺陷或者得到不正确的分析结论的技术缺陷。
为解决上述技术问题,本发明提供了一种时变时滞多智能体的分布式脉冲一致性控制方法,包括以下步骤:
构建含有多个独立智能体的时变时滞非线性多智能模型,基于时变时滞非线性多智能模型获得目标状态模型;
基于脉冲牵制控制方法构建分布式脉冲牵制控制器;
根据智能体的时变时滞非线性多智能模型和分布式脉冲牵制控制器,定义误差向量并建立对应的误差系统模型;
基于误差系统模型构建李雅普诺夫函数并证明其导数值有上界;
求得脉冲发生时刻李雅普诺夫函数值的关系;
根据李雅普诺夫函数的导数值有上界和脉冲发生时刻李雅普诺夫函数值的关系,利用脉冲比较原理和参数变分方法进行上界分析,得出带有时变时滞的非线性多智能体系统能够实现一致性。
作为优选的,所述构建含有多个独立智能体的时变时滞非线性多智能模型,具体包括:
Figure BDA0003117367280000021
其中,
Figure BDA0003117367280000022
是第i个智能体的状态向量;系统中矩阵A,B,D为常数矩阵;
系统中的非线性函数
Figure BDA0003117367280000031
并且满足条件
Figure BDA0003117367280000032
Figure BDA0003117367280000033
Figure BDA0003117367280000034
τ(t)为系统中的时变时滞;
Figure BDA0003117367280000035
均满足李普希兹条件。
作为优选的,所述基于时变时滞非线性多智能模型获得目标状态模型,具体包括:
目标状态s(t)的模型:
Figure BDA0003117367280000036
作为优选的,所述基于脉冲牵制控制方法构建分布式脉冲牵制控制器,包括:
分布式脉冲牵制控制器
Figure BDA0003117367280000037
其中,c是系统的耦合强度;
Figure BDA0003117367280000038
为系统的牵制增益;
矩阵L=(lij)N×N为拉普拉斯矩阵,满足耗散条件;函数δ为狄拉克函数;
当且仅当具有目标状态的智能体与第i个智能体存在直接交换信息的路径时,
Figure BDA0003117367280000039
否则
Figure BDA00031173672800000310
智能体i称为被牵制智能体或者被控智能体;
脉冲序列
Figure BDA00031173672800000311
满足:
Figure BDA00031173672800000312
作为优选的,所述根据智能体的时变时滞非线性多智能模型和分布式脉冲牵制控制器,定义误差向量并建立对应的误差系统模型,具体包括:
定义误差向量ei(t)=xi(t)-s(t);
所述误差模型对应的误差系统:
Figure BDA0003117367280000041
其中,非线性函数
Figure BDA0003117367280000042
满足李普希兹条件;
矩阵
Figure BDA0003117367280000043
为牵制控制矩阵,代表着牵制控制策略,当第i个智能体被牵制时,则对应的
Figure BDA0003117367280000044
否则
Figure BDA0003117367280000045
e(t)在t=tk右连续,即
Figure BDA0003117367280000046
误差系统的初值ei(t)=Φi(t),-τ*≤t≤0,i=1,2,…,N,
其中,
Figure BDA0003117367280000047
是从[-τ*,0]至
Figure BDA0003117367280000048
的连续函数集。
作为优选的,所述基于误差系统模型构建李雅普诺夫函数并证明其导数值有上界,包括:
构建如下的李雅普诺夫函数为:
Figure BDA0003117367280000049
根据迪尼导数的定义,可以得出:
Figure BDA0003117367280000051
两项相加,可以得出:
Figure BDA0003117367280000052
将误差系统模型带入上式得:
Figure BDA0003117367280000053
对上式中的每一项进行处理,可以得出:
Figure BDA0003117367280000054
针对函数f1利用李普希兹条件,可得:
Figure BDA0003117367280000055
可以得出:
Figure BDA0003117367280000056
Figure BDA0003117367280000057
则得出:
Figure BDA0003117367280000058
针对函数f2利用李普希兹条件,得出:
Figure BDA0003117367280000061
Figure BDA0003117367280000062
转变为
Figure BDA0003117367280000063
得出:
Figure BDA0003117367280000064
根据基本不等式可以得出:
Figure BDA0003117367280000065
Figure BDA0003117367280000066
Figure BDA0003117367280000067
将李雅普诺夫函数V(t)带入到上式中,得到:
Figure BDA0003117367280000068
作为优选的,所述求得脉冲发生时刻李雅普诺夫函数值的关系,包括:
当t=tk时,则:
Figure BDA0003117367280000069
其中,
Figure BDA00031173672800000610
Figure BDA00031173672800000611
若下式成立:
Figure BDA0003117367280000071
则得出在脉冲时刻时的关系如下:
Figure BDA0003117367280000072
即:
Figure BDA0003117367280000073
作为优选的,所述根据李雅普诺夫函数的导数值有上界和脉冲发生时刻李雅普诺夫函数值的关系,利用脉冲比较原理和参数变分方法进行上界分析,得出带有时变时滞的非线性多智能体系统能够实现一致性,包括:
设存在常数μ>0,对称矩阵P,
Figure BDA0003117367280000074
平均脉冲间隔小于Na
令不等式满足:
Figure BDA0003117367280000075
矩阵满足:
Figure BDA0003117367280000076
则当0<μ≤1时,时变时滞非线性多智能模型以如下全局指数的形式达到一致性,且收敛率为λ:
Figure BDA0003117367280000077
当μ>1时,通过所设计的分布式脉冲牵制控制器,时变时滞非线性多智能模型以如下全局指数的形式达到一致性,且收敛率为
Figure BDA0003117367280000081
Figure BDA0003117367280000082
其中,参数λ与
Figure BDA0003117367280000083
分别为方程
Figure BDA0003117367280000084
以及方程
Figure BDA0003117367280000085
的特解;
参数
Figure BDA0003117367280000086
参数
Figure BDA0003117367280000087
l1与l2为李普希兹条件中的李普希兹常数;L为拉普拉斯矩阵;In为n阶单位矩阵;INn为Nn阶单位矩阵;sup表示上确界。
本发明公开了一种时变时滞非线性多智能体的分布式脉冲一致性控制系统,包括:
时变时滞非线性模型构建模块,所述时变时滞非线性模型构建模块用于构建含有多个独立智能体的时变时滞非线性多智能模型,基于时变时滞非线性多智能模型获得目标状态模型;
脉冲牵制控制器构建模块,所述脉冲牵制控制器构建模块用于基于脉冲牵制控制方法构建分布式脉冲牵制控制器;
误差系统模型构建模块,所述误差系统模型构建模块用于根据智能体的时变时滞非线性多智能模型和分布式脉冲牵制控制器,定义误差向量并建立对应的误差系统模型;
函数构建模块,所述函数构建模块基于误差系统模型构建李雅普诺夫函数并证明其导数值有上界;
关系求取模块,所述关系求取模块用于求得脉冲发生时刻李雅普诺夫函数值的关系;
系统一致性模块,所述系统一致性模块根据李雅普诺夫函数的导数值有上界和脉冲发生时刻李雅普诺夫函数值的关系,利用脉冲比较原理和参数变分方法进行上界分析,得出带有时变时滞的非线性多智能体系统能够实现一致性。
本发明公开了一种时变时滞非线性多智能系统,包括上述的时变时滞非线性多智能体的分布式脉冲一致性控制系统。
本发明的上述技术方案相比现有技术具有以下优点:
1、本发明采用分布式脉冲牵制控制策略,基于设计的分布式脉冲牵制控制器,应用蔡氏电路模型,李雅普诺夫法,脉冲比较原理,参数变分法以及平均脉冲间隔的定义,基于脉冲的不同作用效果,提出了一种使得带有时变时滞非线性多智能体系统达到一致性的新型方法并分别给出了系统达到一致性时需要满足的条件以及具体的收敛速度。
2、本发明中,带有时变时滞的非线性多智能体系统的一致性与耦合强度、牵制增益以及平均脉冲间隔等参数有关,在引入设计的分布式脉冲牵制控制器之后,带有时变时滞的非线性多智能体系统的一致性误差将会以指数形式收敛到零。
附图说明
图1为带有时变时滞的非线性多智能体系统实现一致性的流程图;
图2为仿真对象的拓扑结构图;
图3为目标状态示意图;
图4为脉冲示意图;
图5为系统的一致性误差演化曲线;
图6为误差状态||ei||(i=1,2,3,4,5,6)演化曲线;
图7为智能体状态||xi||(i=1,2,3,4,5,6)演化曲线。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好地理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
参照图1-图7所示,本发明公开了一种时变时滞多智能体的分布式脉冲一致性控制方法,包括以下步骤:
步骤一、构建含有多个独立智能体的时变时滞非线性多智能模型,基于时变时滞非线性多智能模型获得目标状态模型,具体包括:
Figure BDA0003117367280000101
其中,
Figure BDA0003117367280000102
是第i个智能体的状态向量;系统中矩阵A,B,D为常数矩阵;
系统中的非线性函数
Figure BDA0003117367280000103
并且满足条件
Figure BDA0003117367280000104
Figure BDA0003117367280000105
Figure BDA0003117367280000106
τ(t)为系统中的时变时滞;
Figure BDA0003117367280000107
均满足李普希兹条件。
目标状态s(t)的模型:
Figure BDA0003117367280000108
步骤二、基于脉冲牵制控制方法构建分布式脉冲牵制控制器,包括:
分布式脉冲牵制控制器
Figure BDA0003117367280000111
其中,c是系统的耦合强度;
Figure BDA0003117367280000112
为系统的牵制增益;
矩阵L=(lij)N×N为拉普拉斯矩阵,满足耗散条件;函数δ为狄拉克函数;
当且仅当具有目标状态的智能体与第i个智能体存在直接交换信息的路径时,
Figure BDA0003117367280000113
否则
Figure BDA0003117367280000114
智能体i称为被牵制智能体或者被控智能体;
脉冲序列
Figure BDA0003117367280000115
满足:
Figure BDA0003117367280000116
步骤三、根据智能体的时变时滞非线性多智能模型和分布式脉冲牵制控制器,定义误差向量并建立对应的误差系统模型,具体包括:
定义误差向量ei(t)=xi(t)-s(t);
所述误差模型对应的误差系统:
Figure BDA0003117367280000117
其中,非线性函数
Figure BDA0003117367280000118
满足李普希兹条件;
矩阵
Figure BDA0003117367280000119
为牵制控制矩阵,代表着牵制控制策略,当第i个智能体被牵制时,则对应的
Figure BDA00031173672800001110
否则
Figure BDA00031173672800001111
e(t)在t=tk右连续,即
Figure BDA00031173672800001112
误差系统的初值ei(t)=Φi(t),-τ*≤t≤0,i=1,2,...,N,
其中,
Figure BDA0003117367280000121
是从[-τ*,0]至
Figure BDA0003117367280000122
的连续函数集。
步骤四、基于误差系统模型构建李雅普诺夫函数并证明其导数值有上界,包括:
构建如下的李雅普诺夫函数为:
Figure BDA0003117367280000123
根据迪尼导数的定义,可以得出:
Figure BDA0003117367280000124
两项相加,可以得出:
Figure BDA0003117367280000125
将误差系统模型带入上式得:
Figure BDA0003117367280000126
对上式中的每一项进行处理,可以得出:
Figure BDA0003117367280000127
针对函数f1利用李普希兹条件,可得:
Figure BDA0003117367280000128
可以得出:
Figure BDA0003117367280000131
Figure BDA0003117367280000132
则得出:
Figure BDA0003117367280000133
针对函数f2利用李普希兹条件,得出:
Figure BDA0003117367280000134
Figure BDA0003117367280000135
转变为
Figure BDA0003117367280000136
得出:
Figure BDA0003117367280000137
根据基本不等式可以得出:
Figure BDA0003117367280000138
Figure BDA0003117367280000139
Figure BDA00031173672800001310
将李雅普诺夫函数V(t)带入到上式中,得到:
Figure BDA0003117367280000141
步骤五、求得脉冲发生时刻李雅普诺夫函数值的关系,包括:
包括:
当t=tk时,则:
Figure BDA0003117367280000142
其中,
Figure BDA0003117367280000143
Figure BDA0003117367280000144
若下式成立:
Figure BDA0003117367280000145
则得出在脉冲时刻时的关系如下:
Figure BDA0003117367280000146
即:
Figure BDA0003117367280000147
步骤六、根据李雅普诺夫函数的导数值有上界和脉冲发生时刻李雅普诺夫函数值的关系,利用脉冲比较原理和参数变分方法进行上界分析,得出带有时变时滞的非线性多智能体系统能够实现一致性,包括:
设存在常数μ>0,对称矩阵P,
Figure BDA0003117367280000148
平均脉冲间隔小于Na
令不等式满足:
Figure BDA0003117367280000151
矩阵满足:
Figure BDA0003117367280000152
则当0<μ≤1时,时变时滞非线性多智能模型以如下全局指数的形式达到一致性,且收敛率为λ:
Figure BDA0003117367280000153
当μ>1时,通过所设计的分布式脉冲牵制控制器,时变时滞非线性多智能模型以如下全局指数的形式达到一致性,且收敛率为
Figure BDA0003117367280000154
Figure BDA0003117367280000155
其中,参数λ与
Figure BDA0003117367280000156
分别为方程
Figure BDA0003117367280000157
以及方程
Figure BDA0003117367280000158
的特解;
参数
Figure BDA0003117367280000159
参数
Figure BDA00031173672800001510
l1与l2为李普希兹条件中的李普希兹常数;L为拉普拉斯矩阵;In为n阶单位矩阵;INn为Nn阶单位矩阵;sup表示上确界。
本发明公开了一种时变时滞非线性多智能体的分布式脉冲一致性控制系统,包括时变时滞非线性模型构建模块、脉冲牵制控制器构建模块、误差系统模型构建模块、函数构建模块、关系求取模块和系统一致性模块。
所述时变时滞非线性模型构建模块用于构建含有多个独立智能体的时变时滞非线性多智能模型,基于时变时滞非线性多智能模型获得目标状态模型。
所述脉冲牵制控制器构建模块用于基于脉冲牵制控制方法构建分布式脉冲牵制控制器。
所述误差系统模型构建模块用于根据智能体的时变时滞非线性多智能模型和分布式脉冲牵制控制器,定义误差向量并建立对应的误差系统模型。
所述函数构建模块基于误差系统模型构建李雅普诺夫函数并证明其导数值有上界。
所述关系求取模块用于求得脉冲发生时刻李雅普诺夫函数值的关系;
所述系统一致性模块根据李雅普诺夫函数的导数值有上界和脉冲发生时刻李雅普诺夫函数值的关系,利用脉冲比较原理和参数变分方法进行上界分析,得出带有时变时滞的非线性多智能体系统能够实现一致性。
本发明还公开了一种时变时滞非线性多智能系统,包括上述的时变时滞非线性多智能体的分布式脉冲一致性控制系统。
下面,结合具体实施例,对本发明的技术方案做进一步说明。
首先考虑如下的单一智能体模型:
Figure BDA0003117367280000161
其中,
Figure BDA0003117367280000162
是第i个独立智能体的状态向量;系统中矩阵A,B,D为常数矩阵;系统中的非线性函数
Figure BDA0003117367280000163
并且满足条件为
Figure BDA0003117367280000164
Figure BDA0003117367280000165
并且满足条件为
Figure BDA0003117367280000166
τ(t)为系统中的时变时滞;
Figure BDA0003117367280000171
均满足李普希兹条件。
接下来,对于含有N个多智能体的系统,根据单一的智能体模型(1)可以得出带有时变时滞的非线性多智能体模型:
Figure BDA0003117367280000172
其中,
Figure BDA0003117367280000173
是第i个智能体的状态向量;系统中矩阵A,B,D为常数矩阵;系统中的非线性函数
Figure BDA0003117367280000174
并且满足条件
Figure BDA0003117367280000175
Figure BDA0003117367280000176
Figure BDA0003117367280000177
τ(t)为系统中的时变时滞;
Figure BDA0003117367280000178
均满足李普希兹条件。
以带有时变时滞的非线性多智能体(4-3)为基础,给出带有时变时滞系统的目标状态表达式:
Figure BDA0003117367280000179
随后,为了使得带有时变时滞的非线性多智能体模型(2)达到一致性,本发明设计了一种分布式脉冲牵制控制器如下:
Figure BDA00031173672800001710
其中,c是系统的耦合强度;
Figure BDA00031173672800001711
为系统的牵制增益;矩阵L=(lij)N×N为拉普拉斯矩阵,满足耗散条件;函数δ为狄拉克函数。值得说明的是,当且仅当具有目标状态的智能体与第i个智能体存在直接交换信息的路径时,
Figure BDA00031173672800001712
否则
Figure BDA00031173672800001713
系统中智能体i称为被牵制智能体或者被控智能体。
脉冲序列
Figure BDA00031173672800001714
满足式(5):
Figure BDA0003117367280000181
式(5)表明该脉冲序列是一个脉冲时刻严格单调递增的脉冲序列,可以用来作为系统中的不连续信号输入。通过引入分布式脉冲牵制控制器(5),带有时变时滞的非线性多智能体系统(2)将仅仅在离散脉冲时刻采集智能体本身及相邻智能体的状态信息,而后通过控制器产生的控制信号进行分布式脉冲牵制控制。该控制方法可以很好地减小控制成本,达到理想的控制效果。
随后,通过收集系统中每一个独立智能体的信息以及智能体之间相互交换的信息,定义误差ei(t)=xi(t)-s(t)。在加入设计的分布式脉冲牵制控制器(4)之后,可以得出带有时变时滞的非线性多智能体的误差系统:
Figure BDA0003117367280000182
其中,非线性函数
Figure BDA0003117367280000183
满足李普希兹条件;e(t)在t=tk右连续,即
Figure BDA0003117367280000184
系统(6)的初值为
ei(t)=Φi(t),-τ*≤t≤0,i=1,2,...,N (7)
其中,
Figure BDA0003117367280000185
是从[-τ*,0]至
Figure BDA0003117367280000186
的连续函数集。与无时滞情况类似,由于本文假设e(t)在t=tk右连续,因此对于
Figure BDA0003117367280000187
(7)式的解是在t=tk不连续时刻的分段右连续函数。
为了方便后续得出系统实现一致性的条件,可以在脉冲时刻t=tk利用克罗内克积将误差系统(6)写成下列矩阵的形式:
Figure BDA0003117367280000191
其中,k=1,2,...,;误差向量
Figure BDA0003117367280000192
矩阵L为拉普拉斯矩阵;矩阵
Figure BDA0003117367280000193
为牵制控制矩阵,代表着牵制控制策略:第i个智能体被牵制,则对应的
Figure BDA0003117367280000194
否则
Figure BDA0003117367280000195
定义:对于非线性多智能体系统中的任何初值Φi(i=1,2,...,N),存在λ>0,T0>0以及θ>0,如果||ei(t)||=||xi(t)-s(t)||≤θe-λt,t>T0,i=1,2,...,N成立,则带有时变时滞的非线性多智能体系统(2)将会以全局指数的形式达到目标状态s(t)从而实现一致性。
上述定义是全局指数一致的基本概念,在本发明中,它对于判断系统是否达到一致性是十分关键的。
接下来,主要探讨应用本发明所提出的方法,在引入设计的分布式脉冲牵制控制器(4)之后,带有时变时滞的非线性多智能体系统(2)实现一致性需要满足的条件。首先构建如下的李雅普诺夫函数为:
Figure BDA0003117367280000196
根据迪尼导数的定义,可以得出:
Figure BDA0003117367280000197
两项相加,可以得出:
Figure BDA0003117367280000198
将式(6)带入式(10)得:
Figure BDA0003117367280000201
接下来对于式(11)中的每一项进行处理,可以得出:
Figure BDA0003117367280000202
针对函数f1利用李普希兹条件,可得:
Figure BDA0003117367280000203
可以得出:
Figure BDA0003117367280000204
Figure BDA0003117367280000205
则可以得出:
Figure BDA0003117367280000206
针对函数f2利用李普希兹条件,可以得出:
Figure BDA0003117367280000207
Figure BDA0003117367280000208
转变为
Figure BDA0003117367280000209
可以得出:
Figure BDA00031173672800002010
根据基本不等式可以得出:
Figure BDA0003117367280000211
Figure BDA0003117367280000212
将式(12)与式(13)带入式(11)中可得:
Figure BDA0003117367280000213
将李雅普诺夫函数V(t)带入到上式中可得:
Figure BDA0003117367280000214
接下来考虑脉冲时刻的情况,当t=tk时,由式(8)可得:
Figure BDA0003117367280000215
其中
Figure BDA0003117367280000216
因此可以得出:
Figure BDA0003117367280000217
若下式成立:
Figure BDA0003117367280000218
则可以得出在脉冲时刻时的关系如下:
Figure BDA0003117367280000219
即:
Figure BDA00031173672800002110
基于式(14)与(15),给出如下带有时变时滞的脉冲比较系统,其中ε>0,υ(t)是如下脉冲系统的特解:
Figure BDA0003117367280000221
所以对于任何的t≥0,υ(t)≥V(t)≥0成立。
接下来应用参数变分法,可以得出如下对于υ(t)的积分等式:
Figure BDA0003117367280000222
其中,W(t,s)(t>s≥0)是如下线性脉冲系统的柯西矩阵:
Figure BDA0003117367280000223
情形1:当0<μ≤1,此时脉冲扰动有利于实现一致性,柯西矩阵的右侧可以通过考虑平均脉冲间隔的定义进行计算进而得出:
Figure BDA0003117367280000224
因此可以得出:
Figure BDA0003117367280000225
定义
Figure BDA0003117367280000231
因此可以得出:
Figure BDA0003117367280000232
将θ,β带入上式得:
Figure BDA0003117367280000233
定义函数
Figure BDA0003117367280000234
Figure BDA0003117367280000235
可以得出
Figure BDA0003117367280000236
并且得出ζ(∞)>0。由于
Figure BDA0003117367280000237
即该函数为单调递增函数,因此存在特解λ>0使得等式ζ(λ)=0成立,此时
Figure BDA0003117367280000238
基于上述分析,可以得出:
Figure BDA0003117367280000239
接下来,利用反证法证明在t>0时式(22)成立。假设t>0时式(22)不成立,则存在t*>0使得:
Figure BDA00031173672800002310
以及
Figure BDA00031173672800002311
然而存在t*>0,利用方程
Figure BDA0003117367280000241
以及式(21)等,可以得出:
Figure BDA0003117367280000242
在上式加入
Figure BDA0003117367280000243
可以得出:
Figure BDA0003117367280000244
进而可以得出:
Figure BDA0003117367280000245
将上式进行积分运算,可以得出:
Figure BDA0003117367280000246
将上式每一项进行整合并且根据
Figure BDA0003117367280000247
可以得出:
Figure BDA0003117367280000248
可以看出,式(23)成立,与假设相矛盾,则假设错误,式(22)成立。因此,基于上述推导与分析,在时间t≥0时,式(23)中ε→0,可以得出:
υ(t)≤θe-λt
利用脉冲比较原理,可以推导出:
Figure BDA0003117367280000251
情形2:当μ>1,此时脉冲扰动不利于实现一致性,柯西矩阵的左侧可以通过考虑平均脉冲间隔的定义进行计算进而得出:
Figure BDA0003117367280000252
与情形1的过程相类似,利用脉冲比较原理,参数变分等数学方法,可以得出:
Figure BDA0003117367280000253
可以看出,式(24)与式(26)符合之前定义的全局指数一致的基本概念,即通过采用本发明的方法,在引入设计的分布式脉冲牵制控制器(4)之后,带有时变时滞的非线性多智能体能够以全局指数的形式实现一致性并且当0<μ≤1时收敛率为λ;当μ>1时收敛率为
Figure BDA0003117367280000254
结论:通过上述的推导,给出了系统实现一致性时需要满足的条件。即:
假设存在常数μ>0,对称矩阵P,
Figure BDA0003117367280000255
平均脉冲间隔小于Na。若不等式满足
Figure BDA0003117367280000256
矩阵满足
Figure BDA0003117367280000257
则当0<μ≤1时,通过所设计的分布式脉冲牵制控制器,带有时变时滞的非线性多智能体系统将会以如下全局指数的形式达到一致性,且收敛率为λ:
Figure BDA0003117367280000261
当μ>1时,通过所设计的分布式脉冲牵制控制器,带有时变时滞的非线性多智能体系统将会以如下全局指数的形式达到一致性,且收敛率为
Figure BDA0003117367280000262
Figure BDA0003117367280000263
其中,参数λ与
Figure BDA0003117367280000264
分别为方程
Figure BDA0003117367280000265
以及方程
Figure BDA0003117367280000266
的特解;参数
Figure BDA0003117367280000267
参数
Figure BDA0003117367280000268
l1与l2为李普希兹条件中的李普希兹常数;L为拉普拉斯矩阵;In为n阶单位矩阵;INn为Nn阶单位矩阵;sup表示上确界。
注:在上述推导的过程中构造的李雅普诺夫函数结构较为简单从而简化了推导步骤。通过应用参数变分法与脉冲比较原理等数学方法,引入中间变量,给出了通过应用本发明的方法,在引入所设计的分布式脉冲牵制控制器(2)之后,带有时变时滞的非线性多智能体系统(4)需要达到一致性的具体条件。
步骤1:模型与参数的选取。
考虑如下包含6个智能体的时变时滞非线性多智能体系统:
Figure BDA0003117367280000269
以及目标状态:
Figure BDA00031173672800002610
其中,系统的拓扑结构见图2;智能体初始值在[0,1.1]随机选取;非线性函数
Figure BDA0003117367280000271
非线性函数
Figure BDA0003117367280000272
矩阵:
Figure BDA0003117367280000273
时变时滞τ(t)为:
Figure BDA0003117367280000274
满足耗散条件的拉普拉斯矩阵L为:
Figure BDA0003117367280000275
耦合强度为c=0.2,牵制矩阵
Figure BDA0003117367280000276
即非线性多智能体中智能体1,2,3,5,6被牵制。取常数μ=0.1,Na=0.02。
步骤2:Simulink仿真并得出结论。
针对带有时变时滞的非线性多智能体进行数值仿真,观测系统的一致性误差E,误差状态ei(i=1,2,3,4,5,6),智能体状态xi(t)依次见图(3),图(4)以及图(5)。从图中可以很明显地看出,通过应用本发明提出的方法,引入分布式脉冲牵制控制器(4)之后,带有时变时滞的非线性多智能体系统达到了一致性,验证了发明方法的正确与可靠性。
图2为带有时变时滞的非线性多智能体系统拓扑结构图;图3为系统需要达到的目标状态示意图;图4为脉冲信号示意图;图5为系统的一致性误差曲线,可以看出0.6s之后系统一致性误差指数收敛到零,系统实现一致性;图6为每一个智能体的误差状态,可以看出0.6s之后每一个智能体误差为零,系统实现一致性;图7为每一个智能体的状态,可以看出0.6s之后每一个智能体的状态均保持一致,系统实现一致性。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然,上述实施例仅仅是为清楚地说明所作的举例,并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (10)

1.一种时变时滞多智能体的分布式脉冲一致性控制方法,其特征在于,包括以下步骤:
构建含有多个独立智能体的时变时滞非线性多智能模型,基于时变时滞非线性多智能模型获得目标状态模型;
基于脉冲牵制控制方法构建分布式脉冲牵制控制器;
根据智能体的时变时滞非线性多智能模型和分布式脉冲牵制控制器,定义误差向量并建立对应的误差系统模型;
基于误差系统模型构建李雅普诺夫函数并证明其导数值有上界;
求得脉冲发生时刻李雅普诺夫函数值的关系;
根据李雅普诺夫函数的导数值有上界和脉冲发生时刻李雅普诺夫函数值的关系,利用脉冲比较原理和参数变分方法进行上界分析,得出带有时变时滞的非线性多智能体系统能够实现一致性。
2.根据权利要求1所述的时变时滞多智能体的分布式脉冲一致性控制方法,其特征在于,所述构建含有多个独立智能体的时变时滞非线性多智能模型,具体包括:
Figure FDA0003117367270000011
其中,
Figure FDA0003117367270000012
是第i个智能体的状态向量;系统中矩阵A,B,D为常数矩阵;
系统中的非线性函数
Figure FDA0003117367270000013
并且满足条件
Figure FDA0003117367270000021
Figure FDA0003117367270000022
τ(t)为系统中的时变时滞;
Figure FDA0003117367270000023
均满足李普希兹条件。
3.根据权利要求2所述的时变时滞多智能体的分布式脉冲一致性控制方法,其特征在于,所述基于时变时滞非线性多智能模型获得目标状态模型,具体包括:
目标状态s(t)的模型:
Figure FDA0003117367270000024
4.根据权利要求3所述的时变时滞多智能体的分布式脉冲一致性控制方法,其特征在于,所述基于脉冲牵制控制方法构建分布式脉冲牵制控制器,包括:
分布式脉冲牵制控制器
Figure FDA0003117367270000025
其中,c是系统的耦合强度;
Figure FDA0003117367270000026
为系统的牵制增益;
矩阵L=(lij)N×N为拉普拉斯矩阵,满足耗散条件;函数δ为狄拉克函数;
当且仅当具有目标状态的智能体与第i个智能体存在直接交换信息的路径时,
Figure FDA0003117367270000027
否则
Figure FDA0003117367270000028
智能体i称为被牵制智能体或者被控智能体;
脉冲序列
Figure FDA0003117367270000029
满足:
Figure FDA00031173672700000210
5.根据权利要求4所述的时变时滞多智能体的分布式脉冲一致性控制方法,其特征在于,所述根据智能体的时变时滞非线性多智能模型和分布式脉冲牵制控制器,定义误差向量并建立对应的误差系统模型,具体包括:
定义误差向量ei(t)=xi(t)-s(t);
所述误差模型对应的误差系统:
Figure FDA0003117367270000031
其中,非线性函数
Figure FDA0003117367270000032
满足李普希兹条件;
矩阵
Figure FDA0003117367270000033
为牵制控制矩阵,代表着牵制控制策略,当第i个智能体被牵制时,则对应的
Figure FDA0003117367270000034
否则
Figure FDA0003117367270000035
e(t)在t=tk右连续,即
Figure FDA0003117367270000036
误差系统的初值ei(t)=Φi(t),-τ*≤t≤0,i=1,2,…,N,
其中,
Figure FDA0003117367270000037
是从[-τ*,0]至
Figure FDA0003117367270000038
的连续函数集。
6.根据权利要求5所述的时变时滞多智能体的分布式脉冲一致性控制方法,其特征在于,所述基于误差系统模型构建李雅普诺夫函数并证明其导数值有上界,包括:
构建如下的李雅普诺夫函数为:
Figure FDA0003117367270000039
根据迪尼导数的定义,可以得出:
Figure FDA0003117367270000041
两项相加,可以得出:
Figure FDA0003117367270000042
将误差系统模型带入上式得:
Figure FDA0003117367270000043
对上式中的每一项进行处理,可以得出:
Figure FDA0003117367270000044
针对函数f1利用李普希兹条件,可得:
Figure FDA0003117367270000045
可以得出:
Figure FDA0003117367270000046
Figure FDA0003117367270000047
则得出:
Figure FDA0003117367270000048
针对函数f2利用李普希兹条件,得出:
Figure FDA0003117367270000051
Figure FDA0003117367270000052
转变为
Figure FDA0003117367270000053
得出:
Figure FDA0003117367270000054
根据基本不等式可以得出:
Figure FDA0003117367270000055
Figure FDA0003117367270000056
Figure FDA0003117367270000057
将李雅普诺夫函数V(t)带入到上式中,得到:
Figure FDA0003117367270000058
7.根据权利要求6所述的时变时滞多智能体的分布式脉冲一致性控制方法,其特征在于,所述求得脉冲发生时刻李雅普诺夫函数值的关系,包括:
当t=tk时,则:
Figure FDA0003117367270000059
其中,
Figure FDA00031173672700000510
Figure FDA00031173672700000511
若下式成立:
Figure FDA0003117367270000061
则得出在脉冲时刻时的关系如下:
Figure FDA0003117367270000062
即:
Figure FDA0003117367270000063
8.根据权利要求1所述的时变时滞多智能体的分布式脉冲一致性控制方法,其特征在于,所述根据李雅普诺夫函数的导数值有上界和脉冲发生时刻李雅普诺夫函数值的关系,利用脉冲比较原理和参数变分方法进行上界分析,得出带有时变时滞的非线性多智能体系统能够实现一致性,包括:
设存在常数μ>0,对称矩阵P,
Figure FDA0003117367270000064
平均脉冲间隔小于Na
令不等式满足:
Figure FDA0003117367270000065
矩阵满足:
Figure FDA0003117367270000066
则当0<μ≤1时,时变时滞非线性多智能模型以如下全局指数的形式达到一致性,且收敛率为λ:
Figure FDA0003117367270000071
当μ>1时,通过所设计的分布式脉冲牵制控制器,时变时滞非线性多智能模型以如下全局指数的形式达到一致性,且收敛率为λ:
Figure FDA0003117367270000072
其中,参数λ与
Figure FDA0003117367270000073
分别为方程
Figure FDA0003117367270000074
以及方程
Figure FDA0003117367270000075
的特解;
参数
Figure FDA0003117367270000076
参数
Figure FDA0003117367270000077
l1与l2为李普希兹条件中的李普希兹常数;L为拉普拉斯矩阵;In为n阶单位矩阵;INn为Nn阶单位矩阵;sup表示上确界。
9.一种时变时滞非线性多智能体的分布式脉冲一致性控制系统,其特征在于,包括:
时变时滞非线性模型构建模块,所述时变时滞非线性模型构建模块用于构建含有多个独立智能体的时变时滞非线性多智能模型,基于时变时滞非线性多智能模型获得目标状态模型;
脉冲牵制控制器构建模块,所述脉冲牵制控制器构建模块用于基于脉冲牵制控制方法构建分布式脉冲牵制控制器;
误差系统模型构建模块,所述误差系统模型构建模块用于根据智能体的时变时滞非线性多智能模型和分布式脉冲牵制控制器,定义误差向量并建立对应的误差系统模型;
函数构建模块,所述函数构建模块基于误差系统模型构建李雅普诺夫函数并证明其导数值有上界;
关系求取模块,所述关系求取模块用于求得脉冲发生时刻李雅普诺夫函数值的关系;
系统一致性模块,所述系统一致性模块根据李雅普诺夫函数的导数值有上界和脉冲发生时刻李雅普诺夫函数值的关系,利用脉冲比较原理和参数变分方法进行上界分析,得出带有时变时滞的非线性多智能体系统能够实现一致性。
10.一种时变时滞非线性多智能系统,其特征在于,包括权利要求9所述的时变时滞非线性多智能体的分布式脉冲一致性控制系统。
CN202110667285.4A 2021-06-16 2021-06-16 时变时滞多智能体的分布式脉冲一致性控制方法及系统 Active CN113325719B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110667285.4A CN113325719B (zh) 2021-06-16 2021-06-16 时变时滞多智能体的分布式脉冲一致性控制方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110667285.4A CN113325719B (zh) 2021-06-16 2021-06-16 时变时滞多智能体的分布式脉冲一致性控制方法及系统

Publications (2)

Publication Number Publication Date
CN113325719A true CN113325719A (zh) 2021-08-31
CN113325719B CN113325719B (zh) 2022-04-01

Family

ID=77421209

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110667285.4A Active CN113325719B (zh) 2021-06-16 2021-06-16 时变时滞多智能体的分布式脉冲一致性控制方法及系统

Country Status (1)

Country Link
CN (1) CN113325719B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114415517A (zh) * 2022-01-20 2022-04-29 北京航空航天大学 一种基于时间同步稳定的多智能体一致性控制方法
CN114721269A (zh) * 2022-04-11 2022-07-08 江南大学 基于脉冲窗的受扰非线性多智能体准一致性方法及系统
CN114967460A (zh) * 2022-05-31 2022-08-30 江南大学 时滞非线性多智能体系统的分布式协同控制方法
CN115562037A (zh) * 2022-11-01 2023-01-03 江南大学 一种非线性多智能体系统控制方法、装置、设备及应用

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109031958A (zh) * 2018-10-16 2018-12-18 廊坊师范学院 分数阶多智能体追踪一致性的迭代学习控制方法
CN109683477A (zh) * 2018-12-12 2019-04-26 中国地质大学(武汉) 一种随机多智能体系统有限时间控制器的设计方法及系统
CN111221311A (zh) * 2020-02-26 2020-06-02 江南大学 基于参数变分法的复杂网络分布式脉冲同步方法及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109031958A (zh) * 2018-10-16 2018-12-18 廊坊师范学院 分数阶多智能体追踪一致性的迭代学习控制方法
CN109683477A (zh) * 2018-12-12 2019-04-26 中国地质大学(武汉) 一种随机多智能体系统有限时间控制器的设计方法及系统
CN111221311A (zh) * 2020-02-26 2020-06-02 江南大学 基于参数变分法的复杂网络分布式脉冲同步方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨坤: "时变时滞连续切换奇异系统的一致有限时间稳定分析", 《控制与决策》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114415517A (zh) * 2022-01-20 2022-04-29 北京航空航天大学 一种基于时间同步稳定的多智能体一致性控制方法
CN114415517B (zh) * 2022-01-20 2024-02-02 北京航空航天大学 一种基于时间同步稳定的多智能体一致性控制方法
CN114721269A (zh) * 2022-04-11 2022-07-08 江南大学 基于脉冲窗的受扰非线性多智能体准一致性方法及系统
CN114967460A (zh) * 2022-05-31 2022-08-30 江南大学 时滞非线性多智能体系统的分布式协同控制方法
CN115562037A (zh) * 2022-11-01 2023-01-03 江南大学 一种非线性多智能体系统控制方法、装置、设备及应用
CN115562037B (zh) * 2022-11-01 2023-04-25 江南大学 一种非线性多智能体系统控制方法、装置、设备及应用

Also Published As

Publication number Publication date
CN113325719B (zh) 2022-04-01

Similar Documents

Publication Publication Date Title
CN113325719B (zh) 时变时滞多智能体的分布式脉冲一致性控制方法及系统
Liang et al. A deep reinforcement learning network for traffic light cycle control
CN112884136B (zh) 耦合神经网络有界聚类投影同步调节控制方法及系统
CN114967460B (zh) 时滞非线性多智能体系统的分布式协同控制方法
CN111523648A (zh) 含有聚类拓扑耦合的神经网络脉冲同步方法及系统
CN115860096A (zh) 一种混合时变时滞的惯性神经网络的指数同步控制方法
CN113485124B (zh) 一种考虑通信时滞的异构车辆队列稳定性控制方法及系统
CN114237041A (zh) 一种基于预设性能的空-地协同固定时间容错控制方法
CN115562037B (zh) 一种非线性多智能体系统控制方法、装置、设备及应用
CN111638648A (zh) 一种具有比例延迟复杂动态网络的分布式脉冲准同步方法
CN117055605A (zh) 多无人机姿态控制方法及系统
CN116604532A (zh) 一种上肢康复机器人智能控制方法
CN112757299B (zh) 质心轨迹的获取方法、装置、机器人及存储介质
CN110610140A (zh) 人脸识别模型的训练方法、装置、设备及可读存储介质
CN111694272B (zh) 基于模糊逻辑系统的非线性多智能体的自适应控制方法及装置
CN111221311B (zh) 基于参数变分法的复杂网络分布式脉冲同步方法及系统
CN113359711A (zh) 一种具有未知信息的多智能车系统分布式自触发控制方法
CN114721269B (zh) 基于脉冲窗的受扰非线性多智能体准一致性方法及系统
Zhu et al. Simultaneous stability of large-scale systems via distributed control network with partial information exchange
Yue et al. Distributed Cooperative Control and Communication for Multi-agent Systems
Engelmann Distributed Optimization with Application to Power Systems and Control
CN106647411A (zh) 信息物理融合系统CPS的Agent构造方法及装置
Zheng et al. Event-triggered distributed optimization for model-free multi-agent systems
CN109921924A (zh) 一种智能汽车多层立体网联控制系统建模方法
CN111814333B (zh) 奇异Lur’e网络聚类同步的牵制节点选择方法

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
GR01 Patent grant
GR01 Patent grant