(3) 对于卫生部门所采取的措施做出评论，如：提前或延后5天采取严格的隔离措施，对疫情传播所造成的影响做出估计。（附件2提供的数据供参考。）
3、收集SARS对经济某个方面影响的数据，建立相应的数学模型并进行预测。（附件3提供的数据供参考。）
4、给当地报刊写一篇通俗短文，说明建立传染病数学模型的重要性。
二、 对早期模型的评价
附件1的模型主要采用“数据拟合”和“借鉴参数”的方法对北京疫情走势进行预测。
在数据拟合方面，该模型中有两个疑点：
1、感染期限L的确定。由于被严格隔离、治愈、死亡等原因，感染者在某一时段后不再具有对易感人群的传染力，故对病毒的传染加上感染期限是合理的。但在对该参数的确定上，作者为了较好地拟合各阶段的数据 ,通过人为调试来确定L的取值，缺乏医学上的支持，使模型的说服力减弱，合理性和可靠性大大降低。
2、文中认为“K代表某种环境下一个人传染他人的平均概率”。但从模型的公式中可以看出，参数K的实际意义是一个病人平均每天传染其他人的个数。两者之间有实质的区别，文中的说法显然不妥。
从预测思想来看，该模型是借鉴先发地区——广东、香港的有关参数对北京的疫情进行预测的。由于广东、香港的疫情和控制都在北京之前，已经过了高峰期，到5月8日为止每日新增病例已降至10来例，基本处于后期控制阶段。而当时北京的疫情刚过了高峰期，正处于社会剧烈调整时期，数据较为凌乱，略有下降趋势，但不明显。可见在当时，采取这种借鉴是无奈之举。
但是由于城市之间的政策，风俗习惯等不同，城市之间的可比性不强，借鉴存在很大的局限性。如在香港，由于对传播机制认识不足，中途又出现高度感染的特殊情况。另外使用借鉴法无法对首发城市进行预测。
三、 传播模型
（一）问题的分析
在SARS爆发的初期，由于潜伏期的存在，人们对病毒传播速度和危害程度的认识不够，未能及时识别这一传染病的存在。但当病患数不断增加，政府开始采取应对措施对其进行控制，同时社会舆论加大宣传力度，人们的警觉性提高，病毒的传播速度下降。因此，我们通常把传染病的传播模式近似分为两个阶段：
第一、自由传播阶段（即控前阶段）：在采取切实有效的控制措施之前的一段时间。
第二、控后阶段：介入人为因素之后的一段时间。
由于SARS的传播涉及的因素很多，如潜伏期、人群的迁入迁出，感病者的数量、易感者的数量、传染率和治愈率的大小等。而且在以上因素中，潜伏期的大小、传染率和治愈率的大小因人而易，具有一定的随机性。不可能一开始就把所有的因素全部考虑在内建立模型。对此，我们将作出相应的假设进行简化。
分析附件2所给出的数据，发现6月1日至15日，已确症病例累计数为2522人，其中夹杂3天累计数为2523人。但6月16日后累计数降至2521人，认为
2
累计数的减少可能是误诊引起的。由于误诊的可能性很低，故在这里忽略不计。
（二）基本假设
1、国家卫生部提供的北京疫情统计真实可信。（误诊数仅为1，可忽略不计）。 2、由于非典的主要传播途径是近距离接触，通过受感染者咳嗽或打喷嚏时产生的飞沫传播，这里将所有传播途径都视为与病源的直接接触。
3、不考虑出生与自然死亡的过程和人群的迁入迁出（或认为迁入和迁出基本平衡），认为疾病传播期间所考察地区的总人数为常数。
4、根据国家卫生部资料可知处于潜伏期的SARS病人不具有传染性。
5、目前尚不清楚康复患者是否具有免疫力，但据国家卫生部资料可知康复后的病患无一例复发。故假设康复患者退出传染系统。
6、根据资料显示，SARS病毒的潜伏期一般为2~7天，平均约为5天。（这一条件将在后期的模型中有所改动）
（三）常用基本模型
目前常用的传染病模型，通常将传染病流行范围内的人群分为三类： S类：易感者，指未得病者，但与感病者接触后容易受到感染。 I类：感病者，指染上传染病的人。
R类：移出者，指因患病而被隔离，或因病愈而具有免疫力的人，他们即非感病者，也非易感者，实际上他们已经退出了传染病系统。
并通过三类之间的互相转化关系建立微分方程组进行求解：
?dS?dt??kIS?dI???kIS?hI (1) ?dt?dR?hI?dt???S?I?R?N变量和符号说明：
k ——传染率：每个病人平均每天有效接触（足以使被解除者感染）的人
数。
h ——退出率：单位时间内治愈和死亡人数占感病者人数的百分数。
S(t)——易感人群的总数。
I(t)——感病者总数。 R(t)——退出者总数。 N——一个城市总人口数。
观察附件二中给出的数据，我们发现截至6月23日，感病者累计为2521人，远远小于北京城市的总人口数150万人，故认为感病者和退出者对易感人群的总数影响不大，易感者总人数I为一常数。原方程变形为：
3
?dI?kIN?hI??dt （2） ??dR?hI??dt注意到退出者不是我们研究的范围，故方程组（2）实际上是一个常微分方程
dI?kNI?hI??I （3） dt其中??kN?h，
不难用分离变量法解出：
I(t)?I0e??t （4）
其中I0为初始值。
根据以上分析我们可以看出，常微分方程的传染病模型只适用于病例数与总人口数具有可比性的情况。当病例数远小于总人口数时，常微分方程模型的实质与附件1的模型相同，感病人数将随时间以指数增长。
考虑这一特点，我们用计算机跟踪病毒的个体传播情况，建立了模拟模型。
（四）计算机模拟模型：
在该模型中，我们将传染系统中的人分为五类：
自由携带者（f[t]）——身上携带病毒并均匀散布在人群中的患者，根据基本假设自由携带者在潜伏期内不具有传染力，
日增患者（x[t]）——每天被医疗部门发现并加以隔离的感病者
被隔离者（y[t]）——因曾与自由携带者接触而被怀疑携带SARS病毒的人 有效接触者（z1[t]）——每日与自由携带者接触并感染上病毒的人 无效接触者（z2[t]）——每日与自由携带者接触但未染上病毒的人 并作出如下假设：
1、由于传染性SARS最初（1~2天）的症状通常为发热（?38o），发热通常为高热[1]。症状明显，易于辨认，故可认为自由携带者发病后当天或第二天就立即入院治疗，入院后不会再参与疾病的传播。
2、根据实际情况，假设SARS病人被发现的三天内，有关部门将采取措施，将部分与病源有效接触者隔离，这部分人即使发病后也不会参与疾病的传播。
3、与病源有效接触者必然发病。根据基本假设，潜伏期一般为2至7天，这里取为5天。（这一假设在改进模型中有进一步的讨论。）
另外，对模拟模型中出现的符号变量说明如下：
k1 ——有效接触率，表示一个自由携带者平均每天有效接触的人数。
4
k2 ——无效接触率，表示一个自由携带者平均每天无效接触的人数。
（包括有效接触和无效接触）的人群中可以控制?——与自由携带者接触后
的人数所占的百分比。
模拟模型中个体传播情况如图1所示：自由携带者1~5天处于潜伏期，不具有传染能力；5天后发病，发病后每天有效接触k1人，2天后（第7天）被隔离，在隔离前每天无效接触k2人。与病源接触后的可控人群（占接触者总人数的?）在3天后被视为疑似病人。疑似病人中的有效接触者在接触病源的第7天被发现确认为日增患者，而有效接触者中其他人作为自由携带者留在人群中，继续这之前的个体传播。
有效接触者 自由携带者 第6，7天 3天后 自由携带者 7天中 无效接触者 被隔离者 7 天后 日增患者 3天后 4天后 日增患者 被隔离者
图1：个体传播示意图
用数学模型描述各个变量之间的关系如下：
?z1[i?5]?k1f[i]?z[i?6]?kf[i]11???z2[i?j]?k2f[i](j?0,1,2,?6) （5） ?f[i]?(1??)z[i]1???x[i]?f[i?6]??z1[i?6]由于北京在4月20日才开始建立每日疫情报道制度，故认为政府采取严格
的隔离措施开始于4月20日，以这一天为分界线，之前属于自由传播阶段。
根据这一模型，用计算机模拟北京5月10日之前SARS的传播情况，并对5月10日以后的传播情况进行预测。
5
图2：5月10日以前数据拟合图
图3：5月10日以后的预测曲线
通过图2两条曲线的拟合，得到控前的有效接触率（表征病毒的传染力）k=1.35>1,可控率（表征政府的控制力度）??0.5；控后k1=0.8， ?=0.7。根据这两个参数作出5月10日的预测曲线见图3。，根据预测，北京将在第97天（6月下旬）实现零增长，累计病例数2448人。
经过分析，以上确定型的模拟模型存在以下两点问题：
第一、SARS病毒的潜伏期一般是2~7天，模型将潜伏期确定为5天，从感染到被发现的时间确定为7天，可以明显地看到以7天为传播周期的曲线变动，
6