plot - Julia 语言微分方程示例

Question

我试着玩弄 this Julia 文档中的示例。我的尝试是将细胞分成两部分，每部分含有一半的蛋白质，所以我设置 Theta=0.5。然而，情节看起来像这样:

很明显，每次达到目标蛋白质量时，细胞数量都会同时翻倍，因为它们是相等的。我怎么能策划这个？我也不明白为什么在下面的例子中单元格的数量会停在 3 个。

Answer 1

绘制每个细胞中的蛋白质含量并考虑您创建的模型。第一次除法后，两个单元格具有相同的值。所以在同一时间，你有一个事件火灾。 “最大值”(以较低的索引为准，因此 1)将 split ，而 2 将继续增长到 1 以上。但是现在 u[2] > 1，求根条件 1-maximum (u) 将永远不会再次归零，因此不会再发生 split 。这意味着您将总共有两个拆分，即 3 个单元格。

请记住，程序会完全按照您的吩咐去做。我假设您的意思是，作为您的效果，拆分大于或等于 1 的 任何 单元格。如果那是您想要的 affect!，那么您d 必须写它:

function affect!(integrator)
  u = integrator.u
  idxs = findall(x->x>=1-eps(eltype(u)),u)
  resize!(integrator,length(u)+length(idxs))
  u[idxs] ./ 2
  u[end-idxs:end] = 0.5
  nothing
end

这只是一种方法，当然还有很多其他方法。