问题标签 [runge-kutta]

我想知道是否有人可以帮助我解决我在 excel VBA 中的打印功能遇到的问题。我正在尝试使用 RK4 打印 ODE 的输出，并且我想知道如何循环 printvalue 函数，以便每个循环都会发生单元格值的增加。

还想知道为什么 PrintValues (x.ToString()) 行说 x 是无效的限定符？如果您觉得代码中有任何其他错误，请告诉我，尽管我相当确定一切正常 任何帮助将不胜感激:)

我想使用deSolve R 包中的显式 Runge-Kutta 方法ode45（别名 rk45dp7）来解决可变步长的 ODE 问题。

根据 deSolve 文档，可以使用 ode45 方法对rk求解器函数使用自适应或可变时间步长，而不是等距时间步长，但我不知道如何做到这一点。

rk 函数是这样调用的：

其中times是需要明确估计 y 的时间向量。

对于距离为 0.01 的等距时间步，我可以将时间写为

假设我想求解从 0 到 100 的区间的方程，如何在不给出步长的情况下定义时间？

任何帮助将不胜感激。

I've been working on a 4th Order Runge-Kutta solver, and am having some difficulties. I've written the solver based on the article on gafferongames, but when I run a small contained example, the errors I'm getting are far far worse than what I was getting with simple Euler Integration, even for a simple gravity force. I've tidied it up into a self contained example (~60 lines of code including printing), but it requires GLM to run.
It shows my issue in its entirety. Line 55 prints out the difference between the analytic solution and the RK4 solution. This should be relatively small, but it explodes even after the 10 odd steps it takes.

If anyone could help me out here, I'm at the end of my tether with this.

当我需要在 MATLAB 上通过 ode45 求解 ODE 系统时，我可以使用这种形式吗

我的意思是我的表达正确吗？

例如我该如何解决这个问题：

初始条件是：

如果有人可以帮助解决以下问题，我将不胜感激。我有以下 ODE：

我以两种不同的方式解决了（1）。通过Runge-Kutta 方法（四阶）和ode45在 Matlab 中。我将这两个结果与解析解进行了比较，解析解由下式给出：

当我绘制每种方法相对于精确解的绝对误差时，我得到以下信息：

对于 RK 方法，我的代码是：

对于ode45：

我的问题是，为什么我在使用时会出现振荡ode45？（我指的是绝对错误）。两种解决方案都是准确的 ( 1e-9)，但ode45在这种情况下会发生什么？

当我计算 RK 方法的绝对误差时，为什么它看起来更好？

我正在尝试用 Runge-Kutta 逼近 Bessel 函数。我在这里使用 RK4，因为我不知道 RK2 对于 ODE 系统的方程是什么.....基本上它不起作用，或者至少它不是很好，我不知道为什么。

无论如何，它是用python编写的。这是代码：

my problem is following: I need to integrate Jiles-Atherton model(it's described in answer)

I used the same parameters as in answer but I got bad results ;( First in my program I increase H from 0 to 2500, then decrease from 2500 to -2500 and then increase again from -2500 to 2500.

My code:

What's the problem?

我正在尝试在 Haskell 中实现 4 阶 Runge-Kutta，但我发现很难将 Haskell 类型系统用于此任务。有人可以帮忙吗？我希望在以下代码中将“State”和“DState”类型更改为类型类：

似乎很难制定类型类，因为 State 和 DState 在某种意义上是交织在一起的，即 State 的特定实例需要 DState 的特定实例。或者是否有其他我看不到的方法？

我有以下代码：

我发现 s1 和 s2 给出不同的结果，这本质上是方程（mRNA和蛋白质）的不同顺序。

我是否可以进行任何修改以使订购不会产生影响？

提前致谢。

我正在尝试通过 Runge-Kutta 四阶数值求解双摆。

要求解的微分方程在下一页上：

http://www.myphysicslab.com/dbl_pendulum.html

甚至还有一个动画显示角度 1 和 2。我想通过 Runge-Kutta 方法求解出现在页面下方的方程。

将这些结果打印到文件中，然后使用 gnuplot 绘制它们。但是，我的程序无法运行，找不到原因，希望您能帮帮我。

我在 Fortran 中的代码：

这是该方法的一个很好的描述：http: //mathworld.wolfram.com/Runge-KuttaMethod.html

运行程序时出现的错误是：