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.

Comments and corrections are welcome!

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

  1. Mirko von Leipzig says:

    I remember being super excited when I found this paper a year or so ago. I couldn’t follow most of the assumptions/maths behind it, but I attributed that to my poor understanding of cyclostationarity. Glad to see it it wasn’t entirely my fault.

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )


Connecting to %s