我有一个拟合函数,其形式为:
def fit_func(x_data, a, b, c, N)
其中,a、B、c是长度为N的列表,其中的每个条目都是要在scipy.optimiz.curve_fit()中优化的可变参数,并且N是用于循环索引控制的固定数。
在this question之后,我想我能够修复N,但是我现在调用curve_fit如下:
params_0 = [a_init, b_init, c_init]
popt, pcov = curve_fit(lambda x, a, b, c: fit_func(x, a, b, c, N), x_data, y_data, p0=params_0)
我收到一个错误:lambda()正好接受Q个参数(给定P)
Q和P的变化取决于我的设置方式。
那么:对于初学者来说,这是可能的吗?我可以将列表作为参数传递给curve_fit,并具有我所希望的行为,即它将列表元素作为单独的参数处理吗?假设答案是肯定的,我在函数调用中做错了什么?
4条答案
按热度按时间vaqhlq811#
这里的解决方案是编写一个 Package 器函数,它接受你的参数列表,并将其转换为fit函数理解的变量。这实际上是唯一必要的,因为我正在处理别人的代码,在一个更直接的应用程序中,这将不需要 Package 器层。基本上
要固定N,必须在curve_fit中调用它,如下所示:
其中
jexiocij2#
我可以用稍微不同的方法解决同样的问题。我使用scip.optimize.least_squares来解决,而不是curv_fit。我已经在链接-https://stackoverflow.com/a/60409667/11253983下讨论了我的解决方案
unftdfkk3#
请看这篇文章https://stackoverflow.com/a/73951825/20160627,我建议使用scipy.optimize.curve_fit,其中包含任意数量和位置的参数,以便在列表中进行拟合或修复
q0qdq0h24#
可能不是一个最优的解决方案,但我遇到了类似的问题。我解决它的方法是硬编码“传递”函数,字面上输出我的广义函数,但需要特定数量的参数。