MACDONALD DISTRIBUTED DELAY MODEL

3 次查看（过去 30 天）
Chukwuma Ugoala
Chukwuma Ugoala 2019-10-17
编辑： Chukwuma Ugoala 2019-10-23
Can any one help me with the plot of the system below. The code I have is taking forever to plot.
Here are the scripts for the above system
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function M = macdistime(t,x)
global a b c mu m r alp
dxdt = a*b*m*x(2)*(1-x(1))-r*x(1);
dydt = a*c*x(4)-mu*x(2);
dzdt = a*c*x(1)*(1-x(2)-x(3))-a*c*x(4)-mu*x(3);
dgdt = alp*x(1)*(1-x(2)-x(3))- (alp+mu)*x(4);
M = [dxdt;dydt;dzdt;dgdt];
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
mac_parameters.m
global a b c mu m r alp
a = 0.25; %Man biting rate
b = 0.09; %Proportion of bites that produce infection in human
c = 0.47; %Proportion of bites by which one suceptible mosquito becomes infected.
mu = 0.12; %Per capita rate of mosquito mortality
m = 40.0; %Ratio of number of female mosquitoes to that of humans
r = 0.0035;%Average recovery rate of human
alp = 10; % latency period of mosquitoes set to 10 days
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
run mac_parameters.m
for reps = 1:10;
xstart = [0:0:0:0]; % sets each of our
%initial variable values to a number between 0 and 100 and
[t,x]=ode45(@macdist,0:1:200,xstart); % integrate the system
%forward in time 200 days, with interval of 1 day
for j=1:4;
s(j) = subplot(3,3,j); plot(t,x(:,j)); xlabel('Time(days)'); ylabel('Number of infected humans');
hold on
end
title(s(1), 'x'); title(s(2), 'y'); title(s(3), 'z'); title(s(4), 'g')
end

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

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by