hzi-notebook
6/7/2019 - 4:26 PM

Periodogramm

from scipy import signal

f, Pxx_dens = signal.welch(df.dT, nperseg=60)
plt.plot(f, Pxx_dens)
plt.xlabel('Frequenz [1/min]')
plt.ylabel('power spectral density')
plt.grid()
plt.axvline(.067, c='r', alpha=0.3)
plt.axvline(.14 , c='r', alpha=0.3)
plt.show()