Python如何代入给定三角函数方程求解x？sympy方程数值求解方法咨询

这个方程属于超越方程，不存在通用的闭式解析解，所以sympy的solve()无法返回符号解，你可以通过以下两种常用方法代入C、D的固定值求解数值解：

方法1：使用sympy内置的nsolve求解 #

直接在现有符号体系下代入参数、指定初始猜测值即可求解，适合快速验证：

from sympy import symbols, sin, sqrt, nsolve

x = symbols('x')
D = symbols('d')
C = symbols('c')
# 注意原代码中的sin^2(x)需要修正为Python支持的sin(x)**2写法
eq1 = 2*sin(x)+sqrt(1-sin(x)**2)+D*sin(x/2)-sin(2*x)-1-C

# 代入固定参数示例：C=2，D=3，初始猜测值设为1弧度（可根据实际场景调整）
C_val = 2
D_val = 3
sol = nsolve(eq1.subs({C:C_val, D:D_val}), x, 1)
print(f"数值解x={sol.evalf()}（弧度），换算角度为{sol.evalf()*180/3.1415926}°")

注意：初始猜测值会影响最终返回的解，因为超越方程可能存在多个根，你可以根据实际业务场景调整初始值，或者限定求解区间。

方法2：使用scipy的数值求解工具（更灵活） #

如果需要更灵活的区间限定、求解算法选择，可以把符号方程转为数值函数，用scipy的优化工具求解：

import numpy as np
from scipy.optimize import root_scalar

# 定义数值化的方程，C、D为入参
def equation(x, C, D):
    return 2*np.sin(x) + np.sqrt(1 - np.sin(x)**2) + D*np.sin(x/2) - np.sin(2*x) - 1 - C

# 代入参数示例：C=2，D=3，限定求解区间为0到π（可根据实际场景调整）
C_val = 2
D_val = 3
res = root_scalar(equation, args=(C_val, D_val), bracket=[0, np.pi], method='bisect')
if res.converged:
    print(f"数值解x={res.root}（弧度），换算角度为{res.root*180/np.pi}°")
else:
    print("求解未收敛，请调整求解区间或参数")

注意事项 #

• 公式中的sqrt(1-sin²(x))等价于abs(cos(x))，如果你的实际应用场景中cos(x)的符号是确定的，可以直接替换为cos(x)或者-cos(x)，避免求解时出现分支跳转导致的误差。
• 方程存在多个解的可能性很高，需要结合你的实际物理场景限定x的合理取值范围，避免得到不符合业务逻辑的解。

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