0

我喜欢将语法保持在最通用的形式，通过iminuit使用以下参数给出它的参数*：

import iminuit
import numpy as np

x_data = np.array([0,1,2,3,4,5,6,7,8,9])
y_data = np.array([0,1,2,3,4,5,4,3,2,1])

def fit_function(x, *p):
    return p[0]*np.exp(-((x-p[1])**2/(2*p[2]**2)))

def minimize_me(*p):
    return sum((fit_function(x, *p) - y)**2 for x, y in zip(x_data, y_data))

p=[4.5, 5, 0.4]
print(minimize_me(*p)) # works!  --> gives: 57.1645229329

m = iminuit.Minuit(minimize_me, *p)
m.migrad() # fails!

失败并出现错误：

AttributeError: 'float' object has no attribute 'print_banner'

任何想法我做错了什么？谢谢。

ps这个例子是基于这个SO帖子：https ://stackoverflow.com/a/22540079/5177935

4

1 回答 1

1

有什么问题？

您正在调用Minuit初始化程序，如下所示：

p = [4.5, 5, 0.4]
Minuit(minimize_me, *p)

这相当于：

Minuit(minimize_me, 4.5, 5, 0.4)

即在 Python 中，星号导致参数列表解包，在这种情况下，将浮点数作为不应该是浮点数的参数的位置参数传递：

Minuit(fcn=minimize_me, throw_nan=4.5, pedantic=5, frontend=0.4)

调用Minuit()不正确应该立即失败并给出一个好的错误信息。它目前没有，因为在初始化程序中没有实现输入验证。感谢您在https://github.com/iminuit/iminuit/issues/189报告此问题。

怎么做？

在您的情况下，您不关心参数名称。但是 Minuit 需要为每个参数命名。这是内部数据结构的一部分，例如用于报告拟合结果。

这是处理此问题的通用方法：

p_vals = [4.5, 5, 0.4]
p_names = ['par_{}'.format(_) for _ in range(len(p_vals))]

m = iminuit.Minuit(
    fcn=minimize_me,
    forced_parameters=p_names,
    **dict(zip(p_names, p_vals))
)
print(m.parameters)
m.migrad()  # works!
于 2016-11-02T21:10:44.073 回答