我正在尝试使用耦合的 DE 对 Lotka-Volterra Predator-Prey 系统进行建模：
dy(1)/dt = rx(1-x/k) - ay(1)y(2) % 猎物种群
dy(2)/dt = aby(1)y(2) - dy(2) % 捕食者种群
这是我的代码：
% Solves equations using numerical ODE solver 45 (nonstiff runge kutta)
runtime = 1000; % Duration time of simulation in seconds.
%Parameter values used in simulation
r = 0.5; % exponentional growth rate of prey in absence of predator
a = 0.01;% conversion efficiency of predator
b = 0.02; % attack rate
d = 0.10; % death rate
k = 750; % carrying capacity
y0 = [10, 10] % initial conditions y(1)= 10, y(2) = 10
deq1=@(t,y) [r.*y(1)*(1-(y(1)./k))- a.*y(1)*y(2); a.*b.*y(1).*y(2)- d.*y(2)];
[t,sol] = ode45(deq1,[0 runtime],y0);
如何更改我的代码以使时间序列图（y(1) 与时间）显示振荡，而不是稳定状态？在集成步骤中似乎出了点问题，因为这些图的结构是我想要的，但函数的行为不是我所期望的。
稳态时间序列图
相图达到平衡
