0

如何使用 Minimizer 对象来最小化函数调用的数量。我从日志中得到这条消息：

函数调用过多（最大值设置为 %i）！使用：minimize(func, params, ..., maxfev=NNN) 或设置 leastsq_kws[\'maxfev\'] 来增加这个最大值。

from numpy import sqrt, pi, exp, loadtxt
from lmfit import  Model
from lmfit import Minimizer
import matplotlib.pyplot as plt

data = loadtxt('data/model1d_gauss.dat')
x = data[:, 0]
y = data[:, 1]

def gaussian(x, amp, cen, wid):
    "1-d gaussian: gaussian(x, amp, cen, wid)"
    return (amp/(sqrt(2*pi)*wid)) * exp(-(x-cen)**2 /(2*wid**2))


gmodel = Model(gaussian)

result = gmodel.fit(y, x=x, amp=5, cen=1, wid=1,method= 'least_squares')

print(result.fit_report())

r = result.fit_report()

plt.plot(x, y,         'bo')
plt.plot(x, result.init_fit, 'k--')
plt.plot(x, result.best_fit, 'r-')
plt.show()

lmfit 文档链接

Github 链接

我想通过最小化函数调用来加快拟合速度（无论如何，很多数据都是嘈杂的废话）

4

1 回答 1

1

不确定我是否会推荐这个，但是当您在内部使用 scipy 的 least_squares时，您可以只传递所需的选项，例如：

result = gmodel.fit(y, x=x, amp=5, cen=1, wid=1,method='least_squares',
                    fit_kws={'max_nfev': 100})

当然，这是假设 lmfit 在观察状态后会接受maxiter reached状态。但是从您的问题听起来，这仅被视为警告。

请记住，此优化器基于容差标准（请参阅文档），并且当不收敛给定x数量的步骤时，它实际上认为它仍然可以改进最小化！

如评论中所问：是的，您也可以更改这些标准，例如：

result = gmodel.fit(y, x=x, amp=5, cen=1, wid=1,method='least_squares',
                fit_kws={'ftol': 1-07,   # default 1e-08
                         'xtol': 1-07,   #         1e-08
                         'gtol': 1-07})  #         1e-08
于 2017-12-13T20:04:03.107 回答