(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) 的初始条件,并…