如何向solve_ivp传递带关键字参数的自定义微分方程函数？

核心问题 #

scipy.integrate.solve_ivp要求微分方程函数的签名为func(t, y, *args)，它只会自动传递t（时间）和y（状态向量）两个参数，自定义的关键字参数需要通过特定方式传递，不能直接混入solve_ivp的调用参数中。

正确传参方式 #

假设你的微分方程函数和辅助函数定义如下：

import numpy as np
from scipy.integrate import solve_ivp

def sineCurrent(t, A=1, freq=1):
    return A * np.sin(2 * np.pi * freq * t)

def integrate(t, y, A=1, freq=1):
    i = sineCurrent(t, A, freq)
    # 示例微分方程：dy/dt = i - y
    dy_dt = i - y
    return dy_dt

方式1：用solve_ivp的args参数传递

把自定义关键字参数以位置参数的形式打包成元组，传给args：

t_span = (0, 10)
y0 = [0]
# 对应integrate的A=2, freq=0.5
sol = solve_ivp(integrate, t_span, y0, args=(2, 0.5))

如果需要明确参数对应关系，确保args里的顺序和integrate函数中除t、y外的参数顺序一致即可。

方式2：用lambda函数包裹（不修改原函数）

直接在lambda中接收t和y，调用integrate时传入关键字参数：

sol = solve_ivp(lambda t, y: integrate(t, y, A=2, freq=0.5), t_span, y0)

这种方式完全保留原函数的独立性，符合课程要求。

方式3：用functools.partial绑定参数

通过partial固定关键字参数，生成一个新函数传递给solve_ivp：

from functools import partial

# 绑定关键字参数
integrate_with_params = partial(integrate, A=2, freq=0.5)
sol = solve_ivp(integrate_with_params, t_span, y0)

之前错误的原因 #

1. 基础语法调用报错：你可能直接在solve_ivp中传入了自定义关键字参数（如solve_ivp(integrate, t_span, y0, A=2)），但solve_ivp的参数列表中没有这些自定义参数，导致语法错误。
2. lambda语法错误：大概率是lambda的签名写错了，比如多余定义了参数（lambda t, y, A=2: ...），但solve_ivp只会传递t和y，引发参数不匹配。
3. 关键字移至末尾报错：同样是因为solve_ivp不识别自定义关键字参数，即使放在末尾也会被判定为无效语法。

内容的提问来源于stack exchange，提问作者Anna Jacobsen