回复：小推广的Runge-Kutta方法正确吗？

我以前用 Runge-Kutta 方法算过 chaos 的问题，记得不太清楚了，只能说个大概的印象。

你的天体轨道应该是二维平面问题， 有两个 coupled differential equations， 时间是变量参数，the increment （step size h） 应该总是 > 0。

或者，如果你能用 polar coordinate 表达轨道方程， 以角度为参数，那么角度的 increment 也是大于零的。

没有时间细看你的程序，不一定说得对，仅供参考。

 

所有跟帖： 

• 谢谢！我的问题非常简单，没有你想象的那么复杂， -9$- ♂ (68 bytes) () 12/12/2012 postreply 10:54:50

• 呵呵，比我想象的复杂。 -孩子长了翅膀- ♀ (21 bytes) () 12/12/2012 postreply 11:00:49

• 在平面上处处（！）给定 dy/dx， 你就能积分得出平行曲线簇， -9$- ♂ (31 bytes) () 12/12/2012 postreply 11:10:04

• 在2D问题中，如果你只有包含x变量，y，y dot，y double ldot 等叫做常微分方程， -孩子长了翅膀- ♀ (68 bytes) () 12/12/2012 postreply 11:19:45

• 我研究的是最简单PDE： partial y/partial x = f(x,y) -9$- ♂ (0 bytes) () 12/12/2012 postreply 11:28:19

• 这叫常微分方程。 -孩子长了翅膀- ♀ (214 bytes) () 12/12/2012 postreply 11:34:05