0 个投票

I have a problem simulating the first equations of this article https://www.sciencedirect.com/science/article/pii/S2468042717300039#bib48
I need to solve it using the given parameters below the first graphics (fig.2).
I'been trying to solve the equations with ode23 with the same parameters and nothing.
Where should i start...

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

 采纳的回答

Ameer Hamza
Ameer Hamza 2020-11-6
编辑：Ameer Hamza 2020-11-6

1 个投票

You can use ode45(). Following code plot one set of results from the paper
IC = [1e6; 0; 1e-6];
tspan = [0 30];
[t, X] = ode45(@odefun, tspan, IC);
T = X(:,1);
I = X(:,2);
V = X(:,3);
figure();
subplot(2,1,1);
semilogy(t, V, 'r');
ylim([1e2 1e8])
subplot(2,1,2);
semilogy(t, I, 'b--');
hold on
semilogy(t, T, 'b-');
ylim([1e2 1e6])
function dXdt = odefun(t, X)
s = 10000;
d = 0.01;
beta = 1e-6;
rho = 1;
p = 1200;
c = 23;
T = X(1);
I = X(2);
V = X(3);
dTdt = s-d*T-beta*T*V;
dIdt = beta*T*V-rho*I;
dVdt = p*I-c*V;
dXdt = [dTdt; dIdt; dVdt];
end

更多回答（0 个）

类别

在 帮助中心 和 File 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