1、首先是SIR传裔瞍爿嘞染病模型,也就是三个微分方程,如图。其中中的S表示易感者,I表示感染者,R表示移出者。a表示易感者的感染概率,b表示感染者的恢复概率,也就是移出率。
2、然后是编一个模型的代码,保存为sir.m。a和b的值可根据自己的情况设定。
3、再编写一个运行的m文件,可命名为sirrun.m,代码如下。[0,40]是t的范围,[0.04 0.95 0.01]分别是I、S和R的初始值。
4、运行的结果如下 ,其中红绿线表示S,蓝线表示I,线表示R。
5、运行完也会显示结果数据,由于数据较长,这里只截图一部分,如下: