Computando um espectro de potência

Gostaria de calcular um espectro de potência usando Python3. De outro tópico sobre esse tópico, obtive os ingredientes básicos. Eu acho que deveria ser algo como: ps = np.abs(np.fft.fft(x))**2 timeres = t[1]-t[0] freqs = ...