用SOLVE-IVP求根并找到ODE的零交叉

时间:2019-04-28 20:57:41

标签: python root ode solver

我一直在寻找一种方法来找到微分方程改变其符号的点,它可以是零交叉或求根。

有人告诉我使用Brent's Method,但这并不准确。例如,图形在0.12处切零,但是该方法获得了0.19之类的值或更多,这是不可接受的。

现在,我发现了该算法scipy.integrate.solve_ivp,在这种方法中,它表示为The solver looks for a sign change over each step, so if multiple zero crossings occur within one step, events may be missed.

我不知道该怎么做。

考虑我们有一个微分方程

dE/dz = f(E,z)

我们知道dE/dz的答案,我们知道起始点E0,但是我们想知道在哪个z中答案为dE/dz。因此我们应该获取z的值,从而得出dE/dz

的答案

请帮忙吗?

谢谢

0 个答案:

没有答案