求matlab大神把 Lingo 代码转化为matlab代码

常见的传染病模型按照具体的传染病的特点可分为 SI、SIS、SIR、SIRS、SEIR 模型。其中“S”“E”“I”“R”的现实含义如下:

S (Susceptible),易感者,指缺乏免疫能力健康人,与感染者接触后容易受到感染;

E (Exposed),暴露者 ,指接触过感染者但暂无传染性的人,可用于存在潜伏期的传染病;

I (Infectious),患病者,指有传染性的病人,可以传播给 S,将其变为 E 或 I ;

R (Recovered),康复者,指病愈后具有免疫力的人,如是终身免疫性传染病,则不可被重新变为 S 、E 或 I ,如果免疫期有限,就可以重新变为 S 类,进而被感染。

SIS传染病模型假设:
1.在疾病传播期内,研究区域内总人数N不变,既不考虑生死,也不考虑迁移。将该区域内的人群分为:易感染者、感病者和病愈免疫的移出者三类,以下简称易感染者,病人和移出者。在t时刻将这三类人群的数量分别记为S(t)、I(t)和R(t),t时刻这三类人在总人数中所占的比例分别记为s(t)、i(t)和r(t)。
2.每个病人每天有效接触的平均人数是常数λ,称为日接触率。当病人与健康者有效接触时,会使健康者受感染变为病人。
3.每天被治愈的病人占病人总数的比例为常数μ,称为日治愈率。病人被治愈后仍有可能被感染为病人,故假设1/μ是该疾病的平均传染期。
4.在初始时刻,只有少数个体处于感染状态,其他都是易感染状态。
5.假设疾病的时间尺度远小于个体生命周期,从而不考虑个体的出生和自然死亡对数据统计造成的影响。
6.完全混合:每一个个体与其他个体接触的机会均等。
根据以上假设可得:每个病人每天可使λ·s(t)个健康者变为病人,每天有μ·Ni个病人被治愈。
由该式中,接触数σ是病人在平均传染期内有效接触的人数。

SIS传染病模型构建:
由假设1可得,整个群体由易感染者,感染者,移出者构成,三者之间的关系如下:
将数据带入上述公式,求解微分方程即可得到各个变量随时间变化的情况。

}

身份认证 购VIP最低享 7 折!

基于MATLAB和LINGO的数学实验_3_肖华勇编著_西北工业大学出版社 绝对的稀缺资源

}

我要回帖

更多关于 matlab工具箱生成代码 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信