基于元胞自动机SARS传播模型.pdfVIP

第38卷第3期 数学的实践与认识 V01．38No．3 Feb．，2008 2008年2月MATHEMATICSINPRACTICEANDTHEORY ，’’’’’’’’’’’’、 ；应 用i ～‘ltl‘l‘·tl‘l‘， 基于元胞自动机的SARS传播模型 贺明峰1’2， 邓成瑞2 (1．大连理工大学应用数学系，大连116024) (2．大连理工大学大学生刨新院，大连116024) 摘要： 根据sARs的局部性、潜伏性特征以及传染病个体传播的特点，人群接魅通过近距离邻居发生，利 用元胞自动机模型时间矗!问上离散，以局部规则为基础，以同步更新为前提讨论整体性质的机制，将人群分 为易感者，带菌者，病人，免疫者之后，以易感者对sARs的抵抗能力、病人或带菌者对SARS的传染能力、人 群的大小为参效，建立了基于元胞自动机的sARs传播模型。得到的结果可为从传染动力学角度控制SARS 的传播提供依据． 关键词： 系统模型；元胞自动机；模拟；sARs；传染病 1 引 言 SARS作为一种局部区域性传播的传染病，具有明显的局部性与潜伏性特征[·_3]．在现 实中，人群的接触传播途径体现于近距离的邻居之间，而传染病(包括SARS)的传播有显著 的个体传播特点[4]．对其进行传播动力学研究，一方面可借助于传统的基于微分方程的数学 方法，如SIR模型，SEIR模型[53等．但微分方程模型对局部特点的处理不够灵活[钉．而元胞 自动机模型具有时间、空间的离散性特点，且以局部性规则为基础以同步更新为前提来讨论 整体性质口]．所以用元胞自动机模型来讨论传染病的传播规律，有望通过局部规则的设定而 使其更接近于实际情况．Ahmed等定义了一类基于元胞自动机的传染病模型，讨论了传播 于社会及生态领域问题的分析当中． 本文结合SARs的特点，将一个局部区域的人群分为易感者，带菌者(处于潜伏期)，病 人和免疫者(指治愈或死亡，本文不考虑天然的免疫者)四种情形，以易感者对SARS的抵抗 能力，病人和带菌者对SARS的传染能力以及人群的大小为参数，建立了一个新的基于元胞 自动机的SARS传播模型．讨论了抵抗能力，传染能力，人群的大小等对传染规模和传染时 间的影响，得到的定性结果可加深对SARS传播的了解，为从传染动力学角度控制SARS的 传播提供了依据． 2模 型 考虑一L×L的网格．每一格点(i，歹)处有一人．记。(1．如)为f时刻(i，歹)处的人的状态， 并规定：z“如)一O为易感者状态，z(f埘一1为非典潜伏期状态，z瓴如)一2为非典发病状态， 收稿日期：2004一07一06 万方数据 万方数据 3期 贺明峰。等：基于元胞自动机的SARS传播模型 43 g(1，J，件1)2gf．』，t)+1． ‘ (iii)如果z(f’如)一2，则 z“，j，件1)。=2，g(；，J，f+1)=。g(f，J，t)+1，若g(f，，，1)瓦； z(¨。件1)==3，g(z，J，t十1)==0，否贝Ⅱ． 一 (iv)如果z(州f)=3，则 O． z“．，，件1)23，g(“件1)2 3结 果 数，从病发期进入免疫状态所需的最长时间步孔=100．主要讨论凡，P”只及人群规模￡ 对SARS传播的影响． 3．1 n匕 此时，对任何规模的人群，只要有一个带菌者，则整个群体终将全部染病．但只，R的大 小对SARS的传播速度有明显的影响，见图2． (a) L=lOO (b)