首页
所有分类
TAGS
文字工具 EMOJI BIBLE
编程 Python

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实现病毒仿真器的方法示例(附demo)

- Author -

HL是限量版

声明：登载此文出于传递更多信息之目的，并不意味着赞同其观点或证实其描述。

Python 相关文章推荐
Python抓取淘宝下拉框关键词的方法
Jul 08 Python
Python中shutil模块的常用文件操作函数用法示例
Jul 05 Python
Python实现的中国剩余定理算法示例
Aug 05 Python
对python字典元素的添加与修改方法详解
Jul 06 Python
python绘制多个曲线的折线图
Mar 23 Python
Python定时发送天气预报邮件代码实例
Sep 09 Python
Python拆分大型CSV文件代码实例
Oct 07 Python
python定义类self用法实例解析
Jan 22 Python
pytorch 模型的train模式与eval模式实例
Feb 20 Python
Python利用FFT进行简单滤波的实现
Feb 26 Python
python+opencv实现车道线检测
Feb 19 Python
python图像处理 PIL Image操作实例
Apr 09 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
redis shake(1) 速龙(1) 分析统计(1) ipad(1) 排名(1) 窗口函数(3) Tomcat(18) win11(40) 数据同步(1) i5 10400(1)
You might like
PHP4实际应用经验篇(8)
2006/10/09 PHP
PHP生成静态页
2006/11/25 PHP
隐性调用php程序的方法
2009/03/09 PHP
Yii数据模型中rules类验证器用法分析
2016/07/15 PHP
Laravel框架Request、Response及Session操作示例
2019/05/06 PHP
用PHP做了一个领取优惠券活动的示例代码
2019/07/05 PHP
PHP中类与对象功能、用法实例解读
2020/03/27 PHP
JavaScript 编程引入命名空间的方法与代码
2007/08/13 Javascript
JS 分号引起的一段调试问题
2009/06/18 Javascript
理解JavaScript变量作用域更轻松
2009/10/25 Javascript
分享十五个最佳jQuery 幻灯插件和教程
2010/03/27 Javascript
判断ie的两种简单方法
2013/08/12 Javascript
HTML页面登录时的JS验证方法
2014/05/28 Javascript
JavaScript字符串对象replace方法实例（用于字符串替换或正则替换）
2014/10/16 Javascript
JavaScript实现将xml转换成html table表格的方法
2015/04/17 Javascript
Node.js中Request模块处理HTTP协议请求的基本使用教程
2016/03/31 Javascript
Ext JS框架中日期函数的用法及日期选择控件的实现
2016/05/21 Javascript
详解Vue.js动态绑定class
2016/12/20 Javascript
JS实现动态添加DOM节点和事件的方法示例
2017/04/28 Javascript
javascript删除数组元素的七个方法示例
2019/09/09 Javascript
vue组件 keep-alive 和 transition 使用详解
2019/10/11 Javascript
vue中echarts引入中国地图的案例
2020/07/28 Javascript
python实现在每个独立进程中运行一个函数的方法
2015/04/23 Python
关于Django外键赋值问题详解
2017/08/13 Python
python 矩阵增加一行或一列的实例
2018/04/04 Python
Selenium 模拟浏览器动态加载页面的实现方法
2018/05/16 Python
Python实现基于C/S架构的聊天室功能详解
2018/07/07 Python
Pyqt5实现英文学习词典
2019/06/24 Python
Python实现平行坐标图的两种方法小结
2019/07/04 Python
100行Python代码实现每天不同时间段定时给女友发消息
2019/09/27 Python
Python 抓取数据存储到Redis中的操作
2020/07/16 Python
python3爬虫中多线程进行解锁操作实例
2020/11/25 Python
python 制作简单的音乐播放器
2020/11/25 Python
2015年计生协会工作总结
2015/04/24 职场文书
运动会加油稿
2015/07/22 职场文书
2015年度对口支援工作总结
2015/07/22 职场文书