由于是非专业人士,直接找到一个例子,copy如下
首先自定义一个函数
function z=weifen_bulie(t,y)
z(1,:)=y(1)-0.05*y(1).*y(2)
z(2,:)=-y(2)+0.03*y(1).*y(2)
end
之后直接运行
% 捕食者与被捕食者的数量变化
% 捕食者与被捕食者的数量相位变化图
Y0=[100,20]
[t,Y]=ode23('weifen_bulie',[0,20],Y0);
x=Y(:,1);
y=Y(:,2);
figure(1)
plot(t,x,'b',t,y,'r')
title( '捕食者与被捕食者的数量变化');
xlabel( '演化代数');
ylabel( '数量');
legend('兔子','狐狸')
figure(2);