使用SymPy求解常微分方程时积分常数C1计算错误如何解决

错误排查 #

• 核心错误为变量名拼写不匹配：你定义微分方程时用的变量名为eq，但调用sym.dsolve()求解时传入的是未定义的变量eqn，变量名不一致导致输入异常，最终得到错误的常数计算结果。
• 额外优化点：定义符号时可以明确指定为实数类型，避免sympy默认返回复数域的多解，减少结果歧义。

修正后代码 #

import sympy as sym

sym.init_printing()
# 明确指定符号为实数，过滤无效复数解
t,m,k,g,v0,psi = sym.symbols("t m k g v_0 psi", real=True)

vz = sym.Function('v_z', real=True)(t)
a = vz.diff(t)

eq = sym.Eq(-k*vz**2 - m*a - m*g, 0)
# 修正变量名，把eqn改为正确的eq
sol = sym.simplify(sym.dsolve(eq)).rhs

constant = sym.solve([sol.subs(t, 0) - v0*sym.sin(psi)], dict=True)

可选简化写法 #

你也可以直接在dsolve时传入初始条件，不需要单独求解常数：

sol_final = sym.dsolve(eq, ics={vz.subs(t, 0): v0*sym.sin(psi)})

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