def __run(): N = 3328 d = test_data() p = Periodogram(d, 10000., N, window='HFT116D', overlap=0.5) print("FFT is " + _fft_name) print("Params: " + repr(p.get_parameters())) #plot(f, sqrt(ps_to_psd(ps, params_dict))) figure() title("linear spectrum") ylabel("Vrms") plot(p.f, p.ls()) figure() title("linear spectral density") ylabel("Vrms/sqrt(Hz)") plot(p.f, p.lsd()) yscale('log') show()