legion model equations i want to make this graph from these equation using ode45

1 次查看(过去 30 天)
  8 个评论
David Wilson
David Wilson 2019-6-19
Here's my effort: Your above effort needs to have two differential eqautions, and it would help to know what sort of parameters you are interested in.
%% Oscillatory phase plot
I = 0.2;
e = 0.1; % very small
a = 6; b = 0.1; % guessed from the plot
xdot = @(t,x) [ ...
3*x(1)-x(1).^3+2-x(2)+I; ...
e*(a*(1+tanh(x(1)/b))-x(2))];
tspan = [0,500];
x0 = [0,3]'; % initial conditions in the middle
[t,x] = ode45(xdot, tspan, x0);
% plot(t,x)
plot(x(:,1), x(:,2),'.-')
xlabel('x'); ylabel('y')

请先登录,再进行评论。

回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Ordinary Differential Equations 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by