In this post I provide plots of the spectral correlation for a variety of simulated textbook signals and several captured communication signals. The plots show the variety of cycle-frequency patterns that arise from the disparate approaches to digital communication signaling. The distinguishability of these patterns, combined with the inability to distinguish based on the power spectrum, leads to a powerful set of classification (modulation recognition) features (My Papers [16, 25, 26, 28]).
In all cases, the cycle frequencies are blindly estimated by the strip spectral correlation analyzer (The Literature [R3, R4]) and the estimates used by the FSM to compute the spectral correlation function. MATLAB is then used to plot the magnitude of the spectral correlation and conjugate spectral correlation, as specified by the determined non-conjugate and conjugate cycle frequencies.
There are three categories of signal types in this gallery: textbook signals, captured signals, and feature-rich signals. The latter comprises some captured signals (e.g., LTE) and some simulated radar signals. For the first two signal categories, the three-dimensional surface plots I’ve been using will suffice for illustrating the cycle-frequency patterns and the behavior of the spectral correlation function over frequency. But for the last category, the number of cycle frequencies is so large that the three-dimensional surface is difficult to interpret–it is a visual mess. For these signals, I’ll plot the maximum spectral correlation magnitude over spectral frequency versus the detected cycle frequency (as in this post).
Textbook Signals (Simulated)
For the textbook signals the sampling rate is set to . Multiply the and axes by any sampling rate of interest to see the actual spectral widths of the features. For those signals with a symbol rate, it is set to , and for all signals the carrier frequency is set to . The FSM is used to estimate the spectral correlation function for each of the obtained cycle frequencies. The number of samples (block length) used by the FSM is and the width of the frequency smoothing window is (% of the sampling rate).
A key aspect of the textbook gallery is that some distinct signal types give rise to identical spectral correlation functions. In particular, the spectral correlation functions for rectangular-pulse QPSK, 8PSK, 16QAM, and -DQPSK are identical. Similarly, the spectral correlation functions for square-root raised-cosine (SRRC) QPSK, 8PSK, 16QAM, and -DQPSK are identical, because they all use the same roll-off factor of . The higher-order statistics for these signals are, fortunately, not identical.
The measured spectral correlation functions for four captured-over-the-air communication signals are shown here:
The cycle frequencies for WCDMA (downlink) include the chip rate MHz as well as kHz and kHz. For CDMA EVDO (downlink), we have the CDMA chip rate of MHz as well as several harmonics of MHz. For CDMA (downlink) we have the chip rate MHz, and for ATSC Digital TV, we have the two conjugate cycle frequencies equal to . See My Papers  for more details on the cyclostationarity of these and other real-world signals. Notice that the PSDs for the four captured signals show the effects of a multipath propagation channel between the transmitter and my lab. The multipath channel distorts both the spectrum and the spectral correlation function, but does not modify the values of the cycle frequencies.
Finally, the following sequence of plots shows the spectral correlation for a few signal types possessing a great many cycle frequencies. The usual way of plotting the spectral correlation function is not helpful here as our eyes cannot untangle the many individual features in the hidden-line display.
The two radar signals are simulated. The LFM radar signal is simply a train of identical linear FM pulses (chirps), whereas the phase-coded radar signal is a train of identical PSK-like signals employing a binary minimum peak sidelobe sequence. The pulse repetition interval for both is samples (sampling rate is for the radar signals). The duty cycles are about . These signals are periodic, and so are in a sense trivially second-order cyclostationary signals. This is because they consist of a sum of finite-strength sine-wave components, and so naturally there are finite-strength sine-wave components in second- and higher-order nonlinear transformations of the signal.
The LTE signals produce nearly as many cycle frequencies as the radar signals, but their cyclostationarity is not trivial. That is, the LTE signals are not periodic. But they are fiendishly complex, and this complexity leads to a complex cyclostationarity. Finally, the GSM signal is also quite complex relative to its base modulation of GMSK, so we see the appearance of many cycle frequencies separated by a harmonic of the frame rate, which is Hz.
The signal files I used to create these estimates are available on my Downloads page.