多智能体一致性 matlab 仿真.zip

clear all clc tspan=[0,40]; y0=[7 6 2 1 0.5 -0.3 -0.5 -3.5 -5 -7.2]'; lags=[0.268 0.134 0.188 0.3]; sol=dde23('System_timedelays',lags,y0,tspan); %yy1=yy.^2; %将矩阵的每一个元素平方 %kase=sum(yy1,2);%对每行元素求和 L=[3 -1 0 0 0 0 0 0 -1 -1; -1 3 -1 0 0 0 0 0 -1 0; 0 -1 3 -1 0 0 0 -1 0 0; 0 0 -1 3 -1 0 -1 0 0 0; 0 0 0 -1 3 -1 -1 0 0 0; 0 0 0 0 -1 2 -1 0 0 0; 0 0 0 -1 -1 -1 4 -1 0 0; 0 0 -1 0 0 0 -1 3 -1 0; -1 -1 0 0 0 0 0 -1 4 -1; -1 0 0 0 0 0 0 0 -1 2] %计算第二小特征值是注意要改变系统的状态方程 landa=eig(L); landa2=real(landa(2)); %subplot(2,1,1); figure(1); plot(sol.x,sol.y(1,:)); hold on; plot(sol.x,sol.y(2,:)); hold on; plot(sol.x,sol.y(3,:)); hold on; plot(sol.x,sol.y(4,:)); hold on; plot(sol.x,sol.y(5,:)); hold on; plot(sol.x,sol.y(6,:)); hold on; plot(sol.x,sol.y(7,:)); hold on; plot(sol.x,sol.y(8,:)); hold on; plot(sol.x,sol.y(9,:)); hold on; plot(sol.x,sol.y(10,:)); axis([0 50 -10 10]); xlabel('Time(sec)'); ylabel('State'); grid; title('Time-Delay=0.268'); % subplot(2,1,2); % plot(tt,kase); % axis([0 40 0 200]); % xlabel('time(sec)'); % ylabel('disagreement'); % grid;