使用SymPy仅能解析求解一个微分方程，求另一个的可行解法

问题背景 #

我在解决问题时遇到两个微分方程，尝试用SymPy进行解析求解，仅成功解出其中一个，想找到求解另一个方程的可行方法。尝试的代码如下：

from sympy import symbols, Function, dsolve

x = symbols('x')
y = Function('y')(x)

# 已成功求解的方程
eq1 = y.diff(x) + y - x
sol1 = dsolve(eq1)
print("第一个方程的解：", sol1)

# 未成功求解的方程（示例为y’² + y² = 1）
eq2 = y.diff(x)**2 + y**2 - 1
sol2 = dsolve(eq2)
print("第二个方程的解：", sol2)

可行求解方法 #

• 手动分离变量结合SymPy积分：对于像y’² + y² = 1这类可分离变量的非线性方程，先手动变形为dy/dx = ±√(1-y²)，分离变量后调用SymPy的积分功能计算：

  from sympy import integrate, Eq, sqrt
  
  # 处理正号分支
  left_pos = integrate(1/sqrt(1 - y**2), y)
  right_pos = integrate(1, x)
  sol_pos = Eq(left_pos, right_pos + symbols('C1'))
  print("正号分支解：", sol_pos)
  
  # 处理负号分支
  left_neg = integrate(1/sqrt(1 - y**2), y)
  right_neg = integrate(-1, x)
  sol_neg = Eq(left_neg, right_neg + symbols('C2'))
  print("负号分支解：", sol_neg)
  

  最终可整理得到y = sin(x + C)或y = sin(-x + C)形式的解析解。

• 指定SymPy求解方法：调用dsolve时指定method参数，针对非线性方程可尝试'separable'或'lie_group'方法：

  sol2 = dsolve(eq2, method='separable')
  print(sol2)
  

• 转换为标准方程形式：检查方程是否属于特殊类型（如黎卡提、伯努利方程等），通过变量替换转化为SymPy可直接求解的形式后再尝试。

• 数值求解补充验证：若确认无解析解，可使用SciPy的数值求解工具（如solve_ivp）得到数值解，辅助验证问题趋势：

  from scipy.integrate import solve_ivp
  import numpy as np
  
  def ode_func(x, y):
      return np.sqrt(1 - y**2)  # 选取正号分支计算
  
  x_interval = [0, 10]
  initial_val = [0]
  num_sol = solve_ivp(ode_func, x_interval, initial_val, t_eval=np.linspace(0, 10, 100))
  # 可通过matplotlib绘制数值解曲线
  

内容的提问来源于stack exchange，提问作者Kazi Abu Rousan