问题标签 [jitcode-jitcdde-jitcsde]

我一直在尝试使用 JiTCDDE 求解以下延迟微分方程组：

模型方程

这是我的代码：

问题是，对于解决方案的第二个数组（我正在尝试访问），当我指定 t<0 时，前几个值是负值，该值为 0。我已经尝试过这个相同的模型使用ddeint它给出了一个单调递增的正值曲线，这是我所期望的。不过我想jitcdde工作，因为即使没有延迟项，这个模型也应该运行。第一个数组看起来不错，我检查了我的模型，看看是否有任何错别字，但对我来说一切都很好。我也尝试过使用adjust_diffand integrate_blindly，但问题仍然存在。

我正在尝试使用 JiTCODE 解决动态食物网。该网络的一个方面是将经历阈值的人口设置为零。所以我得到了一个不可微分的方程。有没有办法在 JiTCODE 中实现它？另一个类似的问题是网络的 Heaviside 依赖性。

示例代码：

我正在使用 JiTCDDE 模块来求解延迟微分方程。我的工作要求我让模型发展一段时间，然后对其进行扰动。为此，我尝试使用jitcdd.purge_past()后跟jitcdde.add_past_points(). 问题是，后者要求我不仅提供输入时的值，还需要提供导数的值。

有没有办法从jitcdde实例本身获取这些，还是我需要手动笨拙地计算它们？

编辑：更多信息

我的系统由两个耦合在一起的非线性振荡器组成（它们不是相位振荡器）。我让系统发展一段时间，直到它达到稳定状态，然后通过稍微移动两个振荡器中的一个来扰乱它。I shift it 的量计算为振荡周期的百分比，即有效相移。

我的系统是 12 维的，我正在努力获得一个最小的工作示例，所以这是代码的非工作模型：

这样做之后，我想执行应该是的扰动，比如说一个周期的 10%，假设我知道静止周期长度是T。我要做的是用来get_state获取系统和导数的当前状态，以及状态和导数perturbation时间单位之前的状态。然后，我可以构造一组新的anchors，我可以将其传递给它DDE.add_past_points以从那里开始集成。