Python实现病毒仿真器的方法示例(附demo)


Posted in Python onFebruary 19, 2020

最近新冠在神州大陆横行,全国上下一心抗击疫情。作为一枚程序员,我也希望可以为抗击疫情做出自己的贡献,钟院士一直劝说大家不要出门,减少人口间的流动。对此,我特意做了一个病毒仿真器,探询冠状病毒传播。

1. 仿真效果

仿真开始,一开始只有5个发病者,传播率为0.8,潜伏期为14天

Python实现病毒仿真器的方法示例(附demo)

由于人口的流动,以及医院床位的隔离,一开始病毒扩撒不是很速度

Python实现病毒仿真器的方法示例(附demo)

随着医院床位满了,隔离失败,加上人口的流动,病患数开始几何式的增加

Python实现病毒仿真器的方法示例(附demo)

2. 什么是仿真器

仿真器(emulator)以某一系统复现另一系统的功能。与计算机模拟系统(Computer Simulation)的区别在于,仿真器致力于模仿系统的外在表现、行为,而不是模拟系统的抽象模型。

3. python如何实现仿真器

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

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

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

# 每次Drawing刷新,都会调用该方法
  def paintEvent(self,event):
    qp = QPainter()
    qp.begin(self)
    self.drawing(qp) # 绘制城市的各种状态的市民
    qp.end()

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

仿真器设置一个人员池,来进行人员的流动性

class Persons(metaclass=Singleton):
  def __init__(self):
    self.persons = []          # 保存所有的人员
    self.latency_persons = []      # 保存处于潜伏期的人员
 
    city = City(Params.city_center_x,Params.city_center_y)
 
    for value in range(0, Params.city_person_count):
      x = Params.person_position_scale * next_gaussian() + city.center_x
      y = Params.person_position_scale * next_gaussian() + city.center_y
      if x > Params.city_width:
        x = Params.city_width
      if y > Params.city_height:
        y = Params.city_height
      self.persons.append(Person(city,x,y))
  # 获取特定人群的数量
  def get_person_size(self,state):
    if state == -1:
      return len(self.persons)
    count = 0
    for person in self.persons:
      if person.state == state:
        count += 1
    return count

仿真代码

from PyQt5.QtWidgets import *
from socket import *
 
class Transmission:
  def __init__(self,ui):
    self.ui = ui
    self.host = 'localhost'
    self.port = 5678
    self.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 = 0
 
    data = 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 setup(self):
    self.ui.horizontalSliderBedCount.valueChanged.connect(self.bed_count_value_change)
    self.ui.pushButtonUpdateBedCount.clicked.connect(self.update_bed_count)
 
    self.ui.horizontalSliderFlowIntention.valueChanged.connect(self.flow_intention_value_change)
    self.ui.pushButtonFlowIntention.clicked.connect(self.update_flow_intention)
 
    self.ui.horizontalSliderBroadRate.valueChanged.connect(self.broad_rate_value_change)
    self.ui.pushButtonBroadRate.clicked.connect(self.update_broad_rate)
 
    self.ui.horizontalSliderLatency.valueChanged.connect(self.latency_value_change)
    self.ui.pushButtonLatency.clicked.connect(self.update_latency)
 
    self.ui.pushButtonClose.clicked.connect(self.close_virus_simulation)
 
  def bed_count_value_change(self):
    self.ui.labelBedIncrement.setText(
      f'<html><head/><body><p><span style=\" font-size:24pt; color:#0000ff;\">{self.ui.horizontalSliderBedCount.value()}</span></p></body></html>')
 
  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)
 
  def flow_intention_value_change(self):
    self.ui.labelFlowIntention.setText(
      f'<html><head/><body><p><span style=\" font-size:24pt; color:#fc02ff;\">{self.ui.horizontalSliderFlowIntention.value() / 100}</span></p></body></html>')
 
  def update_flow_intention(self):
    result = self.send_command('set_flow_intention', self.ui.horizontalSliderFlowIntention.value())
    if result == 'ok':
      QMessageBox.information(self.ui.centralwidget, '消息',
                  f'成功设置流动意向为{self.ui.horizontalSliderFlowIntention.value() / 100}', QMessageBox.Ok)
 
  def broad_rate_value_change(self):
    self.ui.labelBroadRate.setText(
      f'<html><head/><body><p><span style=\" font-size:24pt; color:#fc0107;\">{self.ui.horizontalSliderBroadRate.value() / 100}</span></p></body></html>')
 
  def update_broad_rate(self):
    result = self.send_command('set_broad_rate', self.ui.horizontalSliderBroadRate.value())
    if result == 'ok':
      QMessageBox.information(self.ui.centralwidget, '消息',
                  f'成功设置传播率为{self.ui.horizontalSliderBroadRate.value() / 100}', QMessageBox.Ok)
 
  def latency_value_change(self):
    self.ui.labelLatency.setText(
      f'<html><head/><body><p><span style=\" font-size:24pt; color:#ffff0a;\">{self.ui.horizontalSliderLatency.value()}</span></p></body></html>')
 
  def update_latency(self):
    result = self.send_command('set_latency', self.ui.horizontalSliderLatency.value())
    if result == 'ok':
      QMessageBox.information(self.ui.centralwidget, '消息', f'成功设置传播率为{self.ui.horizontalSliderLatency.value()}',
                  QMessageBox.Ok)
 
  def close_virus_simulation(self):
    reply = QMessageBox.information(self.ui.centralwidget, "请问", "是否真的要关闭病毒扩散仿真器?",
                    QMessageBox.Yes | QMessageBox.No, QMessageBox.No)
    if reply == QMessageBox.Yes:
      result = self.send_command('close')
      if result == 'ok':
        QMessageBox.information(self.ui.centralwidget, '消息', '已经成功关闭病毒扩散仿真器', QMessageBox.Ok)

设置变量

class Params:
  success = False
 
  # 初始感染人数
  original_infected_count = 50
 
  # 病毒传播率
  broad_rate = 0.8
 
  # 病毒潜伏期,14天
  virus_latency = 140
 
  # 医院收治响应时间
  hospital_receive_time = 10
 
  # 医院床位
  hospital_bed_count = 100
 
  # 安全距离
  safe_distance = 2
 
  # 平均流动意向[-3,3]  值越大,流动意向越强
  average_flow_intention = 3
 
  # 城市总人口数量
  city_person_count = 5000
 
  # 病死率
  fatality_rate = 0.02
 
  # 死亡时间
  dead_time = 30
  # 死亡时间方差
  dead_variance = 30
 
  # 城市宽度
  city_width = 1100
 
  # 城市高度
  city_height = 800
 
  # 医院宽度(需要计算获得)
  hospial_width = 0
 
  # 城市中心x坐标
  city_center_x = 550
 
  # 城市中心y坐标
  city_center_y = 400
 
  # 用于计算城市中每个人随机位置的scale(用于正态分布)
  person_position_scale = 200
 
  current_time = 1 # 当前时间

我们可以根据这些修改仿真页面

Python实现病毒仿真器的方法示例(附demo) 

以上就是Python实现病毒仿真器的方法示例的详细内容,更多关于请关注三水点靠木其它相关文章!

Python 相关文章推荐
Python中的Numpy入门教程
Apr 26 Python
用PyQt进行Python图形界面的程序的开发的入门指引
Apr 14 Python
Python中集合的内建函数和内建方法学习教程
Aug 19 Python
Python下载网络文本数据到本地内存的四种实现方法示例
Feb 05 Python
Python生成短uuid的方法实例详解
May 29 Python
自定义django admin model表单提交的例子
Aug 23 Python
python rsa实现数据加密和解密、签名加密和验签功能
Sep 18 Python
pytorch实现对输入超过三通道的数据进行训练
Jan 15 Python
python去除删除数据中\u0000\u0001等unicode字符串的代码
Mar 06 Python
python两种注释用法的示例
Oct 09 Python
Pytest中skip和skipif的具体使用方法
Jun 30 Python
使用Python拟合函数曲线
Apr 14 Python
python内打印变量之%和f的实例
Feb 19 #Python
Python 实现日志同时输出到屏幕和文件
Feb 19 #Python
python 控制台单行刷新,多行刷新实例
Feb 19 #Python
python tqdm 实现滚动条不上下滚动代码(保持一行内滚动)
Feb 19 #Python
python 解决tqdm模块不能单行显示的问题
Feb 19 #Python
python 实现在shell窗口中编写print不向屏幕输出
Feb 19 #Python
Python换行与不换行的输出实例
Feb 19 #Python
You might like
PHP自定义函数实现assign()数组分配到模板及extract()变量分配到模板功能示例
2018/05/23 PHP
用js怎么把&amp;字符换成&quot;&amp;amp:&quot;
2006/10/19 Javascript
js实现连个数字相加而不是拼接的方法
2014/02/23 Javascript
jQuery插件开发详细教程
2014/06/06 Javascript
Node.js插件的正确编写方式
2014/08/03 Javascript
javascript学习笔记(三)BOM和DOM详解
2014/09/30 Javascript
jQuery实现锚点scoll效果实例分析
2015/03/10 Javascript
JS实现3D图片旋转展示效果代码
2015/09/22 Javascript
javascript实现自动填写表单实例简析
2015/12/02 Javascript
详细分析Javascript中创建对象的四种方式
2016/08/17 Javascript
vuejs2.0实现一个简单的分页示例
2017/02/22 Javascript
利用Node.js对文件进行重命名
2017/03/12 Javascript
浅谈vue同一页面中拥有两个表单时,的验证问题
2018/09/18 Javascript
通过实例学习React中事件节流防抖
2019/06/17 Javascript
JavaScript 作用域实例分析
2019/10/02 Javascript
JS将指定的某个字符全部转换为其他字符实例代码
2020/10/13 Javascript
[01:00:17]DOTA2-DPC中国联赛 正赛 SAG vs Dynasty BO3 第二场 1月25日
2021/03/11 DOTA
Python解析网页源代码中的115网盘链接实例
2014/09/30 Python
python实现简单点对点(p2p)聊天
2017/09/13 Python
利用Python如何批量更新服务器文件
2018/07/29 Python
python实现梯度下降算法
2020/03/24 Python
Python中的单下划线和双下划线使用场景详解
2019/09/09 Python
django框架使用views.py的函数对表进行增删改查内容操作详解【models.py中表的创建、views.py中函数的使用,基于对象的跨表查询】
2019/12/12 Python
Python基于Dlib的人脸识别系统的实现
2020/02/26 Python
python使用PIL剪切和拼接图片
2020/03/23 Python
Python使用内置函数setattr设置对象的属性值
2020/10/16 Python
美国零售商店:Blue&Cream
2017/04/07 全球购物
微软巴西官方网站:Microsoft Brasil
2019/09/26 全球购物
Oral-B荷兰:牙医最推荐的品牌
2020/02/25 全球购物
介绍一下Cookie和Session及他们之间的区别
2012/11/20 面试题
在Java开发中如何选择使用哪种集合类
2016/08/09 面试题
党员组织关系介绍信
2014/02/13 职场文书
团代会闭幕词
2015/01/28 职场文书
求职自我评价怎么写
2015/03/09 职场文书
多属性、多分类MySQL模式设计
2021/04/05 MySQL
Oracle数据库中通用的函数实例详解
2022/03/25 Oracle