SymPy求解ODE遇NotImplementedError：初始条件导致常数解过多

问题背景 #

尝试求解微分方程 $\dot{x}(t) = \sqrt{1 - x(t)^2}$（初始条件 $x(0)=0$），编写SymPy代码后出现报错：

NotImplementedError: Initial conditions produced too many solutions for constants

已知该方程的解为分段函数，需要修正代码以正确求解。

原代码 #

import sympy as sp

t = sp.symbols('t', real=True)
x = sp.Function('x', real=True)

diffeq = sp.Eq(x(t).diff(t), sp.sqrt(1 - x(t)**2)) 
res = sp.dsolve(diffeq, ics={x(0): 0})

报错原因 #

该微分方程的通解包含多分支（三角函数周期性解、常数解$x=1$/$x=-1$），初始条件$x(0)=0$对应多个常数解，SymPy无法自动选择合适分支并生成分段解，因此抛出错误。

修正方案 #

手动求解通解、确定常数后构造分段函数，代码如下：

import sympy as sp

t = sp.symbols('t', real=True)
x = sp.Function('x', real=True)

# 定义微分方程
diffeq = sp.Eq(x(t).diff(t), sp.sqrt(1 - x(t)**2)) 

# 求解通解
general_sol = sp.dsolve(diffeq)
print("通解：", general_sol)

# 代入初始条件确定常数（以通解x(t)=sin(t+C)为例）
C = sp.symbols('C')
sol_expr = sp.sin(t + C)
eq = sp.Eq(sol_expr.subs(t, 0), 0)
C_solutions = sp.solve(eq, C)
print("常数C的可能解：", C_solutions)

# 构造符合物理意义的分段解
# 因方程右边平方根非负，导数≥0，解在[-π/2, π/2]为sin(t)，超出区间后保持±1
piecewise_sol = sp.Piecewise(
    (-1, t < -sp.pi/2),
    (sp.sin(t), sp.And(t >= -sp.pi/2, t <= sp.pi/2)),
    (1, t > sp.pi/2)
)

print("\n最终分段解：")
sp.pprint(piecewise_sol)

# 验证解的正确性
print("\n验证初始条件x(0)=0：", piecewise_sol.subs(t, 0))
print("验证微分方程是否成立：", sp.simplify(sp.Eq(piecewise_sol.diff(t), sp.sqrt(1 - piecewise_sol**2))))

说明 #

• 通解中的常数$C$有无限多个解（$C=k\pi, k\in\mathbb{Z}$），选择$C=0$得到单调递增的$\sin(t)$分支，结合方程导数非负的特性，扩展为分段函数。
• 分段函数在各区间均满足原微分方程：当$|t|>\pi/2$时，$x(t)=\pm1$，导数为0，$\sqrt{1-x^{2}=0$，等式成立；当$|t|\leq\pi/2$时，$\sin(t)$的导数为$\cos(t)$，而$\sqrt{1-\sin}2(t)}=|\cos(t)|=\cos(t)$（该区间内$\cos(t)\geq0$），等式成立。

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