MATLAB 中的延迟微分方程 (dde23 )
我一直在尝试创建一个每周注射药物的生物模型，该模型会导致少数免疫细胞类型仅在 7-14 天发生反应。该模型由一个带有周数索引的 for 循环和一个用于…
使用 odeint 求解谐波振荡器与 matlab 解不匹配
我正在使用 odeint 和 matlab 求解谐波振荡器。我试图求解的方程是x'' = -x + 0.15*x'。 odeint代码可以在这里找到 确切的解决方案是（每个 matlab）…
在 MATLAB 上使用 Crank-Nicolson 方法计算热方程
我正在尝试在matlab中实现这个方程的曲柄尼科尔森方法： du/dt-d²u/dx²=f(x,t) u(0,t)=u(L,t)=0 u(x,0)=u0(x) 与： - f(x,t)=20*exp(-50(x-1/2)²) if …
求解一阶方程的7个方程 ode45 Matlab
我正在尝试求解具有七个微分方程的系统。我很难掌握 ode45 求解器。 这些是等式： ω2_dot = -0.75 ω1 ω3 ω1_dot = 0.75 ω2 ω3 + 0.2 ω3_dot = …
如何指定求解微分方程的最终值（而不是初始值）
我想求解 R 中的微分方程（使用 deSolve？），我没有初始条件，只有状态变量的最终条件。这怎么能做到呢？ 典型的代码是：ode(times, y,parameters, f…
Python：求解自治二阶 ODE
我有以下自治常微分方程系统，其中x和v是向量（初始条件x0和 >v0) dx/dt = v dv/dt = F(x, v) 等价地，可以将其编写为二阶 ODE 系统 d^2 x/ dt^2 = F(…
转换 dsolve 输出以求解 sympy 中的值
我 import sympy as sm x = sm.symbols('x', cls=sm.Function) t = sm.symbols('t') expr = x(t).diff(t) + 0.05*x(t) sol = sm.dsolve(expr,x(t), ic…
在不使用根的情况下满足特定条件时提前停止 R 的 deSolve 中的积分
我知道 R 的 deSolve 中的提前终止可以通过使用根函数而不提供事件函数来实现，这将导致在找到根时终止积分。然而，通过使用此过程，我们仅限于应用具…
如何使用 mathematica 绘制斜率场？
我正在尝试使用 mathematica 绘制一些微分方程的斜率场，但无法弄清楚。假设我有方程 y' = y(t) y(t) = C * E^t 如何绘制斜率场？ 我找到了一个例子，…
python求解具有复变量的微分方程
伙计们， 是否可以在 python 中使用复杂变量求解 ODE？我的方程具有以下形式， dx/dt = -a x -i y(t) 其中 y(t) 是已知函数，a 是已知数，i 是 -1 的…
如何在函数定义中使用 Maxima 的 diff？
我想在另一个函数中使用一个函数的导数。在 Maxima 中应该如何完成此操作？ 例如： f(x) := 2*x^4; g(x) := diff(f(x),x)-8; 现在 g(x) 按预期产生 8x…
matlab中的微分
我需要找到物体的加速度，文本中给出的公式是 a = d^2(L)/d(T)^2 ，其中 L= 长度，T= 时间 我通过使用这个方程在matlab中计算了这个 a = (1/(T3-T1))*…
求解 ODE，直到 Octave 中发生事件
在 Matlab 中，我知道如何求解 ODE 直到偶数发生（例如，求解 ODE 直到 y 坐标为零）： http://www.mathworks.com/help/techdoc/ref /odeset.html#f92…
