时间戳

首页 > 解决方案 > 如何在使用 lmfit 进行最小化时修复“函数返回的数组在调用之间改变了大小”？

python - 如何在使用 lmfit 进行最小化时修复“函数返回的数组在调用之间改变了大小”？

问题描述

如何在使用lmfit最小化时修复错误“函数返回的数组在调用之间更改大小”的代码？

请在下面找到我的代码：

import numpy as np
import pandas as pd
import lmfit as lf

#model needs to be fitted
x0 = 75
def func(params, x, Tsky):
    A = params['amp']
    w = params['width']
    t = params['thickness']
    v0 = params['mid_freq']
    b0 = params['b0']
    b1 = params['b1']
    b2 = params['b2']
    b3 = params['b3']
    b4 = params['b4']
    B = (4 * (x - v0)**2. / w**2.) * np.log(-1./t * np.log((1 + np.exp(-t))/2))
    T21 = -A * (1 - np.exp(-t * np.exp(B)))/(1 - np.exp(-t))
    model = T21 + b0 * ((x/x0)**(-2.5 + b1 + b2 * np.log(x/x0))) * np.exp(-b3*(x/x0)**-2.) + b4 *  (x/x0)**-2.
    return (Tsky-model)

#read the data
df = pd.read_csv('figure1_plotdata.csv')
data_list = df.T.values.tolist()
xdata = np.array(data_list[0])
Tsky = np.array(data_list[2])

#initial value of the parameters
params = lf.Parameters()
params.add('amp', value=0.2)
params.add('width', value=10)
params.add('thickness', value=5)
params.add('mid_freq', value=70)
params.add('b0', value=500)
params.add('b1', value=-0.5)
params.add('b2', value=-0.5)
params.add('b3', value=-0.5)
params.add('b4', value=500)

#minimize the function
out = lf.minimize(func, params, args=(xdata, Tsky), method='leastsq', kws=None, iter_cb=None, scale_covar=True, nan_policy='omit', calc_covar=True)

print(lf.fit_report(out))

这是错误消息：

File "/home/ankita/Dropbox/Python/Bowman_work/min.py", line 81, in <module>
    out = lf.minimize(func, params, args=(xdata, Tsky), method='leastsq', kws=None, iter_cb=None, scale_covar=True, nan_policy='omit', calc_covar=True)

  File "/home/ankita/anaconda3/lib/python3.7/site-packages/lmfit-0.9.13-py3.7.egg/lmfit/minimizer.py", line 2300, in minimize
    return fitter.minimize(method=method)

  File "/home/ankita/anaconda3/lib/python3.7/site-packages/lmfit-0.9.13-py3.7.egg/lmfit/minimizer.py", line 1949, in minimize
    return function(**kwargs)

  File "/home/ankita/anaconda3/lib/python3.7/site-packages/lmfit-0.9.13-py3.7.egg/lmfit/minimizer.py", line 1492, in leastsq
    lsout = scipy_leastsq(self.__residual, variables, **lskws)

  File "/home/ankita/anaconda3/lib/python3.7/site-packages/scipy/optimize/minpack.py", line 394, in leastsq
    gtol, maxfev, epsfcn, factor, diag)

**ValueError: The array returned by a function changed size between calls**

标签： pythonanacondaminimizationlmfit

解决方案

如果你用过

out = lf.minimize(func, params,...,nan_policy='raise')

你会看到一个异常被抛出，告诉你有 NaN。使用 时nan_policy='omit'，模型生成的任何此类 NaN 都会从残差数组中删除，因此数组的大小在调用之间会发生变化。配件无法处理 NaN 或数组大小的更改——您必须消除它们。

特别np.log(x)是 是 NaN 时x<0。你有np.log()一个复杂的参数，它取决于参数的值t。如果对于 的某个值，该参数低于 0 t，则该模型具有 NaN 并且没有任何意义。您必须确保此参数不能低于 0。可能是使用

params.add('thickness', value=5, min=0)

足够了。但是您应该更详细地检查您的模型并确定这是否有意义。

你的模型对我来说看起来很复杂。我猜不出这样的模型是从哪里来的。取多个np.exp()andnp.log()有点要求数值不稳定性。所以，我不知道仅仅强迫t自己积极会很好地适应，但它可能会为你指明正确的方向。

推荐阅读