Time domain averaging and correlation-based improved spectrum sensing method for cognitive radio

Based on the combination of time domain averaging and correlation, we propose an effective time domain averaging and correlation-based spectrum sensing (TDA-C-SS) method used in very low signal-to-noise ratio (SNR) environments. With the assumption that the received signals from the primary users are deterministic, the proposed TDA-C-SS method processes the received samples by a time averaging operation to improve the SNR. Correlation operation is then performed with a correlation matrix to determine the existence of the primary signal in the received samples. The TDA-C-SS method does not need any prior information on the received samples and the associated noise power to achieve improved sensing performance. Simulation results are presented to show the effectiveness of the proposed TDA-C-SS method.


Introduction
Cognitive radio (CR) networks allow unlicensed (or secondary) users to opportunistically exploit the underutilized spectrum bandwidth of the licensed (or primary) users. Spectrum sensing is a key operation performed by the CR networks to determine the spectrum holes of the spectrum allocated to a primary user. In the literature, several kinds of typical spectrum sensing methods have been reported, including energy detection methods [1,2], matched filter detection methods [3,4], and cyclostationarity feature detection methods [5]. The energy detection methods need the prior knowledge of noise power and are vulnerable to the noise uncertainty. The matched filter detection methods need to know the waveform of the primary user in advance. The methods of cyclostationarity feature detection require the information on the cyclic frequencies of the primary user. The prior knowledge requirements of these methods often limit their realistic applicability. To avoid these undesirable requirements, *Correspondence: shli@sjtu.edu.cn 1 Department of Electronic Engineering, Shanghai Jiao Tong University, Dongchuan Road, Shanghai 200240, China Full list of author information is available at the end of the article some new statistical covariance-based methods for spectrum sensing in additive i.i.d white noise environments have been proposed [6][7][8]. These reported methods do not need the prior information about the signal and the noise power to achieve a sensing performance in a low signal-to-noise ratio (SNR) environment of -22 dB. However, these methods assume that the primary user's signal is a stationary random process and that all the received data samples contain the primary user's signal if it exists. In fact, most existing spectrum sensing methods require the last assumption. In fact, it is very possible that only a part of the received samples in practice contain the signal from the primary user or a part of or all the received samples are from some signals that are not stationary. For example, the initial transient signal, known as fingerprint of a wireless device [6], from the primary user is not stationary. If such samples are used, the methods reported in [7][8][9] will lose the supporting prerequisite and become unusable. Even in the cases that these required assumptions are valid, it is difficult for the methods in [7][8][9] to obtain the theoretical sensing performance in lower SNR environments when the number of the received samples available is limited. http://asp.eurasipjournals.com/content/2014/1/46 To overcome these problems and support more effective sensing in low SNR environments, we propose an effective time domain averaging and correlation-based spectrum sensing (TDA-C-SS) method based on time domain averaging and correlation for spectrum sensing in additive i.i.d white noise environments. Similar to the methods reported in [7][8][9], the proposed TDA-C-SS method does not need any aforementioned prior information and achieve desirable sensing performance in very low SNR environments. The TDA-C-SS method assumes the signal from the primary user to be deterministic. By making use of time domain averaging, the SNR of the received samples is increased. Then, the task of spectrum sensing is obtained by performing the correlation operation. Simulation results from various environments are presented to show the effectiveness of the proposed TDA-C-SS method.
The rest of this paper is organized as follows. Section 2 describes the system model and sampling operations to obtain the input sample sequence. Section 3 gives the proposed TDA-C-SS spectrum sensing method. Simulation results and discussions are presented in Section 4. Finally, conclusion is drawn in Section 5.

System model and sampling description
In a duration, [ t 1 , t 2 ], of spectrum sensing, let z(t) be the continuous-time signal received by the secondary user. The process of the spectrum sensing is regarded as a twohypotheses-test problem at a time duration, [ t 1 , t 2 ], i.e., H 0 : the primary signal is absent and H 1 : the primary signal is present, described by and s pri−c (t) is any received primary signal, and w c (t) is the i.i.d white noise with a zero mean and a variance represented by δ 2 . Let us be interested in the frequency band with central frequency f 0 and a narrow bandwidth W. The signal z c (t) defined in (1) is sampled at a frequency f s ≤ L mul (f 0 + W /2), where, L mul is a large constant, for example, L mul = 100. The choice of such values of f s and L mul is for the purpose of time domain averaging in our proposed method, which is to be described in Section 3. The two hypotheses in (1) are then expressed as It is worth noting that if f 0 + W /2 is very high, it is difficult to implement the sampling process because the available high speed samplers can support up to several tens of GHz [10][11][12]. Therefore, we should suitably select the value of f 0 + W /2 to allow an appropriate large value of L mul possible. It is also possible to use down conversion process [12] to generate the received samples when the value of f 0 + W /2 is very high.

The proposed spectrum sensing method
The proposed TDA-C-SS method treats the samples of the primary users to be deterministic. It firstly averages the received data samples in the time domain to increase SNR (or reduce the noise) then obtains the spectrum sensing by means of correlation operations.

Sample time domain averaging
Time domain averaging is an effective method to decrease noise power for periodic signal detection [13] (and therein [13]). For the obtained data sample sequence z(n), n ∈[ 1, N 0 ] , the time average operation is defined by where M is a small positive integer, for example, M=5, and If the samples from the primary users are treated as deterministic ones, f s is L mul times of f 0 + W /2 and M is small, it means that with the assumption that w c (t) is i.i.d white noise, we easily obtain the following from (2): denotes the variance of a random variable. It can be easily proved that the SNR of the sequence produced by the time averaging process is increased by 10lg(2M + 1) dB. If M = 5, for example, the SNR of the time averaged output is increased by about 10 dB. This gain in SNR is very valuable for effective spectrum sensing in the environment of strong noise.

Correlation operation and sensing decision
Similar to the method in [7], let us perform the correlation operation on the input samples and make a decision on the signal presence of the primary user based on a constructed correlation matrix. Being different from that in [7], however, the correlation used here is to be discussed from the view point of deterministic signal samples.
Let us consider the correlation defined by r e (i) = 1 N N n=1 e(n)e(n − i). Then, under H 0 , we have Under H 1 , we have It is noted that the samples from the i.i.d white noise, ω c (t) , which has a zero mean and a variance of δ 2 2M+1 , the product of s(nL mul )ω(n − i) are the samples from mutually independent random variables with a zero mean and a variance of s 2 (nL mul )δ 2 2M+1 . Therefore, based on the central limit theorem, it is concluded that the term  (6) when N → ∞. In addition, according to statistics theory, it is also noted that when N → ∞, (5) or the last term in (6) is the expectation operation. Finally, the first term in (6) Let us construct an L × L correlation matrix in the same form as that in [6], i.e., Note that we have tested P d versus L by experiments, and the experiment results show that when L ≥ 10, Pd is not affected greatly by L, as a result, we choose the integer L = 10 for the purpose of lower computation complexity. Based on the previous analysis of r e (i) values, apparently, the sensing test statistics and decision policies of the computer-assisted videodensitometry (CAV) method and the GCBA method in [7] are still applicable for our method. For the sake of simplicity, here, we use the test statistics and decision policy of the CAV method in [7], i.e., the test statistics is defined by [7] where r nm is the element located in the nth row and mth column of the matrix R e . With a positive threshold value λ, the sensing decision for the primary signal is made by

The TDA-C-SS Method
Based on the discussion above, we propose the an improved spectrum sensing method as follows.
Step 1: sample the received signal at frequency f s to obtain the discrete samples z(n), n = 1, 2, . . . , N 0 .
Step 4: calculate T in (8), and properly choose the value of λ for the sensing decision.
Similar to the approach reported in [7,9], let us use the computer simulation approach based on the given probability of false alarm, P f , to choose the threshold λ. That is, first, a P f value is given and white noise is generated as the input, and then with a number of simulation results http://asp.eurasipjournals.com/content/2014/1/46 of T in (8), the threshold value, λ, is selected to meet the requirement of P f .

Performance analysis of the proposed TDA-C-SS
From the description above, it is seen that the proposed TDA-C-SS method does not need any prior information about the waveform and the cyclic frequencies of the primary user's signal and the noise power, which is similar to those methods in [7][8][9]. In contrast, our proposed method is also valid when only part of the received sample sequence contain the primary user's signal, which is the main different from those in [7][8][9] that assume that the primary user's signal is a stationary random process and that all the received samples must contain the primary user's signal. Therefore, the proposed method is more general and flexible. By the averaging operation in the time domain, the proposed TDA-C-SS method is able to achieve an SNR improvement by 10lg(2M + 1) dB compared with the CAV method in [7] for the same values of N and L. Therefore, the proposed one is expected to improve the sensing performance substantially particularly in low SNR environments.

Simulations
In this section, simulation results are reported for the following three signal settings in the AWGN environments to verify the effectiveness of the proposed method.
Case I: the signal from the primary user is stationary. During spectrum sensing, all the received samples contain the signal of the primary user. Case II: the signal from the primary user is stationary. During spectrum sensing, only a part of the received samples contain the signal of the primary user. Case III: the signal from the primary user is not stationary.
For comparison, we also present the simulation results of the CAV method [6] applicable for case I. For all the simulations, the values of threshold λ are chosen by the computer simulation approach described previously, and 1,000 Monte Carlo runs are carried out to estimate the value of λ.

Simulation for case I
We use a wireless microphone signal generated by the method in [14] with the following parameters: central frequency f 0 = 100 MHz, bandwidth W = 36.8 KHz. Based on the discussion in Section 2, the sampling frequencies used for the proposed TDA-C-SS method and the CAV method are f s = 10.3 GHz ≥ L mul (f 0 + W /2), where, L mul = 100, and f s = 103 MHz ≥ f 0 + W /2, respectively. Figure 1 illustrates the detection probability P d versus SNRs with P f = 0.1, L = 10, N = 50, 000. Figure 2 presents the detection probability P d versus the false alarm probability, P f , with SNR = −30 dB, L = 10, N = 50, 000, for the proposed TDA-C-SS method, with different values of M, and the CAV method. From Figures  1 and 2, it is seen that under the same conditions, the detection performance of the proposed method is better than that of the CAV and can be improved further with increase of M. In addition, Figure 1 shows that, compared with the CAV method, the TDA-C-SS method achieves the improvement on the sensing performance by about 10lg(2M +1)dB in SNR. To achieve P d = 0.985, the values of the SNR needed by the CAV method and the TDA-C-SS with M = 2, 5, and 10 are −18.5, −26, −29.3, and 32.2 dB, respectively.

Simulation for case II and case III
In case II, the received sample sequence is 45, 000 zero-valued signal samples followed by 5,000 wireless microphone signal samples that are generated in case I. The segment of zero-valued signal samples is included for the frequently observed situation in which the primary user begins to occupy its channel after the secondary user has monitored the channel for certain period of time. The sampling frequency for the TDA-C-SS method is f s = 10.3 GHz. For case III, we use a fingerprint signal, which is generated by using the unit-step response of a second-order RLC parallel circuit described by a differential equation: It is not reasonable to use the CAV method because the signals used in case II and case III are not stationary. Therefore, we only verify the effectiveness of the proposed TDA-C-SS method. Figures 3 and 4 give, respectively, detection probability versus SNRs for the two cases with P f = 0.1, L = 10, N = 50, 000, and the detection probability P d versus the false alarm probability P f with SNR = −30 dB, L = 10, N = 50, 000. From these figures, it is seen that the TDA-C-SS method is still effective because this method does not consider the properties of the received signal. Similarly, the detection performance can be improved further with the increased value of M.

Conclusion
In this paper, an improved spectrum sensing method, TDA-C-SS, based on time domain averaging and correlation  has been proposed. The time domain averaging process has been typically used to decrease the noise effects, and correlation matrix has been constructed to decide the existence of the primary user's signal. In comparing with other reported method, such as CAV method, the proposed one can sense a primary user's signal in the white noise environment in very low SNR environments without requiring any prior knowledge about the signal and noise power. In particular, the proposed method is flexible to effectively sense the signals that are not stationary. Our simulation results have shown the desirable advantages of the proposed methods.