Advanced
A Comparison of Spectrum-Sensing Algorithms Based on Eigenvalues
A Comparison of Spectrum-Sensing Algorithms Based on Eigenvalues
Journal of information and communication convergence engineering. 2015. Dec, 13(4): 241-247
Copyright © 2015, The Korean Institute of Information and Commucation Engineering
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/li­censes/by-nc/3.0/)which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : November 05, 2015
  • Accepted : December 01, 2015
  • Published : December 31, 2015
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Syed Sajjad Ali
Jialong Liu
Chang Liu
Minglu Jin
mljin@dlut.edu.cn

Abstract
Cognitive radio has been attracting increased attention as an effective approach to improving spectrum efficiency. One component of cognitive radio, spectrum sensing, has an important relationship with the performance of cognitive radio. In this paper, after a summary and analysis of the existing spectrum-sensing algorithms, we report that the existing eigenvalue-based semi-blind detection algorithm and blind detection algorithm have not made full use of the eigenvalues of the received signals. Applying multi-antenna systems to cognitive users, we design a variety of spectrum-sensing algorithms based on the joint distribution of the eigenvalues of the received signal. Simulation results validate that the proposed algorithms in this paper are able to detect whether the signal of the primary user exists or not with high probability of detection in an environment with a low signal-to-noise ratio. Compared with traditional algorithms, the new algorithms have the advantages of high detection performance and strong robustness
Keywords
I. INTRODUCTION
Recently, the issue of the limited supply of radio spectrum resources has become a growing concern due to the enormous demand for wireless communication services. At the same time, the current system of fixed allocation of the spectrum leads to low usage for most of the spectrum, according to the research of the US Federal Communi-cations Commission (FCC) [1] . In the face of this so-called “spectral crisis”, Mitola [2] proposed the concept of ‘cognitive radio (CR)’.
Wireless communication technology has entered the era of 4G/5G. The demands for higher data transmission speed through wireless communication systems increase with each passing day. How to improve the reliability and bandwidth efficiency of the system has become the vital target of wireless communication technology design in the next generation and beyond. Transmitting data streams using multiple antennas, known as multiple-input/multiple-output (MIMO), can increase the capacity and bandwidth efficiency of communication systems significantly without greatly increasing the system bandwidth at the same time. This is considered one of the key technologies of modern wireless communication [1] . However, MIMO has some drawbacks [2] : high inter-antenna synchronization is required among transmitting antennas to meet the demand of simultaneous data transmissions; simultaneous data transmissions by multiple antennas can cause high inter-channel interference (ICI), which increases the difficulty of decoding, as well as the complexity of the system; and multiple RF chains are needed when multiple antennas work at the same time, increasing system costs.
A CR system is an intelligent radio communication system that is able to sense the surrounding radio envi-ronment and analyze the sensed result. Then the radio operating parameters will be dynamically adjusted in real time on the basis of that result [3] . The drawbacks of the present spectrum allocation scheme can be resolved effect-tively by the secondary utilization of spectrum resources. According to the IEEE 802.22.1 standard [4] , the secondary user (SU) should detect TV and broadcast signal and idle channels within 2 seconds, which requires cognitive users to access and exit the frequency band authorized by the primary user (PU) promptly, on the premise of guaranteed high detection probability.
The existing spectrum sensing algorithms can be roughly classified into three categories:
  • 1. Regular detecting algorithms requiring power infor-mation on signals and noise, such as the likelihood ratio test[5]and matched filtering detection[5,6]. This kind of algorithm has the best performance of the known algorithms. However, in a real-time situation, these two types of needed information usually cannot be acquired precisely.
  • 2. Half-blind detecting algorithms requiring power infor-mation on noise only, such as energy detection[7], which has the merits of reduced calculation and high detection performance, whereas it can be easily influenced by noise uncertainty, and its sensing performance is closely related to the threshold as well.
  • 3. Blind algorithms with no need for any additional information, such as the maximum-minimum eigenvalue algorithm[8]. No prior information is needed in this algorithm, yet it adopts the limit value of the maximum-minimum eigenvalue to calculate the threshold beyond which it cannot follow changes in the environment, which degrades detection performance.
On account of the problems with the existing algorithms, the concept of using the generalized likelihood ratio has been raised [9] . In a multi-antenna system, when there is only one PU in the space, the structural properties of the received signal can be utilized to design an algorithm. In such a case, the application environment is more realistic because only a small number of sample points are involved. This algorithm has a higher detection performance under a low signal-to-noise ratio (SNR). However, this algorithm is only appropriate for a situation with one PU, while real situations usually involve several PUs at the same time. The detection performance of this algorithm can drop off severely when several PUs exist.
So far, most eigenvalue-based spectrum sensing algo-rithms have been designed based on the analysis of detection performance using random matrix theory (RMT) [10] , and significant branches of RMT can be classified into asymptotic theory and non-asymptotic theory. Asymptotic theory studies the convergence property of the spectrum distribution function of the random matrix in infinite dimensional space. These concepts have now been quite thoroughly explored, and asymptotic theory has already been widely applied to all kinds of proposed eigenvalue algorithms. Meanwhile, non-asymptotic theory [11] is the latest achievement within RMT: the convergence property of the spectrum distribution function of the random matrix in finite-dimensional space. In an actual situation, the number of sampling points is usually small; thus non-asymptotic theory is more widely applicable.
Meanwhile, existing eigenvalue algorithms only use the characteristics of some particular eigenvalues, for example, the maximum eigenvalue, the minimum eigenvalue, and the ratio of the maximum and minimum eigenvalue. Other eigenvalues have not been used. Therefore, the purpose of this paper is to excavate potential uses of these neglected eigenvalues and design eigenvalue detection algorithms with joint distribution of eigenvalues, in order to improve the performance of spectrum sensing.
II. SYSTEM MODEL
A typical multi-antenna spectrum sensing scenario is as shown in Fig. 1 . In the space, the number of PUs is D and the number of SUs is K . The multi-antenna system is applied to the SUs, the receiving equipment has Nr receiving antennas, and each antenna has N sampling signals. There are two hypotheses on whether a PU’s signal exists. H 0 represents the situation in which the PU does not exist, and there is only noise, while H 1 represents the situation in which both the PU and noise exist. The spectrum sensing model of a multi-antenna CR system can be described as a binary hypothesis model:
PPT Slide
Lager Image
where i =1,2,3… Nr denotes i receiving antennas of the SU; yi ( n ) denotes the signal received from the PU at receiving antenna i ; xi ( n ) is the transmitting signal of the j th PU; hij ( n ) is the channel gain between the ith receiving antenna and the j th PU; and ui ( n ) is the Gaussian white noise with its mean value as 0 and variance as
PPT Slide
Lager Image
.
PPT Slide
Lager Image
Multi-antenna spectrum sensing scene. yi(n)
The statistical matrix of the received signal at the SU can be expressed as follows:
PPT Slide
Lager Image
where yi ( n ) , xi ( n ) , and ui ( n ) are an N ×1 vector. According to the definition of matrix, (1) can be expressed as
PPT Slide
Lager Image
The covariance matrix of the received signal is thus
PPT Slide
Lager Image
From (3) and (4),
PPT Slide
Lager Image
In an actual algorithm, to approximate an observed signal with finite sampling points N , (2) can be written in the form of an Nr × N estimate covariance matrix:
PPT Slide
Lager Image
With (4) recalculated,
PPT Slide
Lager Image
Suppose that the eigenvalues of the sampling covariance matrix
PPT Slide
Lager Image
are λ 1 = λ max λ 2 ≥…≥ λNr = λ min , and the eigenvalues of the covariance matrix
PPT Slide
Lager Image
are ρ 1 = ρ max ρ 2 …≥ ρNr = ρ min . From Eq. (7),
PPT Slide
Lager Image
If there is no PU signal, ρ max = ρ min =0 and thus
PPT Slide
Lager Image
Based on the analysis of the system model, we can use eigenvalues of the received signal to detect whether there is a PU. For the detection threshold γ of the detection algorithm, if the detection statistic T < γ , a PU does not exist, indicating that H 0 is confirmed; if T > γ , a PU exists, confirming H 1 .
III. SPECTRUM-SENSING ALGORITHM
Many spectrum-sensing algorithms have been developed based on eigenvalues. Relatively speaking, the maximum-minimum eigenvalue (MME) algorithm based on RMT is a classic one. In light of the drawbacks of its performance, taking the pre-difference, Lagrange’s interpolation, and the double threshold method, for example, is put forward. What is undeniable is that most of the algorithms have only used the properties of some eigenvalues. The algorithms proposed in this paper were developed with the goal of finding effective detection algorithms that use the joint distribution properties of the received signal’s eigenvalues. The existing algorithms and proposed algorithms can be classified as follows.
- A. Based On the Limit Distribution Character of the Maximum Eigenvalue
- 1) Maximum-Minimum Eigenvalue
When H 0 is true, the ratio of the maximum and minimum eigenvalues is 1, and when H 1 is true, the ratio of the maximum and minimum eigenvalues is
PPT Slide
Lager Image
According to the difference of the ratios under the two hypotheses, the MME algorithm uses λ max / λ min as the test statistic to decide whether the PU exists or not.
PPT Slide
Lager Image
If TMME > γMME , then there are PUs in the space. On the other hand, if TMME < γMME , then the space contains no PUs. γMME represents the detection threshold of the MME algorithm.
- 2) Maximum Eigenvalue Detection
Under a binary hypothesis, H 0 indicates a lack of PUs, and thus
PPT Slide
Lager Image
while H 1 indicates the presence of a PU and
PPT Slide
Lager Image
Therefore, the ratio of the maximum eigenvalue and noise variance under the two hypotheses can be used as the detection statistic, written as follows:
PPT Slide
Lager Image
If TMAX > γMAX , then there are PUs in the space, and vice versa. γMAX represents the detection threshold of the maximum eigenvalue detection (MED) algorithm.
- 3) Generalized Likelihood Ratio Test
Under a binary hypothesis regarding the CR system, the general model of the generalized likelihood ratio can be written as follows:
PPT Slide
Lager Image
In accordance with the maximum likelihood estimation, the detection statistic can be determined by
PPT Slide
Lager Image
If TGLRT > γGLRT , then there are PUs in the space, and vice versa. γGLRT represents the detection threshold of the generalized likelihood ratio test (GLRT) algorithm.
- B. Joint Distribution Based on Eigenvalues
- 4) Arithmetic-to-Geometric Mean
This algorithm uses the ratio of the arithmetic mean value and geometric mean value of the eigenvalues of the received signal covariance matrix as the detection statistic, which is indicated as follows:
PPT Slide
Lager Image
If TAGM > γAGM , then there are PUs in the space, and vice versa. γAGM represents the detection threshold of the arithmetic-to-geometric (AGM) algorithm.
- 5) Eigenvalue Energy Summation
The eigenvalues of the received signal covariance matrix can represent the energy of the signal. We perform energy summation of the eigenvalues and extract a root to find its amplitude, written as:
PPT Slide
Lager Image
If TEES > γEES , then there are PUs in the space, and vice versa. γEES represents the detection threshold of the eigenvalue energy summation (EES) algorithm.
- 6) Double Maximum Eigenvalue Detection
Under a binary hypothesis, H 0 indicates that no PUs exist, and thus
PPT Slide
Lager Image
while H 1 indicates the opposite and thus
PPT Slide
Lager Image
This algorithm is similar to the MED algorithm, but instead of the maximum eigenvalue, the sum of the two largest eigenvalues is used. The detection statistic can then be written as follows:
PPT Slide
Lager Image
If TDMED > γDMED , then PUs exist in the space, and vice versa. γAGM represents the detection threshold of the double maximum eigenvalue detection (DMED) algorithm.
- 7) Eigenvalue Weight
The eigenvalues of the received signal covariance matrix reflect its projected length in the direction of the eigenvectors to some extent, which is also the energy of the signal in this direction. Thus we consider making full use of this property of the eigenvalues, taking a proportion of the energy of each signal in the received signal as a weighted value, and the weighted energy as the detection threshold:
PPT Slide
Lager Image
If TEW > γEW , PUs exist in the space, and vice versa. γEW represents the detection threshold of the eigenvalue weight (EW) algorithm.
- 8) Multiple Eigenvalue Weight
For the eigenvalues of the received signal, we can conclude that
PPT Slide
Lager Image
where the left part of the inequality is the detection statistic of the EW algorithm. If the right part of the inequality is shifted to the left, and the new detection statistic is expressed as
PPT Slide
Lager Image
If TMEW > γMEW , PUs exist in the space, and vice versa. γMEW represents the detection threshold of the multiple eigenvalue weight (MEW) algorithm.
IV. SIMULATION RESULTS
In this paper, we suppose that the SU’s receiver has 4 antennas, D —the number of PUs in the space—is smaller than 4, and the transmitting signal of the PUs is binary phase shift keyed (BPSK). Meanwhile, we assume that the channels between the PUs and the SUs are Rayleigh channels. We initialize the false alarm probability as 0.01. For the number of sampling points, we set 3 values, N =4, N =8, and N =100, and compare the performance of each algorithm with each number of sampling points. All the results are averaged over 10,000 Monte Carlo realizations.
V. RESULTS OF ANALYSIS
From the simulation results shown in Figs. 2 4 , we can conclude the following.
PPT Slide
Lager Image
The detection probability under different SNRs (Pfa=0.01, PU=1, N=4, M=4).
PPT Slide
Lager Image
The detection probability under different SNRs (Pfa=0.01, PU=1, N=8, M=4).
PPT Slide
Lager Image
The detection probability under different SNRs (Pfa=0.01 PU=1 N=100 M=4).
When there is only one PU in the environment, as the number of sampling points N increases, the detection performance of each algorithm improves as well. The performance of EW, MED, DMED, and EES are better; GLRT, MEW, and AGM are ordinary; and MME performs the worst. Because of the existence of only one PU, eigenvalues are noise variance except for the maximum eigenvalue. The weight of the maximum eigenvalue in EW is approximately 1, and the second eigenvalue of DMED is noise variance, which can be ignored, as can be the energy of the noise component in EES. Thus the performance of these three algorithms are virtually the same as MED. MED and DMED are half-blind algorithms, since they both require noise information, while EW and EES are blind algorithms since they need no noise information and are not affected by noise. To sum up, EW is more practical in this case.
From the simulation results shown in Figs. 5 7 , we can conclude the following.
PPT Slide
Lager Image
The detection probability under different SNRs (Pfa=0.01, PU=2, N=4, M=4).
PPT Slide
Lager Image
The detection probability under different SNRs (Pfa=0.01, PU=2, N=8, M=4).
PPT Slide
Lager Image
The detection probability under different SNRs (Pfa=0.01, PU=2, N=100, M=4).
When there are two PUs in the space, as the number of sampling points N increases, the detection performance of each algorithm improves as well. The performance of EW, MED, DMED, and EES are better; AGM is ordinary; and MME performs the worst. When N is 4, GLRT and MEW perform badly. The detection probability can only reach about 0.65. When N increases, the performance of these two algorithms improves.
Regardless of whether one or two PUs are in the environment, the EW algorithm has favorable detection performance under a low SNR. For instance, when the detection probability comes to 0.6 and N =4, the perfor-mance gain of the EW algorithm is 21 dB higher than that of the MME algorithm, and 9 dB higher than that of the AGM. This algorithm is suitable for a situation in which the PU’s signals are unknown, the signal energy is inferior and the existence of the PUs needs to be quickly detected.
From the simulation result shown in Fig. 8 , we can conclude the following.
PPT Slide
Lager Image
The detection probability of the EW algorithm with different numbers of sample points (Pfa=0.01, PU=1, M=4).
For the EW algorithm, as the number of sampling points N increases, the detection performance improves as well. When the detection probability is 0.6 and the number of sampling points is smaller than 10, as the detection performance increases by one decibel, the number of sampling points increases by only 2 points. When the number of sampling points is larger than 500, as the detection performance increases by a decibel, the number of sampling points increases by about 500 points, and this demand will greatly increase as the number of sampling points increases. Thus, the EW algorithm can meet the requirements consistently and effectively under the conditions in which the signal and noise information is unknown, high system detection performance is needed, and the detection has to be implemented quickly.
VI. CONCLUSION
This paper has proposed various new spectrum-sensing algorithms based on the joint distribution of the eigenvalues in a multiple-user cognitive radio environment, which make full use of the impact of the eigenvalue component during detection, acquire favorable detection performance under a low SNR, and ensure the robustness of detection with a small number of sampling points. In a study to follow, we plan to analyze the proposed algorithms theoretically on the basis of RMT.
BIO
Syed Sajjad Ali
received the B.Sc. degree from the University of Karachi, Karachi, Pakistan, and the M.S. degree from the Mohammed Ali Jinnah University, Karachi, Pakistan. He worked as a lecturer in the Institute of Business and Technology, Karachi, Pakistan. Before, he also worked in Lucky Cement as a Network and Communication Engineer. He is currently pursuing the Ph.D. degree from School of Information and Communication Engineering, Dalian University of Technology, Dalian, P.R. China under the supervision of Prof. Minglu JIN. His research interests are in wireless communications and signal processing. He is particularly interested in Signal detection in Cognitive radio.
Jialong Liu
is currently studying for his M.S. in Electronic and Communication Engineering at Dalian University of Technology. He received his B.S. in Electronic and Information Engineering from the Department of Information and Communication Engineering, Dalian University of Technology, in 2012. His research interests include Spectrum Sensing in Cognitive Radio.
Chang Liu
received the B.Eng. degree in Electronic Information Engineering from Dalian Maritime University, Dalian, China. He is currently a Ph.D. candidate in School of Information and Communication Engineering, Dalian University of Technology, China. He is also a visiting scholar in department of Electrical Engineering & Computer Science at University of Tennessee, Knoxville, USA from 2015 to 2016. His research interests include Spectrum Sensing in Cognitive Radio, Statistical Signal Processing, Random Matrix Theory and Beamforming Technology.
MingluJin
is a professor at the School of Information and Communication Engineering, Dalian University of Technology, Dalian, China. He received his Ph.D. and M.Sc. degrees from Beihang University, Beijing, China, and his B.Eng. degree from University of Science & Technology, Hefei, China. He was a visiting scholar at the Arimoto Lab., Osaka University, Osaka, Japan, from 1987 to 1988. Further, he was Research Fellow at the Radio & Broadcasting Research Lab at Electronics Telecommunications Research Institute (ETRI), Korea, from 2001 to 2004. His research interests lie in the general areas of signal processing and communications systems. His specific current interests are cognitive radio, multiple input multiple output (MIMO) radio antenna design and wireless sensor networks.
References
Federal Communications Commission 2002 “Spectrum Policy Task Force Report,” Washington, DC ET Docket No. 02-135
Mitola J. , Ph.D. dissertation 2000 “Cognitive radio: an integrated agent architecture for software defined radio,” Royal Institute of Technology (KTH) Stockholm, Sweden Ph.D. dissertation
Yucek T. , Arslan H. 2009 “A survey of spectrum sensing algorithms for cognitive radio applications,” IEEE Communications Surveys & Tutorials 11 (1) 116 - 130    DOI : 10.1109/SURV.2009.090109
IEEE 802.22.1-2010 Standard for the Enhanced Interference Protection of the Licensed Devices [Internet] Available:
Kay S. M. 1998 Fundamentals of Statistical Signal Processing: Detection Theory. PTR Prentice-Hall Englewood Cliffs, NJ
Sahai A. , Cabric D. “Spectrum sensing: fundamental limits and practical challenges,” in Proceedings of IEEE International Symposium on New Frontiers in Dynamic Spectrum Access Networks (DySPAN) Maui, HI 2005
Urkowitz H. 1967 Energy detection of unknown deterministic signals,” Proceedings of the IEEE 55 (4) 523 - 531    DOI : 10.1109/PROC.1967.5573
Zeng Y. , Liang Y. C. “Maximum-minimum eigenvalue detection for cognitive radio,” in Proceedings of IEEE 18th International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC2007) Athens, Greece 2007 1 - 5
Wang P. , Fang J. , Han N. , Li H. 2010 “Multiantenna-assisted spectrum sensing for cognitive radio,” IEEE Transactions on Vehicular Technology 59 (4) 1791 - 1800    DOI : 10.1109/TVT.2009.2037912
Tulino A. M. , Verdu S. 2004 Random Matrix Theory and Wireless Communications. Now Publishers Inc. Hanover, MA
Rudelson M. , Vershynin R. “Non-asymptotic theory of random matrices: extreme singular values,” in Proceedings of the International Congress of Mathematicians (ICM’10) Hyderabad, India 2010