Python中SymPy求解复杂方程无响应问题求助

核心原因 #

• 超越方程无解析解，SymPysolve无法处理：你的方程包含三角函数、指数函数、平方根等复杂非线性运算，属于超越方程——这类方程几乎不存在解析解。SymPy的solve函数主要针对多项式或简单超越方程设计，遇到这类复杂方程会持续尝试推导解析解，最终陷入无限循环导致程序挂起。
• 函数参数不匹配：定义的fun_v2仅接受4个参数（pi, f, T_pi, gamma），但调用时传入了5个参数（pi, f, T_pi, gamma, tau），这会引发参数错误，若未提前报错，可能导致后续计算逻辑混乱，加剧程序异常。
• 混用NumPy与SymPy常量：你使用NumPy的pi和e参与SymPy的符号计算，两类库的常量类型不同，会增加符号计算的复杂度，可能导致内部处理逻辑卡顿。
• 工具选择错误：solve适合求解析解，而工程中这类复杂方程更适合用数值求解方法，而非符号求解。

解决建议 #

1. 修正函数参数：调整fun_v2的定义，加入tau参数，确保调用与定义一致。
2. 统一符号计算常量：替换NumPy的pi和e为SymPy自带的sympy.pi和sympy.E，避免类型冲突。
3. 改用数值求解：放弃SymPy的solve，使用SciPy的数值求解器（如root_scalar）或SymPy的nsolve来寻找数值解。

修正后的示例代码 #

from scipy.optimize import root_scalar
import sympy as sp

# 使用SymPy原生常量
pi = sp.pi
E = sp.E

# 修正函数参数，加入tau
def fun_v2(pi, f, T_pi, gamma, tau):
    """v2计算公式参考NEN-EN-IEC 60865-1:2012的A7部分"""
    part1 = 1 - (sp.sin(4*pi*f*T_pi-2*gamma)+sp.sin(2*gamma))/(4*pi*f*T_pi)
    part2 = (f*tau)/(f*T_pi)*(1-E**-(2*f*T_pi/(f*tau)))*(sp.sin(gamma))**2
    part3 = 8*pi*f*tau*sp.sin(gamma)/(1+(2*pi*f*tau)**2)
    part4 = 2*pi*f*tau*sp.cos(2*pi*f*T_pi-gamma)/(2*pi*f*T_pi)
    part5 = sp.sin(2*pi*f*T_pi-gamma)/(2*pi*f*T_pi)
    part6 = E**(-f*T_pi/(f*tau))
    part7 = (sp.sin(gamma)-2*pi*f*tau*sp.cos(gamma))/(2*pi*f*T_pi)
    return part1 + part2 - part3 * ((part4 + part5) * part6 + part7)

# 计算已知参数值
I_p   = 50000
mu_0  = 4*pi*10**-7
m_s   = 2.402
n_cond= 2
f     = 50
a_s   = 0.1
d     = 0.03624
kappa = 1.999999999
kappa = max(1.1, kappa)
tau = (-2*pi*f/3  *  sp.log((kappa-1.02)/0.98))**-1
gamma = sp.atan(2*pi*f*tau)
v1 = f/(sp.sin(pi/n_cond)) * sp.sqrt((a_s-d)*m_s/(mu_0/(2*pi) * (I_p/n_cond)**2 * (n_cond-1)/a_s))

# 转换为数值类型，用于数值求解
v1_num = float(v1)
tau_num = float(tau)
gamma_num = float(gamma)

# 定义待求解的数值方程
def equation(v2):
    T_pi = v1_num/(f*sp.sqrt(v2))
    v2_eval = fun_v2(pi, f, T_pi, gamma_num, tau_num)
    return float(f*T_pi*sp.sqrt(v2_eval) - v1_num)

# 使用Brentq方法数值求解，需给出合理的根的范围
result = root_scalar(equation, bracket=[0.01, 100], method='brentq')
print(f"v2的解为: {result.root}")

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