fun.m

来自「智能控制中的新热点问题」· M 代码 · 共 30 行

clc
clear

global m n a b c d1 d2 ;

m=5;%“鱼”的数量
n=3;
a=1;
b=20;
c=0.2;
d1=0.1;
d2=1;
S=50;%活动范围
x0=[S*(rand(m*n,1));0;-200;-200];%产生一个随机的矩阵12*2，所以每次画出来的图是不一样的
[T,X]=ode45(@ode1,[0 15],x0);%这就是求解这个微分方程，用的是ode45函数,dx定义用的是ode1函数,0-100是时间范围
for i=1:3:3*(m)
    plot3(X(:,i),X(:,i+1),X(:,i+2));
    hold on
end;
    plot3(X(:,3*m+1),X(:,3*m+2),X(:,3*m+3),'g');

for i=1:3:3*(m+1)
   comet3(X(:,i),X(:,i+1),X(:,i+2),0.01);
hold on
end;