Is there an easy way to make numerical simulations of the ODE of the form dx1/dt=x1(2-x1-x2), dx2/dt=x2(3-x1-x2-x3) for any xn?

Question

Editor 2018-4-5

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/392827-is-there-an-easy-way-to-make-numerical-simulations-of-the-ode-of-the-form-dx1-dt-x1-2-x1-x2-dx2-dt

评论： Editor 2018-4-6

I have tried doing for 'n=5' and here is my code;

dx = @(t,x) [x(1); 
    x(1)*(2-x(1)-x(2));
    x(2); x(2)*(3-x(1)-x(2)-x(3));
    x(3); x(3)*(3-x(2)-x(3)-x(4));
    x(4); x(4)*(3-x(3)-x(4)-x(5));
    x(5); x(5)*(2-x(4)-x(5))];
tspan=[0 15];
x0=[0 2 0 3 0 4 0 6 0 8];
[t,x] = ode45(@(t,x) dx(t,x), tspan, x0);
figure(1)
plot(t, x)

Can someone please check if this code is correct, if not, how can I improve if and for any 'n'?. Thanks in advance.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Abraham Boayue 2018-4-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/392827-is-there-an-easy-way-to-make-numerical-simulations-of-the-ode-of-the-form-dx1-dt-x1-2-x1-x2-dx2-dt#answer_313637

在 MATLAB Online 中打开

 See code comments for answers.  
clear variables
close all;
% dx must have 5 terms for n = 5
% the number terms in system equations
% must be the same as the number of 
% elements in your dx vector.
% Remember your dot before *
dx = @(t,x) [x(1);x(1).*(2-x(1)-x(2)); 
             x(2).*(3-x(1)-x(2)-x(3)); 
             x(3).*(3-x(2)-x(3)-x(4)); 
             x(4).*(3-x(3)-x(4)-x(5))]; 
tspan=[0 15]; 
x0=[0 2 0 3 0]; % same for the initial vector.
[t,x] = ode45(dx, tspan, x0); 
figure(1)
plot(t, x);

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Editor 2018-4-6

Thank you very much Abraham Boayue for your response and your code too does what I want. I have noted the mistakes I had made in my earlier code. However, is there a way in which this code can be be generalised to compute for any 'xn', say 'n=10'?.

请先登录，再进行评论。