本文介绍了在solveset()中使用erf()函数时，Sysmy和mpath给出&q；类型错误：无法创建MPF&q；的处理方法，对大家解决问题具有一定的参考价值，需要的朋友们下面随着小编来一起学习吧！
问题描述
我有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
我想使用mpath和sendy来求解该方程，但它给出了以下错误消息：
类型错误：无法从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是由mpath创建的浮点类型。
我想我已经将问题范围缩小到erf()函数，它返回
EmptySet()
当我运行时
solveset(log(Xmax) - (log(mode) + s ** 2), s, domain=S.Reals)
我想不出下一步该尝试什么，如果有任何帮助，我将不胜感激！
我最初认为这是数学上的问题，但在MatLab中成功地解出了这个方程，所以问题很可能是来自sendy或mpath。
推荐答案
丢失mpmatherf的导入并使用症状版本可以解决您的错误。
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)
另请注意：
- 您不必从
mpmath导入log和sqrt。从sympy获取它们不会对您的结果产生影响 - 您可以在变量s上指定Real域，从而省去了调用
solveset时的麻烦。
如果需要，请参阅包测试here中的进一步用法示例。
这篇关于在solveset()中使用erf()函数时，Sysmy和mpath给出&q；类型错误：无法创建MPF&q；的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持编程学习网！
本站部分内容来源互联网,如果有图片或者内容侵犯您的权益请联系我们删除！