如何使用FiniteSet的元素进行以下计算？

Question

我试图用solveset（方程，变量）来解一个方程，这将得到变量的值作为结果。结果得到一个FiniteSet{value}。这个FiniteSet可以访问吗？这样我就可以在下一次计算中使用这个值（可能是float类型的）？你知道吗

问题如下： 我定义了三个SymPy符号（drift1, drift2, drift3 = symbols('drift1 drift2 drift3')），它们应该是下面计算的占位符。他们没有价值观。它们是两个方程的一部分（equation1和equation2，都是drift1, drift2和drift3的函数）。 我的目标是编程一个迭代循环，通过求解这两个方程来优化所有三个符号drift1、drift2和drift3，并将结果用于下一个迭代步骤。为此，我定义了三个Python变量x, y和z，a和y的初始值是0.0188和0.0099。z是x和y的函数，每次迭代都会计算新的值。 在第一步中，z是用x和y（以及我在代码中定义的其他变量）计算的，我用等式1中的y和z替换了drift2和drift3（等式现在是drift1的函数）。通过对drift1激活equation1=0，我使用了solveset(equation1.subs(drift2, y).subs(drift3, z), drift1)，得到了一个FiniteSet，它包含了结果。我将用这个结果通过solveset(equation2.subs(drift1, result from FiniteSet).subs(drift3, z), drift2)来解第二个方程。这不起作用，因为结果不是float变量，而是sympy.sets.sets.FiniteSet类型。你知道吗

有人能告诉我如何在下一次计算中使用FiniteSet（'inside'the FiniteSet；{value}）的值吗？我用Pyhton 3和SymPy。你知道吗

import numpy as np
from sympy import *
import scipy as sp

drift1, drift2, drift3 = symbols('drift1 drift2 drift3')
init_printin(use_unicode=True)

equation1 =  0.005*(drift1*(-23.88*drift2*(15.16*drift3 + 1.15)-20.31*drift3
                       +00.47*drift2*(15.16*drift3+1.15)+1.03*drift3+0.05))/drift1

equation2 = 0.005*(drift1*(35.73*drift2*(-13.79*drift3 + 0.86)+7.91*drift3+2.09)
                       +1.65*drift2*(-13.79*drift3+0.86)+0.75*drift3+0.07)/drift1

x, y = 0.0188, 0.0099
gesamtlaenge = 1.676
laenge1 = 0.04
laenge2 = 0.02

# the loop for the iteration is not included in this code
# the only thin I wanna know is how to substitute drift1 by the solution of 
# the first calculation sol1

z = gesamtlaenge-laenge1-laenge2-x-y
expr1 = equation1.subs(drift2, y).subs(drift3, z)
sol1 = solveset(equation, drift1)
# up to this point, the code is running, but now, I get a FiniteSet {0.047} as
# result for sol1 and I wanna use this value to substitute drift1 in the next calculation
expr2 = equation2.subs(drift1, sol1).subs(drift3, z)
sol2 = solveset(equation, drift2)
print(drift1, drift2, drift3)

当我执行expr2 = equation2.subs(drift1, sol1).subs(drift3, z)时，我得到等式2作为drift1和drift2的函数，因为它不能使用FiniteSet{0047}中的值（0.047是我需要的值）。我明白为什么会这样，但我没有解决办法。你知道吗

我很感谢你的帮助！干杯

Answer 1

您可以通过转换为列表来访问FiniteSet的元素。在这种情况下，如果有一种解决方案，list(sol1).pop()将起作用：

>>> sol1 = FiniteSet(1)
>>> x.subs(x, list(sol1).pop())
1