The Wiener–Khinchin theorem states that the power spectral density of a wide-sense stationary stochastic process can be obtained through the Fourier transform of the autocorrelation of the signal i.e. $$ S(\omega) = \int_{-\infty}^{\infty} R(\tau) e^{-j \omega \tau} d\tau $$
I would like to find out what happens when I use this theorem for non-WSS signals?
To be specific, I want to assess the validity, for non-WSS, of a derivation which uses this theorem to arrive to the conclusion that when two signals with phase noise are multiplied, the output phase noise $S_{\phi_{IF}}$ is related to the input phase noise $S_{\phi}$ by $$ S_{\phi_{IF}}(f) = S_{\phi}(f) 4 \sin^2 (\frac{\alpha f}{2}) $$
I am aware of the more generally applicable approach for computing PSD through Fourier decomposition but it's a little difficult to apply to the above mentioned derivation.
Thank you.