病毒传播SIS模型研究1.docVIP

摘要

问题一要求我们做出节点连接度和其出现频率的关系。我们先是对数据进行处理，得到A、B类节点连接度、及其出现概率、的表格。为了发现它们之间服从什么分布，我们先是用Matlab软件中的cftool指令拟合出表格中的数据服从的曲线，然后再对这条曲线进行验证，是否可以作为概率分布的曲线，最终得出它们服从的分布。最后得出A类和B类节点连接度的概率分布的函数关系均是。即A、B两类节点连接度服从度为2.17的幂律分布。

问题二要求先生成一个小型的二部图，并用一定的生成机制，得到满足数据文件特点的网络。我们做出了BA无标度网络模型，在新生成的小型二部网络的扩充中，在A类中加入一个新节点并向B类中的节点连边，该边指向B类中号节点的概率正比于号节点当前的连接度。这里的概率我们也取，并把此时的概率分布在0-1之间，而所对应的区间为。最终新加的点是否连接到那个节点，由随机数落入的区间来考虑。最后我们得出A、B两类节点的连接度及其出现频率是符合度为2.31的幂律分布，其具体的函数表达是为。

问题三要求我们在这类网络上考虑“易感－感染－易感”（SIS）模型，得到较平稳时期的得病数量以及A类和B类的得病比例。我们在问题二的生成机制上，通过产生随机数来判断各个节点是否患病且能否传染，利用Matlab软件对该产生机制进行编程并通过循环。并对该模拟程序进行10次，我们得到一个模拟的值，取平均值后得到：当、、时，A、B两类个体的得病率分别为1.14%，1.22%、1.1%，1.44%、

1.36%，1.54%。

问题四要求对我们的模型进行理论的分析，看看是否和我们的模拟结果一致。最后我们得出，两者在感染率上升时的A、B两类节点的感染比例变化上存在较小的差异。但在感染节点的个数上都是较少的。这说明模型和模拟得出的结果基本相同。

关键词：幂律分布，BA无标度网络模型，二部图，SIS模型，

问题重述

病毒传播问题的研究由来已久，而一再的病毒流行使得这一领域长期以来吸引着人们的注意。在对病毒传播过程的描述各种模型中，“易感－感染－易感”（SIS）模型是研究者经常的选择。关于SIS模型，可以简单的描述为：一个易感的个体在和一个具有传染性的个体的接触中，在单位时间以一定的概率（β）被感染，同时，已感染的个体以概率（γ）被治愈又重新成为健康（易感）的个体。

实际中大量的问题可以利用网络（图）进行描述，比如在传染病问题的描述中，个体（人、动物、计算机等）可以看作网络的节点，当个体之间有可以导致病毒传播的接触时在两个个体之间连边。比如，对于接触性传染病，个体存在两种状态，健康的（易感的）和已感染的；将这些个体作为网络的节点，由于两个个体之间的亲密接触可能导致病毒的传播，因此可在两者之间进行连边。一个个体所接触的其它个体数量称为该节点的度（边数）。所谓二部网络（图），是网络中的节点可分成两类（比如男性和女性，雄性和雌性等），边仅仅存在于两类节点之间。

在经典的传染病学模型中，总是假定病毒赖以传播的网络具有匀质性，即网络中节点有基本相同的度，但一些研究表明，这一假设远远背离实际情况。因此，发现实际网络的一些特性，并研究这样的网络上的病毒传播问题具有理论和实际意义。

本题我们主要研究二部网络上的病毒传播问题，根据附件提供的一个二部网络（由10000个A类节点和10000个B类节点构成）的节点度的数据，完成以下任务：

1．根据“附件”提供的数据data.xls，选择适当的坐标，作出节点连接度和其出现频率的图形，观察这种类型的连接度数据大致服从什么分布？

2．生成上述网络，可以采用如下的机制：先生成一个小型的二部图，随后在A类中加入一个新节点并向B类中的节点连边，该边指向B类中号节点的概率正比于号节点当前的连接度，而后在B类中产生新节点，以同样的方式向A类连边，当这两个步骤进行足够多次之后即可得到满足数据文件特点的网络。根据这里所提供的生成机制，发现节点连接度分布的表达式。

3．在这类网络上考虑“易感－感染－易感”（SIS）模型，得到较平稳时期的得病数量以及A类和B类的得病比例。(参数γ＝0.1,考虑到两类个体的感染率可以不同，分析中假定A类个体的感染率为B个体感染率的2倍，即=2，并分别取B类个体的感染率=0.01,0.02,0.03)。由于考虑PC机的计算速度，模拟时网络规模不要太大，可选择500＋500的二部网络。

4．对我们的模型进行理论的分析，看看是否和我们的模拟结果一致。

问题分析

问题背景的分析：

随着卫生设施的改善，医疗水平的提高以及人类文明的不断改善，诸如霍乱、天花等曾经肆虐全球的传染性疾病已经得到有效的控制。但是，一些新的、不断变异着的传染病毒却悄悄向人们袭来。20世纪80年代十分险恶的艾滋病毒突袭人间，至今仍在蔓延；随后SAS病毒、H1