如何在python中实现随机ODE(SDE)系统？

import numpy as np
from scipy.integrate import odeint

def random_ode_system(y, t):
    # 定义随机ODE系统的参数
    a = 0.5
    b = 0.1
    sigma = 0.2

    # 定义随机项
    dW = np.random.normal(0, sigma, size=len(y))

    # 定义ODE系统的方程
    dydt = a * y + b + dW

    return dydt

# 定义初始条件和时间点
y0 = [0.0]
t = np.linspace(0, 10, 100)

# 求解随机ODE系统
sol = odeint(random_ode_system, y0, t)

# 打印结果
print(sol)