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内置的字符串处理函数详细整理(覆盖日常所用)
Aug 19 Python
Python随机生成信用卡卡号的实现方法
May 14 Python
详解python进行mp3格式判断
Dec 23 Python
浅析python递归函数和河内塔问题
Apr 18 Python
python使用mysql数据库示例代码
May 21 Python
Python基于生成器迭代实现的八皇后问题示例
May 23 Python
解决Python找不到ssl模块问题 No module named _ssl的方法
Apr 29 Python
opencv转换颜色空间更改图片背景
Aug 20 Python
python对象转字典的两种实现方式示例
Nov 07 Python
Python MOCK SERVER moco模拟接口测试过程解析
Apr 13 Python
Python爬虫爬取糗事百科段子实例分享
Jul 31 Python
pip/anaconda修改镜像源,加快python模块安装速度的操作
Mar 04 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读取txt文件的内容并赋值给数组的代码
2011/11/03 PHP
利用PHPStorm如何开发Laravel应用详解
2017/08/30 PHP
PHP中如何使用Redis接管文件存储Session详解
2018/11/28 PHP
PHP赋值的内部是如何跑的详解
2019/01/13 PHP
让div层随鼠标移动的实现代码 ie ff
2009/12/18 Javascript
原生javascript获取元素样式
2014/12/31 Javascript
jQuery实现鼠标划过展示大图的方法
2015/03/09 Javascript
js实现的四级左侧网站分类菜单实例
2015/05/06 Javascript
jQuery表格插件datatables用法汇总
2016/03/29 Javascript
AngularJS bootstrap启动详解及实例代码
2016/09/14 Javascript
js实现网页定位导航功能
2017/03/07 Javascript
浅谈JavaScript的innerWidth与innerHeight
2017/10/12 Javascript
javascript+html5+css3自定义弹出窗口效果
2017/10/26 Javascript
node.js基础知识小结
2018/02/26 Javascript
JS 中可以提升幸福度的小技巧(可以识别更多另类写法)
2018/07/28 Javascript
解决layer弹出层中表单不起作用的问题
2019/09/09 Javascript
vue中使用element ui的弹窗与echarts之间的问题详解
2019/10/25 Javascript
ES6 Symbol在对象中的作用实例分析
2020/06/06 Javascript
[06:04]DOTA2英雄梦之声Vol19卓尔游侠
2014/06/20 DOTA
[00:34]TI7不朽珍藏III——地穴编织者不朽展示
2017/07/15 DOTA
浅谈python中的正则表达式(re模块)
2017/10/17 Python
Python File(文件) 方法整理
2019/02/18 Python
python 两种方法修改文件的创建时间、修改时间、访问时间
2020/09/26 Python
施华洛世奇巴西官网:SWAROVSKI巴西
2019/12/03 全球购物
瑞士最大的图书贸易公司:Orell Füssli
2019/12/28 全球购物
建筑个人求职信范文
2014/01/25 职场文书
护理专业自我鉴定
2014/01/30 职场文书
售后服务承诺书范文
2014/03/26 职场文书
骨干教师培训方案
2014/05/06 职场文书
投资入股合作协议书
2014/10/28 职场文书
自主招生自荐信怎么写
2015/03/24 职场文书
优质服务标语口号
2015/12/26 职场文书
干部理论学习心得体会
2016/01/21 职场文书
幼儿园教师教学反思
2016/03/02 职场文书
幽默导游词应该怎么写?
2019/08/26 职场文书
frg-100简单操作(设置)说明
2022/04/05 无线电