使用scipy的odeint求解运动方程时报错无法对float对象解包是什么原因

你遇到的float对象无法解包报错，核心原因是**scipy.integrate.odeint对传入的微分方程函数参数顺序有默认要求**：

• 默认状态下，odeint要求被调用的微分方程函数参数顺序为f(状态量数组, 时间t)
• 你当前定义的微分方程函数为def f(t,xs)，顺序是先时间后状态量，和odeint的默认要求完全相反
• 运行时odeint会把状态量数组传入你定义的t形参，把单浮点型的时间值传入xs形参，后续执行r,theta,phi,V,gamma,psi = xs时，xs是单个时间浮点数，自然无法解包为6个变量，直接触发报错。

有两种任选的修复方式，二选一即可：

• 方式1：调整微分方程的参数顺序，将def f(t,xs)修改为def f(xs, t)
• 方式2：调用odeint时增加tfirst=True参数，明确告知odeint你的微分方程是时间在前、状态量在后，将y = odeint(f,y0,t)修改为y = odeint(f,y0,t, tfirst=True)

修复参数顺序后你还需要调整代码里的其他逻辑错误，否则仍会运行异常：

• 大气密度函数rhos(h)的输入应为海拔高度，你当前传入的r是地心距，需要改为rhos(r - Ro)
• 升力系数clnew、阻力系数cdnew的输入应为攻角，你当前传入的是航迹倾角gamma，二者物理定义不同，需要确认对应关系是否正确
• 运动方程存在多处括号位置错误，例如np.sin(gamma/r**2)应为np.sin(gamma)/r**2、np.cos(gamma/(V*r))应为np.cos(gamma)/(V*r)、np.tan(phi/r)应为np.tan(phi)/r，请对照原始运动方程核对修正

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