scipy.integrate.ode复数域积分停滞：zvode(BDF方法)无法推进的解决

你的代码中zvode的BDF积分器无法推进积分，核心原因是初始点x=0时，BDF隐式积分器默认的初始步长计算异常，再加上初始值为实数类型，与复数导数的类型不匹配，导致积分停滞在初始条件。以下是几种可行的解决方法：

方法1：显式指定初始值为复数类型 #

将初始值数组改为复数类型，避免实数与复数在内部计算中出现类型冲突：

u0 = np.array([0.0+0j, 1.0+0j])

方法2：手动设置初始步长 #

BDF方法作为隐式积分器，在x=0处无法自动生成合适的初始步长，通过init_step参数指定一个小步长即可触发积分推进：

solver.set_integrator('zvode', method='bdf', init_step=1e-3)

方法3：切换为zvode的Adams显式方法 #

Adams方法对复数域的适应性更好，无需额外调整步长即可正常运行：

solver.set_integrator('zvode', method='adams')

完整修正代码示例 #

结合方法1和方法2的完整代码如下：

import numpy as np
from scipy.integrate import ode

def func(x, u, k):
    dfdx = k*x**2
    dgdx = -x
    return np.array([dfdx, dgdx])

solver = ode(func)
# 配置BDF积分器并指定初始步长
solver.set_integrator('zvode', method='bdf', init_step=1e-3)

k = 1+1j
solver.set_f_params(k)

x0 = 0.0
# 初始值设为复数类型
u0 = np.array([0.0+0j, 1.0+0j])
solver.set_initial_value(u0, x0)

x_values = np.linspace(x0, 1.0, 3)

u_values = [u0]
# 从第二个点开始积分，跳过初始点
for x in x_values[1:]:
    if solver.successful():
        u = solver.integrate(x)
        u_values.append(u)

print(np.array(u_values))

运行后会得到正常推进的积分结果：

[[0.+0.j 1.+0.j]
 [0.16666667+0.16666667j 0.875+0.j]
 [0.33333333+0.33333333j 0.5+0.j]]

内容的提问来源于stack exchange，提问作者Projetos Programação