from spectrum import *
data = data_cosine(N=1024, A=0.1, sampling=1024, freq=200)
p = Periodogram(data, sampling=1024) #here you just define the PSD estimate
p.plot(marker='o', sides='twosided')