%% 无中心领导的一致性算法(按照论文中的图模型计算)topology=[1,2;2,1;1,3;3,2;3,4;4,3;3,7;7,3;4,5;5,6;6,5;7,8;8,7;9,7;6,9;9,6;9,10;10,9];Link=zeros(10);for i=1:18 Link(topology(i,1),topology(i,2))=1;endw=zeros(10);%一致性系数v=zeros(10);%反馈项系数for i=1:10 for j=1:10 if i==j w(i,j)=0.5; v(i,j)=0.5; end if Link(j,i)==1;%j是i的父节点 w(i,j)=0.5/sum(Link(:,i)); end if Link(i,j)==1%i是j的父节点 v(i,j)=0.5/sum(Link(i,:)); end endendE_last=zeros(10,1);%上一轮次的共享计划xigma=zeros(10,1);%反馈项lagrant_data=[lagrant];E_data=[];ticfor t=1:1000 D1=max(min((a1-lagrant)./b1,Dmax),0); S1=max(min((lagrant-a2)./b2,Smax),0); E=D1-S1; lagrant=w*lagrant+0.0001*xigma;%一致性变量的更新