Berechnung eines Leistungsspektrums

Ich möchte mit Python3 ein Leistungsspektrum berechnen. Von einem anderen Thread zu diesem Thema habe ich die Grundzutaten bekommen. Ich denke, es sollte so etwas wie sein: ps = np.abs(np.fft.fft(x))**2 timeres = t[1]-t[0] freqs = ...