冠状病毒传播仿真器的原理和实现（Python版)

本文主要是介绍冠状病毒传播仿真器的原理和实现（Python版)，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

李宁老师免费视频已经发布了，扫码看免费视频
在这里插入图片描述
或点击：免费视频
本文摘要：

本文首先会解释一下到底什么是"冠状病毒"，以及杀死"冠状病毒"的方法。然后会利用Python实现一个"冠状病毒"传播仿真器，来演示一下为何“不出门“ +“疯狂建医院”会间接杀死病毒（动态模拟了从发生疫情，到疫情结束的整个过程）。以及如果控制不好，会有什么后果（一定是很严重的）。

相信大家现在一定和我一样，无比郁闷。好不容易盼到过年了，买了一大堆好吃的，结果哪也去不了，只能在家消耗这些美食了。而且很多人宅在家里的每一天都会做同一个伟大的计划！！

在这里插入图片描述
那么我们为什么不能出门呢？答案大家都已经知道了，这就是一种被称为“冠状病毒”的东西在捣乱。据钟南山院士和其他一些专家推测，这种“冠状病毒”的潜伏期是14天左右，所以需要尽可能避免在14天内接触更多的人。而宅在家里是最好的选择。不过这是头一次放这么长的假，而且还必须呆在家里，也是相当郁闷，恍惚间竟想起了猴哥，500年你是怎么熬过来的呢！
在这里插入图片描述
下面就先来了解下到底什么是冠状病毒！

1. 什么是冠状病毒

如果详细解释什么是冠状病毒，先要从病毒讲起，估计要写一本书了，不过完全没有必要。大家只要知道一点就好。病毒不是生物（类生物，有生物的部分特征），不能独立生存，需要依赖于宿主，也就是人体的细胞。换句话说，病毒需要攻击人体的细胞才能生存。病毒可以利用细胞中的成分复制自身，从而达到繁殖的目的。如果受到感染的细胞过多，人体就会生病。

那么"冠状病毒"呢？当然也属于病毒的一种，不过它不是普通的病毒。首先看看"冠状病毒"的样子，病毒表面有很多凸起的冠（看起来还挺漂亮的），所以称为"冠状病毒"。

在这里插入图片描述
"冠状病毒"除了外形上与普通病毒有差异，在本质上与普通病毒也是不同的。普通病毒通常都是DNA病毒，也就是双螺旋结构。而"冠状病毒"属于RNA病毒，也就是单螺旋结构。双螺旋结构的DNA病毒更稳定，不容易变异。但RNA病毒只有一个链，非常不稳定，容易变异。这才是"冠状病毒"的杀手锏，通过变异抵抗人类的反击。目前"冠状病毒"是否变异，还很难说。

在这里插入图片描述

2. 如何杀死"冠状病毒"

现在的问题是，如何杀死"冠状病毒"。当然，最直接的方法是将"冠状病毒"赶出细胞，或直接从细胞中清除，不过目前人类的技术还做不到这一点。这应该属于比细胞医学更高端的医学：分子医学或原子医学。目前人类的科技水平甚至还没完全达到细胞医学的高度（最多0.7个细胞医学），因为癌症等细胞疾病还无法完全治愈。

既然病毒目前还无法直接杀死，那么这么多痊愈的病人是如何做到的呢？其实可以采用如下两种方法间接杀死"冠状病毒"，也包括其他大多数病毒。

利用人体内的抗体细胞（如T细胞），将被感染的细胞连同"冠状病毒"一起消灭
干扰"冠状病毒"在细胞内的复制过程（例如，RNA干扰），让"冠状病毒"无法正常复制自身，也就是让"冠状病毒"没有后代，这样"冠状病毒"就会由于自身的生命终结而死亡。

目前绝大多数痊愈的病人属于第一种情况，也就是通过自身的抗体细胞（如T细胞）检测人体内被感染的细胞，然后通知这些被感染的细胞启动自毁程序，当这些细胞被销毁后，那么病毒也就被消灭了。

不过可能有的同学会问，既然病毒可以被抗体细胞消灭，那么这么多医护人员和医疗设备岂不是多余了。其实并不是多于的，而是非常必要。这是因为病毒在攻击正常细胞的同时，也让人体的抵抗力开始下降，人体内的抗体细胞是有限的，这些抗体细胞会到处救火，顾不上对付"冠状病毒"了。正是由于这些医疗器械（如呼吸机），尽可能让人体机能恢复到接近正常人的水平，这样各种抗体细胞就可以集中力量对付"冠状病毒"了。所以说，这些痊愈的病人其实是通过自身的抗体细胞消灭了"冠状病毒"，而医疗器械、各种药物，医务人员的护理，其实是抗体细胞的援军。另外，由于不同人的体质不同。极个别的个体，抗体细胞非常强大，造成了"冠状病毒"无法快速复制（但仍然在复制，只是增速变慢），这也是为什么有的个体的潜伏期会超过14天的原因，但这毕竟是极少数。潜伏期是被病毒感染的细胞数量达到足以致病的时间（人体内如果只有少量的被病毒感染的细胞，是不会表现出任何症状的）。

还有就是，为什么"冠状病毒"感染者在发病时大多会有肺炎的症状，而且伴随着发烧呢？其实这就是抗体细胞在和"冠状病毒"进行较量呢！进医院隔离，是为了增强抗体细胞对抗"冠状病毒"的筹码。不过由于一些人的抗体细胞的战斗力实在太差（可能有部分人是因为年龄太大的原因），所以就算抗体细胞等来了援军，也于事无补，这些就是已经死亡的被感染者（大多是60岁以上的老年人）。所以等疫情结束后，好好锻炼身体吧，拥有强壮的身体，不能保证你不得病，但至少可以增加得病后活下来的几率。

在这里插入图片描述对于目前正在研发的抗"冠状病毒"的药物主要是通过第2种方式消灭病毒的。也就是干扰病毒的复制过程，不过很遗憾，到现在为止，还没有被证明非常有效的药物可以做到这一点。

3. 对付"冠状病毒"的手段

从生物学角度，我们已经了解了"冠状病毒"的发病原理，但在现实中，如何操作呢？

其实对付"冠状病毒"以及其他大多数病毒，基本上就是基于12个字：有症状赶快治，没症状要隔离。这也是国际上通用的原则。

前6个字容易理解，有症状了，就直接进医院了。如果没症状呢？没症状有两种情况：疑似和正常人。疑似主要是指与被感染者近距离接触，或从外地来本地的人员（输入者），由于"冠状病毒"的潜伏期是14天左右，所以这些疑似者至少需要被隔离14天才会确定是否真的被感染。而正常人只要没和被感染者近距离接触，就不会被感染。这些人之所以也需要隔离，是因为怕被别人感染。不过这里的隔离通常是在自己的家中，不出门。通过隔离，可以大幅度减少病毒感染新的人群，也就是让存量不再增加或少增加。而还有很多被感染者，这些人就需要在医院里接受治疗了。不过由于被感染者太多，所以武汉等地区快速建起了很多临时医院。这是用来减少存量的。当存量不但不会增加、而且在不断减少，直到被感染者为0，疑似者为0时，疫情才会彻底结束，这也是本文要介绍的病毒扩散仿真器的基本原理。

4. 用病毒扩散仿真器来演示病毒扩散和疫情结束的全过程

在实现这个仿真器之前，先来演示下这个仿真器。

仿真器可以对多个数据进行模拟，包括健康者人数、潜伏期人数、发病者人数、已经隔离的人数、已经死亡的人数、空余床位、继续床位、病毒传播率、病毒潜伏期、医院收治响应时间、医院当前床位、安全距离、平均流动意向。

启动程序，会利用初始值进行模拟，初始发病人数为50人，市民总数为5000人。如下图所示。
在这里插入图片描述中间区域的若干个点表示各种状态的市民。白色的表示健康市民、黄色表示潜伏期市民、红色表示发病市民、黑色表示死亡的市民。右侧的竖条表示医院的床位，初始值是100。如果用参数值进行模拟，100张床位很快就会被填满，然后病毒在人群中就会大爆发，很快红点就会遍布人群，如下图所示

在这里插入图片描述当前天数已经显示过了31天（耽误了一个月），感染者已经接近1000了，这时政府开始采取紧急措施。主要有两个：封城和关闭娱乐场所、增加医院的床位。前者是为了避免感染更多的人，后者是为了消耗被感染者的存量。所以通过下面的设置来调整参数。例如，将流动意向调整为-1.71。并且增加床位323个。

在这里插入图片描述这时总床位数变成了423。这里的流动意向在-3和3之间，如果是3，表示市民的活动意愿非常强烈，例如，正好是春节时期，市民逛商场，聚会非常频繁。流动意愿越小，流动意愿就越低。这里调成-1.76，表示市民不能参加聚会、不能出城、出门需要戴口罩，但市民仍然可以在市内流动。流动意愿远低于春节正常的值。不过尽管政府采取了一定的措施，但由于是在疫情开始后一个月才采取了紧急措施，所以病毒已经扩散了，因此，疫情并没有得到非常明显的缓解。如下图所示。主要表现为市民仍然可以自由活动（尽管不能参加聚会），仍然存在一定的感染风险。而且医院床位明显不足。

在这里插入图片描述为了更进一步控制疫情，政府开始封闭小区，更进一步限制人员的活动，以及军方开始干预，派出了大量的医护人员以及各种医疗设备，并且建立的多个方舱医院。医院床位得到了很大的缓解。这里将参数设置成最大值来模拟这一过程，增加床位1200个，流动意向设置为-3.0，也就是说基本上市民不流动了。如下图所示。

在这里插入图片描述这时看到床位已经增加到了1623，比急需的床位多了不少，而且人员趋于不流动，发病人数不断减少（都被送进了医院），而且潜伏期人数逐渐转换为发病人数，也被送进了医院，最终，潜伏期人数和发病者人数都是0，疫情结束，如下图所示。共耗费了60天。当然，实际情况没这么顺利。仿真器可以立刻增加医院床位数，可以立刻隔离人员，但在实际操作中，建立医院需要时间，隔离也需要协调，尤其是上千万人的大城市。

在这里插入图片描述不过只要能做到隔离和及时就医，冠状病毒疫情结束也只是时间问题。当然，这要在这两点做的比较好的情况下，如果处理失当，那么仿真器就会呈现下图的状态，完全失控，人类将面临一场浩劫。

在这里插入图片描述

5. 病毒传播仿真器的实现

现在来谈谈仿真器实现的原理。仿真器使用Python和PyQt5实现。PyQt5是封装了Qt library的跨平台GUI开发库，基于Python语言。

这里主要涉及到仿真器效果绘制，以及如何模拟多个参数。先来说一下绘制市民的状态。绘制的工作通过drawing.py文件的Drawing类来完成。该类是QWidget的子类，这也就意味着Drawing类本身是PyQt5的一个组件。与按钮、标签类似。只是并不需要往Drawing上放置任何子组件。只要在Drawing上绘制各种图形即可。

在PyQt5中，任何一个QWidget的子类，都可以实现一个paintEvent方法，当组件每次刷新时，就会调用paintEvent方法重新绘制组件的内容。Drawing类中paintEvent方法的代码如下：

    def paintEvent(self, event):qp = QPainter()qp.begin(self)# 绘制城市的各种状态的市民self.drawing(qp)qp.end()

在绘制图像前，需要创建QPainter对象，然后调用QPainter对象的begin方法，结束绘制后，需要调用QPainter对象的end方法。上面代码中的drawing方法用于完成具体的绘制工作。

仿真器可以模拟5000个市民的状态，所以需要用5000个小矩形来表示这5000个市民。也就是在drawing方法中需要绘制这5000个表示市民的小矩形。代码如下：

def drawing(self, event):... ...# 绘制代表市民的小矩形persons = Persons().personsif persons == None:returnnormal_person_count = 0latency_person_count = 0confirmed_person_count = 0freeze_person_count = 0death_person_count = 0# 扫描内一个人的状态for person in persons:if person.state == NORMAL:# 健康人qp.setPen(Qt.white)normal_person_count += 1elif person.state == LATENCY:# 潜伏期感染者qp.setPen(QColor(255,238,0))latency_person_count += 1elif person.state == CONFIRMED:# 确诊患者qp.setPen(Qt.red)confirmed_person_count += 1elif person.state == FREEZE:# 已隔离者qp.setPen(QColor(72, 255, 252))freeze_person_count += 1elif person.state == DEATH:# 死亡患者qp.setPen(Qt.black)death_person_count += 1person.update()   # 更新每一个人的状态bed_half_size = Hospital().bed_size // 2rect = QRect(person.x - bed_half_size, person.y - bed_half_size,Hospital().bed_size//2, Hospital().bed_size//2)brush = QBrush(Qt.SolidPattern)brush.setColor(qp.pen().color())qp.setBrush(brush)qp.drawRect(rect)... ...

在上面的代码中，通过 Persons对象的persons属性获取表示市民的对象（Person对象）列表。并在循环中根据Person对象的状态设置小矩形的颜色，以及分别统计不同人群的数量，这些数量会显示在仿真器右侧的组件中。最后，使用drawRect方法绘制表示每一个市民的小矩形。这样就绘制了当前状态的5000个市民。

当然，这些状态要不断更新。这里使用线程每100毫秒刷新一次，这些功能在refresh.py文件的Refresh类中，代码如下：

from PyQt5.QtCore import *
from params import  *class Refresh(QThread):def __init__(self, drawing):super(Refresh, self).__init__()self.drawing = drawingdef run(self):while not Params.success:try:QThread.msleep(100)# 刷新Drawingself.drawing.update()Params.current_time += 1except:pass

每次刷新Drawing，需要调用update方法，调用该方法后，Drawing就会调用自身的paintEvent方法重新绘制整个组件的内容。

在paintEvent方法中，还调用了Person对象的update方法，该方法是我们自己编写的，用于不断更新每一个人的状态，这些状态会根据多个参数进行协调。该方法属于Person类，代码如下：

    def update(self):# 如果已经隔离或者死亡了，就不需要处理了if self.state == FREEZE or self.state == DEATH:return# 处理已经确诊的感染者（即患者）if self.state == CONFIRMED and self.dead_time == 0:destiny = random.randrange(1,10001)  # 幸运数字，[1,10000]随机数if destiny >= 1 and destiny <= int(Params.fatality_rate * 10000):# 幸运数字落在死亡区间dt = int(sp.random.normal(Params.dead_time,Params.dead_variance))self.dead_time = self.confirmed_time + self.dead_timeelse:self.dead_time = -1   # 逃过了死神的魔爪if self.state == CONFIRMED and Params.current_time - self.confirmed_time >= Params.hospital_receive_time:# 如果患者已经确诊，且（世界时刻-确诊时刻）大于医院响应时间，即医院准备好病床了，可以抬走了bed = Hospital().pick_bed()  # 查找空床位if bed == None:# 没有空床位，报告需求床位数if not self.need_bed:Hospital().need_bed_count += 1self.need_bed = Trueelse:# 安置病人self.used_bed = bedself.state = FREEZEself.x = bed.x + Hospital().bed_size // 2self.y = bed.y +  Hospital().bed_size // 2if self.need_bed and Hospital().need_bed_count > 0:Hospital().need_bed_count -= 1bed.is_empty = False# 处理病死者if (self.state == CONFIRMED or self.state == FREEZE) and Params.current_time >= self.dead_time and self.dead_time > 0:self.state = DEATH                          # 患者死亡personpool.Persons().latency_persons.remove(self)      # 已经死亡，无法传染别人，需要从确诊者中删除Hospital().empty_bed(self.used_bed)         # 腾出床位if Hospital().need_bed_count > 0:Hospital().need_bed_count -= 1# 增加一个正态分布用于潜伏期内随机发病时间latency_symptom_time = sp.random.normal(Params.virus_latency / 2,25)# 处理发病的潜伏期感染者if Params.current_time - self.infected_time > latency_symptom_time and self.state == LATENCY:self.state = CONFIRMED                                  # 潜伏者发病self.confirmed_time = Params.current_time              # 刷新确诊时间# 处理未隔离者的移动问题self.action()# 处理健康人被感染的问题persons = personpool.Persons().persons# 不是健康人，返回if self.state >= LATENCY:return# 通过一个随机幸运值和安全距离决定感染其他人latency_persons = personpool.Persons().latency_persons.copy()for person in latency_persons:random_value = random.random()if random_value < Params.broad_rate and self.distance(person) < Params.safe_distance:self.be_infected()break

update方法主要就是根据在params.py中的各种参数变量，以及随机值，计算下一次状态中潜伏期人数、感染人数、被隔离人数等数据，并且在每次刷新页面时更新这些数据。

以上的描述就是如何绘制表示5000个市民的状态。右侧各种数据并不是绘制在页面上的，而是通过QtDesigner设计的右侧的界面，然后将Drawing对象作为标准的组件放在了主界面的左侧。设计界面如下图所示：

在这里插入图片描述然后通过pyuic将.ui文件生成.py文件，在程序中调用即可。这些组件的更新同样是在前面给出的drawing方法中。

另外，这个仿真器还提供了动态设置参数的功能。这是通过另外一个程序实现的，两个程序通过socket通讯。这个设置程序同样是通过QtDesigner设计的，设计界面如下图所示。

在这里插入图片描述在设置程序中，通过Transmission类的send_command方法向仿真器发布命令，例如，更新床位数的代码如下：

from PyQt5.QtWidgets import *
from socket import *
class Transmission:def __init__(self,ui):self.ui = uiself.host = 'localhost'self.port = 5678self.addr = (self.host, self.port)# 向仿真器发布命令def send_command(self, command, value = None):tcp_client_socket = socket(AF_INET, SOCK_STREAM)tcp_client_socket.connect(self.addr)if value == None:value = 0data = command + ':' + str(value)tcp_client_socket.send(('%s\r\n' % data).encode(encoding='utf-8'))data = tcp_client_socket.recv(1024)result = data.decode('utf-8').strip()tcp_client_socket.close()return result# 更新床位数def update_bed_count(self):print(self.ui.horizontalSliderBedCount.value())result = self.send_command('add_bed_count',self.ui.horizontalSliderBedCount.value())if result == 'ok':QMessageBox.information(self.ui.centralwidget, '消息', f'成功添加了{self.ui.horizontalSliderBedCount.value()}张床位', QMessageBox.Ok)

在仿真器端，通过Receiver以及TCPServer来接收设置程序发过来的命令，如果成功设置，返回ok。Receiver类以及相关的代码如下：

from socketserver import (TCPServer as TCP,StreamRequestHandler as SRH)
from common import *
from params import *
from hospital import *
from PyQt5.QtCore import *
import sys
# 响应客户端请求事件的类
class MyRequestHandler(SRH):def handle(self):# 读取客户端发送的数据data = str(self.rfile.readline(),'utf-8')index = data.find(':')command = data[:index]value = data[index + 1:]value = int(value)# 执行具体的命令if command == 'add_bed_count':Params.hospital_bed_count +=  valueHospital().free_bed_count = Hospital().free_bed_count + valueHospital().compute(value)elif command == 'set_flow_intention':Params.average_flow_intention = value / 100elif command == 'set_broad_rate':Params.broad_rate = value / 100elif command == 'set_latency':Params.virus_latency = value * 10elif command == 'close':Params.app.quit()self.wfile.write(b'ok\r\n')# 在线程中监听客户端的请求
class Receiver(QThread):tcp_server = Nonedef __init__(self):super(Receiver,self).__init__()self.host = ''self.port = 5678self.addr = (self.host,self.port)Receiver.tcp_server = TCP(self.addr, MyRequestHandler)def run(self):Receiver.tcp_server.serve_forever()