ich hab ein Zeitsignal was keine dieskreten abstände hab und möchte und möchte eine Frequenz analyse mit multi-taper methode durchführen (mit 4 tapern). Ist das in Python möglich. Also die Lücken können auch mal größer sein. Ich hab es bis jetzt immer mit lombscargle gemacht
Hier ein beispiel Signal:
Und hier die Lonbscargel Lösung:A = 2.
w = 1.
phi = 0.5 * np.pi
nin = 1000
nout = 100000
frac_points = 0.9 # Fraction of points to select
r = np.random.rand(nin)
x = np.linspace(0.01, 10*np.pi, nin)
x = x[r >= frac_points]
y = A * np.sin(w*x+phi)
Kann mir jemand weiter helfen?import scipy.signal as signal
pgram = signal.lombscargle(x, y, f, normalize=True)
plt.subplot(2, 1, 1)
plt.plot(x, y, 'b+')
plt.subplot(2, 1, 2)
plt.plot(f, pgram)
plt.show()