Double-Threshold Cooperative Spectrum Sensing Algorithm Based on Sevcik Fractal Dimension

Spectrum sensing is of great importance in the cognitive radio (CR) networks. Compared with individual spectrum sensing, cooperative spectrum sensing (CSS) has been shown to greatly improve the accuracy of the detection. However, the existing CSS algorithms are sensitive to noise uncertainty and are inaccurate in low signal-to-noise ratio (SNR) detection. To address this, we propose a double-threshold CSS algorithm based on Sevcik fractal dimension (SFD) in this paper. The main idea of the presented scheme is to sense the presence of primary users in the local spectrum sensing by analyzing different characteristics of the SFD between signals and noise. Considering the stochastic fluctuation characteristic of the noise SFD in a certain range, we adopt the double-threshold method in the multi-cognitive user CSS so as to improve the detection accuracy, where thresholds are set according to the maximum and minimum values of the noise SFD. After obtaining the detection results, the cognitive user sends local detection results to the fusion center for reliability fusion. Simulation results demonstrate that the proposed method is insensitive to noise uncertainty. Simulations also show that the algorithm presented in this paper can achieve high detection performance at the low SNR region.


Introduction
With the rapidly increasing data rate requirements of mobile users on multimedia applications, the accessible radio spectrum is becoming critically scarce.Under traditionally fixed spectrum allocation policy, some licensed bands are overcrowded, while many others are underutilized.Cognitive radio (CR) arises as a tempting solution to the spectrum congestion problem by enabling secondary users (cognitive users or SUs) to have opportunistic access to underutilized licensed bands that are lightly occupied by primary users (non-cognitive users or PUs) [1].It is crucial to sense the spectrum holes efficiently and determine the presence of primary users correctly in a CR network.Therefore, spectrum sensing is considered as the premise and foundation of cognitive radio systems [2].
Energy detection [3] possesses low computational costs and is easily implemented.However, it is sensitive to noise uncertainty and performs poorly in the low signal-to-noise ratio (SNR) region.In [4], a low-complexity heuristic decision-making system for the selection of methods was proposed aimed at detecting spectrum availability.An optimization approach for multi-band linear collaborative spectrum sensing based on harmony search heuristics was provided by [5].The problem of maximizing the ratio of the transmission duration over the entire sensing cycle was studied in [6]. Lee et al. proposed an algorithm for the detection of unknown low-power constant-modulus signals in noise of unknown variance for cognitive radio networks [7].An effective method for application to cooperative sensing for positioning in the presence of correlated observations, was provided by [8].The test in [8] combines the OR and the majority voting rules at once, recovering all those unlucky detection cases that are not sensed due to strong shadowing.The authors in [9] provided a fair comparison between the soft and hard schemes in terms of energy efficiency.The cyclostationary feature detection method [10] detects the signal in terms of the characteristic of the cyclic frequency in the signal spectral correlation function.Therefore, it needs a larger amount of computation due to the fact that it is subject to the signal spectral correlation function.
Recently, fractal geometry has been widely used in signal recognition, digital image processing, and other fields.A spectrum sensing algorithm based on fractal dimension obtains the detection information of spectrum sensing without the prior information of the primary user by using the frequency domain fractal dimension of the signal.Furthermore, it has strong anti-noise ability and low computational complexity [11,12].Owing to the above advantages, spectrum sensing algorithms based on fractal dimension has become a research focus.The authors in [13] studied the impact of the range of fractal values, ambient noise levels, and the length of the observation window on the detection result, and concluded that the environment and conditions should be considered when using the fractal dimension method.A spectrum sensing algorithm based on two-dimensional detection was proposed in [14].Though its computational complexity was low, the detection effect in low SNR environments could not meet expectations.A dynamic cooperative spectrum sensing (CSS) based on box dimension was presented by [15].However, this method had poor detection accuracy in the low SNR region due to limitations of box dimension.In [16], a CSS method based on the correlation box dimension was proposed, which eliminated the effect of noise by using the signal autocorrelation operation in the local detection and thus achieved better detection effect.Nevertheless, the autocorrelation operation increased its complexity.The cognitive ultra-wideband (CWUB) spectrum sensing was studied in [17], where the fractal box dimension was used for double-threshold CSS and the information dimension was used to identify the signal.However, it is difficult to identify certain signals by information dimension, and the use of two fractal dimensions increases the computational complexity.The Sevcik fractal dimension (SFD) was used in [18], with which the accuracy of local spectrum sensing was improved and the drawback that the box dimension cannot sense the signal of some modulation modes was overcome.
Noise power changes with time and physical position in the practical communication environment; moreover, there exists a certain deviation for receivers to estimate noise power.Thus, noise uncertainty can seriously affect the accuracy of spectrum sensing.The proposed algorithm senses the presence of primary users based on analyzing different characteristics of the SFD between signals and noise.The proposed scheme first transforms the received PU's signal into a frequency signal by discrete Fourier transform (DFT) and then calculates the SFD of the frequency signal.The judgment is conducted by comparing the above SFD value with the threshold, where the noise SFD is selected as the decision threshold.Since the noise SFD changes randomly in a certain range, the single threshold method cannot meet the requirement of detection precision when the SNR is low or the noise is uncertain.To address this, we resort to the double-threshold method in local sensing.Furthermore, multi-user CSS method is adopted to promote the reliability of sensing results and the detection accuracy.The simulation results demonstrate that the proposed method achieves better detection performance at low SNR.Simulation results also show that the scheme presented in this paper has the advantage of insensitivity to noise uncertainty and is less affected by signal modulation parameters.
The rest of this paper is organized as follows.Section 2 introduces the double-threshold cooperative spectrum sensing algorithm based on Sevcik fractal dimension.Simulation results are given in Section 3. The paper concludes in Section 4.

The Noise Sevcik Fractal Dimension
The SFD of the modulated signal and that of noise have different randomness, so the spectrum sensing model can be constructed using this difference.In the following, we first discuss the calculation of SFD.
Assume that the signal waveform whose length is N consists of the point (x i , y i ), 0 ≤ i ≤ N − 1.We first normalize the waveform. x where x max is the maximum value of x i , and y max and y min respectively denote the maximum and minimum value of y i .The total length of the curve is calculated as The communication signal sequence is a continuous positive integer and 1) and ( 3) are equivalent to The SFD of the signal is The SFDs of different kinds of signals are simulated in Figure 1, where the signals include: 2 amplitude shift keying (2ASK), 2 frequency shift keying (2FSK), 2 phase shift keying (2PSK), 4ASK, 4FSK, quadrature phase shift keying (QPSK) and white Gaussian noise signal.The simulation parameters are set as follows: bandwidth B = 20 kHz, carrier frequency f c = 10 kHz, sampling frequency f s = 40 kHz, baseband symbol rate R b = 1000 bit/s, the frequency deviation of the 2FSK signal is 5000 Hz, and the frequency deviation of the 4FSK signal is 2500 Hz.The noise is additive white Gaussian noise (AWGN).Observing Figure 1, the noise SFD fluctuates randomly between 1.79 and 1.83 with the SNR increasing, whereas the SFDs of the six kinds of modulation signals are decreasing with SNR.From the figure, the SFDs of the modulation signals are gradually separated from the noise SFD for SNRs greater than −25 dB.It can be inferred that the noise SFD can be set as the threshold when using SFD in spectrum sensing, and thus the signal can be distinguished from noise when the SNR is greater than a certain value.
The fluctuation of noise power could affect the performance of spectrum sensing, and thus result in the detection probability decreasing or the false alarm probability increasing.In order to further prove that the calculation of SFD is insensitive to Gaussian noise, we analyze the influence of noise power on SFD.
The received signal follows the Gaussian distribution y(n) ∼ N(0, σ 2 ) when there is no primary user.Applying N-point discrete Fourier transform (DFT) to y(n): Since Ỹ(k) is a linear combination of y(n), it also follows the Gaussian distribution.The mean and variance of Ỹ(k) are expressed as where E[•] and D[•] respectively denote mean and variance.Assume Ỹ(k) = Y r (k) + jY i (k), both the real and imaginary parts of the received signal Ỹ(k) follow the Gaussian distribution N (0, Nσ 2 2 ).Let Y(k) denote the modulus of Ỹ(k), which is calculated by Thus, We perform normalization to Y(k) according to (1) and ( 2).Y max is regarded as a constant.Therefore, the mean and variance of Y(k) after normalization are expressed as According to previous analysis, because x * i+1 − x * i = 1/N − 1 , we only discuss the expectations of (y * i+1 − y * i ) 2 in Formula (3).
Y max has been regarded as a constant so far.However, Y max is related to σ 1 under the specific distribution: where C is a constant which is calculated as follows: where η = Y max /σ 1 .According to ( 5), ( 6), (16), and ( 17), we have 4 As seen from ( 19), E[D] has no correlation with noise power σ 2 0 .Thus, SFD is insensitive to noise uncertainty.
Figure 2 calculates the SFD of white Gaussian noise versus different power.Sequence length N = 4000.From the figure, the mean change range of Gaussian white noise SFD is smaller, which always fluctuates near 1.811 when noise power changes within the range −30-20,dBm.Simulation results validate that the noise SFD values have no correlation with noise power, which means that SFD can overcome the influence of noise uncertainty to a certain extent.

System Model of the Proposed Algorithm
Local spectrum sensing on the basis of SFD can be formulated as the following binary hypothesis problem: where ω(n) is the independently and identically distributed AWGN with mean 0 and variance ω 2 .y(n) is the received signal of the cognitive user.Hypothesis H 0 and H 1 represent the absence and presence of the PU, respectively.Primary signal s(n) can either be a deterministic signal (accounting for AWGN channel) or a stochastic signal (corresponding to channel characteristics like fading) with mean µ and variance σ s 2 .s(n) and ω(n) are mutually independent.y(n) follows the following distribution: The SFD of the signal fluctuates little with the SNR, and it can be distinguished from the SFD of Gaussian white noise.Therefore, SFD can be used to distinguish the signal and noise.The block diagram of spectrum sensing system based on SFD in the frequency domain is shown in Figure 3. Firstly, the received primary user signal is transformed into the frequency domain by DFT, and then SFD is calculated.In the following, we illustrate the spectrum sensing scheme based on SFD in detail.Firstly, applying DFT to (20), we have the frequency domain form of the received signal: Y(k), W(k), and S(k) respectively represent the complex spectrum of y(n), ω(n), and s(n).We perform normalization to Y(k).
Since Y(k) is the DFT transform of y(n), we set Y min = 0. Equation ( 23) can be simplified into: Calculate the SFD of the signal by using (6), and judge whether the PU exists or not through (25).
Equation ( 25) is the local spectrum sensing decision expression based on SFD.If the SFD of the received signal calculated by the SU is less at than the decision threshold λ, the judgment that is the PU exists; otherwise, the judgment is that the PU does not exist.
As mentioned earlier, we can select a typical noise SFD value as the decision threshold λ.However, the noise SFD randomly changes in a certain range, leading to better detection performance at high SNR and poor performance when the SNR is low or noise is uncertain.Therefore, the drawback would be revealed if we use the single threshold method shown in (25): when the SFD of the signal falls within the fluctuation range of the noise SFD, the signal and noise would be easily confused, leading to the inaccurate judgement of the PU's presence, and thus increasing the false alarm probability or even misjudgment.The single threshold cannot meet the requirement of detection precision due to channel environment change and noise uncertainty in practical application.Therefore, based on the existing double-threshold spectrum sensing theory, we propose a double-threshold CSS strategy based on SFD.
Figure 4 shows the double-threshold decision.We select the maximum and minimum values of the noise SFD as the two thresholds, where λ 1 = min{D}, λ 2 = max{D}, and λ 1 < λ 2 .Considering that noise power in the wireless environment has a random change, we resort to the dynamic detection method in each sensing cycle, where the best decision thresholds are decided by the real-time calculation of the noise SFD.
this region is called "decision buffer" and the decision is the uncertain state Ω due to the fact that the SFD of the received signal falls within the random fluctuation range of the noise SFD.
SUs participating in the CSS send the local decision results H 1 and H 0 as hard decision information to the fusion center (FC).If the decision is uncertain state, the test statistic D is sent as the soft decision information to the FC for decision.
Assume the FC receives Z local decision results sent by SUs, in which the number of the hard decision information is K and correspondingly the number of the soft decision information is Z − K.
The fusion rule for the soft decision information at the FC is shown as follows: where α i represents one of the soft decision information and Γ is the threshold for the reliability fusion at the FC.In practice, Γ can be calculated by where λ = (λ 1 + λ 2 )/2 , ∆ is the correction, ∆ ∈ (−1, 1).Since the two thresholds are set in terms of the maximum and minimum values of the noise SFD, their range of variation is related to the current spectral environment and the noise background, and thus ∆ needs to be adjusted according to the actual situation.
The OR rule is adopted as the final decision so as to maximize the global detection probability.The FC makes the fusion decision based on the hard decision information and the soft decision information.The final decision is expressed as: where Λ is the final decision result at the FC.

Simulation
Simulation experiments are under AWGN environment without channel fading effects.The signal of the PU is the modulated QPSK signal.Signal parameters are set the same as that of Section 2.1.Bandwidth B = 20 kHz and carrier frequency f c = 10 kHz.The sampling frequency at the cognitive receiver is set as f s = 40 kHz.Baseband symbol rate R b = 1000 bit/s.Assume that the number of collaborative SUs is 10.Simulation results are averaged over 1000 channels.
In Figure 5, we compare the detection probability of the proposed method with the Katz fractal dimension scheme (Katz), the box dimension method (Box dimension) and the energy detection algorithm.According to Figure 2, the noise fractal dimension is basically maintained, fluctuate around a mean value under different noise powers.Hence, the decision threshold can be determined by the value of the noise fractal dimension, and correspondingly the false alarm probability is a constant which is set as P f = 0.03.The values of the threshold and the correction in different schemes are shown in Table 1.Observing Figure 5, the detection probabilities of the three fractal dimension algorithms all outperform energy detection.The proposed scheme achieves 100% detection probability P d at −20 dB, whereas Katz fractal dimension and box dimension attains 100% detection probability at −17 dB and −13 dB, respectively.It can be inferred that the double-threshold CSS algorithm based on SFD improves the detection performance in low SNR environment.
Table 1.The values of the threshold and the correction in different schemes.

Sevcik
Figure 6 investigates the influence of the number of SUs on the detection performance when applying the proposed algorithm.We consider M = 1, 5, 10 and let the average SNR at each secondary receiver be the same.As seen from the figure, the detection performance of the proposed scheme gets better when the number of cognitive users involved in collaborative sensing increases.We compare the detection performance of the proposed algorithm with that of the single threshold CSS algorithm based on SFD in Figure 7, where the threshold in the single threshold detection equals to the average noise SFD λ S = 1.765.The double-threshold method achieves 100% detection probability at −20 dB, which is 1 dB higher than that of the single threshold method.It can be demonstrated that the proposed approach benefits from higher detection performance compared with the single threshold CSS scheme based on SFD.In Figure 8, the detection performance of the proposed algorithm is plotted when the noise uncertainty is 1 dB, 2 dB, and 5 dB, respectively.As expected, noise uncertainty had little impact on the detection performance of the proposed method, which means the proposed scheme is insensitive to noise uncertainty.Furthermore, even in the case of 5 dB noise uncertainty, the proposed scheme still obtains 100% detection probability when the SNR is greater than −20 dB, and it can achieve more than 90% detection probability when the SNR is greater than −21 dB.Therefore, the algorithm presented by this paper can not only be insensitive to noise uncertainty, but also achieve high detection performance at low SNR.

Conclusions
According to the different SFD characteristics of the signal and noise, a double-threshold CSS algorithm based on SFD is proposed in this paper.Theoretical analysis and simulation show that the proposed scheme can sense some modulation signals that cannot be sensed by the box dimension; meanwhile, it is insensitive to noise uncertainty.Considering the characteristic of the noise SFD fluctuating within a certain range, the double-threshold is adopted in the multi-user cooperation sensing.Simulation experiments are under AWGN environment without channel fading effects.The results show that the SFD algorithm can achieve high detection probability at low SNR region and the double-threshold cooperative sensing effectively improves the accuracy of CSS.The proposed algorithm achieves a stable detection result in the case of noise uncertainty.In conclusion, the double-threshold CSS algorithm based on SFD is suitable for the fast blind spectrum sensing with low SNR and noise uncertainty, at the same time obtaining better detection results.
k) follows a Rayleigh distribution with parameter σ 2 1 = Nσ 2 0 /2.f Y (Y) and F Y (Y) respectively represent the probability density function and cumulative distribution function.

Figure 2 .
Figure 2. The relationship between noise power and the noise Sevcik fractal dimension.

Figure 3 .
Figure 3.The block diagram of spectrum sensing system based on Sevcik fractal dimension (SFD) in the frequency domain.DFT: discrete Fourier transform.

Figure 5 .
Figure 5. Detection probability comparison of the proposed method with the other three methods.SNR: signal-to-noise ratio.

Figure 6 .
Figure 6.Detection probability comparison of different numbers of secondary users (SUs).

Figure 7 .Figure 8 .
Figure 7.Detection performance comparison of the proposed algorithm with the single threshold detection.