问题标签 [jitcode-jitcdde-jitcsde]

我是一名使用 Python 编码的学生，并试图解决以下延迟微分方程：

我想使用JiTCDDE，但没有成功找到适应系统的方法，即使在研究了模块文档中的示例之后也是如此。我遇到的主要问题是我不明白如何处理同时包含y和v的第二个方程。

目标是绘制系统的分岔图（v作为τ的函数）。我是否使用了错误的工具？或者有没有办法在我的情况下使用 JiTCDDE？

我正在使用 JiTCDDE 在修改后的 Oregantor 模型上解决 DDE。我遇到的问题在分岔点附近，它将返回负值。虽然我知道这些是数学上有效的解决方案，但 Oregantor 代表了一个化学系统。因此，否定答案对于真实系统来说是不现实的。有没有办法设置代码以在 <= 0 时返回变量的最小值。以下是我到目前为止的代码的主要部分。

我想使用 python 3.6 求解一个基于网络的微分方程系统。方程组如下：

我编写了下面的代码，试图解决空间网络中的微分方程组。

我期待得到四个微分方程的解和一些模拟来表示方程。我收到的错误消息是“RuntimeError: Not Implemented”

我想为在 Python 3.6 中连接到网络中的个人绘制微分方程组的图表。方程组如下：

我已经为连接在网络中的n 个个体编写了对病毒入侵的免疫反应的代码。

该模型代表了一个耦合方程系统，表示连接个体网络中的抗体和病毒种群。

我期望得到随着时间 t 的抗体和病毒种群的图表。但是，我收到以下错误消息。

ValueError：x 和 y 必须具有相同的第一维，但具有 (50,) 和 (1,) 形状

我正在通过神经网络获取一个函数（一个学习的动态系统），并希望将其传递给JiTCODE以计算轨迹、Lyapunov 指数等。根据 JiTCODE 文档，该函数f必须是一个符号函数。有什么办法可以改变这一点，因为最终 JiTCODE 将要对符号函数进行羔羊化处理？

基本上，这就是我现在正在做的事情：

我想使用 JiTCDDE 来解决 DDE 系统，该系统在 RHS 上具有过去导数ẏ</em>( t −<em>τ )。

我应该使用该get_state功能还是有更合适的方法？

我正在尝试使用 Python 研究以下延迟微分方程的行为：

其中f是一个截止函数，当其参数的绝对值在 1 和 10 之间时，它基本上等于恒等式，否则等于 0（见图 1），并且Nd和τ是T常数。

为此，我使用包 JiTCDDE。这为上述方程提供了一个合理的解。然而，当我尝试在等式右侧添加噪声时，我得到的解在几次振荡后稳定到非零常数。这不是方程的数学解（唯一可能的常数解等于零）。我不明白为什么会出现这个问题以及是否有可能解决它。

我在下面重现我的代码。在这里，为了简单起见，我用高频余弦代替了噪声，它被引入方程系统作为虚拟变量的初始条件（余弦可以直接引入系统，但对于一般噪音，这似乎是不可能的）。为了进一步简化问题，我还删除了涉及f函数的术语，因为没有它也会出现问题。图 2 显示了代码给出的函数图。

顺便说一句，我注意到即使没有噪音，解决方案在零点似乎也是不连续的（对于负时间，y 设置为零），我不明白为什么。

我正在使用 JiTCODE 来计算 Lorenz 振荡器的 Lyapunov 指数。

这是文档之后的简单脚本：

我收到以下错误。

我认为方程式不是问题，因为我已经解决了jitcode但jitcode_lyap无法解决。

我正在尝试用一些参数求解一个二维延迟微分方程。问题是我无法获得正确的解决方案（我知道）并且我怀疑它来自集成步骤，但我不确定并且我不太了解 JiTCDDE 的工作原理。

这是 DDE：

这是我的模型：

我只对 y(1) 解决方案感兴趣

和参数：

这是我对该模型和参数的绘图：

这是（蓝线）正确的解决方案（有人给我情节而不是数据）

我正在积分以下微分方程

使用 JiTCDDE 使用此代码：

并使用此参数

我的麻烦是，当我绘制通过

我最多得到一个非常不平滑的配置文件，如下所示：

当我将集成步骤更改为 0.1 时，我得到了这个

但那个情节去 300 而不是 30 这是我想要的

问题是：有什么方法可以将方程积分到 t=30 ，但像 t=300 一样平滑？我可以通过仅更改参数来做到这一点吗？