查找插值函数的积分何时等于特定值（python）

Question

我有x和y点的数组t_array和dMdt_array。我们打电话给M = trapz(dMdt_array, t_array)。我想找出t的什么值，dM / dt与t的积分等于某个值 - 比如0.05*M。在python中，有一种很好的方法吗？

我在想F = interp1d(t_array, dMdt_array)之类的东西。然后某种根找到F的积分等于0.05 * M的位置。我可以在python中这样做吗？

Answer 1

一种简单的方法是使用CubicSpline类。然后是CubicSpline(x, y).antiderivative().solve(0.05*M)或其左右。