Question

我有4个输入变量（浮点数）：

的Xmax
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）/ S

我的代码如下：

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。

Answer 1

丢失mpmath的{{1}}导入，并使用sympy版本可以解决您的错误。

erf

另请注意：

您无需从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
您可以在变量s上指定Real域，保存您在sympy调用中执行此操作。

如果您需要，可以在包测试here中进一步使用。

sympy和mpmath给出＆＃34; TypeError：无法创建mpf＆＃34;在solveset（）中使用erf（）函数时

1 个答案: