What does the frequency axis of a Power Spectral Density mean?

Question

I have never really understood what the frequency axis meant when we plot the Power Spectral Density(PSD).

Does it correspond to frequency as we get after we take the Fourier Transform of a time domain signal. This doesn't make sense to me because the frequency here(when we do FT of a time domain signal) corresponds to the time variable "t".

To find the PSD we take the Fourier transform of the Autocorrelation function which is a function of either time difference $$\tau = t_2 - t_1$$ i.e. $$R_{xx}(\tau) $$ (for the Wide Sense Stationary case) or $$R_{xx}(t_1,t_2)$$ for other cases.

Here the autocorrelation is a function of either time difference or two time variables. If this is the case then how can we say that the PSD gives the power distributed in the frequency domain?

How do I understand this. If someone could give both a mathematical explanation and an intuitive explanation I would be really grateful.

Thank you.

Answer 1

The fact that the frequency variable of a power spectral density (PDS) equals the one of a Fourier transform of a "normal" time-domain signal can be seen more easily by considering the following definition of the PDS of a wide-sense stationary (WSS) random process $x(t)$:

$$S_x(\omega)=\lim_{T\rightarrow\infty}E\left\{ \frac{1}{T}\left| \int_{-T/2}^{T/2}x(t)e^{-j\omega t}dt \right|^2 \right\}\tag{1}$$

Note that $(1)$ is just the limit of the expectation of the scaled squared magnitude of the Fourier transform of a truncated version of $x(t)$. So the frequency variable $\omega$ is just the frequency variable of the (squared magnitude of the) Fourier transform.

The Wiener-Khinchin theorem says that the autocorrelation function and the PDS of a WSS random process form a Fourier transform pair.

This answer shows that the inverse Fourier transform of $(1)$ indeed equals the autocorrelation function of $x(t)$.