Can you simulate the following article?(Only the first formula: malaria epidemic dynamics) article name: Control of malaria outbreak using a nonlinear robust strategy with adaptive gains

1 次查看（过去 30 天）
af af
af af 2020-12-30
编辑： af af 2021-5-19
Can you simulate the following article?(Only the first formula: malaria epidemic dynamics)
article name: Control of malaria outbreak using a nonlinear robust strategy with adaptive gains
  3 个评论
显示 1更早的评论
Walter Roberson
Walter Roberson 2020-12-30
Some things like still look to be missing, and it has not been clarified whether the derivative is with respect to t or not.
is still not defined. It needs to be defined before we can be certain that it is possible to use ode45.
Knowing the exact values of the constants such as is not important to answer the question asked, which was whether I could simulate that formula. I offerered an extra answer as well, namely that you could simulate the system using ode45() . I recommend you read the documentation on that.

请先登录，再进行评论。

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

回答（1 个）

Walter Roberson
Walter Roberson 2021-1-1
Ran in:
P.R_h = rand; %unknown value, unknown unit
P.t = randn() * 10; %unknown value, unknown unit
P.lambda_h = 0.2; %per day
P.kappa = 0.7902; %per day
P.Beta_m = 0.8333; %probability
P.u_h = 0.00004; %per day
tspan = [0 10];
Sh0 = rand; %unknown initial condition
[x, Sh] = ode45(@(x,Sh)Sh_fun(x,Sh,P), tspan, Sh0 );
plot(x, Sh);
function dS_h = Sh_fun(x, S_h, P)
dS_h = P.lambda_h + P.kappa .* P.R_h - (1-u_1(P.t)).*P.Beta_m .* S_h - P.u_h .* S_h;
end
function u1t = u_1(t)
u1t = sin(t*7*pi).^2 - exp(-t); %function u_1 is not defined so make something up
end
Notice that the t being passed to u_1 has nothing to do with the integration: it is just some parameter whose value you did not specify. t is not the variable of integration of S_h -- if it were then there would have had to have been a in the definition of
  4 个评论
显示 2更早的评论
af af
af af 2021-5-19
编辑：af af 2021-5-19
How to implement Ehd in Simulink as shown below?
Eh0=200, Ehf=0, a=0.055, t=100,
Ehd=(Eh0-Ehf)*exp(-at)+Ehf
Where is the problem in Simulink below?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Multibody Modeling 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by