differential-equations

0热度

1回答

我的问题是：为什么rhsexpr有效，但rhsfunc没有。谢谢。 :) m = l = 1.; \[Mu] = 100; beltv[t_] = .1; spring[x_] = 1000. (l - x); viscous[v_] := -30. (v - beltv[t]); coulomb[v_] := -25. Sign[v - beltv[t]

0热度

2回答

拟合微分方程：如何适应一组数据到差分方程中的R

与数据集： conc <- data.frame(time = c(0.16, 0.5, 1.0, 1.5, 2.0, 2.5, 3), concentration = c(170, 122, 74, 45, 28, 17, 10)) 和我想此数据代入下面的微分等式： dC/dt= -kC 其中C将是数据集中的浓度和时间t。这也会给出k的结果。有谁能给我一个线索如何在R中做到这一点？谢谢。

4热度

1回答

如何在matlab中绘制斜场

我正在寻找一种在Matlab中绘制斜场的方法。这里就是我要找：我有一个公式 dy/dx = f(x,y) 或 dx/dt = f(x,y) dy/dt = g(x,y) ，我希望把它画在一个不错的方式因为关于它的唯一answer没有回答我的问题，它花了我一些时间如何做到这一点。因为这不是我一直在MATLAB中做的事情（很可能直到下一次我需要它，我会忘记它）我正在为我创建一个备忘录如

1热度

1回答

如何绘制并求解枫叶12的微分方程？

考虑初始条件为$y\left(0\right)=1$的微分方程$y^{\prime}=y-2$。一）使用欧拉法有4个步骤尺寸0.2的估计$y\left(0.8\right)$ 我知道如何做手工这一点;然而，我已经安装了枫树12，并试图找出如何用枫树做到这一点，然后制作一个图形显示该功能的每一步。有什么建议么。我试图在枫树上寻找，但它一直指向我的功能新版本的枫树软件，我没有。我发布了这个问题作为

1热度

2回答

实现龙格库塔系统的DE

我试图在MATLAB中实现Runge-Kutta Method for Systems of DEs。我没有得到correct answers，我不确定代码中是否有错误或者我用来运行它的命令。这里是我的代码： function RKSystems(a, b, m, N, alpha, f) h = (b - a)/N; t = a; w = zeros(1, m);

0热度

1回答

使用龙格库塔解决贝塞尔函数

我正在为我的一个类的任务，我应该写一个代码使用我选择的程序（我选择了Matlab）来解决贝塞尔函数微分方程使用4阶Runge-Kutta方法。作为参考，贝塞尔函数DE是： x^2 *（J_n）''+ x *（J_n）'+（x^2-n^2）* J_n = 0。我能够通过此分离成两个耦合的第一顺序的DE：（J_n） '= Z_n和（Z_n）' +（1/X）* Z_n + [（X^2-正^ 2）/

0热度

2回答

如何在Matlab中表示和求解这个微分方程组？

的方程如下： DXI（T）/ dt的= -c * XI（T）*义（吨） DYI（T）/ dt的= A *Σ{i = 1至其中a，b和c是任何你想要的常数值，例如a = 1，其中a，b和c是任意常数值，例如a = 1， b = 2，c = 3。（t）* y1（t））表示从i = 1到n的和，例如n = 3：x1（t）* y1（t）+ x2（x1（t）（t）* y2（t）+ x3（t）* y3（t

0热度

1回答

线性双方程系统，两个变量二阶导数都在

原谅我，如果这被认为是重新发布，但我已被建议我可能已经给出了一个不好的格式。我试图解决两个线性二阶微分方程。我想将它们分解为单阶方程，但我无法看到两个变量在两个问题中都有二阶导数。 (m*a)u” + (I + m*a^2)θ” + (d*a^2)θ’ + (K - m*g*a)θ = 0 (M + m)u” + (m*a)θ” = -F 我已经使用了一些既dsolve和ODE45 这是

4热度

1回答

在Python中绘制轨道轨迹

如何在python中设置三个主体问题？如何定义函数来解决ODE？三个方程是 x'' = -mu/np.sqrt(x ** 2 + y ** 2 + z ** 2) * x， y'' = -mu/np.sqrt(x ** 2 + y ** 2 + z ** 2) * y，并 z'' = -mu/np.sqrt(x ** 2 + y ** 2 + z ** 2) * z。写为6阶，我们有 x' =

1热度

1回答

python执行3d图时会冻结

从plotting orbital trajectories开始，我们有以下代码。价值已经改变为已知的IC工作。如果此代码是正确的（它不能是虽然），它会产生运行这段代码只是冻结我的电脑或输出绝对错误的地块。有人可以帮我找到解决这个问题的方法吗？ import numpy as np from scipy.integrate import odeint import matplotlib.p