在 lmfit 中访问 params 属性

Question

我正在尝试使用文档网站中的 lmfit 完成一个示例：

import sys

import numpy as np
import matplotlib.pyplot as plt
from matplotlib.patches import Rectangle
from matplotlib.backends.backend_pdf import PdfPages

from lmfit import minimize, Parameters, Parameter, report_fit, Minimizer, conf_interval, conf_interval2d, printfuncs
from lmfit import Model

def main():
    def decay(t, N, tau):
        return N*np.exp(-t/tau)

    decay_model = Model(decay)
    print decay_model.independent_vars
    for pname, par in decay_model.params.items():
        print pname, par

if ___name___ == "___main___":
    main()

当我执行最后一条命令时，出现以下错误：

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
AttributeError: 'Model' object has no attribute 'params'

params好像不是Model的属性。有人可以解释为什么这段代码不起作用吗？我在多台计算机和 Python 版本上试过这个。运行 Python 2.7.9 32 位，Windows 7.

Answer 1

可能会导致问题的是 - 假设这确实是您的整个代码 - 您没有使用任何数据进行拟合，也没有在任何地方初始化参数字典。下面是一个最小的例子，它适用于 0.8.3 版本（见下面的输出）；如果你使用 0.9.x 你将不得不稍微调整一部分（检查 here 从 0.8.3 到 0.9.x 做了哪些更改）。

import numpy as np
import matplotlib.pyplot as plt
from lmfit import minimize, Parameters, Parameter, report_fit

# generate some data with noise
# replace xData and yData with your data
xData = np.linspace(0., 100., 50.)
Nf = 5.
tauf = 6.5
yData = Nf * np.exp(-xData / tauf) + np.random.normal(0, 0.5, len(xData))
# plt.plot(xData, yData, 'bo')
# plt.show()

def decay(params, x, data):

    N = params['N'].value
    tau = params['tau'].value

    model = N * np.exp(-x/tau)
    return model - data # that's what you want to minimize

# create a set of Parameters
params = Parameters()
params.add('N', value=10) # value is the initial value
params.add('tau', value=8.)

# do fit, here with leastsq model
result = minimize(decay, params, args=(xData, yData))

# calculate final result
final = yData + result.residual

# write error report
report_fit(params)

# plot the data
plt.plot(xData, yData, 'bo')
plt.plot(xData, final, 'r')
plt.show()

首先，我使用您的衰减函数创建了一些数据并向其添加了一些噪声；这部分应该替换为您的数据。剩下的部分很简单，你会得到以下输出（由于噪声被添加到数据中，确切的值可能会有所不同）：

[[Variables]]
    N:     5.15685000 +/- 0.419115 (8.13%) (init= 10)
    tau:   6.58557758 +/- 0.877337 (13.32%) (init= 8)

如您所见，确定的参数与我为数据生成选择的参数非常接近。绘制数据和拟合结果如下：

lmfit 的好处是您可以为参数值使用边界。您可以查看示例 here.

希望这能让你入门，如果你有任何其他问题，请告诉我。

在 lmfit 中访问 params 属性

Accessing params attribute in lmfit

python

attributes

curve-fitting

python-2.7

lmfit