ncov2019模拟

巫师养成中,看看我是否适合做占卜……

  • 建立新型冠状病毒肺炎预警模型,对温州市和浙江省的确诊人数进行了模拟
  • 相关报告已经在2月10日递交
  • 数据最后更新:2020年2月17日
  • 请用电脑打开以下链接,点击右上角右数第二个按钮compare data on hover以方便观看模拟数据
  • 不保证准确性
  • 个人作图水平还在提高中……枯

新版链接(链接地址会更新,所以请从本网页获取最新链接QVQ)

ncov温州市模拟

ncov浙江省模拟

日志

  • 2月9日:确立模型并模拟
  • 2月10日:上交模型报告(10页)
  • 2月16日:矫正参数
    浙江省和温州市的趋势均明显变缓,尤其是浙江省。尽量保守估计,调整参数使预测曲线尽量在“上面贴近”实际曲线
  • 2月17日:上报了一段话……其实就接下来的新增人数来讲,更倾向于这是一个随机数了,模型已经不再适用。

近日,浙江省和温州市的趋势均明显变缓,尤其是浙江省。我们对此进行了参数调整,但采取尽量保守估计的方式,使预测曲线尽量在“上面贴近”实际曲线。
…给出链接
以上是最新的模拟结果,根据保守估计,温州市和浙江省将在2月底达到平稳,并估计上限分别不超过540人和1230人。

  • 2月19日:我的任务已经结束了。浙江省近两日各一例,分别来自温州和嘉兴。个人认为强隔离政策在其中发挥了巨大的作用,以至于可以在模型后期阶段增加人数骤减,这是模型没有考虑到的地方。当然这本身也难以光从数据上做出预测,但是可以在数据发生变化时及时判断这一点的来临。

    强隔离政策是以损失经济为代价的,是否有采取强隔离的必要以及何时开始采取强隔离政策是值得思考的问题。另一方面,强隔离政策其实也有不人性的一面,这一点我是从自己到社区做防疫志愿者的时候体会到的,一些老父老母在医院里需要子女照看送饭,却被强制两天出一次门;一些在此期间工作的人晚上10点后没回到社区就只能在外过夜等等。当然也许这和奋战在一线的人来说都谈不上什么,但是家家都有本难念的经,愿这样的疫情不再发生,愿人类与自然和谐共处。

    接下来浙江省开始复工,强隔离解除,人员流动增加,因此相应政策的跟进和个人的警惕和防范是必不可少的。希望浙江省的疫情能顺利结束。


下面简单地轻松地说一下我的方法……

我主要是用微分方程组进行了建模,然后通过计算机模拟的方法进行模拟和预测,我将我的模型命名为SICRLD模型。以温州市作为浙江省典型代表,对模型进行了训(调)练(参),做了一个温州市(我的家乡太了不起太不容易了,为家乡点赞)预警模型然后直接放到浙江省进行模拟,最后做了一个治愈情况的模(估)拟(计)。

模型的数学描述

利用微分方程组刻画传染病发展过程,从而模拟并预测疫情人数。根据实际情况,考虑如下变量,总人数,易感者,感染未隔离者,确诊感染并隔离者,死亡/治愈者,武汉及周边流入者,未确诊的感染隔离者. 这六个变量均为累计值且为的函数。

的单位是天,六个变量的单位是人,这六个变量的关系如下:

  • 对于易感者,若感染未隔离者每人每天与人接触,并使其中的易感者受到感染,则每天新增感染数量为,另外有人从武汉及周边流入易感者群体。并且由于人民意识逐渐提高以及政策的控制,有未感染的易感者群体选择自行居家隔离。
  • 对于感染未隔离者,有人被感染但未隔离归入,有武汉及周边流入的人感染但未隔离,同时有人被确诊感染并隔离。
  • 对于未确诊的感染隔离者,在即时被感染后中有人被隔离但未确诊归入,另外中有人感染未确诊直接被隔离进入,同时人被确诊归入.
  • 对于确诊感染并隔离者,日新增为上述另四类群体中确诊的和。
  • 对于死亡/治愈人数,其新增人数为确诊人数与已经治愈的人数的差值即在院人数的倍。

容易列出这个微分方程组,略。

计算机模拟

那个看起来好像还能装点门面的微分方程组的显式解是不易求得的,可以用计算机模拟的方法对它进行离散时间的模拟。事实上它可以看作微分形式的状态转移函数,置模拟时间的初始时刻和初始值,给出当前模拟时间输入变量的值后,根据状态转移函数确定下一时刻各变量的值。

更具体地说,例如对于方程组第一行式子,有



就是这样简单,如此如此就能用每一天的数据模拟第二天的数据辽。