CN113359439B - 多智能体系统有限时间多目标协同追踪方法和存储介质 - Google Patents

Abstract

本发明公开了多智能体系统有限时间多目标协同追踪方法和存储介质。本发明对N个多智能体进行运动学和动力学建模，并把这N个智能体分成M个子群，每个子群中都有一个相应的虚拟目标，设计多智能体系统的通讯拓扑图，由此可以设计出相应的有向图，在满足上一步骤所述有向图的前提下，设计多智能体系统的分层协同控制器，结合第一步的运动学和动力学模型，实现多智能体系统的有限时间多目标协同追踪控制。本发明设计的估计器算法可确保每个智能体都能在有限时间内获得其相应目标的准确信息，设计的本地层算法可以实现智能体的一对一有限时间局部跟踪。

Description

多智能体系统有限时间多目标协同追踪方法和存储介质

技术领域

本发明涉及分布式人工智能领域，特别是涉及多智能体系统有限时间多目标协同追踪方法和存储介质。

背景技术

多智能体系统是20世纪末至21世纪初分布式人工智能领域的国际前沿研究课题，其核心支撑理论是人工智能、分布式控制和分布式计算。进入21世纪，人们在解决大型、复杂的工程问题时，发现单个智能体的能力已经无法胜任，需要多个智能体在网络环境下以信息通讯的方式组成多智能体系统协同地解决工程问题。典型的多智能体系统包括多智能体系统，多无人机系统，智能电网和分布式卫星系统等。

但是，目前大多数的多智能体系统的控制方法仅适用于渐进时间内控制单个目标，且局限于研究用单积分器、双积分器动力学模型描述的智能体系统，并忽略了外部扰动、参数不确定性的影响。而在实际应用中，外部扰动、参数不确定性是不可避免的；同时，可能会要求多智能体系统通过目标追踪同时完成多个任务，且运用欧拉-拉格朗日动力学来描述的多智能体系统能更严谨地描述现实生活中的应用，例如：多个机械臂，多个无人智能小车，多个无人舰艇等；另外，很多实际的应用中对完成任务的时间会有需求，从而需要控制算法来实现更短的时间内执行任务。

因此，在考虑外部扰动及动力学模型参数不确定性的多智能体系统的一种多智能体系统有限时间多目标协同追踪方法具有重要意义。

发明内容

本发明要解决的技术问题在于，针对现有技术的无法满足有限时间多目标协同追踪技术缺陷，提供了多智能体系统有限时间多目标协同追踪方法和存储介质。

根据本发明的一个方面，本发明的提供一种多智能体系统有限时间多目标协同追踪方法，包括：

对多智能体系统中N个智能体进行运动学和动力学建模；

将N个智能体分成M个子群，每个子群中确定一个虚拟跟踪目标，其中N和M均为大于等于1的整数；

对智能体之间的交互作用进行建模，获得所述子群的通讯有向拓扑图；

根据所述通讯有向拓扑图，设计基于估计器的多智能体系统的分层协同控制器；

根据所述运动学和动力学模型，利用所述分层协同控制器，实现每个子群的智能体在有限时间收敛至所述每个子群的智能体对应的虚拟追踪目标的轨迹。

进一步地，所述对多智能体系统中N个智能体进行运动学和动力学建模包括：

运动学和动力学模型表达式为：

其中，

以及

分别表示关节的位置、速度和加速度；

惯性矩阵为：

科氏离心矩阵为：

重力转矩为：

其中，

表示输入力矩，d_i(t)是外部扰动且满足

为一个已知正常数，

以及

是动力学不确定项。

进一步地，所述模型表达式还可以表达为：

其中，

进一步地，当第i个智能体属于第l个子群时，所述跟踪目标的数学表达式为：

其中，q_o,l(t)、ω_o,l(t)和

分别表示第l个子群的追踪目标的位置状态、速度和加速度，l∈{1,2,...,M}表示追踪目标所在子群的序号。

进一步地，所述获得所述子群的通讯有向拓扑图包括：

所述通信拓扑图

其中，

为N个智能体的点集，

为有向边，

为邻接矩阵，且有向边e_ij∈ε，即第i个智能体可以接收第j个智能体的信息；

如果第i个智能体和第j个智能体是竞争关系，则a_ij＜0，如果第i个智能

体和第j个智能体是合作关系，则a_ij＞0，没有自循环时，a_ij＝0，且

根据所述邻接矩阵，确定所述有向图的拉普拉斯矩阵

根据所述拉普拉斯矩阵，所述有向图

定义为

且l_ij＝-a_ij，

进一步地，所述将N个智能体分成M个子群包括：

多智能体系统被划分为M(M＞1)个子群，每一个所述子群被描述为含有n_l个智能体的子图：

其中，

且

进一步地，所述获得所述子群的通讯有向拓扑图还包括：

根据智能体与其对应跟踪目标之间的通讯状况，确定有向图的对角权重矩阵B：

其中，N＝(1,2,...,i),若第i个智能体能够直接接收其跟踪目标的信息，则bi为取值大于0的数，如果不能，则bi＝0。

进一步地，所述设计基于估计器的多智能体系统的分层协同控制器还包括：

所述分层协同控制器的数学表达式为：

其中，i∈V_l表示第l个子群里的第i智能体，集合

是用来标记第l个子群中的所有智能体，

表示N个智能体的集合，l∈{1,...,M}表示M个子群；τ_i1和τ_i2分别表示等效控制律和滑模控制律，λ₁,λ₂＞0，K_p和K_d是两个正定对角矩阵，

0＜α₁＜1，α₂＝2α₁/(α₁+1)，

和

分别表示ε_1,i和ε_2,i的估计值，并且分别定义为：

其中，η＞0，β_i＞0，

估计层中，

χ＞0，1＜p＜2，

κ₁＞0，κ₂＞0，

为

的第u个元素α₂＝2α₁/(α₁+1)，μ_γ(x)是一个非线性函数，定义为：

其中，

γ＞0。

进一步地，所述利用所述分层协同控制器，实现每个子群的智能体在有限时间收敛至所述每个子群的智能体对应的虚拟追踪目标的轨迹还包括：

结合估计器的误差形式：

将所述分层协同控制器的数学表达式转化为：

对闭环系统Ξ₁、Ξ₂、Ξ₃分别构造李雅普诺夫函数

和

1＜r＜2；

结合所述闭环系统和稳定性分析，得到多目标协同追踪的误差满足以下条件：

对于所有的

l∈{1,...,M}，其中，e_i,l＝q_i-q_o,l，

T_g为到达时间。

根据本发明的另一个方面，还提供一种存储介质，其特征在于，其上存储有如前任一所述的一种多智能体系统有限时间多目标协同追踪方法。

在本发明所述的一种多智能体系统的有限时间多目标协同追踪控制方法中，提供一种考虑扰动及模型参数不确定性的多智能体系统的有限时间多目标协同追踪控制方法，通过选择合适的控制参数来提高多智能体系统的有限时间多目标协同追踪的速度和准确性。

总体而言，通过本发明提供的一种多智能体系统的有限时间多目标协同追踪控制方法，具有以下有益效果：

1、与传统的集中控制方式比起来，本项目采用的分布式控制方法具有更强的鲁棒性，更少的能耗和更高的执行效率；

2、考虑智能体模型中的动力学参数不确定性以及外部扰动，使得多目标协同追踪控制实用性更强；

3、与控制单一目标的系统来说，本项目是多目标控制，能够在需要同时处理多个任务的网络系统中发挥作用；

4、与渐进收敛时间相比，本项目所设计的控制算法的优势在于能够保证多智能体在有限时间内追踪到相应的目标，具有更好的跟踪性能。

5、本项目采用的分布式控制方法只要求

矩阵是满秩的，降低了对

矩阵的条件。

附图说明

下面将结合附图及实施例对本发明作进一步说明，附图中：

图1是本发明实施例中的多智能体系统的系统流程图；

图2是本发明实施例中的多智能体系统的通讯有向拓扑图；

图3是本发明实施例中的三个子群中位置状态轨迹追踪图和位置状态估计量轨迹追踪图；

图4是本发明实施例中的三个子群中中速度状态轨迹追踪图和速度状态估计量轨迹追踪图；

图5是本发明实施例中的三个子群中位置估计误差图、局部位置追踪误差图以及在总体位置追踪误差图；

图6是本发明实施例中的三个子群中总体速度追踪误差图；

具体实施方式

为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。

实施例一，如图1所述，为多智能体系统的系统流程图，具体包括以下步骤：

步骤S1：选定13个智能体与3个虚拟追踪目标，进行动力学和运动学建模；参考图2，其为多智能体系统的通讯有向拓扑图，其中，点1-13表示智能体，点L1-L3表示其对应的虚拟追踪目标；将编号1-3的智能体定义为第一个子群，编号4-8的智能体定义为第二个子群，编号9-13的智能体定义为第三个子群；参考图3,所述的多智能体具有不同的物理参数，具体见表1以及表2；

表1.多智能体系统中的物理参数

表2.多智能体系统中的物理参数的估计值：

步骤S11：对N个智能体进行动力学和运动学建模，模型表达式为：

其中，

以及

分别表示关节的位置、速度和加速度，惯性矩阵：

科氏离心矩阵：

重力转矩：

表示输入力矩，d_i(t)是外部扰动且满足

为一个已知正常数，

以及

是动力学不确定项。

在一些实施例中，建好的动力学和运动学模型可以改写为:

其中，

同时把这13个智能体分成3个子群，3个子群中的虚拟追踪目标分别对应为子群L1、L2、L3。

步骤S12、根据有限时间多目标协同追踪控制方法，步骤S1中追踪的多目标的数学表达式为：

其中，q_o,l(t)、ω_o,l(t)和

分别表示第l个子群的追踪目标的位置状态、速度和加速度，l∈{1,2,3}表示多个目标子群的序号，当第i个智能体属于第l个子群时，上述表达式即为当前追踪目标；其中，

步骤S2、如图2为多智能体系统的通讯有向拓扑图，首先建立智能体之间的通讯，采用有向图

对智能体之间的交互作用进行建模，其中

表示13个智能体,

表示有向边,

表示邻接矩阵；

有向边e_ij∈ε，这意味着第i个智能体可以直接接收第j个智能体信息，如果第i个智能体和第j个智能体是竞争关系；则a_ij＜0，如果是合作关系，则a_ij＞0，否则a_ij＝0，

假设没有自循环，即a_ij＝0，

此外，相应的拉普拉斯矩阵

关于有向图

定义为

且有l_ij＝-a_ij对于所有的

所设计的多智能体系统被划分为3个子群，根据前面定义的邻接矩阵

确定有向图的拉普拉斯矩阵

对于每个智能体，至少存在一条路径，使得被考虑的智能体可以接收到领导者的信息，由此确定有向图的对角权重矩阵

所述的

和B设置为：

B＝diag(7,7,0,0,8,0,0,0,0,7,9,0,8)

步骤S3、在满足有向图的前提下，设计多智能体系统的基于估计器的分布式控制器，实现每个子群的智能体在有限时间收敛至其对应的虚拟目标轨迹，完成多智能体系统的有限时间多目标协同追踪控制；其中，基于估计器的分布式控制器中的控制参数的选取具体见表3：

表3分布式控制器中的控制参数

步骤S31、通过步骤S3中已经设计好的有限时间多目标追踪控制器，对构造李雅普诺夫函数

求导进行稳定性分析证明，可以得到三个闭环系统都是有限时间内稳定的。

步骤S32、结合步骤S31得到的稳定性分析结果，从而得到多智能体系统的跟踪误差满足的条件：

当t→T_g时，q_i→q_i,1，ω_i→ω_i,1，i∈{1,2,3}；q_i→q_i,2，ω_i→ω_i,2，i∈{4,5,6,7,8}；q_i→q_i,3，ω_i→ω_i,3，i∈{9,10,11,12,13}即得到多智能体系统的有限时间多目标追踪。

如图3所示，(a)和(b)表示三个子群的多智能体位置状态q_i1和q_i2的追踪轨迹图，(c)和(d)表示三个子群的多智能体位置状态估计量

的跟踪轨迹图，在所设计的状态估计器的作用下，三个子群中的智能体的位置状态估计量都在有限时间内收敛于其对应的跟踪目标轨迹。

如图4所示，(a)和(b)表示三个子群的多智能体速度状态ω_i1，ω_i2的追踪轨迹图，(c)和(d)表示三个子群的多智能体速度状态估计量

的跟踪轨迹图，在所设计的状态估计器的作用下，三个子群中的智能体的速度状态估计量都在有限时间内收敛于其对应的追踪目标速度。

如图5所示，(a)和(b)表示三个子群中位置估计误差

图，位置估计误差都在有限时间内收敛到[-0.005,0.005]内；(c)和(d)表示三个子群的局部追踪位置误差

图，(e)和(f)表示三个子群的总体位置追踪误差e_i1，e_i2图，在所设计的控制器的作用下，三个子群中的智能体的位置状态都在有限时间内收敛到其对应的跟踪目标轨迹。

如图6所示，(a)和(b)表示的是三个子群的总体速度追踪误差

图，(c)和(d)表示的是除颤后的三个子群的总体速度追踪误差

图，表明三个子群中的智能体的速度状态都在有限时间内收敛到其对应的跟踪目标速度，且除颤后的速度状态颤动明显减小。

以上结合附图对本发明的实施方式进行了说明，但是本发明的保护范围并不局限于上述的具体实施方式，任何熟悉本技术领域的技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多变化或替换形式，这些都应该涵盖在本发明的保护范围之内。

Claims (9)

1.一种多智能体系统有限时间多目标协同追踪方法，其特征在于，包括：

对多智能体系统中N个智能体进行运动学和动力学建模；

将N个智能体分成M个子群，每个子群中确定一个虚拟跟踪目标，其中N和M均为大于等于1的整数；

对智能体之间的交互作用进行建模，获得所述子群的通讯有向拓扑图；

根据所述通讯有向拓扑图，设计基于估计器的多智能体系统的分层协同控制器；

所述分层协同控制器的数学表达式为：

其中，i∈V_l表示第l个子群里的第i智能体，集合

是用来标记第l个子群中的所有智能体，

表示N个智能体的集合，l∈{1,...,M}表示M个子群；τ_i1和τ_i2分别表示等效控制律和滑模控制律，λ₁,λ₂＞0，K_p和K_d是两个正定对角矩阵，

0＜α1＜1，α₂＝2α₁/(α₁+1)，

和

分别表示ε_1,i和ε_2,i的估计值，并且分别定义为：

其中，η＞0，β_i＞0，

估计层中，

χ＞0，1＜p＜2，

κ₁＞0，κ₂＞0，

为

的第u个元素α₂＝2α₁/(α₁+1)，μ_γ(x)是一个非线性函数，定义为：

其中，

根据所述运动学和动力学模型，利用所述分层协同控制器，实现每个子群的智能体在有限时间收敛至所述每个子群的智能体对应的虚拟追踪目标的轨迹。

2.如权利要求1所述的一种多智能体系统有限时间多目标协同追踪方法，其特征在于，所述对多智能体系统中N个智能体进行运动学和动力学建模包括：

运动学和动力学模型表达式为：

其中，q_i、

以及

分别表示关节的位置、速度和加速度；

惯性矩阵为：

科氏离心矩阵为：

重力转矩为：

其中，

表示输入力矩，d_i(t)是外部扰动且满足

为一个已知正常数，

以及

是动力学不确定项。

3.如权利要求2所述的一种多智能体系统有限时间多目标协同追踪方法，其特征在于，所述模型表达式还可以表达为：

其中，

4.如权利要求1所述的一种多智能体系统有限时间多目标协同追踪方法，其特征在于，当第i个智能体属于第l个子群时，所述虚拟跟踪目标的数学表达式为：

其中，q_o,l(t)、ω_o,l(t)和

分别表示第l个子群的追踪目标的位置状态、速度和加速度，l∈{1,2,...,M}表示追踪目标所在子群的序号。

5.如权利要求1所述的一种多智能体系统有限时间多目标协同追踪方法，其特征在于，所述获得所述子群的通讯有向拓扑图包括：

通信拓扑图

其中，

为N个智能体的点集，

为有向边，

为邻接矩阵，且有向边e_ij∈ε，即第i个智能体可以接收第j个智能体的信息；

如果第i个智能体和第j个智能体是竞争关系，则a_ij＜0，如果第i个智能体和第j个智能体是合作关系，则a_ij＞0，没有自循环时，a_ij＝0，且

根据所述邻接矩阵，确定所述有向图的拉普拉斯矩阵

根据所述拉普拉斯矩阵，有向图

定义为

且l_ij＝-a_ij，

6.如权利要求5所述的一种多智能体系统有限时间多目标协同追踪方法，其特征在于，所述将N个智能体分成M个子群包括：

多智能体系统被划分为M， M＞1， 个子群，每一个子群被描述为含有n_l个智能体的子图：

其中，

且

7.如权利要求1所述的一种多智能体系统有限时间多目标协同追踪方法，其特征在于，所述获得所述子群的通讯有向拓扑图还包括：

根据智能体与其对应跟踪目标之间的通讯状况，确定有向图的对角权重矩阵B：

其中，N＝1,2,...,i,若第i个智能体能够直接接收其跟踪目标的信息，则bi为取值大于0的数，如果不能，则bi＝0。

8.如权利要求1所述的一种多智能体系统有限时间多目标协同追踪方法，其特征在于，所述利用所述分层协同控制器，实现每个子群的智能体在有限时间收敛至所述每个子群的智能体对应的虚拟追踪目标的轨迹还包括：

结合估计器的误差形式：

和

将所述分层协同控制器的数学表达式转化为：

对闭环系统Ξ₁、Ξ₂、Ξ₃分别构造李雅普诺夫函数

和

1＜r＜2,进行稳定性分析；

结合所述闭环系统和稳定性分析结果，得到多目标协同追踪的误差满足以下条件：

对于所有的

l∈{1,...,M}，其中，e_i,l＝q_i-q_o,l，

T_g为到达时间。

9.一种存储介质，其特征在于，其上存储有如权利要求1-8任一所述的一种多智能体系统有限时间多目标协同追踪方法。

Images