<div dir="ltr">good point!</div><br><div class="gmail_quote"><div dir="ltr">On Fri, Aug 11, 2017 at 4:55 PM Nathaniel Smith <<a href="mailto:njs@pobox.com">njs@pobox.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On Wed, Aug 9, 2017 at 10:19 PM, Phillip Feldman<br>
<<a href="mailto:phillip.m.feldman@gmail.com" target="_blank">phillip.m.feldman@gmail.com</a>> wrote:<br>
> Firstly, white noise means only that the power spectral density is flat, or<br>
> equivalently, that the autocorrelation function is zero everywhere except at<br>
> lag zero.<br>
><br>
> One can generate Gaussian noise with an arbitrary power spectral density by<br>
> filtering white Gaussian noise with a suitable filter. The output of the<br>
> filter is the convolution of the impulse response of the filter with the<br>
> autocorrelation function of white noise, which gives us the impulse response<br>
> of the filter. So, from the convolution theorem and the definition of the<br>
> PSD, the PSD of the output is the squared magnitude of the frequency<br>
> response of filter.<br>
<br>
It's not actually possible to generate true 1/f noise this way --<br>
technically 1/f noise is non-stationary and doesn't have a PSD. (You<br>
can run a PSD estimator on any finite sample of 1/f noise, and get<br>
some answer, but as your samples get larger your estimate won't<br>
converge, because you keep discovering more and more power at lower<br>
and lower frequencies.) So there are specialized methods for<br>
generating 1/f noise, involving things like fractional differencing or<br>
wavelets.<br>
<br>
-n<br>
<br>
--<br>
Nathaniel J. Smith -- <a href="https://vorpus.org" rel="noreferrer" target="_blank">https://vorpus.org</a><br>
_______________________________________________<br>
SciPy-Dev mailing list<br>
<a href="mailto:SciPy-Dev@python.org" target="_blank">SciPy-Dev@python.org</a><br>
<a href="https://mail.python.org/mailman/listinfo/scipy-dev" rel="noreferrer" target="_blank">https://mail.python.org/mailman/listinfo/scipy-dev</a><br>
</blockquote></div>