colored_noise_2regimes (pyleoclim.utils.tsmodel.colored_noise_2regimes)

pyleoclim.utils.tsmodel.colored_noise_2regimes(alpha1, alpha2, f_break, t, f0=None, m=None, seed=None)

Generate a colored noise timeseries with two regimes

Parameters

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

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

• f_break (float) – the frequency where the scaling breaks

• 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).