我相信你们现在都已经看过华盛顿邮报的文章了COVID-19的模拟，谈到了在面对高度传染性疾病时，社会隔离对公众健康的好处。如果你没有，这篇文章:https://www.washingtonpost.com/graphics/2020/world/corona-simulator/
出于对隔离的厌倦，我在MATLAB中重新创建了这个仿真，有更多的旋钮可以转动，有更多的参数可以使用。例如，如果一些携带者可能死亡怎么办?如果是埃博拉病毒(~50%的死亡率)怎么办?
此功能模拟疾病传播的一组N个载波，在密闭空间内，与载体的比例一些社会隔离本身之中。一个简单的多体物理模型（两个相等质量的粒子之间弹性碰撞）确定载体的轨迹。
希望这能让你从社交孤立的郁闷中得到一个有趣的(虽然短暂的)喘息!
引用作为
约书亚Gafford (2020)。Simulitis（冠状模拟）（//www.tatmou.com/matlabcentral/fileexchange/74610-simulitis-a-coronavirus-simulation），MATLAB中央文件交换。检索。
评论和评分（21)
更新
| 1.0.7 | 删除不必要的刷新数据调用以加速执行。 |
|
| 1.0.6 | 修改仿真所以后面检疫墙只有运营商可以开始感染。 |
|
| 1.0.5 | 校正后的输入解析检疫栅极功能 |
|
| 1.0.4 | 添加隔离功能。 |
|
| 1.0.3 | 更新标签 |
|
| 1.0.2中 | 修正了矢量大小不一致的问题和参数检查改进 |
|
| 1.0.1 | 删除注释掉的代码 |
MATLAB版本兼容性
创建与R2018b
兼容R2018b的任何版本
哈维尔（查看资料）
约书亚Gafford（查看资料）
@ James -你是对的!我认为refreshdata调用是将新数据推送到坐标轴句柄所必需的，但是我想drawnow会像您所说的那样处理这个问题。现在跑得快多了!
詹姆斯•史密斯（查看资料）
太酷了!
你需要那两个数据更新电话吗?他们吃了很多时间。我把它们注释掉了，看起来运行得很好，速度也快了很多。我认为“drawnow”负责更新数据。
Mostafa Essuri（查看资料）
@Joshua,谢谢。
现在我正在检查它。
Tsotne Marghia（查看资料）
@Joshua：非常感谢您对这种快速反应。只是要注意：它可能是把一个方向的轨迹成一个随机游走是一个好主意，所以在每一个步骤，将选择一个新的方向向量COS（2 * PI *兰特）。
约书亚Gafford（查看资料）
@迈克：看起来不错，完美的我跑了！
约书亚Gafford（查看资料）
你是完全正确的。修改了代码，所以最初的感染率只适用于隔离墙后面(左边)的人。谢谢!
苏菲Kargo（查看资料）
Tsotne Marghia（查看资料）
@Joshua:感谢您对代码的回复和更新。然而，一旦我运行他们，我注意到一些感染是在大门的两侧。我原以为隔离门的目的就是把被感染的人隔离在隔离区的大门后面。
迈克fienen（查看资料）
@Joshua：这里的Python代码。https://github.com/mnfienen/pysimulitis渴望任何反馈。
约书亚Gafford（查看资料）
@迈克:太棒了!如果您碰巧有github存储库，希望看到代码吗
@穆斯塔法：新添加隔离门功能，检查出来。也许不是最强大的/有效的方式，因为我时间有限，但它的作品。可以通过分别w_loc设置的隔离门的位置，并通过w_init和w_speed控制打开时间和打开速度。
迈克fienen（查看资料）
这太酷了!谢谢你的帖子。我冒昧地将其翻译成Python(带有属性)。非常及时的工作和良好的编码。谢谢你这么做!
Mostafa Essuri（查看资料）
谢谢你的工作。
如何在仿真期间使用打开的门来包含一个仿真。
weshenton后的文章中类似
约书亚Gafford（查看资料）
你好，空间限制变量“LIM”只是指在无量纲单位的2D领域的“尺寸”。所以LIM = 200生成200×200的运营商的变量（“N”）的舞台是指颗粒（人）的总数，可以修改最初通过玩弄“p_init”感染的分数。
Tsotne Marghia（查看资料）
你好，非常感谢你对你的工作。
我对如何正确地确定空间限制的数字感兴趣。你说的携带者是指整个人群还是只是一小部分被感染的人?
谢谢你！
约书亚Gafford（查看资料）
谢谢你的报道。我假设在元素操作中存在一个向量大小的问题，因为v是[n x 2]，而隔离是[n x 1](我的MATLAB版本中，如果列大小不同，我猜会自动假设行元素乘法，但可能旧版本没有?)不管怎样，现在应该对你有用了。
longonauta（查看资料）
对不起，我在评论时不小心按了回车键。我有这个错误
错误使用。*
矩阵尺寸必须一致。
误差在simulitis（线131）
pos_new = POS + V *（〜分离物）* DELT。;
约书亚Gafford（查看资料）
嗨，阿比盖尔，你有什么问题吗?任何错误?
阿比盖尔吉梅内斯（查看资料）
它有一些bug…我不能运行它。
E. Cheynet（查看资料）
斯蒂芬Cobeldick（查看资料）
整洁的代码，良好的帮助，易于使用…最及时的:谢谢:)