我不確定事件處理是否scipy.integrate.solve_ivp工作正常。在下面的例子中，我對一個導數進行積分，得到一個三次多項式，它的根在x=-6，x=-2和x=2。我設置了一個事件函數，返回y，在x值處為零。我希望在解決方案的t_events屬性中看到三個條目，但是我只看到一個條目，盡管很明顯解決方案跨越x軸三次。在

我是不是做錯了什么？在

復制代碼示例：def fprime(x, y):

return 3 * x**2 + 12 * x - 4

def event(x, y):

return y

import numpy as np

from scipy.integrate import solve_ivp

sol = solve_ivp(fprime, (-8, 4), np.array([-120]), t_eval=np.linspace(-8, 4, 10), events=[event])

上述代碼導致：

^{pr2}$

```

問題是從sol.y數組中可以看到應該有三個0(有三個符號改變)，但是只記錄了一個事件。在

Scipy/Numpy/Python版本信息：1.0.0 1.13.3 sys.version_info(major=3, minor=6, micro=0, releaselevel='final', serial=0)

[UPDATE]：如果我使用可選的max_step參數來求解_ivp，并使其足夠小，那么我會看到所有三個根。似乎事件函數并不是在t帴eval步驟中調用的，而是僅在內部解算器步驟上調用，后者遠遠少于t帴eval步驟，并且最終會跳過某些根。這似乎不是很有用，因為你必須知道如何設置最大步長以避免丟失根。在

總結

以上是生活随笔為你收集整理的python scipy库函数solve用法_如何在中使用事件scipy.integrate.solve_ivp的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。