多智能体系统一致性程序

多智能体系统一致性的程序，解微分方程

function h1
a11=0;a12=1;a13=1;a14=0;a15=1;
a21=1;a22=0;a23=0;a24=0;a25=0;
a31=0;a32=1;a33=0;a34=1;a35=0;
a41=0;a42=0;a43=0;a44=0;a45=1;
a51=1;a52=1;a53=1;a54=1;a55=0;
% a11=0;a12=0;a13=0;a14=0;a15=1;
% a21=1;a22=0;a23=0;a24=0;a25=0;
% a31=0;a32=1;a33=0;a34=1;a35=0;
% a41=0;a42=0;a43=0;a44=0;a45=1;
% a51=1;a52=0;a53=0;a54=0;a55=0;
b11=0;b12=1;b13=1;b14=0;b15=1;
b21=1;b22=0;b23=0;b24=0;b25=0;
b31=0;b32=1;b33=0;b34=1;b35=0;
b41=0;b42=0;b43=0;b44=0;b45=1;
b51=1;b52=1;b53=1;b54=1;b55=0;
% b11=0;b12=0;b13=0;b14=0;b15=1;
% b21=1;b22=0;b23=0;b24=0;b25=0;
% b31=0;b32=1;b33=0;b34=1;b35=0;
% b41=0;b42=0;b43=0;b44=0;b45=1;
% b51=1;b52=0;b53=0;b54=0;b55=0;
k1=1;
k2=1;
t0=0:0.001:10;
x0=[10;8;6;4;2;9;7;5;3;1];
[tx]=ode45（@consensust0x0）;
figure（1）
box on
xlabel（‘时间‘）
ylabel（‘智能体状态‘）
plot（tx（:1:5）‘-‘‘linewidth‘3.5）
hold on
plot（tx（:6:10）‘--‘‘linewidth‘3.5）
function dx=consensus（tx）
dx=zeros（101）;
dx（1）=x（6）;
dx（2）=x（7）;
dx（3）=x（8）;
dx（4）=x（9）;
dx（5）=x（10）;
dx（6）=k1*（a11*（x（1）-x（1））+a12*（x（2）-x（1））+a13*（x（3）-x（1））+a14*（x（4）-x（1））+a15*（x（5）-x（1）））+k2*（b11*（x（6）-x（6））+b12*（x（7）-x（6））+b13*（x（8）-x（6））+b14*（x（9）-x（6））+b15*（x（10）-x（6）））;
dx（7）=k1*（a21*（x（1）-x（2））+a22*（x（2）-x（2））+a23*（x（3）-x（2））+a24*（x（4）-x（2））+a25*（x（5）-x（2）））+k2*（b21*（x（6）-x（7））+b22*（x（7）-x（7））+b23*（x（8）-x（7））+b24*（x（9）-x（7））+b25*（x（10）-x（7）））;
dx（8）=k1*（a31*（x（1）-x（3））+a32*（x（2）-x（3））+a33*（x（3）-x（3））+a34*（x（4）-x（3））+a35*（x（5）-x（3）））+k2*（b31*（x（6）-x（8））+b32*（x（7）-x（8））+b33*（x（8）-x（8））+b34*（x（9）-x（8））+b35*（x（10）-x（8）））;
dx（9）=k1*（a41*（x（1）-x（4））+a42*（x（2）-x（4））+a43*（x（3）-x（4））+a44*（x（4）-x（4））+a45*（x（5）-x（4）））+k2*（b41*（x（6）-x（9））+b42*（x（7）-x（9））+b43*（x（8）-x（9））+b44*（x（9）-x（9））+b45*（x（10）-x（9）））;
dx（10）=k1*（a51*（x（1）-x（5））+a52*（x（2）-x（5））+a53*（x（3）-x（5））+a54*（x（4）-x（5））+a55*（x（5）-x（5）））+k2*（b51*（x（6）-x（10））+b52*（x（7）-x（10））+b53*（x（8）-x（10））+b54*（x（9）-x（10））+b55*（x（10）-x（10）））;
end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1927  2013-09-09 16:00  h1.m

----------- ---------  ---------- -----  ----

                 1927                    1