def getPowerSpectralDensity(X,fs=1.0): assert fs > 0 f, Pxx_den = signal.periodogram(X, fs,scaling='density',window=None,detrend=False) return (f,Pxx_den)