二阶微分方程常数替换求助：已求常数但不知如何替换C1

看起来你在处理二阶微分方程的初始条件时，犯了一个常见的小错误——你错误地替换了符号变量x，而不是时间变量t。让我帮你修正代码，并且解释每一步的逻辑：

首先，你的微分方程描述的是受重力与阻力作用的运动问题：$\frac{d^2x}{dt2} = g - \frac{K}{m}\frac{dx}{dt}$，初始条件是$t=0$时$x(0)=0$且$v(0)=\frac{dx}{dt}|_{t=0}=0$。

原代码的问题分析： #

• y_res.subs(x, 0)：这里x是你定义的符号，但x(t)是关于t的函数，你需要代入的是t=0，而非x=0
• 求解C2时，你先替换x再求导，逻辑颠倒了——应该先对位移函数求导得到速度表达式，再代入t=0

修正后的代码： #

import sympy as sym

# 定义符号与函数，明确x是关于t的函数
x, g, K, m, t = sym.symbols('x, g, K, m, t', positive=True, real=True)
x_t = sym.Function('x')(t)

# 定义二阶微分方程
dxx = sym.Eq(x_t.diff(t, 2), g - (K/m)*x_t.diff(t))

# 求解微分方程通解
solution = sym.dsolve(dxx)
y_res = solution.args[1]

# 处理初始条件1：t=0时速度为0
v_res = sym.diff(y_res, t)  # 先求导得到速度函数
C2_eq = sym.Eq(v_res.subs(t, 0), 0)
C2_val = sym.solve(C2_eq, sym.symbols('C2'))[0]

# 处理初始条件2：t=0时位移为0，先代入已求得的C2
y_res_with_C2 = y_res.subs(sym.symbols('C2'), C2_val)
C1_eq = sym.Eq(y_res_with_C2.subs(t, 0), 0)
C1_val = sym.solve(C1_eq, sym.symbols('C1'))[0]

# 得到满足初始条件的特解
final_solution = y_res.subs({'C1': C1_val, 'C2': C2_val})

display(C1_val, C2_val, final_solution)

关键修正点： #

1. 用sym.Function('x')(t)明确定义x_t，避免符号与函数的概念混淆
2. 处理初始条件时，始终替换t=0，而非x=0，符合物理问题的初始条件逻辑
3. 按顺序求解常数：先解出C2，再将C2代入位移表达式后求解C1，保证计算逻辑的正确性

这样你就能正确得到C1和C2的具体值，并且替换回原通解中，得到满足初始条件的特解了。

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