四元二次方程组的求解

2023-12-18 13:21:56 浏览数 (1)

如图,由测量可得图中惠斯通电桥任意两个相邻端口之间的电阻,要求4个分立电阻的阻值。这种解方程组的问题可以用 sympy模块。代码如下

代码语言:javascript复制
# 4元2次方程组的计算。
# 应用在惠斯通电桥测电阻后求每个独立电阻的阻值。
import sympy


a = 1.75  # R1//(R2 R3 R4)
b = 4.00  # R2//(R1 R3 R4)
c = 1.75  # R3//(R1 R2 R4)
d = 3.00  # R4//(R1 R2 R3)
x1, x2, x3, x4, x_sum = sympy.symbols("x1, x2,x3,x4,x_sum",positive=True)
eq0 = x1   x2   x3   x4 - x_sum
eq1 = x1 * (x_sum - x1) - a * x_sum
eq2 = x2 * (x_sum - x2) - b * x_sum
eq3 = x3 * (x_sum - x3) - c * x_sum
eq4 = x4 * (x_sum - x4) - d * x_sum


s = sympy.solve([eq0, eq1, eq2, eq3, eq4], [x1, x2, x3, x4, x_sum], dict=True)
print(s)

0 人点赞