wodjkrvd 发表于 2010-6-16 09:09

这是关于sars传播模型的分析与预测和数值模拟并求得未知参数

N是所考察地区的人口总数
H是潜伏期天数
x是未被控制的感染者每天接触并传染易感者的比例系数
y是潜伏者的日发病率
r是日移除率(为死亡率与治愈率之和)
c是控制因子
μ是感染者日自愈率
n是每日未被控制感染者被控制的比例系数
综合上述补充,我们又得到新的SARS传播模型

S`= -xSIu
E`=xSIu-y(1-c)E-μE
Ii`=nIu-rIi

Iu`=y(1-c)E-nIu
R`=rIi+μE
S`+E`+Ii`+Iu`+R`=0

S+E+Ii+Iu+R=N

其中初值S(0)=S0 、E(0)=E0 、Ii(0)=Ii0 、Iu(0)=Iu0 、R(0)=R0 均大于0,t∈,T为传染病流行时间,y(1-c)+μ=1/H,0≤y,μ,r≤1
由于改良模型的参数增多,这样也增加了我们计算的困难,为了计算简便,我们需将参数逐一分析其实际意义,并对有可靠研究的参数的取值采取直接套用:(1)由历史数据可知,2006年北京市居住半年以上人口为1581万人,在籍人口为1197.6万人,为了计算方便,我们假设北京2003年人口总数为1000万。
(2)y(1-c)是潜伏者的日发病率,μ是潜伏者的日自愈率。假定每个潜伏者实际潜伏天数服从区间上的均匀分布,即潜伏者在H天内等可能的发病获自愈,因而得到y(1-c)+μ=1/H。
(3)n是未控制的感染者每天被控制的比例系数,取决于社会的对SARS的重视程度以及政府和医疗卫生机构所采取措施的严格程度。
(4)r是移除率,即已控制感染者每天治愈率和死亡率之和。根据所给的数据,可以直接得到。
(5)由资料可知,潜伏期天数为2-12天,我们取其中位数,假设潜伏期均为7天,即H=7。
(6)根据SARS机理的研究,大约有90%的处于潜伏期的感染者会自愈,只有很少一部分会真正发病,由(2)可知,y(1-c)+μ=1/7。
页: [1]
查看完整版本: 这是关于sars传播模型的分析与预测和数值模拟并求得未知参数