import matplotlib.pyplot as plt
from numpy import log10
from sk_dsp_comm import sigsys as ss
x,b, data = ss.nrz_bits(10000,10)
Px,f = ss.my_psd(x,2**10,10)
plt.plot(f, 10*log10(Px))
plt.ylabel("Power Spectral Density (dB)")
plt.xlabel("Frequency (Hz)")
plt.show()