我试图从表中给定x和y值来模拟一个特定的trig函数。我应该找到的函数参数的形式是$A\cdot sin（x\cdot t + ph）+d$。其中$A$是振幅，$x $是周期，$ph$是相移，$d$是垂直位移。这里$t$是变量。

from pylab import * 
 from scipy.optimize import curve_fit
 time = [1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23]
 waterlevel = [125,199,223,176,101,73,109,181,219,188,122,85]

 def h(t,x, A , phi, d):
      return A*sin(x*t + phi) + d

 x, A, phi, d = curve_fit(h, time ,waterlevel, p0 = [0,71,-57,150]) # p0 are my approximations

 print("A =", round(A,2))
 print("phi =", round(phi,2))
 print("d =", round(d,2))
 print("x =", round(x,2))

 plot(tid, vannstand, "o")
 xlabel("Time (hours)")
 ylabel("water level (cm)")
 t = linspace(0,24,1000)
 plot(t, h(t,x,A,phi,d), "r")
 show()

这里的技巧是确定$x$，相移，没有最初的猜测。我已经尝试了上述代码，但我得到一个消息，没有足够的值在第七行解包（预期4，得到2）。任何想法，我的错误和潜在的步骤，以确定值？