New answers tagged


Here the solution, a function to create a randomized time-series starting from a PSD: def gen_pkfield(nsample=100, alpha=-11.0 / 3, fknee=1, fsample=2.5): # sampling per block ufreq = np.fft.fftfreq(n=nsample+1, d=1./fsample)[1:] # the Discrete Fourier Transform sample frequencies (with negative values), [1:] to solve the DC kfreq = np.sqrt(...

Top 50 recent answers are included