心房颤动计算机模型的Matlab实现最初由Gordon Moe和他的同事发表在[美国心脏杂志1964;67(2): 200 - 220。要运行模拟,只需键入函数名(afib)。它不接受输入或输出参数。通过编辑函数行,可以改变多个模型参数。在函数体前面的注释中对此进行了解释。 由于不应期分散的机制,该功能基本上是一个细胞自动机模型模拟心房扑动/心房颤动。 一簇(4)细胞被迅速刺激,引起周期性波动 激活扩散到模型表面。由于在整个表面的难熔期的变化,波阵面在与仍然难熔的区域碰撞时被打破。快速刺激在80步后结束,但激活波阵面仍在模型表面传播。作者通过寻找不应期较短的相邻细胞来仔细选择一簇细胞进行刺激。我在这里不这么做;相反,我只是不断改变用来标记K值的随机数发生器的初始状态,直到我得到自我维持的活动。 在包含仿真的图形窗口中,绘制了三种离散状态:激活状态或计划状态(蓝色), 难熔(浅蓝色)和静息/易兴奋(白色)。时间步模拟的数量显示在激活模式下面。要调节模拟的速度,请将输入参数更改为pause语句(时间,以秒为单位,在每一步中暂停)。 随着时间的推移,生成第二幅图,显示“伪电图”。这再现了原文中的一个情节。 # v/ w. W: o% |7 H# |1 d$ r
$ u# Q6 a# M2 {
|