?? fun.m
字號:
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;
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -