Sign up ×
MathOverflow is a question and answer site for professional mathematicians. It's 100% free, no registration required.

I want to numerically study the behavior of a system described by a set of differential equations in the presence of colored noise. It seems that the standard procedure is to use the Langevin equation: $ y=fdt+gdX; dX=-aXdt+bdW$, where $W$ is the Wiener process.

The spectral density function of $X$ is $1/f^2$ type. What about the spectral density function of $dX$? If my MATLAB calculations are correct, in this case the power density increases with the frequency. If this is true, how can I introduce $1/f$ noise in my system?

share|cite|improve this question

1 Answer 1

You cannot create $1/f$ noise as solutions to an SDE. The best you can do is to take a vector $\{\xi_k\}_{k=-N}^N$ of i.i.d. random variables, divide $\xi_k$ by $\sqrt k$, and take the Fourier transform of that. This will be your $1/f$ noise $X$...

share|cite|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.