我正在尝试使用SciPy模拟一些动态模型。
我有模型定义:
def model(y, t, control_signal):
dy/dt = some_function_of_time_and_y
return dy
我定义了要模拟模型的时间戳列表:t_list=np.linspace(0, 5, 100)
。我想使用为每个时间戳定义的control_signal
值来模拟模型。我试图通过使用:
controls = [list_of_values]
scipy.integrate.odeint(model, 0, t_list, args=(controls))
但是我得到The size of the array returned by func (5) does not match the size of y0 (1).
似乎我的controls
被解释为模型状态,而不是每个时间戳中的输入。如何传递controls
作为每个时间戳的值?
谢谢!
在scipy.interpolate.interp1d
中,您可以将插值模式定义为“零保持”或0阶样条,即分段常数,类型为“ kind =“ zero”。用它来定义控件的时间依赖性。
contfunc = interp1d(t_list,control, kind="zero");
def model(y, t, control_signal):
u = contfunc(t);
dydt = some_function_of_time_and_y_and_u
return dydt
尺寸错误可能是另一个问题。要进行调试,请使用调试器(如果可用)或为输入和输出的大小/形状添加打印语句。状态和导数应为相同大小的平面数组。
请不要忘记将最大时间步长设置为小于控件的步长,这无关紧要,但是可能会导致奇怪的错误。