How to plot with various initial conditions?

15 次查看（过去 30 天）
marya
marya 2017-6-9
评论： marya 2017-6-9
Edited: I want to plot
using the model with different initial conditions(use whatever you like)
function dydt = model (t,y)
dydt = zeros (size(y));
Ah=0.000051;Av=0.071;b1=0.071;b2=0.091;g=0.0035;
d1=0.0000043;d2=0.04;e1=0.001;w=0.11;
Sh=y(1);
Ih=y(2);
Rh=y(3);
Sv=y(4);
Iv=y(5);
Nh = Sh+Ih+Rh;
%The model
dydt(1) = Ah - b1*Iv*Sh/Nh +w*Rh - d1*Sh;
dydt(2) = b1*Iv*Sh/Nh - (g + e1 + d1)*Ih;
dydt(3) = g*Ih - (w +d1)*Rh;
dydt(4) = Av - b2*Ih*Sv/Nh - d2*Sv;
dydt(5) = b2*Ih*Sv/Nh - d2*Iv;
% plot
tspan = [0 10000];
y0 = [3600 1000 100 9600 400];% Here is the initial condition
[t,y] = ode45(@model,tspan,y0);
plot(t,y(:,2),'r','Linewidth',1)
title('Plot of human population against time')
xlabel('Time(years)')
ylabel('Number of People')
legend('Infectious')
Regards,
  4 个评论
显示 2更早的评论
Torsten
Torsten 2017-6-9
And what should be the result of the three lines of code below ?
Best wishes
Torsten.

请先登录，再进行评论。

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

采纳的回答

KSSV
KSSV 2017-6-9
tspan = [0 10000];
% initilai conditions
figure
hold on
for i = 1:5
% y0 = [3600 1000 100 9600 400];
y0 = rand(1,5)*100*i ;
[t,y] = ode45(@model,tspan,y0);
plot(t,y(:,2),'r','Linewidth',1)
end
title('Plot of human population against time')
xlabel('Time(years)')
ylabel('Number of People')
legend('Infectious')
  3 个评论
显示 1更早的评论
Torsten
Torsten 2017-6-9
Init=[13413000 3350000 3343000 16500000 38000000
13400 3350 3000 165000 38000
1341300 436000 334300 1650000 3800000
1341300 55000 33430 1650000 3800000
13400 3200 2800 16500000 38000000
134000 43350 3000 165000 38000
15130 6000 4300 1650000 3800000
16000 5000 4000 1650000 3800000];
tspan = [0 10000];
% initilai conditions
figure
hold on
for i = 1:8
y0 = Init(i,:) ;
[t,y] = ode45(@model,tspan,y0);
plot(t,y(:,2),'r','Linewidth',1)
end
title('Plot of human population against time')
xlabel('Time(years)')
ylabel('Number of People')
legend('Infectious')

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Genomics and Next Generation Sequencing 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by