使用Python的fsolve求解三元非线性方程组出现语法错误如何解决？

1. 语法错误修复 #

你遇到的SyntaxError根本原因是括号没有成对闭合，同时存在库函数调用错误：

• eq1行末尾缺少2个右括号，需要补))
• eq2行末尾缺少1个右括号，需要补)
• Numpy没有内置np.cot余切函数，需要替换为1/np.tan(z)

2. 求解器适配说明 #

fsolve完全支持三元非线性方程组的求解，属于无约束求解器，本身不支持变量范围限制，如果需要保证x落在1/4~1/3之间，有两种方案：

• 调整初始猜测值，将x的初始值设为区间内的值（比如0.3），让求解过程优先靠近目标范围
• 换用支持边界约束的求解器，比如scipy.optimize.least_squares，显式指定x的取值边界

3. 修复后的可运行代码 #

import numpy as np
from scipy.optimize import fsolve

def equations(vars):
    x, y, z = vars
    # 补全括号、替换np.cot为对应实现
    eq1 = x/(1-x) - (5/(8*np.pi*np.sin(z)**2)) * (1 + 4*np.tan(z))
    eq2 = y/(1-y) - (5/(8*np.pi * np.cos(z))) * (1 - 1/np.tan(z))
    eq3 = np.tan(z) - ((1-x) / 3 * (1+y))
    return [eq1, eq2, eq3]

# 初始值x设为目标区间内的0.3，避免求解到范围外的解
x, y,z =  fsolve(equations, (0.3, 1, 1))

print(f"求解结果：x={x:.4f}, y={y:.4f}, z={z:.4f}")
# 校验x范围合规性
if 1/4 <= x <= 1/3:
    print("x值符合1/4~1/3的范围要求")
else:
    print("x值超出范围，建议调整初始值或换用带约束的求解器")

4. 带x范围约束的求解方案 #

如果必须严格限制x的取值范围，可以换用支持边界约束的least_squares求解器：

from scipy.optimize import least_squares
# 定义边界：x∈[0.25, 1/3]，y、z无限制
bounds = ([0.25, -np.inf, -np.inf], [1/3, np.inf, np.inf])
res = least_squares(equations, x0=(0.3, 1, 1), bounds=bounds)
x, y, z = res.x
print(f"带约束求解结果：x={x:.4f}, y={y:.4f}, z={z:.4f}")

额外注意事项 #

非线性方程组可能存在多个解，若求解结果不符合预期，可以多尝试不同的初始猜测值调整收敛方向；代码中未使用的pandas、matplotlib导入可以直接删除。

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