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

摘 要 问题重述 病毒传播问题的研究由来已久，而一再的病毒流行使得这一领域长期以来吸引着人们的注意。在对病毒传播过程的描述各种模型中，“易感－感染－易感”（SIS）模型是研究者经常的选择。关于SIS模型，可以简单的描述为：一个易感的个体在和一个具有传染性的个体的接触中，在单位时间以一定的概率（β）被感染，同时，已感染的个体以概率（γ）被治愈又重新成为健康（易感）的个体。 实际中大量的问题可以利用网络（图）进行描述，比如在传染病问题的描述中，个体（人、动物、计算机等）可以看作网络的节点，当个体之间有可以导致病毒传播的接触时在两个个体之间连边。比如，对于接触性传染病，个体存在两种状态，健康的（易感的）和已感染的；将这些个体作为网络的节点，由于两个个体之间的亲密接触可能导致病毒的传播，因此可在两者之间进行连边。一个个体所接触的其它个体数量称为该节点的度（边数）。所谓二部网络（图），是网络中的节点可分成两类（比如男性和女性，雄性和雌性等），边仅仅存在于两类节点之间。在经典的传染病学模型中，总是假定病毒赖以传播的网络具有匀质性，即网络中节点有基本相同的度，但一些研究表明，这一假设远远背离实际情况。因此，发现实际网络的一些特性，并研究这样的网络上的病毒传播问题具有理论和实际意义。 本题主要研究二部网络上的病毒传播问题，附件提供一个二部网络（由10000个A类节点和10000个B类节点构成）的节点度的数据，完成以下任务：1．根据“附件”提供的数据data.，选择适当的坐标，作出节点连接度和其出现频率的图形，观察这种类型的连接度数据大致服从什么分布？ 2．生成上述网络，可以采用如下的机制：先生成一个小型的二部图，随后在A类中加入一个新节点并向B类中的节点连边，该边指向B类中号节点的概率正比于号节点当前的连接度，而后在B类中产生新节点，以同样的方式向A类连边，当这两个步骤进行足够多次之后即可得到满足数据文件特点的网络。根据这里所提供的生成机制，发现节点连接度分布的表达式。 3．在这类网络上考虑“易感－感染－易感”（SIS）模型，得到较平稳时期的得病数量以及A类和B类的得病比例。(参数γ＝0.1, 考虑到两类个体的感染率可以不同，分析中假定A类个体的感染率为B个体感染率的2倍，即=2，并分别取B类个体的感染率=0.01,0.02,0.0)。考虑PC机的计算速度，模拟时网络规模不要太大，可选择00＋00的二部网络。 4．对的模型进行理论的分析，看看是否和的模拟结果一致。病毒传播问题的研究由来已久，而一再的病毒流行使得这一领域长期以来吸引着人们的注意。号节点；再从B类节点中出发，连接到A类中的节点。通过不断地重复这个步骤，找出节点连接度分布的表达式。 对于问题三：利用问题二的产生机制，在A类中加入新的节点，先判断它是否患病，再判断与之相连的B类节点是否患病，在相连的基础上，判断它们能否能够传染。再从B类节点出发，连接到A类中的节点。通过不断地重复这个步骤，得到较平稳时期的得病数量以及A类和B类的得病比例、 A类中第号节点的连接度 B类中第号节点的连接度 、 A类中第号节点出现的概率 B类中第号节点出现的概率 治愈率 A个体的被传染率 B个体的被传染率 具有k个连接度的A类个体的感染密度 具有k个连接度的B类个体的感染密度 节点连接度k 任意一条边指向感染的A类节点的概率 任意一条边指向感染的B类节点的概率 节点的平均连接度 A类节点连接度的出现频数 B类节点连接度的出现频数 A类个体第i号节点的状态 B类个体第i号节点的状态 0 易感者 1 感染者 模型建立 问题一：节点连接度和其出现频率的关系 问题一的思想： 因为这是一个数据处理的过程，所以我们先是对数据进行处理，得到A、B类节点连接度、及其出现概率、的表格。为了发现它们之间服从什么分布，我们先是用Matlab软件中的cftool指令拟合出表格中的数据服从的曲线，然后再对这条曲线进行验证，是否可以作为概率分布的曲线，最终得出它们服从的分布。具体步骤如下： 步骤一：数据处理 由附件中的数据，我们可以整理出A、B两类节点连接度、与其出现频率、的数据，如表1和表2 表1 A类节点连接度及其出现频率 A类节点连接度 出现频率 A类节点连接度 出现频率 A类节点连接度 出现频率 1 0.6626 17 0.0010 33 0.0001 2 0.1679 18 0.0009 35 0.0001 3 0.0703 19 0.0009 37 0.0001 4 0.0316 20 0.0003 38 0.0002 5 0.0185 21 0.0007 39 0.0003 6 0.0109 22 0.0002 41 0.0001