# Spectrum Sensing in Cognitive Radios Based on Multiple Cyclic Frequencies

Invited Paper

###### Abstract

Cognitive radios sense the radio spectrum in order to find unused frequency bands and use them in an agile manner. Transmission by the primary user must be detected reliably even in the low signal-to-noise ratio (SNR) regime and in the face of shadowing and fading. Communication signals are typically cyclostationary, and have many periodic statistical properties related to the symbol rate, the coding and modulation schemes as well as the guard periods, for example. These properties can be exploited in designing a detector, and for distinguishing between the primary and secondary users’ signals. In this paper, a generalized likelihood ratio test (GLRT) for detecting the presence of cyclostationarity using multiple cyclic frequencies is proposed. Distributed decision making is employed by combining the quantized local test statistics from many secondary users. User cooperation allows for mitigating the effects of shadowing and provides a larger footprint for the cognitive radio system. Simulation examples demonstrate the resulting performance gains in the low SNR regime and the benefits of cooperative detection.

## I Introduction

Spectrum sensing is needed in cognitive radios in order to find opportunities for agile use of spectrum. Moreover, it is crucial for managing the level of interference caused to primary users (PUs) of the spectrum. Sensing provides awareness of the radio operating environment. A cognitive radio may then adapt its parameters such as carrier frequency, power and waveforms dynamically in order to provide the best available connection and to meet the user’s needs within the constraints on interference.

In wireless communication systems we typically have some knowledge on the waveforms and structural or statistical properties of the signals that the primary user of the spectrum is using. Such knowledge may be related to the modulation scheme, the symbol or chip rate of the signal, the channel coding scheme, training or pilot signals, guard periods, and the power level or correlation properties of the signal, just to mention a few. These properties may be used to design a detector that works in a very low SNR regime and has low complexity and consequently low power consumption. These are very desirable properties especially for cognitive radios in mobile applications. In the absence of any knowledge of the signal, one may have to resort to classical techniques such as energy detection [1]. An energy detector may need to collect data over a long period of time to detect the primary users reliably. Moreover, controlling the false alarm rates in mobile applications is difficult because the statistics of the signals, noise and interference may be time-varying. Another significant drawback is that energy detection has no capability to distinguish among different types of transmissions or to dichotomize between primary and secondary users of the spectrum.

Cyclostationary processes are random processes for which the statistical properties such as the mean and autocorrelation change periodically as functions of time [2]. Many of the signals used in wireless communication and radar systems possess this property. Cyclostationarity may be caused by modulation or coding [2], or it may be also intentionally produced in order to aid channel estimation or synchronization [3]. Cyclostationarity property has been widely used in intercept receivers [2, 4, 5], direction of arrival or time-delay estimation, blind equalization and channel estimation [6] as well as in precoder design in multicarrier communications [3]. In order to exploit cyclic statistics, the signal must be oversampled with respect to the symbol rate, or multiple receivers must be used to observe the signal. The use of cyclostationary statistics is appealing in many ways: noise is rarely cyclostationary and second-order cyclostationary statistics retain also the phase information. Hence, procedures based on cyclostationarity tend to have particularly good performance at the low SNR regime. Moreover, cyclostationarity allows for distinguishing among different transmission types and users if their signals have distinct cyclic frequencies. A comprehensive list of references on cyclostationarity along with a survey of the literature is presented in [7].

The presence of cyclostationary signals may be determined by using hypothesis testing. Many existing tests, such as [8], are able to detect the presence of cyclostationarity at only one cyclic frequency at a time, and they partly ignore the rich information present in the signals. For example, a communication signal may have cyclic frequencies related to the carrier frequency, the symbol rate and its harmonics, the chip rate, guard period, the scrambling code period, and the channel coding scheme. In this paper we propose a method for detecting multiple cyclic frequencies simultaneously. It extends the method of [8] to take into account the rich information present at different cyclic frequencies. This provides improved detector performance over techniques relying only on single cyclic frequency and facilitates dichotomizing among the primary and secondary user signals and different waveforms used.

In cognitive radio systems, there are typically multiple geographically distributed secondary users (SUs) that need to detect if the primary user is transmitting. The distributed sensors may work collaboratively to decide between two hypotheses: is the primary user active, or is the spectrum unused and available for the secondary users? Decentralized processing has a number of advantages for such situations. Obviously, it allows for a larger coverage area. Furthermore, there are gains similar to diversity gains in wireless communications so that the detection becomes less sensitive to demanding propagation conditions such as shadowing by large obstacles, large numbers of scatterers, differences in attenuation, or fast fading caused by mobility. Moreover, distributed sensory systems may require less communication bandwidth, consume less power, be more reliable and cost less as well. In this paper, we propose a simple decentralized decision making approach based on sharing and combining quantized local decision statistics. This approach may be used in both decision making with or without a fusion center.

This paper is organized as follows. In Section II, there is a short review of cyclostationary statistics. A novel detector for multiple cyclic frequencies is derived in Section III. Section IV addresses the problem of collaborative detection of primary user. Simulation results demonstrating the detector’s reliability in the low SNR regime as well as the gains obtained via collaborative operation are presented in Section V. Finally, conclusions are drawn in Section VI.

## Ii Cyclostationarity: a recap

In this section, we provide a brief overview of cyclostationarity in order to make the derivation of the detector in Section III clearer. A continuous-time random process is wide sense second-order cyclostationary if there exists a such that [2]:

(1) |

and

(2) |

is called the period of the cyclostationary process.

Due to the periodicity of the autocorrelation , it has a Fourier-series representation. By denoting and , we obtain the following expression for the Fourier-series [2]:

(3) |

where the Fourier coefficients are

(4) |

and is called the cyclic frequency. The function is called the cyclic autocorrelation function. If the process has zero mean, then this is also the cyclic autocovariance function.

When the autocorrelation function has exactly one period we have the following set of cyclic frequencies

where is the cyclic autocorrelation function and are the set of cyclic frequencies. The cyclic frequencies are harmonics of the fundamental frequency. If the autocorrelation function has several periods , we may express at the limit [2]

(5) |

The process is almost cyclostationary in the wide sense and the set of cyclic frequencies is comprised of a countable number of frequencies that do not need to be harmonics of the fundamental frequency. In general, the process is said to be cyclostationary if there exists an such that for some value of . Typically cyclic frequencies are assumed to be known or may be estimated reliably.

## Iii Detection using multiple cyclic frequencies

Statistical tests for the presence of a single cyclic frequency have been proposed, for example, in [8]. The tests in [8] have asymptotically constant false alarm rate (CFAR) for testing presence of cyclostationarity at a given cyclic frequency. However, the tests do not retain the CFAR property over a set of tested frequencies.

Typical communication signals exhibit cyclostationarity at multiple cyclic frequencies instead of just a single cyclic frequency. That is, for example a signal that is cyclostationary at the symbol frequency is typically cyclostationary at all integer multiples of the symbol frequency as well. There also may be cyclic frequencies related to the coding and guard periods, or adaptive modulation and coding may be used. In such cases the cyclic frequencies present may vary depending on channel quality and the waveform used. If one is testing for the presence of many different signals at a given frequency band, or in case the cyclic frequencies are not known, it would be desirable to retain the CFAR property over the whole set of tested cyclic frequencies. This would be especially desirable in a cognitive radio application where the interest is in finding unoccupied frequency bands. Otherwise the frequency band may unnecessarily be classified as occupied for most of the time.

In the following we extend the test based on second-order cyclic statistics of [8] to multiple cyclic frequencies. To do so we first define all the terms used in the test statistics.

Let denote an optional complex conjugation. The notation allows convenient handling of both cyclic autocorrelation and conjugate cyclic autocorrelation with only one equation. An estimate of the (conjugate) cyclic autocorrelation may be obtained using observations as

(6) | ||||

(7) |

where the latter term is the estimation error. This estimator is consistent, (see [8]) so that the error goes to zero as .

Now we need to construct a test for a number of lags as well as a set of cyclic frequencies of interest. Let denote the set of cyclic frequencies of interest, and

(8) | ||||

denote a vector containing the real and imaginary parts of the estimated cyclic autocorrelations at the cyclic frequency of interest stacked in a single vector.

The covariance matrix of can be computed as [8]

(9) |

where the th entries of the two covariance matrices and are given by

and

(10) |

Here, and denote the unconjugated and conjugated cyclic spectra of , respectively. These spectra can be estimated using frequency smoothed cyclic periodograms as

(11) | ||||

(12) |

where and is a normalized spectral window of odd length .

Now the hypothesis testing problem for testing if is a cyclic frequency can be formulated as [8]

(13) | ||||

Here is the estimation error which is asymptotically normal distributed, i.e., [8]. Hence, using the asymptotic normality of the generalized likelihood ratio (GLR) is given by

(14) | ||||

Finally, by taking the logarithm and multiplying the result by 2, we arrive at the test statistic in [8]

(15) |

Under the null hypothesis is asymptotically distributed.

Now in order to extend the test for the presence of second-order cyclostationarity at any of the cyclic frequencies of interest simultaneously, we formulate the hypothesis testing as follows

(16) | ||||

For this detection problem, we propose the following two test statistics:

(17) | ||||

(18) |

The first test statistic calculates the maximum of the cyclostationary GLRT statistic (15) over the cyclic frequencies of interest while the second calculates the sum. Assuming independence of cyclic autocorrelation estimates for different cyclic frequencies the test statistic is the GLRT statistic. Depending on the signal and the set of tested cyclic frequencies the test statistics may have different performances. This requires further research.

The asymptotic distribution of is under the null hypothesis where is the number of cyclic frequencies in set . This is due to the fact that the sum of independent chi-square random variables is also a chi-square random variable whose degrees of freedom is the sum of the degrees of freedom of the independent random variables.

In the following we derive the asymptotic distribution of the test statistic under the null hypothesis. As stated above, under the null hypothesis is asymptotically distributed. The cumulative distribution function of the chi-square distribution with degrees of freedom is given by

(19) |

where is the lower incomplete gamma function and is the ordinary gamma function. For a positive integer the following identities hold:

(20) | ||||

(21) |

Hence, the cumulative distribution function of the chi-square distribution with degrees of freedom is given by

(22) |

The cumulative distribution function of the maximum of independent and identically distributed random variables is the cumulative distribution function of the individual random variables raised to the power . Thus, the cumulative distribution function of the test statistic is given by

(23) |

The corresponding probability density function is obtained by differentiating the cumulative distribution function, i.e.,

(24) | ||||

Consequently, the null hypothesis is rejected if where is the false alarm rate and is the number of tested cyclic frequencies.

## Iv Cooperative detection

User cooperation may be used to improve the performance and coverage in a cognitive radio network. The users may collaborate in finding unused spectrum and new opportunities. Many of the collaborative detection techniques stem from distributed detection theory; see [10, 11]. In cognitive radio systems, there are typically multiple geographically distributed secondary users that need to detect whether the primary user is active. All the secondary users may sense the entire band of interest, or monitor just a partial band to reduce power consumption. In the latter case each SU senses a certain part of the spectrum, and then shares the acquired information with other users or a fusion center.

The cooperation may then be coordinated by a fusion center (FC), or it may take place in an ad-hoc manner without a dedicated fusion center. Here we assume that a fusion center collects information from all secondary users and makes a decision about whether the spectrum is available or not. We assume that each secondary user sends a quantized version of its local decision statistics (such as the likelihood ratio) to the FC. In the case of very coarse quantization, binary local decision may be sent. To derive a test for the FC, we assume that the sensors are independent conditioned on whether the hypothesis or is true. Then the optimal fusion rule is the likelihood ratio test over the received local likelihood ratios :

(25) |

In case the secondary users send binary decisions, the sum of ones may calculated and compared to a threshold. Here, we consider the simplest way of making the decision using generalized likelihood ratios. Instead of using the product of the generalized likelihood ratios, we can employ the sum of generalized log-likelihood ratios. We propose the following test statistic for the hypothesis testing problem (13)

(26) |

and the following two for the hypothesis testing problem (16)

(27) | ||||

(28) |

where is the cyclostationarity based test statistic (15) from secondary user. Due to the use of generalized likelihood ratios, no optimality properties can be claimed. The GLRT test does, however, perform highly reliably in many applications.

Under the conditional independence assumption the asymptotic distributions of the test statistic and are under the null hypothesis and , respectively. This is again due to the fact that the sum of independent chi-square random variables is also a chi-square random variable whose degrees of freedom is the sum of the degrees of freedom of the independent random variables. The cumulative distribution function of is under the null hypothesis where is again the number of tested cyclic frequencies. The testing is done similarly as in one secondary user case.

Different techniques for reducing the amount of transmitted data, taking into account the relevance of the information provided by secondary users as well as how to deal with communication rate constraints will be addressed in a forthcoming paper.

## V Simulation examples

In this section the performance of the proposed detectors is considered. The test signal is an orthogonal frequency division multiplex (OFDM) signal. The baseband equivalent of a cyclic prefix OFDM signal may be expressed as

(29) |

where is the number of subcarriers, is the symbol length, denotes the rectangular pulse of length , and the ’s denote the data symbols. The symbol length is the sum of the length of the useful symbol data and the length of the cyclic prefix , i.e., .

The above OFDM signal exhibits cyclostationarity (i.e., complex conjugation is used in (6) and the following equations) with cyclic frequencies of , and potentially other frequencies depending on the coding scheme. The cyclic autocorrelation surfaces for peak at [9].

In the following the performance of cyclic detectors based on one and two cyclic frequencies is compared as a function of signal-to-noise ratio (SNR) in an additive white Gaussian noise (AWGN) channel. The SNR is defined as where and are the variances of the signal and the noise, respectively. The cyclic frequencies employed by the detectors are and . The detector based on one cyclic frequency uses the first frequency and the detectors based on two cyclic frequencies use both frequencies. Each detector uses two time lags .

The cyclic spectrum estimates were calculated using a length-2049 Kaiser window with parameter of 10. A Fast-Fourier transform (FFT) was employed for faster computation. The FFT size was 10000 giving a cyclic frequency resolution of 0.0001.

The OFDM signal has 32 subcarriers and the length of the cyclic prefix is 1/4 of the useful symbol data. The subcarrier modulation employed is 16-QAM. The signal length is 100 OFDM symbols.

Fig. 1 depicts the performance of the detectors as a function of the SNR for a constant false alarm rate of 0.05. Fig. 2 shows a zoom of the important area illustrating the differences in performance more clearly. All the curves are averages over 10000 experiments. It can be seen that the detectors based on multiple cyclic frequencies outperform the detector based on single cyclic frequency in the low SNR regime. Furthermore, the multicycle detector calculating the sum over the cyclic statistics of different frequencies has the best performance.

Fig. 3 plots the probability of detection vs. false alarm rate for SNR of -7 dB. The figure show that the detectors have desirable receiver operating characteristics. That is, the probability of detection increases as the false alarm rate parameter is increased.

Next the performance gain from cooperative detection of several secondary users is analyzed. The signal is the same as above. The cooperative detection is based on the data of 5 secondary users. Each secondary user receives the same data with different noise. SNR is the same for each secondary user.

Fig. 4 depicts the performance for 5 secondary users compared to the single secondary user case. Performance gain of roughly 3 dB is obtained from the cooperation of 5 secondary users. Using two cyclic frequencies provides similar performance improvement as in single secondary user case. Fig. 5 shows the probability of detection vs. false alarm rate for SNR of -9 dB.

In the following simplistic example, we illustrate the gains that may be achieved via collaborative detection in the face of shadowing effects. In order to simulate shadowing, the SNR of each user was independently selected randomly from a normal distribution with a mean of -9 dB and standard deviation of 10 dB. That is, the logarithm of the received power level is normally distributed. Fig. 6 depicts the performance of the multicycle detectors for the simple shadowing scenario. Comparison to Fig. 5 reveals that cooperation among secondary users reduces sensitivity to shadowing effects significantly.

## Vi Conclusion

In this paper, a generalized likelihood ratio test for detecting primary transmissions with multiple cyclic frequencies has been proposed, and the asymptotic distribution of the test statistic has been derived. In this test, impairments such as shadowing and fading are mitigated by combining the quantized local likelihood ratios from a number of secondary users under a conditional independence assumption. Simulation examples demonstrating the improved reliability in the detector performance in the low SNR regime as well as significant gains obtained via collaborative decision making have also been presented.

## References

- [1] H. V. Poor, An Introduction to Signal Detection and Estimation, 2nd edition, Springer, New York, 1994.
- [2] W. A. Gardner, Statistical Spectral Analysis: A Nonprobabilistic Theory, Prentice-Hall, Upper Saddle River, NJ, 1987.
- [3] M. Tsatsanis and G.B. Giannakis, “Transmitter Induced Cyclostationarity for Blind Channel Equalization,” IEEE Trans. Signal Processing, Vol. 45, pp. 1785-1794, Jul. 1997.
- [4] T. Koivisto and V. Koivunen, “Blind Despreading of Short-Code CDMA Signals in Asynchronous Multi-User Systems,” Signal Processing, to appear in 2007.
- [5] J. Lundén and V. Koivunen, “Automatic Radar Waveform Recognition,” IEEE Journal of Selected Topics in Signal Processing, Special issue on Adaptive Waveform Design for Agile Sensing and Communication, to appear in 2007.
- [6] L. Tong, G. Xu, and T. Kailath, “Blind Identification and Equalization Based on Second-Order Statistics: A Time Domain Approach,” IEEE Trans. Information Theory, vol. 40, no. 2, pp. 340-349, Mar. 1994.
- [7] W. A. Gardner, A. Napolitano, and L. Paura, “Cyclostationarity: Half a Century of Research,” Signal Processing, Vol. 86, pp. 639-697, Apr. 2006.
- [8] A. V. Dandawaté and G. B. Giannakis, “Statistical Tests for Presence of Cyclostationarity,” IEEE Trans. Signal Processing, vol. 42, no. 9, pp. 2355–2369, Sep. 1994.
- [9] M. Öner and F. Jondral, “Air Interface Recognition for a Software Radio System Exploiting Cyclostationarity,” in Proc. 15th IEEE International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC’04), Barcelona, Spain, Sep. 5-8, 2004, vol. 3, pp. 1947-1951.
- [10] R. Viswanathan, and P. K. Varshney, “Distributed Detection with Multiple Sensors: Part I – Fundamentals,” Proceedings of the IEEE, Vol. 85, No. 1, pp. 54-63, Jan. 1997.
- [11] R. S. Blum, S. A. Kassam, and H. V. Poor, “Distributed Detection with Multiple Sensors: Part II – Advanced Topics,” Proceedings of the IEEE, Vol. 85, No. 1, pp. 64-79, Jan. 1997.