Skip to content

astrogewgaw/PowerLawNoise.jl

Repository files navigation

PowerLawNoise

Generate power law noise, in Julia.


License GitHub Stars Gitmoji Badge


This package implements a novel algorithm, following the work in "On generating power law noise", Timmer, J. and Koenig, M., Astron. Astrophys. 300, 707-710 (1995), which allows generating power law noise with arbitrary dimensions and arbitrary exponents. An exponent of Ξ² = 0 , 1 and 2 corresponds to white, pink and red noise. The figure below demonstrates the difference between them, by plotting the resultant time series, generated via PowerLawNoise.jl:

Plot: 1D Noise

Here is another one, but this time we are simulating noise in two dimensions:

Plot: 2D Noise

We can also go ahead and verify that the algorithm works by performing a linear fit to the logarithmic plot of the power spectrum of the noisy time series generated via PowerLawNoise.jl. Here, we do this for the case of red noise (that is, Ξ² = 2.0 ):

Plot: Fitting Noise

Install it by typing and running:

] add PowerLawNoise

in the Julia REPL.

About

Generate power law noise, in Julia.

Resources

License

Stars

Watchers

Forks

Packages

No packages published

Languages