colored_noise (pyleoclim.utils.tsmodel.colored_noise)

pyleoclim.utils.tsmodel.colored_noise(alpha, t, f0=None, m=None, seed=None)[source]

Generate a colored noise timeseries

Parameters

  • alpha (float) – exponent of the 1/f^alpha noise

  • t (float) – time vector of the generated noise

  • f0 (float) – fundamental frequency

  • m (int) – maximum number of the waves, which determines the highest frequency of the components in the synthetic noise

Returns

y – the generated 1/f^alpha noise

Return type

array

References

Eq. (15) in Kirchner, J. W. Aliasing in 1/f(alpha) noise spectra: origins, consequences, and remedies.

Phys Rev E Stat Nonlin Soft Matter Phys 71, 066110 (2005).