Seleccionamos algunos puntos en función de nuestro muestreo (opcional)
t = t[0:500]
y = y[0:500]
se realiza un ajuste
a * np.cos(w*t)+b*np.sin(w*t) + d
[0.6,0.6,7,3.6]