求解三ode系统并绘制相空间图
1次观看(过去30天)
显示旧的注释
我试图解决这些方程,这是一个修正的捕食者模型,我必须在x, y和Z轴上分别绘制M, N, Z的相空间图,以获得不同的tau值。
dM/dt = r1M(1−M/k1)−α1MN
dN/dt = βNZ(t−τ)−d1N−α2MN
dZ/dt = r2Z(1−Z/k2)−βNZ(t−τ)−d2Z
我尝试过的程序和代码附在这里。有人能检查错误并帮助我尽快修改程序吗?
14日的评论
Aiswarya K
2020年4月6日
编辑:darova
2020年4月6日
先生,我已经添加了这个术语,但还是出现了错误。实际的图应该像图1,我得到的是图2 tau = 35.95。每个tau值都会出现类似的错误。实际上这是一个时滞微分方程。有没有其他的语法来解这样的方程?
Aiswarya K
2020年4月7日
函数dx = cancer_control(t, x) r1 = 0.18;K1 = 5*10^6;A1 = 1.101*10^-7;B = 6.2*10^-9;D1 = 0.0412;A2 = 3.422*10^-10;R2 = 0.0245;K2 = 1*10^7;Tau = 35.95;Dx = [0; 0; 0]; dx(1) = r1*x(1)*(1-x(1)/k1)-a1*x(1)*x(2); dx(2) = b*x(2)*x(3)*(t-tau)-d1*x(2)-a2*x(1)*x(2); dx(3) = r2*x(3)*(1-x(3)/k2)-b*x(2)*x(3)*(t-tau)-d2z;
Aiswarya K
2020年4月8日
编辑:darova
2020年4月8日
函数Dx = cancer_control(t, x)
R1 = 0.18;
K1 = 5*10^6;
A1 = 1.101*10^-7;
B = 6.2*10^-9;
D1 = 0.0412;
D2 = 0.0412;
A2 = 3.422*10^-10;
R2 = 0.0245;
K2 = 1*10^7;
Tau = 35.95;
Dx = [0;0;0);
Dx (1) = r1*x(1)(1-x(1)/k1)-a1*x(1)*x(2);
Dx (2) = b*x(2)*x(3)*(t-tau)-d1*x(2)-a2*x(1)*x(2);
dx (3) = r2 * x (3) * (1 - x (3) / k2) - b * x (2) * (3) * (t-tau) d2 * x (3);
答案(0)
发生错误
由于对页面所做的更改,无法完成操作。重新加载页面以查看其更新状态。
您也可以从以下列表中选择一个网站:
如何获得最佳的网站性能
选择中国站点(中文或英文)以获得最佳站点性能。其他MathWorks国家站点没有针对您所在位置的访问进行优化。