scipy.integrate.odeint(func,y0,t,args =(),Dfun = None,col_deriv = 0,full_output = 0,ml = None,mu = None,rtol = None,atol = None,tcrit = None,h0 = 0.0,hmax = 0.0,hmin = 0.0,ixpr = 0,mxstep = 0,mxhnil = 0,mxordn = 12,mxords = 5,printmessg = 0,tfirst = False)
解决一阶ode-s的刚性或非刚性系统的初值问题:
dy/dt = func(y, t, …) [or func(t, y, …)]
主要参量
func callable(y,t,…)或callable(t,y,…)
计算t处y的导数。如果签名是,则必须设置参数 tfirst。callable(t, y, ...)True
y0 数组
y的初始条件(可以是向量)。
t 数组
求解y的时间点序列。初始值点应该是此序列的第一个元素。此序列必须单调增加或单调减少;允许重复的值。
args 元组,可选
传递给函数的额外参数。
例子:
摆锤的重力角在摩擦力作用下的角度θ的二阶微分方程可以写成: