Calcular un espectro de potencia

Me gustaría calcular un espectro de potencia usando Python3. De otro hilo sobre este tema obtuve los ingredientes básicos. Creo que debería ser algo como: ps = np.abs(np.fft.fft(x))**2 timeres = t[1]-t[0] freqs = np.fft.fftfreq(x.size, timeres) ...