怎样把ode45求解出来的微分方程组中，某一输入对应输出求出具体值

传染病传播模型

先要将第D天的s e i a r输出

% 未采取限制措施
function [s1,e1,i1,a1,r1]=ill(D)
% Day=D*7.2;
betal=0.00105;k=0.5;p=0.14;omiga=0.53;omiga1=0.83;gama=0.23;gama1=0.24;
f=@(t,x)[-betal*x(1)*(x(3)+k*x(4)),betal*x(1)*(x(3)+k*x(4))-(1-p)*omiga*x(2)-p*omiga1*x(2),(1-p)*omiga*x(2)-gama*x(3),p*omiga1*x(2)-gama1*x(4),gama*x(3)+gama1*x(4)]';
[t,y]=ode45(f,[0 D],[1665 0 2 0 0]);
plot(t,y(t,1),t,y(:,2),t,y(:,3),t,y(:,4),t,y(:,5));
% s1=y(Day,1);e1=y(Day,2);i1=y(Day,3);a1=y(Day,4);r1=y(Day,5);
legend('易感人群','潜伏期感染者','显性感染者','隐姓感染者','恢复者');

D = 100;
betal=0.00105;k=0.5;p=0.14;omiga=0.53;omiga1=0.83;gama=0.23;gama1=0.24;
f=@(t,x)[-betal*x(1)*(x(3)+k*x(4)),betal*x(1)*(x(3)+k*x(4))-(1-p)*omiga*x(2)-p*omiga1*x(2),(1-p)*omiga*x(2)-gama*x(3),p*omiga1*x(2)-gama1*x(4),gama*x(3)+gama1*x(4)]';
[t,y]=ode45(f,[0 D],[1665 0 2 0 0]);
plot(t,y(:,1),t,y(:,2),t,y(:,3),t,y(:,4),t,y(:,5));
% s1=y(Day,1);e1=y(Day,2);i1=y(Day,3);a1=y(Day,4);r1=y(Day,5);
legend('易感人群','潜伏期感染者','显性感染者','隐姓感染者','恢复者');

结果

您好，我是有问必答小助手，您的问题已经有小伙伴解答了，您看下是否解决，可以追评进行沟通哦~

如果有您比较满意的答案 / 帮您提供解决思路的答案，可以点击【采纳】按钮，给回答的小伙伴一些鼓励哦～～

ps: 问答会员【8折】购 ，仅需→￥23.2，即可享受5次/月 有问必答服务，了解详情>>>https://t.csdnimg.cn/RW5m