# Comments on “Blind Cyclostationary Spectrum Sensing in Cognitive Radios” by W. M. Jang

I recently came across the 2014 paper in the title of this post. I mentioned it briefly in the post on the periodogram. But I’m going to talk about it a bit more here because this is the kind of thing that makes things a bit harder for people trying to learn about cyclostationarity, which eventually leads to the need for something like the CSP Blog.

The idea behind the paper is that it would be nice to avoid the need for prior knowledge of cycle frequencies when using cycle detectors or the like. If you could just compute the entire spectral correlation function, then collapse it by integrating (summing) over frequency $f$, then you’d have a one-dimensional function of cycle frequency $\alpha$ and you could then process that function inexpensively to perform detection and classification tasks.

The author sets up an expression for one of the terms of the incoherent suboptimal multicycle detector, which has an idealized counterpart that is equal to the sum (integral) of the magnitude-squared spectral correlation function.  Something like this:

$\displaystyle M(\alpha) = \int_{-\infty}^\infty \left| S_x^\alpha(f) \right|^2\, df \hfill (1)$

The author then tries to simplify this kind of expression, and that is where things start to go very wrong. The root of it all is the author’s equation (3),

$\displaystyle S_x^\alpha(f) = \lim_{T\rightarrow\infty} \frac{1}{T} X_T (f+\alpha/2) X_T^*(f-\alpha/2), \hfill (2)$

where $X_T(f)$ is the Fourier transform of the given data $x(t)$ restricted to a time interval centered at $t=0$ and having width $T$. In other words, the author starts off by stating that the spectral correlation function is the limit of the cyclic periodogram as the amount of data increases without bound. Exactly what we show is not true in the periodogram post. The limit does not exist.

The (idealized) multi-cycle detection statistic of interest to the author is

$\displaystyle Y_{ML} = \sum_{\alpha\neq 0} M(\alpha) \hfill (3)$

This can be re-expressed as

$\displaystyle Y_{ML} = \sum_{\alpha\neq 0} \int_{-\infty}^\infty S_x^\alpha(f) S_x^\alpha(f)^* \, df \hfill (4)$

The author then tries to convert the sum over cycle frequency $\alpha$ into an integral, but does not use impulse functions to do so, ending up with

$\displaystyle Y_{ML} = Y^\alpha - Y^0 \hfill (5)$

$\displaystyle Y^\alpha = \int_{-\infty}^\infty \int_{-\infty}^\infty \left| S_x^\alpha(f) \right|^2 \, df \, d\alpha \hfill (6)$

$\displaystyle Y^0 = \int_{-\infty}^\infty \left| S_x^0(f) \right|^2 \, df \hfill (7)$

But we know that the spectral correlation function $S_x^\alpha(f)$ is continuous in spectral frequency $f$ and discrete in cycle frequency $\alpha$. In many cases of practical interest, there are only a finite number of cycle frequencies for which the spectral correlation function is not identically zero (think of square-root raised-cosine PSK, QAM, PAM), so that in (6) the integral over cycle frequency is zero. In other cases, the number of cycle frequencies is infinite, but countable (think rectangular-pulse BPSK), which still leads to a zero-valued integral over cycle frequency. In general, the spectral correlation function is non-zero only on a set of measure zero in the bi-frequency $(f, \alpha)$ plane, and so such integrals as the author writes in (6) above are always zero. So (6) is a mistaken representation of the desired quantity.

A double-integral representation of $Y_{ML}$ is possible by using impulse functions (Dirac delta functions).

Then the author tries to substitute (2) above into (4), which is a mistake because (2) is not true, but he also does it incorrectly. In particular, he writes

$\displaystyle Y^\alpha = \lim_{T\rightarrow\infty} \frac{1}{T} \int_{-\infty}^\infty \int_{-\infty}^\infty X_T(f+\alpha/2) X_T^*(f-\alpha/2) X_T^*(f+\alpha/2)X_T(f-\alpha/2) \, df \, d\alpha, \hfill (8)$

which he then simplifies to

$\displaystyle Y^\alpha = \lim_{T\rightarrow\infty} \frac{1}{T} \int_{-\infty}^\infty \int_{-\infty}^\infty X_T(f+\alpha) X_T(f) X_T^*(f+\alpha)X_T^*(f) \, df \, d\alpha, \hfill (9)$

by using a change of variables for $\alpha$. But this is an incorrect substitution. The author’s (3) (my (2) above) must be substituted for each $S_x^\alpha(f)$ appearing in the expression, leading to two limits (over two time variables $T_1$ and $T_2$). In other words, he should have ended up with a double limit.

In reality, each spectral correlation function $S_x^\alpha(f)$ is itself correctly represented by a double limit involving the data observation interval length $T$ and a smoothing or averaging parameter. In the case of the frequency-smoothing method, it is the width of the smoothing function $g_\Delta(f)$. So, we would end up with four limits in (9) if we did things properly.

But the problems continue. The author then defines $Y_T(f)$ as the (scaled) periodogram,

$\displaystyle Y_T(f) = X_T(f) X_T^*(f) \hfill (10)$

He ends up with

$\displaystyle Y^\alpha = \lim_{T\rightarrow\infty} \frac{1}{T} \int_{-\infty}^\infty Y_T(\alpha) \otimes Y_T(-\alpha)^* \, d\alpha \hfill (11)$

(the conjugation is superfluous since $Y_T(f)$ is non-negative).

So the quantity $Y^\alpha$, which was defined as the integral of the squares of all non-zero spectral correlation functions, is now reduced (incorrectly!) to the convolution of the periodogram with itself. All the cyclostationarity of the signal–no matter how rich or how poor–can be extracted by computing only the periodogram!

There are lots of other problems with the paper, including in the simulations section, but they don’t really matter because the theoretical development is seriously mathematically flawed. And most of the problem arises from trying to take the limit of the cyclic periodogram.