III. BACKGROUND
在本节中,我们简要回顾了我们在整个本文中使用的有关图论和正多边形的许多重要概念。
A. Graph Theory
一个(无向)图G包含顶点V和边ε,其中边是G中无序的不同顶点对。如果x,y∈V并且(x,y)∈ε,这是可以称x,y相邻,或者为邻居,并可以使用符号x∼y表示。每个顶点的邻居的个数数即为其度数。 从顶点x到顶点y的长度为r的路径是一系列r+1个不同的顶点,这些顶点以x开头,以y结束,依次连接相邻的顶点。 如果图G的任意两个顶点之间有一条路径,则称G已连接。
(无向)图G的邻接矩阵A(G)=[aij]是一个对称矩阵,其中行和列由G的顶点索引,因此如果顶点i和顶点j是邻居,则aij=1,否则aij=0。 我们还假设所有i的aii=0。 图G的度矩阵D(G)是对角矩阵,其中行和列由\mathcal{V}索引,其中(i,i)-项是顶点i的度。
定义为的对称奇异矩阵
L(G)=D(G)−A(G)
称为G的拉普拉斯算子。拉普拉斯算子矩阵可以表征图的许多拓扑特性。 拉普拉斯算子L是一个半正定矩阵,其零特征值(即其核的维数)的代数多重性等于图中连通的分量数。 与零特征值关联的n维特征向量是1的向量,1n=[1,...,1]T 。 有关图论的更多信息,请参见[13]。
B. Regular Polygons
令d<n为正整数,并定义p=n/d。 令y1为单位圆上的一个点。 令Rα为顺时针旋转角度α=2π/p。 广义正则多边形可以表示为p由点yi+1=Rαyi组成和或由点i与i+1之间的边组成。
当d=1时,多边形p被称为普通规则多边形,并且其边缘不相交。 如果d>1,并且n和d为互质,则边线相交,并且多边形为星形。 如果n和d的公共因子l>1,则该多边形由具有n/l个顶点和边缘的同一个多边形的l个重复图形组成。 如果d=n,则多边形n/n对应于同一位置的所有点。 如果d=n/2(n个偶数),则多边形由两个端点和它们之间的一条线组成,点的一端具有偶数索引,而另一端具有奇数索引。 有关正则图的更多信息,请参见[7]。
IV. PROBLEM STATEMENT
考虑一组n个单位速度的平面(2维空间中的)代理(agent)。 每个代理都能够感知来自其邻居的信息。 代理i的邻居集,即Ni,是代理i可以“看到”的代理集。 “看见”的确切含义将在后面阐明。 邻域的大小取决于传感器的特性。 代理之间的相邻关系可以通过连接图G=(V,E,W)描述。
定义1(连通图):连通图G=(V,E,W)是由以下几部分构成
1)由一组移动代理组成的一组顶点集合(V;
2)一组边ε=((i,j)∣i,j∈(V,andi〜j;
3)一组均为正值的,与每条边(i,j)相对应的边权重集合。
代理i的邻居可以被定义为:
Ni≐{j∣i∼j}⊆V\{i}
令ri代表代理i的位置,令vi为其速度矢量。 每个单位速度代理的运动学公式为:
r˙i=vi
v˙i=ωivi⊥
v˙i⊥=−ωivi
其中vi⊥是垂直于速度矢量vi的单位矢量(见图1)。 正交对{vi,vi⊥}形成了主体i的主体框架。我们用v=[v1T,...,vnT]T∈R2n×1表示所有速度的叠加矢量。
每个代理的控制输入是角速度ωi。由于假定代理以恒定的单位速度移动,因此施加到每个代理的力必须垂直于其速度矢量,即,作用在每个代理上的力是陀螺力,并且它不会改变其速度(因此,它的动能)。因此,ωi充当每个代理的转向控制[16]。
让我们正式定义我们要考虑的编队形式。
定义2(平行形式):所有代理的航向相同且速度矢量对齐的配置称为平行形式。
请注意,在此定义中,我们不考虑达成一致的速度的价值,而仅考虑达成一致的事实。在平衡状态下,代理的相对距离决定了队形的形状。另一个有趣的编队队形是平衡的圆形编队。
定义3(平衡圆形结构):代理沿着相同的圆形轨迹移动并且代理的几何中心固定的配置称为平衡圆形结构。这种形状可以由适当的规则多边形表示。
在以下各节中,我们将研究每个编队并设计其相应的分布式控制法则。
V. PARALLEL FORMATIONS
我们在本节中的目标是为每个代理设计一个控制法则,以使移动代理的航向达到一致,即,它们的速度矢量对齐,从而形成类似蜂群的模式。 对于任意连通性图G,请考虑拉普拉斯矩阵L。因此,我们按以下方法定义未对准的度量[27],[35]:
ω(v)=21i ~j∑∣vi−vj∣2=21⟨v,Lˉv⟩(2)
其中累加的是所有的(i,j)∈ω,Lˉ=L⨂I2∈R2n×2n,I2是2x2的单位矩阵。则ω(v)关于时间的导数为:
ω˙(v)=i=1∑n⟨v˙i,(Lˉv)i⟩=i=1∑nω˙i⟨v˙i⊥,(Lˉv)i⟩
其中(Lˉv)i∈R2是与第i个代理相关联的Lˉv的子向量。 因此,以下面的梯度下降控制规则保证了ω(v)单调递减:
ωi=k⟨vi⊥,(Lˉv)i⟩=−kj∈Ni∑⟨vi⊥,vij⟩
其中k<0是已知的,vij=vi−vj。
备注1: 令θi表示在固定世界座标系中测得的代理i的航向(见图1)。 单位速度矢量vi及其正交矢量vi⊥由vi=[cosθi,sinθi]T和vi⊥=[−sinθi,cosθi]T给出。 因此,控制输入(3)变为
ωi=kj∈Ni∑sin(θi−θj)k<0(4)
图1.每个代理的轨迹由平面Frenet框架表示。
值得注意的是,所提出的控制器是用于耦合非线性振荡器的Kuramoto模型控制器,该控制器已在数学物理学以及控制界中进行了广泛的研究[15],[19],[36]。相同的模型也已经用于循环机器人系统的相位调节[18]。
我们有以下定理[27],提供了获得平行形式的充分条件。
定理1: 考虑一个具有动力学的n个单位速度智能体的系统(1)。如果基础连通性图保持固定和连接,则通过应用控制输入(4),系统收敛到ω=[ω1,⋅⋅⋅,ωn]T=0的平衡点。此外,如果θi∈(−π/2,π/2),则速度共识集在局部有吸引力。 。
证明1:证明[27]。
速度共识集是所有主体具有相同速度向量的状态集,它对应于定义2中定义的平行形式的队形。请注意,θi∈(−π/2,π/2)∀i=1,...,n是将初始航向限制为半圆的充分条件。结果可以扩展到switching topology,,如[27]所示。