scipy.solve_ivp 支持多事件检测，通过 sol.t_events 和 sol.y_events 的列表索引顺序与传入 events= 参数的函数顺序严格对应，可直接定位触发的是第几个事件。
scipy.solve_ivp 支持多事件检测，通过 `sol.t_events` 和 `sol.y_events` 的列表索引顺序与传入 `events=` 参数的函数顺序严格对应，可直接定位触发的是第几个事件。
在使用 solve_ivp 进行常微分方程数值积分时，事件（event）机制是捕获系统状态突变（如角度翻转、碰撞、过零等）的关键工具。尤其在双摆这类强非线性系统中，我们常需分别监控第一杆或第二杆是否发生“翻转”（例如广义角 θ₁ 或 θ₂ 跨越 ±π），此时定义多个事件函数并明确区分其触发来源至关重要。
solve_ivp 会将每个事件函数的触发时间与状态自动归类存储：
sol.t_events 是一个
列表
，长度等于传入的 events 列表长度；
sol.t_events[i] 是一个 NumPy 数组，包含第 i 个事件函数所有触发时刻（按积分顺序排列）；
同理，sol.y_events[i] 是对应时刻的状态向量数组（形状为 (n_events, n_states)）。
✅ 正确用法示例（双摆翻转检测）：
⚠️ 注意事项：
索引严格对齐
：events[0] ↔ sol.t_events[0]，不可依赖函数名或返回值内容推断；
空数组表示未触发
：若某事件从未满足条件，对应 t_events[i] 为空数组（array([], dtype=float64)），务必用 len(...)>0 或 .size > 0 判断；
方向控制精度
：direction=1（仅上升过零）、-1（仅下降）、0（任意）——对翻转检测建议用 0，避免漏判；
终端事件优先级
：多个 terminal=True 事件同时满足时，solve_ivp 以
列表靠前的事件为准
并立即终止积分，因此事件顺序也隐含优先级。
总结：无需额外标记或回调，solve_ivp 原生通过结构化输出实现事件溯源——牢记“位置即身份”，合理组织 events 列表，即可在双摆、航天器姿态、电路开关等多事件场景中实现清晰、可靠的状态监控。
import numpy as np
from scipy.integrate import solve_ivp
# 假设双摆动力学函数（简化示意）
def double_pendulum(t, y):
theta1, omega1, theta2, omega2 = y
# 此处省略具体导数计算（如基于拉格朗日方程）
dtheta1_dt = omega1
domega1_dt = ... # 实际表达式
dtheta2_dt = omega2
domega2_dt = ... # 实际表达式
return [dtheta1_dt, domega1_dt, dtheta2_dt, domega2_dt]
# 定义两个事件：第一杆翻转（θ₁ ≡ π mod 2π）、第二杆翻转（θ₂ ≡ π mod 2π）
def flip_rod1(t, y):
return np.sin(y[0] / 2) # 零点等价于 θ₁ = 0, ±2π, ±4π…，但翻转通常关注 θ₁ = ±π → sin(θ₁/2)=±1？更稳健做法是：(y[0] + np.pi) % (2*np.pi) - np.pi → 接近 ±π 时趋近 0
# ✅ 更推荐：检测 θ₁ 是否穿越 ±π，使用无符号距离：np.cos(y[0]) + 1 会在 θ₁=±π 时为 0（但有双根）；或直接：y[0] % (2*np.pi) - np.pi → 零点即翻转点
def flip_rod2(t, y):
return np.cos(y[0]) + 1 # 示例：实际应替换为基于 y[2]（θ₂）的合理事件函数，如 np.sin(y[2]/2)
# ✅ 关键：明确指定事件属性（推荐设 terminal=True 防止多次触发）
flip_rod1.terminal = True
flip_rod1.direction = 0 # 检测所有过零（上升+下降）
flip_rod2.terminal = True
flip_rod2.direction = 0
# 传入事件列表，顺序即索引标识
events = [flip_rod1, flip_rod2]
sol = solve_ivp(
double_pendulum,
t_span=[0, 10],
y0=[0.1, 0, 0.2, 0], # 初始小扰动
events=events,
max_step=0.01,
rtol=1e-6
)
# ? 判断哪个事件被触发
if len(sol.t_events[0]) > 0:
print(f"第一杆在 t = {sol.t_events[0][0]:.4f} 时翻转")
if len(sol.t_events[1]) > 0:
print(f"第二杆在 t = {sol.t_events[1][0]:.4f} 时翻转")
# 若需获取完整触发记录（如多次翻转）：
for i, t_list in enumerate(sol.t_events):
if len(t_list) > 0:
print(f"事件 #{i+1} 共触发 {len(t_list)} 次：{t_list.round(4)}")