我有4个输入变量（浮点数）：
最大值
Xmin
百分比
模式
我想为s求解以下（相当长的）方程：
> (1/2+1/2*erf((log(Xmax)-(log(mode)+s**2))/(sqrt(2)*s))-(1/2+1/2*erf((log(Xmin)-(log(mode)+s**2))/(sqrt(2)*s))) - percentage == 0
我想使用mpmath和sympy来求解方程，但是它给出了以下错误消息：
TypeError：无法从0.707106781186547 *创建MPF *（-s ** 2-
0.287682072451781）/秒
我的代码如下：
from mpmath import erf, log, sqrt
from sympy import Symbol, solveset, S
percentage = 0.95
mode = 2
Xmin = 1.
Xmax = 1.5
s = Symbol('s')
eqn = (1/2+1/2*erf((log(Xmax)-(log(mode)+s**2))/(sqrt(2)*s))-(1/2+1/2*erf((log(Xmin)-(log(mode)+s**2))/(sqrt(2)*s))) - percentage)
solveset(eqn, s, domain=S.Reals)
mpf是由mpmath创建的浮点类型。
我想我将问题缩小到erf（）函数，它返回
EmptySet()
当我跑步
solveset(log(Xmax) - (log(mode) + s ** 2), s, domain=S.Reals)
我不知道下一步该怎么做，我们将不胜感激！
我最初以为这是数学问题，但方程式已在Matlab中成功解决，因此问题可能出在sympy或mpmath。
最佳答案
丢失mpmath的erf导入，并使用sympy版本可解决您的错误。
from sympy import Symbol, solveset, S, erf, log, sqrt
percentage = 0.95
mode = 2
Xmin = 1.
Xmax = 1.5
s = Symbol('s', real=True)
eqn = (1/2+1/2*erf((log(Xmax)-(log(mode)+s**2))/(sqrt(2)*s))-(1/2+1/2*erf((log(Xmin)-(log(mode)+s**2))/(sqrt(2)*s))) - percentage)
solveset(eqn, s)
另请注意：
您不必从
log导入sqrt和mpmath。从sympy获取它们对您的结果没有影响您可以在变量s上指定Real域，这样可以省去
solveset调用。如果需要，可在包测试here中使用进一步使用的示例。
关于python - 当在solveset()中使用erf()函数时，sympy和mpmath给出“TypeError:无法创建mpf”，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/37759214/