Abstract

Eigenvalue based spectrum sensing can make detection by catching correlation features in space and time domains, which can not only reduce the effect of noise uncertainty, but also achieve high detection probability. Hence, the eigenvalue based detection is always a hot topic in spectrum sensing area. However, most existing algorithms only consider part of eigenvalues rather than all the eigenvalues, which does not make full use of correlation of eigenvalues. Motivated by this, this paper focuses on multiantenna system and makes all the eigenvalues weighted for detection. Through the analysis of system model, we transfer the eigenvalue weighting issue to an optimal problem and derive the theoretical expression of detection threshold and probability of false alarm and obtain the close form expression of optimal solution. Finally, we propose new weighting schemes to give promotions of the detection performance. Simulations verify the efficiency of the proposed algorithms.

1. Introduction

The rapid development of wireless services leads to the scarcity of the public radio spectrum becoming more and more serious. Traditionally, licensed spectrum is allocated over relatively long time periods and is intended to be used only by legitimate users. Cognitive radio (CR) technology was proposed to handle the contradiction between the shortage of spectrum resource and the underutilization of licensed spectrum [1, 2]. Spectrum sensing which is a fundamental task of CR is aimed at obtaining the awareness of licensed spectrum usage and existence of primary users (PUs) in a specific geographical location [37]. The main function of spectrum sensing is to frequently explore the spectrum holes for the secondary users (SUs) by detecting the presence of primary users so that the SUs can share the licensed spectrum. Therefore, spectrum sensing becomes critical in cognitive radio system.

There have been many discussions and proposed solutions for spectrum sensing [8]. Of these methods, likelihood ratio test (LRT) [9], cyclostationary detection (CSD) [10, 11], and matched filtering (MF) detection [12, 13] can achieve optimal performance while requiring both source signal and noise power information, which is not available in practice. Hence, semiblind methods such as energy detection (ED) [9, 14] and maximum eigenvalue detection (MED) [15] are proposed. Among these, ED is the most commonly chosen scheme for study and implementation due to its relatively low complexity and satisfactory performance under low signal-to-noise ratio (SNR) environment. However, ED heavily relies on the accuracy of the knowledge of noise power which is generally changing over time. This so-called noise uncertainty problem [16] can significantly degrade the performance of ED algorithm.

To overcome these shortcomings, blind detection algorithms which require no information on source signal or noise power have been intensively studied recently [1721]. The classical blind detection algorithms are the eigenvalue based methods. For example, maximum-minimum eigenvalue (MME) detection [22], arithmetic to geometric mean (AGM) detection [23], and signal-subspace eigenvalues (SSE) method [23] can overcome the shortcoming of ED and achieve outstanding performance. On the other hand, eigenvalue based methods have also been studied in new scenarios, such as cooperative adaptive versions [24] and Multiple Primary Transmit Power (MPTP) scenario [25].

However, most algorithms only consider part of eigenvalues, such as maximum, minimum, and mean value, which does not make full use of all the eigenvalues to make detection. Motivated by this, we focus on the problem of eigenvalue weighting in multiantenna system and analyze the related problems. By analyzing the model of eigenvalue weighting, we transfer the weighting problem to an optimal problem. Using the latest random matrix theory (RMT) [26, 27], we derive the close form expression of probability of detection and probability of false alarm and obtain the optimal solution. Finally, we propose new weighting schemes to give promotions of the detection performance. Simulations verify the efficiency of the proposed algorithms. The main contributions of this paper include the following:(i)Different from the traditional eigenvalue based detection, we consider making detection by utilizing all of the eigenvalues in the multiantenna system. By transferring the weighting problem to an optimal problem, we analyze and derive an energy based maximum ratio combination (EN-MRC) method.(ii) Considering the case of correlated signals is common in applications, we use the idea of MRC weighting in EN-MRC method to design an eigenvalue based MRC (EIG-MRC) scheme: signal eigenvalue weighting (SEW) based detection, which needs the a priori information of signals’ covariance matrix and noise power.(iii) To make the detection more practical, we use the maximum likelihood estimation (MLE) approach to design a method of signal eigenvalue approximation weighting (SEAW) based detection, in which only the noise power is needed.

The rest of the paper is organized as follows. Section 2 explains the system model. The eigenvalue weighting based detection is studied in Section 3. Section 4 presents simulation results and conclusion is presented in Section 5. Some notations used in the paper are listed as follows: superscripts and stand for transpose and Hermitian transpose (transpose-conjugate), respectively.

2. System Model

Figure 1 illustrates a classical multiantenna spectrum sensing scenario with randomly distributed primary users (PU in figure) and randomly distributed secondary users (SU in figure). Once PUs begin to communicate, the surrounding SUs can receive the PU signals and then capture the samples to operate the spectrum sensing.

According to Figure 1, the SUs are equipped with receiving antennas and there are PU signals arriving in the antenna array. In this case, the sensing problem in multiantenna cognitive radio system can be written aswhere represents the th receiving antenna and is the th sample. is the sample of the th receiving antenna. is the channel gain between the th PU signal and the th receiving antenna. is the additive white Gaussian noise (AWGN) with 0 mean, variance.

Stacking the samples at the same time we can get the following receiving vector of antenna array:Hence, formula (1) can be rewritten as the matrix form:where and are the antenna receiving matrix and noise matrix, respectively. is the channel gain matrix of the signal matrix .

3. Eigenvalue Weighting Based Detection

3.1. Fundamental of Eigenvalue Weighting Based Detection

Based on (3), the corresponding covariance matrix can be written asHence, we can rewrite asLet and represent the eigenvalues of and , respectively. Obviously, when PUs are present, we can get ; when PUs are absent, that is, , we can have .

Based on the analysis above, we can make detection by weighting the eigenvalues. Considering the number of samples is finite in reality we can get the following test statistic:where is the eigenvalues and is the weighting coefficient. is the samples covariance matrix. Obviously, if ( is the test threshold), then PUs are present; otherwise, PUs are absent.

Finally, we summarize the general eigenvalue weighting algorithm steps as follows.

Eigenvalue Weighting Based Spectrum Sensing Algorithm for Multiantenna Cognitive Radio System

Step 1 (compute the sample covariance matrix of the received signal). Since the number of samples is finite, we can only use the sample covariance matrix .

Step 2 (obtain the eigenvalues of sample covariance matrix). Make eigenvalue decomposition (EVD) of , obtain eigenvalues, and sort them in a descending order: .

Step 3 (calculate the test statistic of the eigenvalue weighting). Let all the eigenvalues be weighted by and compute the sum of them. Thus, we can obtain the test statistic in (6).

Step 4 (decision). If , then signal exists (“yes” decision); otherwise, signal does not exist (“no” decision), where is a threshold.

3.2. Theoretical Analysis of Eigenvalue Weighting Based Detection

Note that how to select weights is of great importance, which can affect the performance of the algorithm directly. Based on the Neyman-Pearson rule, we can express the weighting selection problem as the following optimal problem [28, 29]:where is the weighting coefficient vector; and represent the probability of detection and the probability of false alarm. and are the probability density function of test statistic under and , respectively.

Based on (6), we find that it is possible to analyze the distribution of the test statistic whereas the joint probability density function is rather complex, whose close form expression is not available. However, we can transfer the problem of eigenvalue weighting of the matrix to a problem of the trace of a new matrix and the analysis of distribution of the trace is a simple problem. The detailed analysis is showed in the following.

Let and . Hence,where . When the number of samples tends to infinite, the tends to a diagonal matrix and we can get the following:where represents the diagonal matrix of eigenvalues and the equality holds when the number of samples tends to infinite. Hence, if we make eigenvalue weighting of by and calculate the sum of the eigenvalues after weighting, then it is equivalent to compute the trace of . Since , we can rewrite (6) as the following:where is the th row, th element of . Let and thus the test statistic can be written asFor simplification, we assume the noise variance . When the number of samples is large enough we can get the following expression based on central limit theorem (CRT):where is the power of the PU signals. Therefore, we can obtain the expressions of and , respectively:where . Hence, based on (13) and (14), we can finally get the expression aswhere and . Since the SNR of spectrum sensing is rather low (−20 dB), which leads to , we can get . Hence, (15) can be approximated asTherefore, problem (7) can be rewritten asNote that this problem can be solved by Lagrangian multiplier method and the solution is written asLet and we can finally get the weighting coefficient:Note that this weighting scheme is exactly identical to the maximal ratio combination (MRC) weighting scheme in [28, 29] and we call it energy based MRC (EN-MRC) detection. Hence, by studying the idea of MRC weighting scheme, we apply this idea into eigenvalue weighting and finally develop a kind of energy based MRC algorithm. The test statistic can be written aswhere is the power of the PU signals.

3.3. Eigenvalue Weighting Based Detection

Note that the transformation from eigenvalue to energy in (9) is approximately equivalent and the equality holds when tends to infinite. Hence, the corresponding analysis should be more accurate when the number of samples tends to be very large. On the other hand, the analysis under this case is based on the assumption that the received signals are independent and identically distributed (i.i.d.) for each other, which is not very accurate for the case of highly correlated signals. For example, as for (12), the distribution of under is considered as a linear combination of Gaussian variables with -variance, which is based on the assumption that the received signals under are i.i.d. for each other. However, this assumption is only available under a cooperative spectrum sensing model whose samples are collected from different sensing nodes. Hence, the weighting coefficient in (20) is not an appropriate weighting scheme especially for the case of highly correlated signals and thus it needs to be improved for better catching the signals’ correlation.

Motivated by this, we try to analyze the weighting coefficients from the aspect of eigenvalue directly. Since is the power of the PU signals and , we can then replace the power of the PU signals in (20) with the eigenvalues of signal covariance matrix

In this case, the test statistic can further capture the correlation among signals and may achieve better performance especially when there are highly correlated PU signals. Hence, we propose a signal eigenvalue weighting (SEW) based detection and the test statistic is given aswhere and are the eigenvalues of sample and signals’ covariance matrix, respectively. Although the SEW based detection may perform better performance, it is not available in practice as it needs the a priori information of the channel, signal, and noise. Hence, we try to use the maximum likelihood estimation (MLE) of these parameters to design semiblind detection, in which only noise power is needed. Hence, we will analyze and derive the MLE of eigenvalues of the PU signals’ covariance matrix in the following.

According to the analysis in [23], the MLE of signals’ covariance matrix can be expressed aswhere is the eigenvector of sample covariance and represents the maximum between and . Hence, the MLE of eigenvalues of PU signals’ covariance matrix can be written asSubstituting the MLE of signal eigenvalue into (21) we can obtain the test statistic of signal eigenvalue approximation weighting (SEAW) based detection asAs a summary, we propose three weighting schemes: one is traditional MRC based detection (i.e., EN-MRC) and the other two are improvement eigenvalue weighting schemes, that is, SEW based detection and SEAW based detection. For the convenience of comparison, we summarize these three methods in Table 1.

Remark. Since eigenvalue weighting problem can not be solved directly, we first loose the constraint conditions and assume that the PU signals follow the i.i.d. model and the number of samples is very large. In this case, we can obtain an inaccuracy solution: EN-MRC. Based on the MRC weighting scheme, we then tighten the constraint conditions and modify the assumption to make it satisfy the requirements of the practical system, that is, correlated signal model. Considering the eigenvalues can further capture the correlations of signals, we finally replace the energies with eigenvalues and design the eigenvalue based MRC (EIG-MRC) schemes: SEW and SEAW based detection. The corresponding illustration is shown in Figure 2.

4. Simulations and Discussions

This section provides some simulation results for multiantenna cognitive radio systems in the MATLAB environment. Since this paper focuses the eigenvalue weighting schemes for spectrum sensing, we will compare the proposed EN-MRC, SEW, and SEAW based detection with eigenvalue based methods, including MED, MME, and AGM detection. We assume there is 1 PU or 2 PUs transmitting signal over the Nakagami- () channel in presence of AWGN. The SUs are equipped with 4-element antenna array. The stopping criterion set is at 10 000 iterations and the is set as 0.1 (this has been specified as the maximum allowable by the WRAN 802.22 working group).

The simulation results of detection performance in terms of number of samples with 1 PU and 2 PUs are presented in Figures 3 and 4, respectively. It is shown that when compared with eigenvalue based methods, such as MED, MME, and AGM, the proposed SEW and SEAW based detection perform much higher probability of detection with different SNRs, while the EN-MRC performs a relatively lower detection probability when compared with MED, MME, and AGM detection. It is because algorithms SEW and SEAW are regarded as “EIG-MRC” weighting scheme and MED and AGM belong to the selection combination (SC) and equal gain combination (EGC) weighting schemes for eigenvalues, respectively. As for EN-MRC, it is the energy based weighting coefficients, which can not fully capture the correlations. In addition, the MME is just a kind of partial eigenvalue based nonweighting detection and thus it has limited detection performance. However, since low SNR approximation has been adopted to derive the EN-MRC scheme, the EN-MRC is able to achieve a relatively higher detection probability. For example, the EN-MRC is slightly better than MME and AGM when the SNR is ranging from −35 dB to −13 dB. On the other hand, when the SNR increases, the probability of detection of EN-MRC drops a little and presents a slightly worse performance (since the number of eigenvalues for the simulation is very small, the advantages of making detection by using all the eigenvalues or the energies are not obvious, which means the AGM or EN-MRC may not achieve a better performance than MME). When comparing Figure 3 with Figure 4, we can find that the performance increases with the increasing number of PUs, such as a nearly 30% detection probability improvement in terms of −15 dB.

Similarly, Figures 5 and 6 present the simulation results of probability of detection in terms of number of samples with 1 PU and 2 PUs. Again, the proposed SEW and SEAW methods achieve a higher detection performance and the EN-MRC outperforms MME and AGM under low SNRs. Hence, the simulation results can further verify that it is just the replacement of energy with eigenvalue that leads to the high improvements in terms of detection probability.

In addition, as for the three new methods, we can find that SEW performs the best among these proposed methods, EN-MRC performs the worst, and the performance of SEAW is between these two methods. For example, the probability of detection of SEAW with 2 PUs (i.e., SEAW in Figure 3) is 0.5 in terms of SNR = −15 dB, which is in the middle of of SEW (i.e., 1) and of EN-MRC (i.e., 0.2).

According to Figures 36, a more interesting phenomenon can be found; that is, the SEAW’s performance shifts from the lower area (close to EN-MRC) to a higher area (close to SEW) with the increasing of number of samples and number of PUs, which is like a kind of lower and upper bounds of the performance of SEAW. If we consider the performance-complexity tradeoff, the proposed SEAW can be selected as an alternative for its low complexity and relatively better performance. Hence, the SEAW may be more suitable for the application in reality.

5. Conclusion

This paper focuses on the problem of the eigenvalue weighting based spectrum sensing in multiantenna cognitive radio system. Through the analysis of system model, we transfer the eigenvalue weighting issue to the energy based weighting problem and derive the theoretical expression of detection threshold and probability of false alarm and finally obtain the close form expression. Considering the case of correlated signals is common in applications, we then design the signal eigenvalue based detection methods and they can achieve more higher detection probability. Simulation results verify the efficiency of the proposed algorithms.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This research was supported by the MSIP (Ministry of Science, ICT and Future Planning), Korea, under the ITRC (Information Technology Research Center) support program (IITP-2016-H8501-16-1019) supervised by the IITP (Institute for Information & Communications Technology Promotion).