什么 ODE 求解器使用步进函数中的计算进行插值?
我对具有不同初始条件的 ODE 的多个解进行平均,因此所有解同时具有值非常重要;例如,增量为 0.01。 我一直在使用数值食谱 3 (nr3) 中的 ODE 例…
有没有c++常微分方程 (ODE) 求解器库?
更具体地说,我对八阶 Dormand-Prince 嵌入方法感兴趣,它基于 Runge-Kutta 和刚性方程。 我使用 Numerical Recipes 3,但我经常在编译他们的库时遇到…
假设高斯分布,如何为几个常微分方程添加白噪声过程项?
这个问题已经让我困惑好几天了。虽然我向高年级学生提出了要求,但他们也无法给出答复。 我们有十个 ODE,每个 ODE 都应添加一个噪声项。噪声定义如下…
矢量<双>比双倍快*:为什么?双>
这是我使用 std::vector 和普通旧 double* 尝试过的循环。 对于 1000 万个元素,矢量版本的运行时间始终是 double* 版本的 80% 左右;对于几乎任何 N …
Maple ODE 求解,图形不完整。为什么?
所以我们尝试求解简单的 ODE。我们使用文本书籍(警告 - 俄语)。这是书中示例的代码(第 18 页): > with(DEtools); with(plots); > eq1 := diff(x(…
(Matlab)如何根据参数绘制两个ODE的相图(不是图中的x(1),x(2))
亲爱的大家,我是 matlab 的新用户,正在做一个关于动态系统的项目。我有两个方程,其中 T0、epsilon (E)、a 是固定的。我想查看当 b 和 c 分别在 (0,…
scipy.integrate.ode 具有两个耦合的 ODE?
我目前正在尝试使用 SciPy 的integrate.ode 包来求解一对耦合的一阶 ODE:例如 Lotka-Volterra 捕食者-被捕食者方程。但是,这意味着在集成循环期间,…
寻找三阶 ODE 和直线之间的交点?
如何找到三阶 ODE 的解与直线 y=x 之间的交点? 我的 ODE 代码是 sol=dsolve('D3y-4*D2y+Dy+2*y=0,y(0)=-4,Dy(0)=-6,D2y(0)=-4') x=0:2 y=subs(sol,'t…
Java 中微分方程组的 Runge-Kutta (RK4)
这个方程主要是这个线程的结果:Java 中的微分方程。 基本上,我尝试遵循 Jason S. 的建议,并通过 Runge-Kutta 方法 (RK4) 实现微分方程的数值解。 …
Java 中的微分方程
我正在尝试用java创建一个简单的SIR流行病模型模拟程序。 基本上,SIR 由三个微分方程组定义: S'(t) = - l(t) * S(t) I'(t) = l(t) * S(t) - g(t) * …
在 Mathematica 中编写欧拉方法
我想编写一个函数,其中有一个循环,它执行欧拉方法所需的操作。下面是我糟糕的尝试。 In[15]:= Euler[icx_,icy_,h_,b_,diffeq_] := curx; cury; n=0;…
Mathematica 中非线性 ODE 的初始条件
我正在尝试使用 Mathematica 的 NDSolve[] 使用耦合 ODE 计算沿球体的测地线: x" - (x" . x) x = 0 问题是我只能输入 x(0) 和 x'(0) 的初始条件,并…
- 共 1 页
- 1