SymPy调用solve_poly_system解含sqrt的方程组报PolynomialError

报错现象 #

调用sympy.solvers.solve_poly_system求解包含平方根函数的方程组时，触发报错：
sympy.polys.polyerrors.PolynomialError: sqrt(x1) contains an element of the set of generators

对应复现代码：

# 符号变量声明逻辑（该部分此前用于普通方程求解、导数计算均运行正常）
for x in range(len(coordArr)):
    coordArr[x][0] = symbols(xCoords[x], real = True)
    coordArr[x][1] = symbols(yCoords[x], real = True)

# 代入含sqrt的方程组调用多项式求解器
print(sympy.solvers.solve_poly_system(
    [5*(sympy.sqrt(coordArr[0][0])*2 + 2 - 3), 2*(coordArr[0][1] - 3)],
    coordArr[0][0], coordArr[0][1]
))

已按常见调整方案，将平方根计算替换为原生sympy.sqrt写法，报错仍持续触发。

报错根因 #

• solve_poly_system的设计定位是仅求解严格符合多项式定义的方程组，要求所有表达式中，求解变量仅以非负整数次幂、整数系数的组合形式存在，不支持变量的开平方（即1/2次幂）、倒数、三角函数、对数等非多项式运算。
• 代码中出现的sqrt(coordArr[0][0])本质是求解变量的1/2次幂运算，已经超出多项式范畴。求解器在解析方程组的生成元集合时，识别到平方根表达式包裹了核心求解变量，就会判定当前系统不满足多项式要求，抛出PolynomialError。
• 此前尝试的替换sympy.sqrt写法，只是调整了平方根的实现接口，并没有改变“表达式含变量非整数次幂”的本质，因此无法解决报错。

可落地解决方法 #

方法1：替换为适配非多项式系统的通用求解器 #

对于带根式的非多项式方程组，直接使用通用求解函数sympy.solvers.solve即可，该函数可自动完成带根式方程的化简、求解逻辑，不需要严格满足多项式约束。

调整后代码示例：

# 将solve_poly_system替换为通用solve函数
print(sympy.solvers.solve(
    [5*(sympy.sqrt(coordArr[0][0])*2 + 2 - 3), 2*(coordArr[0][1] - 3)],
    coordArr[0][0], coordArr[0][1]
))

注意：通用求解函数返回的解需要额外校验定义域，排除开平方操作可能引入的增根。

方法2：通过变量代换转化为标准多项式系统 #

如果需要使用solve_poly_system的高效多项式求解逻辑，可通过引入中间变量消去根号，把原方程组转化为纯多项式系统：

1. 引入中间变量u对应原表达式中的sqrt(x1)，补充恒等关系u**2 = x1，将根号运算转化为多项式形式的约束
2. 把原方程中所有sqrt(x1)替换为中间变量u，将额外的u**2 -x1 =0约束加入方程组
3. 对转化后的纯多项式方程组调用solve_poly_system求解，最后剔除中间变量，得到原求解变量的解

对应实现代码示例：

# 声明中间变量，约束其非负避免引入增根
u = symbols('u', real=True, positive=True)
x1 = coordArr[0][0]

# 构造转化后的纯多项式方程组
poly_eqs = [
    5*(u*2 + 2 - 3),  # 原方程中的sqrt(x1)替换为中间变量u
    2*(coordArr[0][1] - 3),
    u**2 - x1  # 补充开平方对应的多项式等价约束
]

# 求解转化后的多项式系统
sol = sympy.solvers.solve_poly_system(poly_eqs, u, x1, coordArr[0][1])
# 移除中间变量，得到原变量的最终解
final_sol = [(item[1], item[2]) for item in sol]

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