curve_fit多变量阵列非线性回归

Question

我正在尝试使用curve_fit拟合多元函数的系数。所有变量均为以下形状的数组：（1000）手动我可以按以下方式拟合曲线。首先，我定义函数，其中变量= [dphi1，dphi2，phi1，phi2，m]，并且系数= [c1，c2，c3，c4，c4，c5，c5，c6，c7，c8，c9]：

    import numpy as np 
    import matplotlib.pyplot as plt
    import tensorflow as tf
    import math
    
    # Function definition
    def ddphi1(dphi1,dphi2,phi1,phi2,M,c1,c2,c3,c4,c5,c6,c7,c8,c9):
        return (-(c1*np.sin(phi1-phi2)*np.cos(phi1-phi2)*dphi1**2)-(c2*np.sin(phi1-phi2)*dphi2**2)+(c3*np.cos(phi1-phi2)*np.sin(phi2))+(c4*np.cos(phi1-phi2)*(dphi2-dphi1))-(c5*np.sin(phi1))+c6*M-c7*dphi1)/(c8-(c9*np.cos(phi1-phi2)*np.cos(phi1-phi2))) 

我的第一个预测系数的第一个预测：

p = [0.5625, 0.375, 27.590625000000003, 0.09375, 55.18125, 62.5, 0.425, 1, 0.5625]

我迭代地计算函数的值。我占用所有变量的长度已经具有相同的大小：

n = len(time1)
y = np.empty(n)

for i in range(n):
    y[i] = ddphi1(dphi11[i],dphi22[i],phi11[i],phi22[i],M[i],p[0],p[1],p[2],p[3],p[4],p[5],p[6],p[7],p[8])

plt.plot(time1, ddphi11)
plt.plot(time1, y, 'r')

现在想法是用curve_fit自动计算系数：** ddphi1 ist我的回调函数和ddphi11我的形状数据（1000）以及其他变量

from scipy.optimize import curve_fit

g = [0.56, 0.37, 27.63, 0.094, 55.18, 62.5, 0.625, 1, 0.56]

c,cov =curve_fit(ddphi1,(dphi1,dphi2,phi1,phi2,M),ddphi11,g)
print(c)

，我收到了这个错误，

--------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-158-e8e42e7b1216> in <module>()
      1 from scipy.optimize import curve_fit
      2 
----> 3 c,cov =curve_fit(ddphi1,(dphi1,dphi2,phi1,phi2,M),ddphi11,g=[0.56, 0.37, 27.63, 0.094, 55.18, 62.5, 0.625, 1, 0.56])
      4 print(c)

1 frames
/usr/local/lib/python3.7/dist-packages/scipy/optimize/minpack.py in curve_fit(f, xdata, ydata, p0, sigma, absolute_sigma, check_finite, bounds, method, jac, **kwargs)
    719         # non-array_like `xdata`.
    720         if check_finite:
--> 721             xdata = np.asarray_chkfinite(xdata, float)
    722         else:
    723             xdata = np.asarray(xdata, float)

/usr/local/lib/python3.7/dist-packages/numpy/lib/function_base.py in asarray_chkfinite(a, dtype, order)
    484 
    485     """
--> 486     a = asarray(a, dtype=dtype, order=order)
    487     if a.dtype.char in typecodes['AllFloat'] and not np.isfinite(a).all():
    488         raise ValueError(

ValueError: setting an array element with a sequence. The requested array has an inhomogeneous shape after 1 dimensions. The detected shape was (5,) + inhomogeneous part.

我已经看到大多数数据，进入curve_fir在表单列表中。也许在处理数组时有一个解决方案？ Espero我希望您能为我提供帮助，因为我是Python的新手。

Answer 1

我终于能够解决它。只有阵列应在全局变量

x = np.column_stack（[DPHI11，DPHI22，PHI11，PHI22，M]）

中描述基于全局变量的模型

def model(X,c1,c2,c3,c4,c5,c6,c7,c8,c9):
  dphi1 = X[:,0]
  dphi2 = X[:,1]
  phi1 = X[:,2]
  phi2 = X[:,3]
  M = X[:,4]

  f = (-(c1*np.sin(phi1-phi2)*np.cos(phi1-phi2)*dphi1**2)-(c2*np.sin(phi1-phi2)*dphi2**2)+(c3*np.cos(phi1-phi2)*np.sin(phi2))+(c4*np.cos(phi1-phi2)*(dphi2-dphi1))-(c5*np.sin(phi1))+c6*M-c7*dphi1)/(c8-(c9*np.cos(phi1-phi2)*np.cos(phi1-phi2)))
  return f

，魔术开始

guesses = [0.56, 0.37, 27.63, 0.094, 55.18, 62.5, 0.625, 1, 0.56]

from scipy.optimize import curve_fit

popt, pcov = curve_fit(model, X, ddphi11, guesses)
print(popt)