Spectrum sensing based on cumulative power spectral density

This paper presents new spectrum sensing algorithms based on the cumulative power spectral density (CPSD). The proposed detectors examine the CPSD of the received signal to make a decision on the absence/presence of the primary user (PU) signal. Those detectors require the whiteness of the noise in the band of interest. The false alarm and detection probabilities are derived analytically and simulated under Gaussian and Rayleigh fading channels. Our proposed detectors present better performance than the energy (ED) or the cyclostationary detectors (CSD). Moreover, in the presence of noise uncertainty (NU), they are shown to provide more robustness than ED, with less performance loss. In order to neglect the NU, we modified our algorithms to be independent from the noise variance.


Introduction
Due to an increasing demand of wireless devices and the limitation of natural spectrum resources, the cognitive radio (CR) has been introduced to optimize the use of the available spectrum [1].In CR network, primary (PU) and secondary (SU) users can share the same bandwidth.PU is a licensed user, whereas SU is an opportunistic user.To avoid any interference with the PU's signal, SU can only be active when PU is idle.Therefore, SU should continuously monitor PU status (active or idle).
Many factors (such as low signal to noise ratio (SNR), shadowing, channel fading, etc.) lead to a situation where SU is no longer able to correctly diagnose the status of the PU.To overcome this problem, new techniques for cooperative spectrum sensing (CSS) have been proposed [2,3].In the CSS context, several SUs cooperate to reach a decision on the presence of the PU signal; a fusion center (FC) should make the final decision.We can distinguish two schemes of CSS: the soft and the hard combining schemes.In the soft combining scheme, each SU sends its own measured test statistic to the FC who performs a linear combination of all received test statistics and compares the result to a threshold.With the hard combining scheme, each SU compares its own test statistic to a predefined threshold, makes his decision on the PU status, and sends it to the FC, where a final decision is made by combining the decisions of all the SUs using a predefined strategy (such as OR, AND, or a majority strategy).
In the literature, many techniques have been proposed to perform test statistics.Most known methods are energy detection, cysclostationary detection, and waveform detection.The last one requires perfect information about the PU signal [3,4].Hence, the use of this method is limited to the case of cooperative relationship between PUs and SUs.
Due to its simplicity and the fact that no prior information about the PU signal is required, the energy detector (ED) is widely used in spectrum sensing.ED consists of measuring the energy of the received signals and comparing it to a predefined threshold based on the known noise variance.Therefore, ED requires an accurate estimation of the noise variance; otherwise, we can face a SNR-wall problem, where the ED becomes incapable to made a robust decision about the PU status, even with a very large observation time [5,6].However, many studies have been proposed to enhance the energy detector (ED) performance and to overcome its limitations [5,[7][8][9].
The cyclostationary features detection (CSD) [10] is also widely used in spectrum sensing.CSD requires the knowledge of the cyclic frequencies of the PU, in addition to a long observation time and computing efforts.In [11], a cooperative CSD is proposed, where each SU performs a single cyclic-frequency detection and sends his decision to a FC to make the final decision.The authors of [12] analyzed theoretically the performance of the algorithm proposed in [11] in multi-path and log-normal channels.In [13], a blind CSD is proposed.This detector performs the spectrum sensing without the knowledge of the cyclic frequencies.
Although energy and cyclostationary detectors are widely used in the field of spectrum sensing, various other methods are also proposed [14][15][16].The goodness of fit (GoF) algorithm introduced in [14] compares the empirical distribution of the received samples to a known distribution of the noise (when PU is idle).If the empirical distribution is not compatible with the known distribution of the noise then the PU signal exists.To enhance the performance, [15] extends the algorithm of [14] by using the square of the received samples instead of the samples themselves.
By assuming the oversampling aspect of the baseband received signal (i.e., number of samples per symbols N s ≥ 2), the autocorrelation for a non-zero lag only vanishes when the PU signal is absent and the channel is only occupied by a white noise [16].The corresponding test statistic combines linearly the autocorrelation measures for different non-zero lags before making a decision on the PU status.The performance of this algorithm increases with N s .
In this paper, we propose new spectrum sensing detectors based mainly on the cumulative sum of the power spectral density (PSD) of the received signal 1 .It is known that the PSD of a white noise is flat.However, PSD losses this property with an oversampled PU signal.If the PU is absent, the cumulative sum of the received signal PSD has a close shape to a straight line.Whereas, a curved shape is obtained when PU exists.
To enhance the robustness of our contribution, hard and soft combining schemes are introduced.In those two schemes, the spectrum is divided into two parts: at first, the negative frequency points are considered while the second part deals with the positive frequency points.Hence, two test statistics are calculated based on the cumulative PSD of each part, and they are then combined according to the considered scheme.
The false alarm and detection probabilities are evaluated analytically under Gaussian and Rayleigh fading channels.Our detectors are compared to ED and CSD [8,16].Our detectors present a better performance than the energy detector, even at N s = 2 samples per symbol, where the CSD detector provides a poor performance relatively to ED.Furthermore, our detectors are less sensitive to the noise uncertainty than ED.In particular, we demonstrate that our detectors can be modified to become independent from the noise variance.This case represents an important advantage in real scenarios.
The rest of this paper is organized as follows.In Section 2, the system model and the spectrum sensing hypothesis are presented, in addition to an overview of the PSD and its estimation.Our proposed detectors based on the cumulative sum of the PSD are discussed in Section 3. Section 4 provides an analytic study on the statistical distributions of the test statistics as well as the calculus of the false alarm and detection probabilities.In Section 5, the probability of detection over Rayleigh flat-fading channel is provided.The numerical results of our detectors will be presented in Section 6.The effects of the noise uncertainty problem on our detectors are shown in Section 7. To overcome the noise uncertainty problem, this section presents modified versions of our detectors which are independent of the noise variance.At the end, a conclusion and perspective section of our work is provided.

System model and generality
The spectrum sensing consists in making a decision on the presence of PU in a bandwidth of interest (BoI).The PU baseband signal, s(n), can be modeled as follows: b m are the symbols to be transmitted; g(n) is the shaping window; N s satisfies the Nyquist criterion; N s = F s B ≥ 2 samples per symbol (sps), where F s is the sampling frequency and s p (n) and s q (n) are respectively the real and imaginary parts of s(n).s(n) has a bandwidth B and an even power spectral density (PSD).s(n) is assumed to be complex-valued zero mean unknown deterministic 2 signal [11,12,[17][18][19].
The presence/absence of PU can be presented in a classic Bayesian detection problem.Under H 0 , the PU is absent, whereas under H 1 PU exists.
where h is the complex channel gain.w(n) is N (0, σ 2 w ), where N (m, V ) stands for a normal distribution of a mean m and a variance V. Further, w(n) = w p (n) + jw q (n) is an i.i.d complex circular symmetric, i.e., E[ w 2 (n)] = 0 and the real part, w p (n), and the imaginary part, w q (n), of w(n) are independent Gaussian processes with equal variance. Where ] and E[ .]stands for the expectation.Without loss of generality, we assume that s(n) is a unit power signal.In this case, the signal to noise ratio (SNR), γ , is defined as follows: A wrong decision about the channel status can affect either the PU transmission or the efficient use of the channel.In fact, a missed detection can cause a harmful interference by the transmission of the SU in the same band of the PU.A false alarm, however, decreases the profit of the opportunity of the channel.Therefore, the probability of detection (p d ) should be increased as much as possible, by keeping the probability of false alarm (p fa ) low.Neyman-Person's detection method consists in a trade-off between a high p d and a low p fa .

Power spectral density
The power spectral density (PSD), P x (k), of a wide sense stationary signal x(n) is the Fourier transform (FT ) of its autocorrelation function, r xx (m) [20]: Thanks to the whiteness of the noise, the autocorrelation function of w(n) becomes where δ(m) is Kronecker's function.Based on Eq. ( 7), the PSD of the white noise, P w (k), becomes real, even and constant over the frequency band [ −B; B], with an amplitude σ 2 w .It is obvious that the cumulative sum of the PSD becomes a straight line, with a slope σ 2 w .According to the model (1), s(n) is a cyclostationary signal characterized by its cyclic spectral density (CSD), S α s (k) [13,18]: where α is a cyclic frequency of s(n), S N (k) = FT {s(n)} is the FT of N received samples of s(n) 3 and R(α, m) is the cyclic autocorrelation function of s(n) and can be defined as [18,19]: According to Eq. ( 10), the PSD, P s (k), of s(n) can be evaluated for null cyclic frequency (i.e., α = 0), Due to the fact that the PU signal is not an i.i.d.signal (i.e.R(0, m) is not a Kronecker's function), then its PSD, P s (k), should be not constant on [ −B; B].Based on this fact, the distinguish between H 0 and H 1 can be realized using the shape of the cumulative sum of PSD.

Estimation of the power spectral density
The PSD of the received signal y(n) can be estimated by its periodogram as follows [20]: where Y (k) is the discrete Fourier transform (DFT) of the signal y(n) with N samples: Therefore, the estimated PSD of the signal is related to the modulus of its DFT.Since w(n) is a circular symmetric process, then its DFT W (k) becomes also a zero mean circular symmetric process [21][22][23].

Cumulative power spectral density-based detector
The cumulative power spectral density (CPSD), ψ y (k), of the received signal y(n) is defined, over a frequency interval I =[ v; l] ∀k ∈ I as follows : where Py (u) is the estimated PSD of y(n) using Eq.(11).Unlike ED, which makes the sum of energy on all frequency points of PSD, CPSD makes the sum of the energy on a given frequency points interval.The sense of variation of CPSD will be tested in order to make a decision on the PU signal existence as detailed below.
The expected value of ψ w (k), v ≤ k ≤ l, of w(n) can be found as follows: The last value shows that ψ w (k) increases linearly on the frequency interval [ v, l] with k till its highest value Due to this normalization, under H 0 , (k) still increase to one.Thus, under H 0 and thanks to the flat PSD of w(n), the shape of (k) becomes closed to a straight line Under H 1 , this constraint is not satisfied, since the PSD of s(n) is not a constant, and (k) has higher values than the one obtained under H 0 due to the additional power of s(n).
Figures 1 and 2 show (k) under H 0 and H 1 for rectangular and raise-cosine (roll-off factor = 0.5) pulse shaping filters respectively, for various values of N s and different SNR.The signal modulation is 16-QAM, and the number of used samples is N = 10000.As shown in Figs. 1 and  2, the gap between the normalized CPSD shape under H 1 and that under H 0 increases with the SNR for both considered pulse shaping filters.In addition, the non-linearity of the CPSD shapes grows with N s .Therefore, we define a test statistic T as the difference between (k) and the reference straight line R(v, l; k).Accordingly, we introduce two detectors: , in this case, the corresponding CPSD, ψ p (k), is defined as follows: Accordingly, the normalized CPSD, p (k) becomes and the reference straight line R 1, N 2 ; k is obtained by T p detector aims at finding the difference between p (k) and the corresponding reference shape R 1, N 2 ; k : 2. T a : This detector is based on the CPSD of all frequencies of y(n Where a (k) can be found as follows:

Proposed combining detectors
In this section, two combining detectors are proposed.The first proposed detector, T or , aims at exploiting all the frequency points of the signal y(n), by applying two test statistics: the first one tests the shape of the CPSD for positive frequency points (which is T p ), and T n tests the CPSD shape of the CPSD of the symmetric of P y (k) part standing for the negative frequency points (i.e., where n (k) in this case is given by Since P s (k) is symmetric and deterministic (as s(n) is deterministic) and the components of W (k) are i.i.d.(as it is shown in Section 4.1.1);therefore T p and T n become independent and have same mean and variance.
Once T p and T n make their own decisions, T or , acting as hard cooperative detector, processes those two decisions using an OR rule.
where D T p and D T n are the detection results of T p and T n respectively.
The second proposed cooperative detector, T av , performs the average, P av , between the positive frequency PSD and the symmetric of the negative frequency PSD of the received signal.CPSD is then performed on P av .The averaging process smooth the PSD, since P s (k) is symmetric, and the components of Pw (k) are independent as shown in Section 4.
T av can be considered as a soft combining detector of T p and T n .
where the first term of the previous equation becomes

Statistical analysis
Distributions of test statistics are essential in order to find analytically the probability of false alarm p fa and the detection probability p d .T p and T a have the same statistical distribution since W (k) is i.i.d. and S(k) is deterministic as s(n) is assumed to be deterministic.In the following, we develop the distribution of T p under H 0 and H 1 over a Gaussian channel, where the channel effect h is assumed to be constant.Similarly, the distribution of T a can be found.T or made its decision by applying the logical OR on the decisions of T p and T n which have the same distribution.

False alarm and detection probabilities of T p
The distribution of T p depends on k=1 p (k) as presented in Eq. (20).A simplification of the term N 2 k=1 p (k) can be obtained as follows (see Appendix 1):

False alarm probability of T p
Under H 0 , the test statistic T p is only related to the noise w(n).Using Eq. (29), T p can be written as follows: Being the discrete Fourier transform of a white noise w(n), W (k) asymptotically follows a normal distribution since it is the sum of independent terms.It is known that two Gaussian variables are independent if they are uncorrelated [22].As Being the sum of independent terms and according to the central limit theorem, the distribution of T p tends towards N (μ 0 , V 0 ) under H 0 .In this case, the probability of false alarm p p fa of T p can be found as follows: where Q(.) is the Q-function 4 , and λ is the threshold of comparison.Since E |W (k)| 2 = Nσ 2 w and based on Eq. (30), we can evaluate μ 0 as follows: In this case, the variance, V 0 , of T p becomes (see Appendix 4)

Probability of detection of T p
Under H 1 , Y (k) = hS(k) + W (k), then Eq. ( 29) becomes as follows: Since S(k) is deterministic and the terms of W (k) are independent, the distribution of T p under H 1 tends also towards N (μ 1 , V 1 ).In this case, the probability of detection p p d , of T p can be found as follows: μ 1 and V 1 should be evaluated in order to find where hS(k) = DFT{hs(n)} and Re{X} is the real part of X.The mean value of T p under H 1 can be found as follows (see Appendix 3): , and γ is the SNR as defined by Eq. ( 4).
Under H 1 , the variance V 1 is given by the following equation (see Appendix 5): where c = 8 T a is based on similar idea to T p , but it covers the N frequency points instead of just positive frequency points ( N 2 points).Since W (k) is i.i.d. and S(k) is deterministic, then by following the same process for p p fa and p p d , the probability of false alarm p a fa and the probability of detection p a d of the detector T a can be found as follows: where

Probabilities of T or
T or applies the OR rule between the decisions of T p and T n , then T or can be considered as a hard cooperative detector of these two detectors.Since T p and T n are independent and have the same statistics as we defined previously, the probability of false alarm p or fa and the probability of detection p or d of T or can be found as follows [2]:

Probabilities of T av
T av can be developed following similar steps as Eq.(30).
Under H 0 , Y (k) = W (k), then T av becomes the sum of independent terms.Based on CLT, T av asymptotically follows N (μ av 0 , The probability of false alarm, p av fa , and detection, p av d , of T av are expressed as follows: Since W (k) is i.i.d. and P s (k) is even, μ av 0 , V av 0 , μ av 1 , V av 1 can be found by following similar steps to μ 0 , V 0 , μ 1 and V 1 , we can find that μ av 0 = μ 0 , V av 0 = V 0 2 , μ av 1 = μ 1 , and The theoretical and the simulated ROC curves of proposed detectors are with good agreement as shown in Figs.3a, b.Simulations were done under following conditions: 16-QAM modulation, γ = −12 dB, N = 1000 samples, and N s = 4 sps.
As shown in Fig. 3a, T av is the most efficient detector for both considered shaping filters.For the simulations of Section 6 under Gaussian channel, only T av and T or are compared to other well-known detectors.The rest of simulations in this paper are done with a rectangular pulse-shaping filter.

Probability of detection over Rayleigh fading channel
In this section, we derive the detection probability over the Rayleigh flat-fading channel.The false alarm probability remains the same since it is independent of the channel gain h.The distribution of the SNR, γ , in a Rayleigh channel is given by [21]: where γ is the average SNR.Over a Rayleigh channel, the probability of detection, p dr , is found by averaging the probability of detection, p dg under Gaussian channel with respect to f γ (γ ).
Concerning our detectors, T p , T a , and T av , their detection probabilities over Rayleigh fading channel can be derived similarly, since they have a similar probability of detection: The above integral does not have an analytic solution.Therefore, Taylor series of the first order are used to approximate the argument of the Q-function, g(γ ) = λ−bγ √ V 0 +cγ , around γ 0 = λ/b as follows (see Appendix 6): where g 1 (γ ) is the first-order Taylor series approximation of g(γ ).The approximation of p p dr , pp dr is given by : . According to [24], we can find Using Eq. ( 51) and the fact that , the integral of Eq. (50) becomes For the approximation pa dr , of p a dr , it can be found by replacing the number of samples N/2 used in T p by N used in T a .Similarly to pp dr ,the approximation pav dr of p av dr can be derived as follows: where γ 0 is not modified in the expression of pav dr , since T p and T av have the same mean under H 1 .
The detection probability of T or under Gaussian channel is a non-linear combination of the probability of detection of T p and T n .Over a Rayleigh channel, the fading coefficient, h, is the same for T p and T n , p or dr of T or becomes As there is no analytic solution of Eq. ( 54), we solve it in numerical way.In Fig. 4, analytical results are very closed to simulation results.T av and T or both achieve the best performance among the proposed detectors, thus they will be compared to the ED and CSD under Rayleigh channel.

Performance evaluation
In this section, we compare the performance of our detectors to that of the energy detector (ED) [8] and the cyclostationary detector (CSD) [10], under Gaussian and Rayleigh channels.Throughout the upcoming simulations, A 16-QAM baseband-modulated PU signal with a rectangular pulse-shaping filter is considered.

Performance analysis over Gaussian channel
Figure 5 presents the ROC curves under a Gaussian channel of various numbers of samples per symbol N s .Simulations are done using N = 1500 samples and γ = −12 dB.According to that figure, the performance increases with N s .For the limiting case, i.e., when N s = 2 sps, CSD shows a poor performance relatively to T av , T or , and ED, while T av outperforms all other detectors.However, T av and T or outperform ED and CSD for different values of N s .
Figure 6 shows the variation of the probability of detection with respect to SNR for a constant p fa = 0.1.The number of samples is fixed N = 1000 samples, and various values are assigned for N s .T av and T or reach higher probabilities of detection than ED and CSD for similar SNR and different values of N s .In addition, increasing N s leads to enhancing the performance of T or , T av , and CSD.For example, T av reaches p d = 0.9 at SNR −7 dB for N s = 2 sps, while the same probability of detection is reached for SNR −9 dB and SNR −10 dB at N s = 4 sps and N s = 8 sps respectively.

Performance analysis over Rayleigh channel
Figure 7 shows the ROC curves under Raleigh fading channel for N s = 4 and N s = 8 sps, with N = 1000 samples and an average SNR of −5 dB.
Over Rayleigh fading channel, T av and T or are still outperforming ED and CSD.The same fading suffered by the negative and the positive frequency parts of PSD affects the performance of our detectors.This fact makes the performances of ED and CSD closed to the proposed detectors, as shown in Fig. 7, where the gap of performance among our detectors, ED and CSD becomes smaller comparing it to a Gaussian channel.

Complexity analysis
To obtain the decision of T or , we have to previously find T p and T n .Since those two test statistics have the same complexity, therefore, the complexity of T or becomes the double of that of T p except for the calculus of the FFT which should be evaluated once.Then the required number of operations of T or , C(T or ), becomes For T av , we can find the corresponding number of operations, C(T or ), similarly to T or : To compare the complexity of our detectors to that of ED and CSD, we should emphasize that ED needs N operations of multiplication and N − 1 for addition: For CSD, the corresponding number of operations, C(CSD), can be derived as follows [25]: where L is an odd number and stands for the length for the unit window used in the detection process of CSD [10].
According to Eqs. ( 56) and ( 57), the complexity T or and T av is independent of the number of samples per symbol N s , contrary to CSD where the complexity depends on N s as shown in Eq. ( 59).As shown previously, increasing the oversampling rate leads to obtain a robust performance for CSD, T or and T av , while the performance of ED is not affected by the oversampling rate.
Figure 8 shows the number of samples and the complexity of the various detectors in terms SNR for a target (p fa ; p d ) = (0.1; 0.9) under a Gaussian channel.As shown in Fig. 8a, T or and T av need a number of samples less than that of ED and CSD, in order to reach the target probabilities.In addition, the required number of samples decreases with an increasing N s .Figure 8b shows the number of performed operations corresponding to the number of required samples given in Fig. 8a.The complexity of T or and T av decreases if N s increases due to the fact that the number of required samples decreases with an increasing N s .Contrariwise, the complexity of CSD grows with N s even if the total number of samples decreases, this is because the complexity of CSD depends on N s .On the other side, our proposed detectors are slightly more complicated than ED.
To summarize, ED is less complicated than the proposed detectors.However, T av and T or require shorter observation time to fulfill the equivalent performance of ED and CSD.

Robustness of our proposed detectors under noise uncertainty
Due to many factors (such as thermal noise, ambient interference, receiver non-linearity, etc.), the variance of the white noise cannot be estimated accurately.This fact introduces the noise uncertainty (NU) problem which prevents the detector from reaching a target (p fa , p d ) even with a large observation time.According to Eq. ( 15), the noise variance should be pre-estimated in order to perform the normalization.That means our proposed detectors are sensitive to the estimation of the noise variance.In this section, the impact of the NU on the robustness of our proposed detectors is evaluated.The estimated noise variance σ 2 w can be bounded as follows: w is the nominal value of the noise variance and r ≥ 1 stands for the NU factor.The distribution of σ 2 w , f σ 2 w ( σ 2 w ), is assumed to be uniform in a logarithmic scale red [26].
In a conventional detection mechanism, a suitable threshold λ is fixed according to a target p fa .In our proposed detectors, the choice of λ depends on setting the estimated value of the noise variance and the observation time (i.e. the number of the received samples).A wrong estimation of σ 2 w may lead to an inappropriate λ.This fact increases the false alarm rate and deteriorates the spectrum sensing performance.
In Fig. 9, the variation of p fa with ρ is illustrated for the proposed detectors and the conventional ED.All considered detectors are assumed to have p fa = 0.1 when no NU exists (i.e.ρ = 0), then the threshold of each detector is chosen according to p fa = 0 with a perfect estimation of the noise variance.Figure 9 shows that the p fa of ED ia higher than the ones of the proposed detectors.Indeed, our detectors exploit both the energy and the CPSD shape of the received signal unlike ED which exploits only the energy.On the other hand, T av and T a have the same p fa as these two detectors are based on linear combination of |W (k)| 2 which is i.i.d.. T or has a superior p fa than T p since these two detectors are related to each other by a non-linear Eq. (41).
To show the effect of the noise uncertainty on the ROC curves, we evaluate numerically the performance loss p d = p d (0) − p d (ρ) where p d (0) stands for the case where there is no NU, and p d (ρ) stands for probability of detection for the case where NU is equal to ρ dB. Figure (10) presents p d of our detectors and ED for p fa = 0.1, N=1000 samples, SNR of -10 dB, and N s = 3 sps.As shown in this figure, our detectors are less sensitive to the noise uncertainty than the energy detector for p fa = 0.1.We observe that T a is most affected among the proposed detectors.T av and T or have the same detection loss since they both are related to both T p and T n , which exhibit the same NU.

Spectrum sensing based on self-normalized CPSD:
As discussed previously, the normalization by the noise variance (Eq.( 15)) leads to a possible NU problem.In order to avoid this problem, a new model for the proposed detectors is introduced.Instead of normalizing the CPSD by the mean value of the last term of ψ w (k), the CPSD is normalized using the last term.The selfnormalized CPSD of the received signal, y(n), is defined as follows: becomes the estimated energy of the received signal.Such a normalization makes η y (k) not vulnerable to the NU as its calculus does not depend on the estimation of the noise variance (as presented in Eq. ( 62)), and its distribution and its statistical parameters (mean and variance) are also independent of the noise variance (see Appendix 7).
This normalization is equivalent to a scaling and does not change the shape form of the CPSD.
Figure 11 shows η y (k) under H 0 and H 1 and a SNR = 0 dB for various values of N s .We have under H 0 a shape like a straight line, but a curved shape under H 1 .As shown in this figure, the difference between η y (k) and the straight line increases as N s increases, which means the detection becomes more reliable with the increasing of N s .
Without taking into account the relative position of η y (k) with respect to the reference straight line, the decision about the presence of the PU signal is made by comparing the η y (k) shape to the reference line shape.For that reason, the dissimilarity between the CPSD and the reference straight line is computed (see Eqs. ( 63)-( 66)).Similarly to the detectors T p , T a , T or , and T av , we define the new detector model based on the selfnormalization detectors: T s p , T s a , T s or , and T s av as follows respectively: where η av y (ν) is given by

Conclusions
In this paper, we proposed spectrum sensing detectors based on the cumulative power spectral density (CPSD).
Our proposed detectors verify the linearity of the CPSD shape of the received signal.
Hard and soft schemes are used to combine the CPSD measures, which are derived based on the two symmetric parts of the power spectral density.False alarm and detection probabilities are derived analytically under both Gaussian and Rayleigh flat-fading channels.Our simulation results show the performance superiority of our detectors comparing to classic detectors such as the energy and the cyclostationary detectors.
In addition, simulation results show that the proposed detectors are less affected by the noise uncertainty than the energy detector.However, to avoid the impact of the noise uncertainty, the measured CPSD is normalized by the estimated energy of the received signal.By doing this, we make our detectors independent from the noise variance.
In a future work, we will enhance the power spectral density estimator and extend our algorithms to deal with multiple antennas spectrum sensing problem.

Endnotes
1 This work was presented in part as a book chapter in Springer book [27] 2 According to [18] and [19], almost-periodic signals can be considered as deterministic in a fraction of time (i.e., for our application, this is the time of sensing).As the linear modulated signal s(n) is almost cyclostationary signal according to [19], then it is almost-periodic [18,19] and consequently it can be considered as deterministic during a period of spectrum sensing.On the other hand, in [28] the author motivates the use of the deterministic assumption of the cyclostationary signals in a fraction of time while the purely stationary noise signals (do not have any cyclic features) are considered random. 3for simplicity, we refer to S N (k) by only S(k) in the forthcoming sections 2 dt [22] Appendix 1 Simplification of Eq. ( 20) . . .

Autocorrelation of the DFT of W(k)
Let us consider the autocorrelation function of W (k), r WW (m).
The mean of T p under H 1

Variance of T p under H 0
As μ 0 = E[ T p ] = 0, using Eq. ( 30), the variance V 0 can be written as follows: The term 4 be simplified as follows: Since W (k) is Gaussian and i.i.d.then: 1) The kurtosis of W (k) is zero: Using Eqs. ( 73) and (74), the Eq.(72) becomes Back to Eq. (71), V 0 becomes

Variance of T p under H 1
The calculation of Eq. (77) in the next page stands for finding the variance, In Eq. (77), the part A 2 = 0 because +h p S q (k)W q (k) + h q S p (k)W q (k) − h q S q (k)W q p(k) Since W p (k) and W q (k) are independent and E W p (k 1 )W 2 p (k 2 ) = E W 2 q (k 1 )W q (k 2 ) = 0 ∀k 1 and k 2 , since W q (k) and W p (k) are Gaussian, then Eq. (78) becomes zeros.
Using the i.i.d. and the circular properties of W (k) and the fact that s(n) is deterministic, Eq. (77) becomes

Approximation of the detection probability
The approximation of the Q-function using a first-order Taylor series was proposed in [16] without justification.
Here, we show by simulation the effectiveness of this approximation.The Taylor series of a function f (t) around t 0 can be developed as follows: where f (n) (t 0 ) is the nth order derivative of f (t) at t 0 .According to Eq. ( 48), we aim at simplifying in order to find the analytic probability of detection under Rayleigh channel.Let us define g(γ ) = λ−bγ √ V 0 +cγ .The first g 1 (γ ) and the second order, g 2 (γ ), Taylor series approximations of g(γ ) around γ 0 can be found as follows: Let γ 0 = λ/b as g(γ 0 ) = 0, and then Q(g(γ 0 )) = 0.5, which is the middle point of the Q-function.Accordingly, g (γ 0 ) and g (γ 0 ) can derived as follows: ) and its approximations Q(g 1 (γ )) and Q(g 2 (γ )) under different conditions.Figure 14a shows the variation of p p d and its approximations in terms of SNR for different p fa .The number of samples is fixed to 1500 and N s = 4 sps, while Fig. 14b shows under SNR of 10 dB and N s = 4 sps, the variation of p p d and its approximations in terms of N for different p fa .As shown, the analytic and the approximated curves are closed to each others under the various conditions.As expected, g 2 (γ ) leads to a more robust approximation, since Q(g 2 (γ )) is almost colinear with Q(g(γ )).Even though, g 1 (γ ) results are very closed to exact ones.As no important loss is obtained when g 1 (γ ) is used, and since it is more simple to deal with the first-order Taylor series, g 1 (γ ) will be considered to approximate the detection probability under Rayleigh fading channel.

1 .Fig. 1 Fig. 2
Fig. 1 The normalized CPSD shapes for N = 10000 samples and several values of N s and SNR: rectangular pulse-shaping filter.a N s = 2 sps.b N s = 4 sps.c N s = 8 sps

Fig. 3
Fig. 3 Simulated vs. analytical results over Gaussian channel for rectangular and raised-cosine pulse-shaping filters.a Rectangular pulse-shaping filter.b Raised-cosine pulse-shaping filter: roll-off factor = 0.5

Fig. 5 Fig. 6 Fig. 7
Fig. 5 ROC curves of our proposed detectors comparing to ED and CSD for various N s , SNR of -12 dB and N = 1500 samples.a N s = 2 sps.b N s = 4 sps.c N s = 8 sps

Fig. 8 a
Fig. 8 aThe number of samples needed by the detectors in order to reach (p fa ; p d ) = (0.1; 0.9) for various SNR and b the number of required operations performed by each detectors corresponding to the number of samples given in a

Fig. 9 Fig. 10
Fig.9 The evolution of p fa in terms of the NU ρ

Figure 12
Figure12shows the performance of proposed detectors under a NU of 0, 0.75, and 1.5 dB with N s = 4 sps, SNR of -10 dB, and N = 1000 samples under Gaussian and Rayleigh flat-fading channels.For a NU = 0 dB, T av and T or outperform T s av and T s or respectively for both Gaussian and Rayleigh channels.Beside that, when the NU grows, the performance of T s av and T s or is not affected, whereas the detectors T av , T or , and ED suffer a performance degradation and become less robust than the self-normalization detectors, as shown in Figs.12b, c and 13b, c.However, T av and T or have a superior performance relative to ED for different values of ρ.

bFig. 14 a
Fig. 14 a The probability of detection p p d and its approximation in terms of SNR for different p fa , N = 1500 samples and N s = 4 sps.b The probability of detection p p d and its approximation in terms of the number of samples N, for different p fa , SNR=−10 dB samples, and N s = 4 sps where δ, V H 0 , and β are constants.Hereinafter, we only derive p