关于捕食者方程在matlab中的应用

由于是非专业人士，直接找到一个例子，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);