^{1}

^{1}

^{1}

^{1}

^{1}

This paper addresses the problem of the opportunistic spectrum access in Cognitive Radio. Indeed, most spectrum sensing algorithms suffer from a high computational cost to achieve the detection process. They need a prior knowledge of signal characteristics and present a bad performance in low Signal to Noise Ratio (SNR) environment. The choice of the optimal detection threshold is another issue for these spectrum sensing algorithms. To overcome the limits of spectrum detectors, we propose in this paper, a blind detection method based on the cyclostationary features of communication signals. Our detector evaluates the level of hidden periodicity contained in the observed signal to make decision on the state of a bandwidth. In order to reduce the computational cost, we take advantage of the FFT Accumulation Method to estimate the cyclic spectrum of the observed signal. Then, we generate the Cyclic Domain Profile of the cyclic spectrum which allows us to evaluate the level of the hidden periodicity in the signal. This level of periodicity is quantified through the crest factor of Cyclic Domain Profile, which represents the decision statistic of the proposed detector. We have established the analytic expression of the optimal threshold of the detection and the probability of detection to evaluate the performance of the proposed detector. Simulation results show that the proposed detector is able to detect the presence of a communication signal on a bandwidth in a very low SNR scenario.

In this paper, we propose the Cyclostationary Features based Detection (CFD) algorithm based on a spectral analysis of the cyclostationarity of digital communications signals. By definition, a cyclostationary process is a random process in which the statistical characteristics are periodic or quasi-periodic [

• First, we adapt the Fast Fourier transform Accumulation Method (FAM) algorithm to estimate the cyclic spectrum of the observed signal.

• Then, we generate the Cyclic Domain Profile (CDP) which allows us to evaluate the level of the hidden periodicity in the signal. This level of periodicity is quantified through the Crest Factor (CF) of CDP, which represents the decision statistic of CFD.

During the detection phase, the CFD compares the CF of the received signal to a predetermined threshold. The simulations which we carried out through the Receiver Operating Characteristic (ROC) curves show the robustness of our detection method.

The rest of this paper is organized into five sections. In Section 2, we briefly present the properties of cyclostationary processes. Section 4 presents the FAM algorithm used to estimate the cyclic spectrum. The principle of the detector is presented in Section 5. We end this paper with Sections 6 and 7 which respectively present the results of our simulations and the conclusion.

The first works on the cyclostationarity of communication signals were published in 1958 by Bennett [

Today, cyclostationarity is exploited in various fields including telecommunications [

The cyclostationarity is defined as an extension of the concept of stationarity or as a special case of non-stationarity. Such a signal is defined by a periodicity at order n of its statistics [

Let X ( t ) be a real-value continuous-time stochastic process. The process X ( t ) is said to be Cyclo-Stationary (CS) with period T_{0} in the strict sense of the order n if and only if its probabilistic characteristics are T_{0}-periodic [

F X ( t ) ( x ) = F X ( t + k T 0 ) ( x ) (1)

When the process X ( t ) contains several periodicities, the process is named a polycyclostationary process. Thus, X ( t ) is said to be polyclostationary in the strict sense with the periods T 0 , T 1 , ⋯ , T P if and only if its statistical properties are periodic in time with the periods T 0 , T 1 , ⋯ , T P [

Let X ( t ) be a real-value continuous-time stochastic process. The process X ( t ) is said to be CS of period T_{0} in the broad sense if and only if its statistics of order 1 (its mean) and of order 2 (its autocorrelation function) vary periodically over time with a period T_{0}.

In the rest of this paper, we will focus on cyclostationarity in the wide sense to develop our CFD detector which only exploits the statistical properties of second order of the digital communications signal.

The theory of second order cyclostationarity is based on the autocorrelation function of the stochastic process. Considering that in practice, the observed signal is only one realization of the cyclostationary process, we will deal in the continuation with the second order cyclostationarity of digital communications signals.

A centered signal x ( t ) is said to be second-order cyclostationary, in the wide sense, if its autocorrelation function R x ( t , τ ) defined by:

R x ( t , τ ) ≜ E [ x ( t + τ 2 ) x * ( t − τ 2 ) ] (2)

is a T_{0}-periodic function for each τ ; which means in other words that: ∀ k ∈ ℤ * ,

R x ( t + k T 0 , τ ) = R x ( t , τ ) (3)

where τ denotes the time offset (correlation variable). T_{0} is called the cyclic period. It is a period hidden in the signal x ( t ) ; E [ . ] denotes the expectation [

Since the auto-correlation function R x ( t , τ ) is T_{0}-periodic in t for each shift τ , then it can therefore be decomposed into Fourier series in the form:

R x ( t , τ ) = ∑ k = 1 K R x k α 0 ( τ ) ⋅ e j 2 π k α 0 t (4)

where α 0 = 1 / T 0 denotes the fundamental cyclic frequency and K the number of harmonics.

The quantity R x k α 0 ( τ ) is called the cyclic correlation function [

R x k α 0 ( τ ) ≜ 1 T 0 ∫ − T 0 2 T 0 2 R x ( t , τ ) ⋅ e − j 2 π k α 0 t d t (5)

When its autocorrelation function contains several cyclic periods (hidden periodicities) T 0 , T 1 , T 2 , ⋯ , the signal x ( t ) is said to be polycyclostationary [

R x ( t , τ ) ≜ ∑ α ∈ A R x α ( τ ) ⋅ e j 2 π α t (6)

where A = { α 0 , α 1 , ⋯ , α N } is a numbered set containing all the possible cyclic frequencies α and the Fourier coefficients R x α ( τ ) .

R x α ( τ ) ≜ l i m T → ∞ 1 T ∫ − T 2 T 2 R x ( t , τ ) ⋅ e − j 2 π α t d t (7)

denote the cyclic autocorrelation function.

By analogy with the Power Spectral Density, the Fourier transform of the cyclic autocorrelation function is called density function of the cyclic correlation or cyclic spectrum [

S x α ( ν ) ≜ ∫ − ∞ + ∞ R x α ( τ ) ⋅ e − j 2 π ν τ d τ (8)

The estimation of the quantities defined by Equations (7) and (8) constitutes the cyclic spectral analysis.

Let w ( t ) be a sliding window of unit area and not zero over the interval [ − T 2 , T 2 ] . The adequate estimator of the cyclic autocorrelation function R x α ( τ )

of the observation x ( t ) for t ∈ [ t 0 − T / 2 , t 0 + T / 2 ] is the cyclic correlogram R x ( T ) ( α , t ) defined by [

R x ( T ) ( α , τ ) ≜ ∫ ℝ w ( t − t 0 ) x ( t + τ ) x ¯ ( t ) e − j 2 π α t d t (9)

x ¯ represents the conjugate of x ( t ) . The cyclic correlogram R x ( T ) ( α , τ ) is a consistent estimator of R x α ( τ ) [

Assuming w ( t ) = r e c t ( t / T ) / T , we show that the cyclic autocorrelation function is equivalent to [

R x ( T ) ( α , τ ) = 1 T ∫ − T 2 + T 2 x ( t + τ 2 ) x ¯ ( t − τ 2 ) ⋅ e − j 2 π α t d t (10)

By introducing e − j π α ( τ 2 − τ 2 ) = 1 into the equation eqref Equation Autocor Cyclic2, we end up with:

R x ( T ) ( α , τ ) = 1 T ∫ − T 2 T 2 x ( t + τ 2 ) e − j 2 π α ( t + τ 2 ) x ¯ ( t − τ 2 ) e − j 2 π α ( t − τ 2 ) d t = 1 T { z ( α , τ ) ∗ z ¯ ( α , − τ ) } (11)

where z ( α , τ ) = x ( t + τ 2 ) e − j 2 π α ( t + τ 2 ) and * denotes the convolution product operator.

The Fourier transform of the cyclic correlogram R x ( T ) ( α , τ ) gives the cyclic periodogram I ( T ) ( α , ν ) [

I ( T ) ( α , ν ; t ) ≜ 1 T T F { R x ( T ) ( α , τ ) } = 1 T T F { z ( α , τ ) ∗ z ¯ ( α , − τ ) } = 1 T X T ( t , ν + α 2 ) X ¯ T ( t , ν − α 2 ) (12)

where X T ( t , ν ) is the Short Term Fourier Transform of x ( t ) defined as follows:

X T ( t , ν ) ≜ ∫ t − T 2 t + T 2 x ( u ) e − j 2 π ν u d u (13)

X T ( t , ν + α 2 ) is named the complex demodulate [

The cyclic periodogram I ( T ) ( α , ν ; t ) is a non-consistent estimator of the cyclic spectrum S x α ( ν ) . As in the case of stationary processes, consistency can be obtained by temporal or frequency averaging of I ( T ) ( α , ν ; t ) . We then obtain two estimators of the cyclic spectrum: the frequency-averaged cyclic periodogram G x ( Δ t , Δ ν ) ( α , ν ) and the time-averaged cyclic periodogram S x ( Δ t , Δ ν ) ( α , ν ) [

G x ( Δ t , Δ ν ) ( α , ν ) = 1 Δ ν ∫ ν − Δ ν / 2 ν + Δ ν / 2 I ( Δ t ) ( α , λ ; t ) d λ (14)

S x ( Δ t , Δ ν ) ( α , ν ) = 1 Δ t ∫ t 0 − Δ t / 2 t 0 + Δ t / 2 I ( 1 / Δ ν ) ( α , ν ; t ) d t (15)

with Δ t the duration of observation of the signal x ( t ) and Δ ν the frequency resolution.

S x ( Δ t , Δ ν ) ( α , ν ) and G x ( Δ t , Δ ν ) ( α , ν ) are defined by a frequency resolution Δ ν and a cyclic frequency resolution Δ α .

Δ t denotes the duration of observation of the signal x ( t ) and T the duration (or width) of the sliding window (apodization window), we have [

Δ α = 1 Δ t (16)

Δ ν = 1 T (17)

To obtain a fine estimate of the cyclic spectrum, the Grenander uncertainty condition must be respected [

Δ t ⋅ Δ ν ≫ 1 (18)

It states that for a better estimate of the cyclic spectrum, the duration T of the apodization window w ( t ) must be much less than the observation duration Δ t of the signal x ( t ) to analyze [

The estimators S x ( Δ t , Δ ν ) ( α , ν ) and G x ( Δ t , Δ ν ) ( α , ν ) of S x α ( ν ) are asymptotically equivalent [

lim Δ ν → 0 lim Δ t → ∞ S x ( Δ t , Δ ν ) ( α , ν ) = lim Δ ν → 0 lim Δ t → ∞ G x ( Δ t , Δ ν ) ( α , ν ) = S x α ( ν ) (19)

However, in [

• FAM [

• Spectral Correlation Algorithm (SSCA) algorithm [

The FAM algorithm is faster and more efficient than the SSCA algorithm. Our detection model therefore uses the FAM algorithm to estimate the cyclic spectrum of the observed signal.

Let x ( n ) , n ∈ ℤ be the discrete signal obtained by uniform sampling of the continuous signal x ( t ) , t ∈ ℝ . Let N be the number of samples observed during the duration Δ t , Δ t = N T e and N P the number of samples contained in the sliding window w ( n ) ; T = N P T e ; T e denotes the sampling period. The cyclic spectrum estimator by time averaging is written [

S x ( N , 1 N P ) ( α , ν ) = 1 N ∑ n = 0 N − 1 1 N P X T ( n , ν + α 2 ) X ¯ ( n , ν − α 2 ) = 1 N ∑ n = 0 N − 1 1 N P T F { w ( k ) x ( n − k ) e − j 2 π ( α 2 ) ( n − k ) T e } T F { [ w ( k ) x ( n − k ) ¯ ] e + j 2 π ( α 2 ) ( n − k ) T e } (20)

The FAM algorithm for Fast Fourier Transform Accumulation Method, based on Equation (20), takes advantage of the advantages of the Fast Fourier Transform (Fast Forier Transform (FFT)) to reduce the computational cost of the spectrum cyclic. It divides the dual-frequency plane ( ν , α ) into sub-blocks

( ν i , α i ), then calculates the estimates of the cyclic spectrum S x ( N , 1 N P ) ( α i , ν i ) of each sub-block using the FFT. It is constitued of three main stages:

1) The computation of complex demodulates X T ( n , ν + α 2 ) and of their conjugates X ¯ T ( n , ν − α 2 ) to obtain the spectral composition of the signal.

2) The two-by-two product between complex demodules and their conjugates corresponds to the calculation of the cyclic correlogram

I ( T ) ( α , ν ; n ) = 1 T X T ( n , ν + α 2 ) X ¯ T ( n , ν − α 2 ) .

3) The computation of FFT of I ( T ) ( α , ν ; n ) to obtain the estimation of cyclic spectrum S x ( N , 1 N P ) ( α , ν ) .

Let s ( t ) be the signal of the primary user (PU) and y ( t ) the signal observed by the secondary user (SU):

y ( t ) = θ ⋅ h ( t ) s ( t ) + b ( t ) (21)

where h ( t ) denotes the impulse response of the channel, θ a Boolean variable and b ( t ) a Gaussian white noise of variance σ b 2 .

The hypotheses H 0 and H 1 are as follows:

• H 0 : no communication signal; θ = 0 ;

• H 1 : communication signal presence; θ = 1 .

The fundamental idea of detection consists in measuring the level of hidden periodicity in the signal y ( t ) .

To measure this hidden periodicity, we first define the Cyclic Domain Profile or CDP D ( α ) which contains the maximum values of the cyclic spectrum S y α ( ν ) :

D ( α ) = max ν ‖ S y α ( ν ) ‖ (22)

Next, we define the crest factor F c of the CDP D ( α ) which will serve as the decision statistic T C F D of the CFD:

T C F D = F c = max α D ( α ) 1 2 N ∑ k = 1 2 N D 2 ( α k ) (23)

The detection criterion is given by:

T C F D ≷ H 0 H 1 λ C F D (24)

λ C F D is the detection threshold predetermined according to the desired P f a [

In order to analytically approach the optimal detection threshold λ C F D , we must estimate the Probability Density Function (PDF) of the crest factor F c under H 0 . Equation (23) shows that the calculation of F c is based on a maximization function; we can therefore approach the PDF of F c both under H 0 and under H 1 by the PDF f Z ( z ) of generalized extreme value or Generalized Extreme Value (GEV) of parameters ( μ , σ , ξ ) defined by [

f Z ( z ) = 1 σ ( 1 + ξ z − μ σ ) ( − 1 ξ − 1 ) exp [ − ( 1 + ξ z − μ σ ) − 1 / ξ ] (25)

with μ location parameter; σ the scale parameter and ξ ∈ ℝ the shape parameter.

P f a = ∫ λ C F D + ∞ f Z / H 0 ( z ) d z (26)

with:

{ μ 0 = μ / H 0 σ 0 = σ / H 0 ξ 0 = ξ / H 0 (27)

P f a becomes as follows:

P f a = ∫ λ C F D + ∞ 1 σ 0 [ 1 + ξ 0 z − μ 0 σ 0 ] − 1 ξ 0 − 1 exp [ − ( 1 + ξ 0 z − μ 0 σ 0 ) − 1 / ξ 0 ] d z (28)

Using the change of variable u = 1 + ξ 0 z − μ 0 σ 0 , we obtain:

P f a = 1 ξ 0 ∫ ( 1 + ξ 0 λ C F D − μ 0 σ 0 ) + ∞ 1 u u ξ 0 exp [ 1 u ξ 0 ] d u = [ exp [ − u − 1 / ξ 0 ] ] ( 1 + ξ 0 λ C F D − μ 0 σ 0 ) + ∞ P f a = 1 − exp [ − ( 1 + ξ 0 λ C F D − μ 0 σ 0 ) − 1 / ξ 0 ] (29)

From Equation (29), we determine the optimal detection threshold λ C F D as follows:

λ C F D = σ 0 ξ 0 [ ( ln ( 1 1 − P f a ) ) − ξ 0 − 1 ] + μ 0 (30)

With:

{ μ 1 = μ / H 1 σ 1 = σ / H 1 ξ 1 = ξ / H 1 (31)

the probability of detection P d can be defined by:

P d = ∫ λ + ∞ f Z / H 1 ( z ) d z P d = 1 − e − [ 1 + ξ 1 λ − μ 1 σ 1 ] − 1 / ξ 1 (32)

Equations (29) and (32) allow to establish the relation between P d et al. P f a :

P d = 1 − exp { − [ 1 + ξ 1 σ 0 ξ 0 [ ( ln ( 1 1 − P f a ) ) − ξ 0 − 1 ] + μ 0 − μ 1 σ 1 ] − 1 / ξ 1 } (33)

A delicate part in this approximation is the estimation of the parameters ξ , σ , k of the probability density function GEV. To do this, we use the Optimal Biased Robust Estimator (OBRE) algorithm proposed in [

In order to compare the theoretical and experimental results, we generate the performance curves P d as a function of Signal Noise Ratio (SNR) and the Reciever Operating Characteristics (ROC) curves. These simulations are performed on different types of signals (16-QAM, QPSK, BPSK, 4-ASK). We group the simulation parameters in

To analyse the evolution of the P d as a function of the P f a , we use the roc curves.

SNR (dB) | −14 | −12 | −10 | −9 | −8 | −7 | −6 |
---|---|---|---|---|---|---|---|

ξ 1 | −0.033 | −0.104 | −0.098 | −0.124 | −0.1016 | −0.0804 | −0.1286 |

σ 1 | 0.439 | 0.528 | 0.557 | 0.627 | 0.647 | 0.663 | 0.718 |

μ 1 | 5.76 | 6.326 | 7.126 | 7.673 | 8.19 | 8.839 | 9.533 |

The detector based on the cyclostationary characteristics of the signal acts in a quasi-blind context. It has no knowledge of the characteristics of the signal emitted by the PU. Therefore, we compare it to the classic detector based on the energy level. We notice that the probability of detection of the CFD is almost identical to that of the Energy-based Detection (ED). For a false alarm probability set at P f a = 0.05 , we obtained (

64-QAM, 16-QAM, BPSK, 4-ASK Signal | ||
---|---|---|

User | Parameters | Values |

Primary | Sampling rate F e | 128 kHz |

Data rate | 16 kHz | |

bandwidth W | 21.6 kHz | |

Secondary | Observation duration Δ t | 16 ms |

Number of samples N | 2048 | |

Frequency resolution | Δ α | 62.5 Hz |

Δ ν | 1 kHz |

Detecting the Primary User (PU) signal, in a very weak SNR condition, is one of the most difficult tasks in cognitive radio. In this paper, we have proposed a non-cooperative detection model which is able to detect the signal of PU in low SNR scenario. This detection model exploits the cyclostationary features of communication signals. This cyclostationarity can be revealed through cyclic spectral analysis. Thus, after the estimation of the cyclic spectrum using the FFT Accumulation Method (FAM) algorithm, we generate the Cyclic Domain Profile (CDP). Then, we quantify this level of periodicity using the CDP crest factor. We use the CDP crest factor as the decision statistic to detect the state of the channel (busy or not). We approached the probability density function of the crest factor, which allowed us to analytically determine the optimal detection threshold, the probabilities of false alarm and detection. The simulations carried out subsequently show a strong similarity between the theoretical results and the experimental results. Cyclostationary Feature based Detection is an efficient detector. In low SNR scenario, it is able to detect a digital communication signal up to SNR = −11 dB. Our work on CFD shows that the CDP makes it possible to visualize the hidden frequencies of the observed signal on the one hand and on the other hand that the CF is a characteristic quantity of the modulation technique. We propose, in future works, to exploit these two tools (the CDP and the Crest Factor) to develop an algorithm, which is able to detect the signal of the PU and also to distinguish the signal of a PU from the signal of a SU, by using the principle of modulation recognition.

The authors declare no conflicts of interest regarding the publication of this paper.

Kadjo, J.-M., Agoua, R., Bamba, A., Konaté, A. and Asseu, O. (2021) Non-Cooperative Spectrum Sensing Based on Cyclostationary Model of Digital Signals in the Context of Cognitive Radio. Engineering, 13, 56-70. https://doi.org/10.4236/eng.2021.131005

CDP: Cyclic Domain Profile

CFD: Cyclostationary Features based Detection

COR: Caractéristique Opérationnelle du Récepteur

CS: Cyclostationnary

ED: Energy based-Detection

FAM: Fast fourier transform Accumulation Method

CF: Crest Factor

FFT: Fast Forier Transform

GEV: Generalized Extreme Value

MLE: Maximum Likewood Estimator

OBRE: Optimal Biased Robust Estimator

OFDM: Strip Spectral Correlation Algorithm

PU: Primary User

ROC: Receiver Operating Characteristic

SNR: Signal Noise Ratio

SU: Secondary User