Evaluating the spectral density of generated noise through the autocovariance

I'm working on generating noise signals $X(t)$ (with $t \in \left[0,T\right]$ with step size $\delta t$) with a prescribed power spectral density $S_{XX}(f)$ and I'm figuring out how well I am generating said noise; I need to know if it has the desired spectral density (in case someone notices, I posted about similar issues a while back. I've since then worked on a different topic and I am now returning to this).

However, estimating the spectral density of a signal is a tedious process. It involves a lot of estimation through periodograms and window functions and such, and it is therefore not so easy to quantify numerically if the spectrum not only has the desired trends but actually implements the desired parameters (in practice I will simply use a spectrum analyzer, but that is a different story).

The auto covariance function $C_{XX}(t') = \left<X(t) X(t+t') \right>$ on the other hand is very easy to calculate and subsequently fit. Moreover, the auto covariance function and the spectral density are related by a Fourier transform: $C_{XX}(t') = \int_0^\infty{S_{XX}(f)\cos{2\pi f}}$ (note that I'm using the cosine transform as this is customary in these settings, due to the evenness of the spectral density).

So my question is, how much do I learn about the spectral density by fitting the auto covariance? Say I am working with the Ornstein Uhlenbeck process, which has auto covariance $C_{XX}(t') = \frac{c \tau}{2}e^{-t'/\tau}$ and spectral density $S_{XX}(f) = \frac{2 c \tau^2}{1+(2\pi f \tau)^2}$. If I fit my numerically generated $X(t)$ to $C_{XX}$ to extract $c$ and $\tau$ and am able to do so with high accuracy, do I then know that my spectral density is also as desired?

Intuitively I would say yes, Fourier pairs contain the exact same information. The reason I am not sure if this is the case is due to the fact that I am unsure of how sampling rates and such come into this. Of course one cannot expect a spectral density to be accurate up to the GHz range if one has an $X(t)$ point every millisecond. But is this equally true for the auto covariance? In other words, if my autocovariance is accurate up to some specific lag time (which means I have a high enough sampling rate I suppose), is the spectral density then also accurate up to some specific frequency? And can I quantify this in some way?

What I mean by that last sentence is, can I tell up to (or down to) what frequencies I can trust my spectral density, given my auto covariance fit, $T$ and $\delta t$?

• I know this question is old, but I wonder what you concluded regarding this issue. I'm in the same situation, I'm generating some noise and I want to test that it satisfies the prescriptions. – Pincopallino Aug 31 '18 at 14:09
• In the end I ended up using periodograms, I believe with Welch's method, to directly look at the power spectral density instead of the autocovariance. – user129412 Sep 7 '18 at 7:51