python 曲线拟合参数能否为数组_Scipy曲线拟合全局参数

我正在尝试用python中的Scipy curve_fit拟合一个星系模型的两个全局参数。我有一个自变量数组和一个因变量数组。数据集的前1/4需要根据两个全局参数和两个局部参数来拟合函数，下一个季度需要拟合另一个函数，这取决于两个全局参数和两个局部变量等

不管怎样，我能写一个函数来调用适当的函数，这个函数有正确的索引和整个数组的全局参数。在

到目前为止，我得到的是：def galaxy_func_inner(time,a,b,c,d):

telescope_inner = lt.station(rot_angle=c,pol_angle=d)

power = telescope_inner.calculate_gpowervslstarray(time)[0]

return a*np.array(power)+b

def galaxy_func_outer(time,a,b,c,d):

telescope_outer = lt.station(rot_angle=c,pol_angle=d)

power = telescope_outer.calculate_gpowervslstarray(time)[0]

return a*np.array(power)+b

def galaxy_func_global(time,R,P,a,b,c,d,e,f,g,h):

for t_index in range(len(time)):

if t_index in range(0,50):

return galaxy_func_outer(t_index,a,b,R,P)

elif t_index in range(50,100):

return galaxy_func_outer(t_index,c,d,R,P)

elif t_index in range(100,150):

return galaxy_func_inner(t_index,e,f,R,P)

elif t_index in range(150,200):

return galaxy_func_inner(t_index,g,h,R,P)

问题是，这种方法只适用于第一次，而不适用于整个时间序列，并且单点只适合于相应的模型点，而不是整个阵列。你能帮我重新制定这个吗？我试图将其重新表述为：

^{pr2}$

但我得到了一个错误：File "galaxy_calibration.py", line 117, in

popt,pcov = curve_fit(galaxy_func_global,xdata,ydata)

File "/Library/Python/2.7/site-packages/scipy-0.14.0.dev_7cefb25-py2.7-macosx-10.9-intel.egg/scipy/optimize/minpack.py", line 555, in curve_fit

res = leastsq(func, p0, args=args, full_output=1, **kw)

File "/Library/Python/2.7/site-packages/scipy-0.14.0.dev_7cefb25-py2.7-macosx-10.9-intel.egg/scipy/optimize/minpack.py", line 369, in leastsq

shape, dtype = _check_func('leastsq', 'func', func, x0, args, n)

File "/Library/Python/2.7/site-packages/scipy-0.14.0.dev_7cefb25-py2.7-macosx-10.9-intel.egg/scipy/optimize/minpack.py", line 20, in _check_func

res = atleast_1d(thefunc(*((x0[:numinputs],) + args)))

File "/Library/Python/2.7/site-packages/scipy-0.14.0.dev_7cefb25-py2.7-macosx-10.9-intel.egg/scipy/optimize/minpack.py", line 445, in _general_function

return function(xdata, *params) - ydata

ValueError: operands could not be broadcast together with shapes (4,) (191,)

任何帮助都将不胜感激。在