MATLAB的答案

我从我的代码有两种不同的图像ODE89和欧拉方法。我需要比较图表。

1视图(30天)
Naveen Krish
Naveen Krish 2022年3月14日
评论道: Naveen Krish2022年3月18日
ODE89
clc
清除所有
A0 = 1;
B0 = 3;
P0 = 0;
K = 5 * 10 ^ 5;
Yb = 1;
Yp = 0.15;
tspan = [0 43200];
[t、Y] = ode89 (@ (t, Y) odefun (t, Y, K, Yb, Yp), tspan, [A0; B0 P0]);
图(1)
情节(t, Y (: 1))
图(2)
情节(t、Y (:, 2))
图(3)
情节(t, Y (:, 3))
函数dYdt = odefun (t, Y, K, Yb, Yp)
dYdt = [(Y - k * Y (1) * (2));
(yb * (K * Y Y (1) * (2)));
(Yp * (K * Y Y (1) * (2))));
结束
欧拉
nsteps = 12;
t = 0 (nsteps, 1);
一个= 0 (nsteps, 1);
B = 0 (nsteps, 1);
P = 0 (nsteps, 1);
(1)= 1;
B (1) = 3;
C (1) = 0;
K = 5 * 10 ^ 5
k = 13
t (k) = (k - 1) + 3600
(k) = (k - 1) + (- k * (k - 1) * B (k - 1)) * 3600;
(k) = B (k - 1) + (yb * (k * (k - 1) * B (k - 1))) * 3600;
P (k) = P (k - 1) + Yp * (k * (k - 1) * B (k - 1)) * 3600;
结束
情节(t)
图(1)
情节(t (: 1))
情节(t, B)
图(2)
情节(t, B (: 1))
情节(t, P)
图(3)
情节(t, P (: 1))
我需要的ODE89图比较,B与欧拉和P A, B和P
2的评论
显示1年长的评论

登录置评。

在回答这个问题。

接受的答案

大卫。马塞罗说道
大卫。马塞罗说道 2022年3月15日
编辑:大卫。马塞罗说道 2022年3月15日
跑:
clc
清晰的所有
% % ODE89
A0 = 1;
B0 = 3;
P0 = 0;
K = 5 * 10 ^ 5;
Yb = 1;
Yp = 0.15;
tspan = [0 43200];
(x, Y) = ode89 (@ (t, Y) odefun (t, Y, K, Yb, Yp), tspan, [A0; B0 P0]);
% %欧拉
nsteps = 12;
t = 0 (nsteps, 1);
一个= 0 (nsteps, 1);
B = 0 (nsteps, 1);
P = 0 (nsteps, 1);
(1)= 1;
B (1) = 3;
C (1) = 0;
K = 5 * 10 ^ 5;
为k = 13
t (k) = (k - 1) + 3600;
(k) = (k - 1) + (- k * (k - 1) * B (k - 1)) * 3600;
(k) = B (k - 1) + (yb * (k * (k - 1) * B (k - 1))) * 3600;
P (k) = P (k - 1) + Yp * (k * (k - 1) * B (k - 1)) * 3600;
结束
图(1)
次要情节(1、3、1)
情节(t, A, x, Y (: 1))
包含(“t”)
ylabel (“一个”)
次要情节(1、3、2)
情节(t, B, x, Y (:, 2))
包含(“t”)
ylabel (“B”)
次要情节(1、3、3)
情节(t、C、x, Y (:, 3))
包含(“t”)
ylabel (“P”)
传奇(“ODE89”,“欧拉”,“位置”,“最佳”)
函数dYdt = odefun (t, Y, K, Yb, Yp)
dYdt = [(Y - k * Y (1) * (2));
(yb * (K * Y Y (1) * (2)));
(Yp * (K * Y Y (1) * (2))));
结束
3评论
显示2年长的评论

登录置评。

更多的答案(0)

类别

找到更多的在数学在帮助中心和文件交换

标签

下载188bet金宝搏


释放

R2021b

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!