我正在尝试复制Tang & Othmer的论文，该论文与盘基网柄菌（一种变形虫物种）中的 G 蛋白模型中的激发和振荡有关。论文中的数学模型如下：

其中是生化物质，其他参数是与盘基网柄菌相关的常数，是分泌率函数。我应该在方形网格中模拟一种情况，其中细胞均匀分布并相互作用。是细胞外cAMP（环状AMP），它从变形虫细胞中扩散出来，在某些条件下在空间中形成螺旋波形。所以我基本上在MATLAB中使用欧拉法模拟了这个。但我无法在模拟中生成螺旋。$w_1,w_2,w_3,w_4,w_5,u_1,u_2,u_4$$sr(.)$$w_5$

有人可以给我一些想法在网格$w_5$

以下是考虑单个细胞且无扩散内容物的结果（5 分钟）。（只是为了表明我已经从论文中复制了一些结果）

一个 2D 解决方案如下所示（来自论文）：

我对这个问题的尝试没有成功。我为启动输出所做的首先是为每个单元格定义一个为等变量设置的初始值，除了这些值的补丁（图中的图）。我从 ODE 图中选择了这些初始值。但这种安排被证明是不成功的。$w_1, w_2 ...$$180°$imagesc$w_5$

更新：现在我改进了我的代码，使输出限制在 2 个值之间。它不再爆炸了。如有必要，我可以提供我的代码。这是我创建的视频输出的链接。