本发明专利技术公开了基于子图处理的多智能体系统一致性问题的优化方法,其特征在于,包括如下步骤:1)构建多智能体系统的图信号模型;2)对待求解问题加上辅助约束条件;3)子图内局部求逆;4)子图间融合求平均;5)迭代消除误差。这种方法能降低运算复杂度、实现一致平均,且收敛速度快,对系统单个节点的计算能力、存储能力要求低。
The optimization method of unified problem of multi-intelligent system based on subgraph processing
【技术实现步骤摘要】
基于子图处理的多智能体系统一致性问题的优化方法
本专利技术涉及多智能体系统
,具体涉及基于子图处理的多智能体系统一致性问题的优化方法。
技术介绍
多智能体系统广泛应用在飞行编队、机器人编队、交通控制系统、传感器网络等诸多领域。多智能体系统中的单个智能体具有有限的传感、通信、计算能力,而整个多智能体系统利用智能体之间的协作来完成高复杂度的任务。多智能体系统作为分布式信息处理的一个载体,在计算机科学、通信工程、生物智能、自动控制等领域均有广泛研究。在多智能体系统的协同控制问题中,达到所需的一致性是智能体进行协调合作的前提条件,因此一致性问题受到更多关注。一致性是指每个智能体与其邻居交换彼此的状态信息,不断协同,最终使得所有智能体的某个状态值达到一致,与一致性问题相对应的算法称为一致性算法。一致性问题中,以智能体初始状态值的平均值为一致性收敛目标的问题称为一致平均问题,其应用较为广泛,如直流微电网作为一种新型的能源网络化供应与管理结构,是一个包含分布式电源、负荷及储能装置的微型系统,在微电网中需要维持母线电压的稳定,文献提出基于离散一致性的微电网自适应下垂控制策略,以全网平均电压差为一致性优化目标,来实现高精度的电压调节。一致性问题被提出以后,在此基础上形成了分布式计算理论,众多学者提出了分布式优化算法。Xiao和Boyd首先提出确定性Gossip算法,该算法中,同一时刻下每个节点都与其邻居节点交互信息,然后计算自己和所有邻居间的加权平均值,当加权系数所构成的权重矩阵是一致平均矩阵时,所有节点的信息会收敛到平均一致;Nedic和Ozdaglar借鉴确定性Gossip算法的思想,提出分布式梯度下降法,该算法令每个节点对自己及邻居节点的信息求加权平均,同时采用梯度下降法极小化本地代价函数来加快一致性收敛速度,此类基于梯度的分布式优化算法实现简单,但收敛速度较慢。为此,有学者提出了基于对偶理论的分布式算法,如Terelius等人提出一种分布式交替方向乘子法(AlternatingDirectionMethodofMultipliers,ADMM),ADMM主要思想是引入对偶变量,将原优化问题分解为两个子问题,每个节点需要更新自身状态并计算对偶变量的更新,相较于一般的梯度算法,分布式ADMM算法提高了收敛速度,但对节点的存储能力、计算能力要求更高。
技术实现思路
本专利技术的目的是针对现有技术的不足,而提供基于子图处理的多智能体系统一致性问题的优化方法。这种方法能降低运算复杂度、实现一致平均,且收敛速度快,对系统单个节点的计算能力、存储能力要求低。实现本专利技术目的的技术方案是:基于子图处理的多智能体系统一致性问题的优化方法,包括如下步骤:1)构建多智能体系统的图信号模型:依据多智能体系统中各智能体之间的通信情况作图信号模型的拓扑结构,智能体的观测数据作为图信号模型中的节点信号,建立多智能体系统的图信号模型G=(V,E,A),其中V={1,2,…,N}表示节点的集合;E={eij}表示边的集合,eij=(i,j)表示节点i和节点j之间有边相连接,A表示全图的邻接矩阵,邻接矩阵A={aij}∈RN×N,当节点i与j之间有边存在时aij=aji=1,定义图的度矩阵D=diag(di),其中图的拉普拉斯矩阵定义为公式(1):L=D-A(1),其特征分解为L=UΛUT,其中特征向量矩阵U=[u1,u2,…,uN],特征值矩阵Λ=diag{λ1,λ2,…,λN},节点i处的信号值为x(i),所有节点的信号值集合为一个列向量x=[x(1),x(2),…,x(N)]T,在研究多智能体系统一致性问题中,一致平均矩阵有重要作用,构造一致平均矩阵为公式(2):其中I为单位矩阵,λmax为图的拉普拉斯矩阵L的最大特征值;2)对待求解问题加上辅助约束条件:一致性优化问题为公式(3):其中f(xt)表示多智能体系统的全局代价函数,fi(xt)是智能体i所观测到的局部代价函数,且其它节点无法获知该代价函数,xt是全局代价函数f(xt)的自变量,当多智能体系统实现一致平均时,该问题为如公式(4)所示的最小二乘问题:其中x0为智能体初始信号,约束条件xt(i)=xt(j),i,j∈V等价于Wxt=xt,因此问题(4)等价于公式(5):又等价于公式(6):其中,α为常数,将等式(6)的约束条件惩罚至目标函数上,则有:将一致平均矩阵代入公式(7),则:对公式(8)求解梯度向量为:令该梯度向量为零,得到该优化问题的最优解:3)子图内局部求逆:节点i从矩阵αI+εL中取出其2r阶邻居对应位置的值组成小矩阵,对小矩阵进行求逆并保留节点r阶邻居对应位置的值,过程为:首先构建子图,令每个节点及其邻居节点所在的区域作为一个子图,将建立的图模型划分为N个子图,对于节点i,取其子图内节点在大矩阵αI+εL中对应位置的值组成小矩阵再进行求逆,其本地优化变量由公式(11)计算:其中表示局部操作、定义为取图中节点i的r阶内邻居节点对应位置的值的局部操作,x(k)为第k步迭代的取值见公式(15),在第一步迭代时,有x(0)=x0,上标表示求矩阵的伪逆、可以等价转化为小规模矩阵的求逆;4)子图间融合求平均:对节点i而言,节点i子图内所有节点j,包括节点i都会计算出本地优化变量在i位置上的值,把所有节点j所计算的求逆结果在i位置的值进行融合平均,然后赋给节点i:进行融合求平均的计算:其中i∈M,mi,r表示节点i的r阶邻居内所有节点个数,vj(i)表示vj在i对应位置的值;5)迭代消除误差:步骤1)-步骤4)将大规模矩阵的求逆问题近似为一系列小规模矩阵的求逆,假设为矩阵αI+εL求逆的近似结果,即但局部求出的解并非全局最优解,存在近似误差,因此采用迭代更新的方式使局部近似解逼近全局最优解,由公式(10),考虑目标函数其梯度矩阵为海森矩阵为采用牛顿法进行求解,步长为单位步长,迭代公式为:将代入公式(13),则:迭代过程等效为:迭代过程可分布式进行,节点i对待求解优化变量xt(i)进行更新对x(k)(i)进行更新:其中di为图的度矩阵中第i行非零元素,aij为图的邻接矩阵第i行第j列的元素,x(k)(i)更新后,返回公式(11)、公式(12)、公式(16)继续迭代,当时,σ为迭代终止条件,迭代终止,输出与现有技术相比,本技术方案基于子图处理的思想,把大规模矩阵的求逆问题近似分解为一系列子图优化问题,每个节点在大矩阵中取出子图内所有节点对应位置的值组成小矩阵,对小矩阵进行求逆并保留节点一阶邻居内所有节点对应位置的值进行融合求平均,极大地降低运算复杂度。这种方法能降低运算复杂度、实现一致平均,且收敛速度快,对系统单个节点的计算能力、存储能力要本文档来自技高网...
【技术保护点】
1.基于子图处理的多智能体系统一致性问题的优化方法,其特征在于,包括如下步骤:/n1)构建多智能体系统的图信号模型:依据多智能体系统中各智能体之间的通信情况作图信号模型的拓扑结构,智能体的观测数据作为图信号模型中的节点信号,建立多智能体系统的图信号模型G=(V,E,A),其中V={1,2,…,N}表示节点的集合;E={e
【技术特征摘要】
1.基于子图处理的多智能体系统一致性问题的优化方法,其特征在于,包括如下步骤:
1)构建多智能体系统的图信号模型:依据多智能体系统中各智能体之间的通信情况作图信号模型的拓扑结构,智能体的观测数据作为图信号模型中的节点信号,建立多智能体系统的图信号模型G=(V,E,A),其中V={1,2,…,N}表示节点的集合;E={eij}表示边的集合,eij=(i,j)表示节点i和节点j之间有边相连接,A表示全图的邻接矩阵,邻接矩阵A={aij}∈RN×N,当节点i与j之间有边存在时aij=aji=1,定义图的度矩阵D=diag(di),其中图的拉普拉斯矩阵定义为公式(1):
L=D-A(1),
其特征分解为L=UΛUT,其中特征向量矩阵U=[u1,u2,…,uN],特征值矩阵Λ=diag{λ1,λ2,…,λN},节点i处的信号值为x(i),所有节点的信号值集合为一个列向量x=[x(1),x(2),…,x(N)]T,
构造一致平均矩阵为公式(2):
其中I为单位矩阵,λmax为图的拉普拉斯矩阵L的最大特征值;
2)对待求解问题加上辅助约束条件:一致性优化问题为公式(3):
其中f(xt)表示多智能体系统的全局代价函数,fi(xt)是智能体i所观测到的局部代价函数,且其它节点无法获知该代价函数,xt是全局代价函数f(xt)的自变量,
当多智能体系统实现一致平均时,该问题为如公式(4)所示的最小二乘问题:
其中x0为智能体初始信号,约束条件xt(i)=xt(j),i,j∈V等价于Wxt=xt,因此问题(4)等价于公式(5):
又等价于公式(6):
其中,α为常数,将公式(6)的约束条件惩罚至目标函数上,则有:
将一致平均矩阵W=I-εL代入公式(7),则:
对公式(8)求...
【专利技术属性】
技术研发人员:蒋俊正,李龙斌,李杨剑,池源,冯海荣,卢军志,黄炟鑫,
申请(专利权)人:桂林电子科技大学,
类型:发明
国别省市:广西;45
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。