标签：sir 0.2 函数 figure 0.8 求解 plot 微分方程

定义函数 sir.m

function y = sir(t,x)
    %UNTITLED 此处显示有关此函数的摘要
    %   此处显示详细说明
    a=0.8;  %感染率0.8
    b=0.2;  %治愈率0.2
    y=[-a*x(1)*x(2),a*x(1)*x(2)-b*x(2)]';
end

运行函数 rum.m

[t,x]=ode45('sir',[0,50],[0.98 0.02]);
% ode45是用来求解常微分函数的方法
[t,x]
% 健康者与病人的图
figure(1);
plot(t,x(:,1),'.',t,x(:,2),'^');
title('SIR模型');
xlabel('时间');
ylabel('比例');
legend('健康者','病人','2');

% 健康者与病人的相轨线
figure(2);
plot(x(:,1),x(:,2));
title('i-s 相轨线');
xlabel('s');
ylabel('i');

运行结果

标签：sir,0.2,函数,figure,0.8,求解,plot,微分方程
来源： https://www.cnblogs.com/mrfanqie/p/matlab001.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。