使用指数减法评估 symfit 模型时的 AttributeError

Question

我已经建立了一个模型来拟合一些指数减法的数据。它适合拟合，但是当我明确评估模型时，我得到了奇怪的结果和一个 AttributeError: 'Mul' object has no attribute 'exp'.

问题似乎源于 lambdify 的工作方式，并在评估包含重要函数（例如增强或指数）的模型时发生。

import symfit as sf

a = sf.Parameter('a',1,0,2) #name, initial value, minimum, maximum
b = sf.Parameter('b',1,0,2)
c=sf.Parameter('c',1,0,2)
x, y = sf.variables('x, y')

model=sf.Model({y: a * (c - sf.exp((x) / b))})
model(1,a,b,c)
Out: Ans(y=a*(c + 1))

model=sf.Model({y: a * (1 - sf.exp((x-c) / b))})
model(1,a,b,c)
Traceback (most recent call last):

  File "<ipython-input-9-faba3b52b923>", line 1, in <module>
    model(1,a,b,c)

  File "C:\PortablePrograms\Python\WPy-3670\python-3.6.7.amd64\lib\site-packages\symfit\core\fit.py", line 334, in __call__
    return Ans(*self.eval_components(**bound_arguments.arguments))

  File "C:\PortablePrograms\Python\WPy-3670\python-3.6.7.amd64\lib\site-packages\symfit\core\fit.py", line 296, in eval_components
    return [expr(*args, **kwargs) for expr in self.numerical_components]

  File "C:\PortablePrograms\Python\WPy-3670\python-3.6.7.amd64\lib\site-packages\symfit\core\fit.py", line 296, in <listcomp>
    return [expr(*args, **kwargs) for expr in self.numerical_components]

  File "C:\PortablePrograms\Python\WPy-3670\python-3.6.7.amd64\lib\site-packages\sympy\utilities\lambdify.py", line 444, in wrapper
    return funcarg(*newargs, **kwargsx)

  File "<string>", line 1, in <lambda>

AttributeError: 'Mul' object has no attribute 'exp'

我希望使用参数的初始值或得到一个象征性的答案，但两者都得不到。

因此，对于第一个模型，结果既错误又不一致，我希望

Ans(y=a*(c-exp(1/b))

或

Ans(y=a*(c-1))

或者只是简单地

Ans(y=0)

对于第二个模型，lambdify 似乎无法解析表达式。

使用sf.sqrt() 或sympy.exp() 时也是如此。我在 python 3.6.7 和 IPython 7.1.1 中工作，使用 Symfit 0.4.6 和 Sympy 1.1.1（因为 Symfit 不适用于根据 pip 的更高版本）。

我可以使用任一模型来拟合数据，并使用此处示例中所示的最佳拟合参数评估模型 (https://pypi.org/project/symfit/)。该行

yfit = model(x=xdata, **fit_result.params)[y]

也不起作用（在 python 2.7、3.5 和 3.6 上试过）除非更改为

yfit = [model(x=x, **fit_result.params) for x in xdata]

Answer 1

在 symfit 中调用 Model 意味着符号表达式被转换为 lambda 函数，即不再是符号的普通 python 函数。

在您的示例中，调用

sf.Model({y: a * (c - sf.exp((x) / b))})

相当于写作

lambda x, a, b, c: a * (c - np.exp((x) / b))

请注意，符号指数现已更改为 numpy。从这个例子可以看出模型不是用符号表达式调用的，而是用数字或数组调用的。

所以你的电话应该改为

model(x=1, a=1, b=1, c=1)

如果您需要在符号级别上工作，请继续直接使用 Model 中的表达式，该调用纯粹用于数值用途。

最后，文档中的示例假设 xdata 是一个数组，但在示例中它是一个列表。那应该（由我）更新。

最后一点，学会喜欢 python 中的关键字参数；)。以下是等价的：

a = sf.Parameter('a',1,0,2) #name, initial value, minimum, maximum
a = sf.Parameter('a', value=1, min=0, max=2)

但只有一个是自我记录的:)。

使用指数减法评估 symfit 模型时的 AttributeError

AttributeError when evaluating symfit model with substraction in exponent

python

sympy

attributeerror

symfit