一种实现方法是改用scipy.optimize.leastsq(curve_fit是对minimumsq的方便包装).
将x数据一维堆叠; y数据的同上. 3个独立数据集的长度甚至无关紧要;我们称它们为n1,n2和n3,这样您的新x和y将具有形状(n1 n2 n3,).
在要优化的功能内,您可以方便地拆分数据.它不是最好的功能,但是可以正常工作:
def function(x, E, T, a, n, m):
return x/E + (a/n+1)*T^(n+1)*x^m
def leastsq_function(params, *args):
a = params[0]
n = params[1]
m = params[2]
x = args[0]
y = args[1]
E = args[2]
T = args[3]
n1, n2 = args[2]
yfit = np.empty(x.shape)
yfit[:n1] = function(x[:n1], E[0], T[0], a, n, m)
yfit[n1:n2] = function(x[n1:n2], E[1], T[1], a, n, m)
yfit[n2:] = function(x[n2:], E[2], T[2], a, n, m)
return y - yfit
params0 = [a0, n0, m0]
args = (x, y, (E0, E1, E2), (T0, T1, T2), (n1, n1+n2))
result = scipy.optimiz