Analyzing Chaos Systems and Fine Spectrum Sensing Using Detrended Fluctuation Analysis Algorithm

1Mathematical Department, Polytechnic University of San Luis Potosi, Urbano Villalón 500, 78363 San Luis Potośı, SLP, Mexico 2Telematics Engineering Department, Polytechnic University of San Luis Potośı, Urbano Villalón 500, 78363 San Luis Potośı, SLP, Mexico 3Centro de Investigación y Desarrollo en Telecomunicaciones Espaciales de Ingenieŕıa Eléctrica, Universidad Autónoma de Zacatecas, 801 Ramón López Velarde, 98000 Zacatecas, ZAC, Mexico


Introduction
Detrended fluctuation analysis (DFA) is an effective numerical tool used to measure the persistency (or antipersistency) of data series with nonstationarities.DFA characterizes their complex behavior by identifying the long-range power law correlations.The DFA algorithm has several advantages over the conventional analysis methods such as power spectrum and correlation analysis owing to that it avoids false detections originated from artifacts of spurious or not precise measured observations from real systems.Besides, DFA is widely used to detect and/or characterize if a time series is random or not.DFA algorithm has been applied in different areas, for example, weather, economy, biology, meteorology, and climate [1,2].According to the new research trends, the DFA algorithm is implemented to analyze the WEB [3], to estimate synchrophasor measurements [4], to analyze oil reservoirs [5], to study volcano seismicity [6], and to study earthquakes [7].
Chaotic systems (CSs) are nonlinear deterministic systems that generate random and unpredictable time series which have statistical properties similar to random processes [8].CSs generate time series with chaotic behavior that can be observed in data from observations [9].It is important to study the effects of the DFA algorithm on the chaotic series because many noisy data observations are extracted from physical, biological, and social systems.These systems exhibit scale-invariant features and contain hidden long-range power law correlations that can be detected using DFA algorithm and can provide interesting and useful information about the structure and evolution of any dynamic system.There are few articles that can be categorized in this research direction, like in [10] where the authors studied and characterized time series from a Chua System.
The DFA algorithm is investigated in many research articles [11,12] and applied to solve some problems in many engineering areas.Furthermore, the DFA algorithm is compared with other well-known and effective algorithms such as power spectral density analysis [13], correlation dimension [14], and wavelet [15].
A very big set of tools have been developed to analyze the deterministic and predictable signals such as periodic, transient, and stochastic signals.Chaotic signals are deterministic signals with unpredictable related time series.In many cases, chaotic signals appear like noise (deterministic but not predictable) and could be analyzed using conventional stochastic tools.In [16], a numerical study of the DFA algorithm response under data series from discrete chaotic systems is presented.It is shown that the DFA algorithm detects, in low scale intervals, the short range correlations in chaotic data series that are very similar to noise.The last characteristic can have important applications in telecommunication systems.The moving DFA (MDFA) algorithm in [17] could detect time-instants of abrupt changes in dynamic structures where the analyzed chaotic system parameters are changed.The MDFA performance is independent of the length of subseries and has a perfect capability to resist the effects of noise.Unfortunately, no direct applications were mentioned in [16,17] for the presented work.Besides, there is no statistical analysis about the robustness of the DFA algorithm under the variation of the parameters and the initial conditions of the chaotic systems.The missed statistical analysis is presented in this paper (the initial conditions of chaotic data series from real systems are usually unknown).
An interesting analysis is presented by [18] in the field of signal processing when the DFA algorithm is employed as a denoising technique.The DFA algorithm in [18] is used to select the band limited intrinsic mode functions (IMFs) of the broken down noisy signal and to define the relevant modes in order to construct the filtered signal.The proposed method is valid only under the use of variational mode decomposition (VMD) and for white Gaussian noise.In [19], the authors proposed a metric based on DFA to define a robust threshold that determines which oscillations (called IMFs) are noise components or noise free signal components under empirical mode decomposition (EMD).The presented results are promising and significant around 0 dB range of signal-tonoise ratio (SNR), limited SNR region, and for the white Gaussian noise as in the case of [18].The multifractal DFA (MF-DFA) algorithm is proposed as a detection approach of ionospheric irregularities for global navigation satellite system (GNSS) receivers [20] (GNSS signal detection).The MF-DFA was combined with complementary ensemble empirical mode decomposition (CEEMD) to analyze nonlinear and nonstationary signals (positive and negative white noises are added to the data).This combination increases the complexity of the detection technique and, in some cases, leads to generate a different number of modes that contain new components not related to the signal at the receiver input (the inclusion of white noise with inappropriate amplitudes).
Since many physical signals are noisy and heterogeneous and exhibit different types of nonstationarities which can affect the correlation properties of these signals and based on the previous review of the related literature, it is evident that the DFA algorithm can be applied in different engineering areas and helps to alleviate some specific problems.
In this paper, with the purpose to extend and observe the DFA algorithm response to data from well-known chaotic systems, the Lorenz, the Rossler, and multistable systems are considered.The main properties of the DFA algorithm that could be used to characterize nonlinear systems time series are presented.The output of this analysis leads us to propose the employment of the DFA algorithm in new engineering area, namely, spectrum sensing in cognitive radio (CR) systems.
It is evidently confirmed that the available licensed spectrum is not efficiently and fully utilized by the related wireless technologies in time and spatial domains.Cognitive radio (CR) introduces a futuristic concept that helps to dynamically use the spectrum in an opportunistic manner.The spectrum sensing (SS) in CR defines accurately the vacant frequency bands (frequency holes) by detecting the existence of the primary user (PU) signals in the spectrum of interest.Up to this day, plenty of SS techniques and approaches were presented under different initial practical conditions such as low signal-to-noise ratio (SNR), spatially correlated multiple antennas, and noise power uncertainty.In [21], a wide review of SS algorithms and some other important topics in CR are discussed.
Owing to its simplicity when no prior knowledge about the PU signal is required, the energy detector (ED) is commonly used as a coarse spectrum sensor [22].The ED suffers from performance degradation at low SNR with noise variance uncertainty (SNR wall problem) and spatial correlation of multiple antennas (between the antenna array elements) [23,24].
The two stages sensing architecture that combines two detectors was proposed by the IEEE 802.22 working group (standards).In this work, we consider two-stage SS scheme for CR systems.The coarse (fast, several tens of micro seconds) sensing stage is performed by the ED and the fine sensing (several tens of milliseconds) stage is based on the DFA algorithm.The topology of two stages SS is proposed with the purpose to maximize the probability of detection under practical imperfections like noise power uncertainty and spatially correlated antenna array elements.The main motivation of using the DFA algorithm for fine SS is that DFA shares with ED the same property of no information about the PU signal is needed and the SS can be effectively performed based on the input data samples.The final decision about the PU signal presence or absence is made by combining the decisions of the two stages using "OR" rule (other types of rules can be applied such as "And" or "Majority").The implementation of DFA algorithm for fine SS demonstrates promising results and gives insight about the possibility of using DFA in CR systems as a new signal detection scheme and also for other applications in telecommunication and sensor systems.All the presented numerical results in this paper are generated using MATLAB.
The remainder of this paper is organized as follows.Section 2 contains the DFA algorithm main steps.In Section 3, the DFA is applied in periodic systems (sinusoidal signals) under different parameter values.In Section 4, various numerical experiments are performed for three different chaotic systems, namely, Lorenz, Rossler, and the multistable system.The DFA algorithm implementation for fine spectrum sensing is presented in Section 5. Finally, the conclusion remarks are given in Section 6.

DFA Main Steps
Define (), the data series to be analyzed by DFA where 1 ≤  ≤  and  is the data series size.The DFA algorithm can be applied using the following standard steps: (i) Calculate the random walk of the data series () using where ⟨⟩ is the mean of the data series () given by (ii) The series () is divided into ⌊/⌋ data intervals of equal size  where the notation ⌊⋅⌋ represents the floor function.
(iii) In each data interval, the data is fitted by using a polynomial function of order  to obtain   () which is called the local trend (order- DFA).In previous studies [25][26][27], the value  = 1 is commonly used, that is, the linear detrended fluctuation of the signal.Nevertheless, any polynomial function of grade  ̸ = 1 could be used for the data fitting step.(iv) The integrated signal () is detrended by subtracting the local trend   () in each data interval generating the detrended fluctuation signal: (v) For an interval size , the root mean square (rms) of the detrended fluctuation signal is calculated: (vi) In order to find the relation between the scaling exponent and the data fluctuation, the previous steps are repeated for different data interval sizes 4 ≤  ≤ /4.The power law relation indicates the scaling presence where  is the scaling or the correlation exponent.Algorithm 1 gives a brief summary about the main DFA steps.In this paper, the polynomial functions of grade  = 1 are used.Thus,  1 () is referred to as ().The correlation properties of a signal are characterized based on the scaling exponent .For the linear case ( = 1), if  = 0.5, then the signal under analysis is uncorrelated (random process).In the case when  < 0.5, the signal is anticorrelated, and for  > 0.5 the signal contains positive correlations.The scaling exponent  can be determined using (5) as follows: where  is a positive integer.The value of  is defined from the slope of the straight line that is adjusted from the points (log  (()), log  ()).For the presented numerical analysis in this paper, the value  = 2 is used.

Periodic (Sinusoidal) Signals
The DFA algorithm response to periodic signals by using sinusoidal sequences is investigated.The periodic signal used in this section is defined as: when the amplitude  = 1, the frequency  = {1, 8, 16}, and the phase  = 0.After analyzing the DFA algorithm for a periodic signal with only one harmonic component, a numerical analysis of the DFA response for signals with more harmonic components is conducted.Figure 2(a) shows the compound signal defined as  1 sin(2 2 )+ 2 sin(2 16 ) case at  1 =  2 = 1,  2 = 2, and  16 = 16 (with periods  2 = 1/2 and  16 = 1/16, resp.)under sampling time   = 2 −12 .Furthermore, the respective DFA response is presented in Figure 2(b).
As shown in the previous case, the dents in Figure 2 are positioned at ( 16 /2)/  = 2 12 /2 5 = 2 7 and ( 2 /2)/  = 2 12 /2 2 = 2 10 that are placed at the instants where each sinusoidal signal crosses the horizontal axes towards the negative values.The case when the amplitudes of the sinusoidal signals are not equal ( 1 ̸ =  2 ) is presented in Figure 2(c) when the component with the higher frequency  = 16 has greater amplitude  2 = 1 in comparison with the amplitude of sinusoidal component with lower frequency.As the amplitude  1 becomes smaller, the DFA algorithm does not detect the component  1 sin(2 2 ).It is inferred that the fluctuation graphs  calculated with the small  1 become horizontal when the period of the component  2 sin(2 16 ) is completed; that is, at the scale value  −1  /16 = 2 12 /2 4 = 2 8 .On the other hand, Figure 2(d) contains the DFA results for a constant  1 = 1 and for  2 <  1 .Analogous to the previous results, the DFA algorithm is not able to detect the component  2 sin(2 16 ).Besides, it is inferred that the deviation at the scale value 2 7 vanishes as  2 becomes smaller.

Chaotic Signals
In this section, we are interested in processing time series of chaotic systems arising from differential equation that can be presented using the following form: where  = [, , ] ∈ R  is the state variable (in our case  = 3), () is a function with domain and range space ⊂ R 3 , (ii) when (⋅) is a piecewise linear function.
The DFA algorithm numerical results are presented for time series of differential systems (8) when the related initial conditions are varied and for various parameter values (such that their orbits conserve a chaotic regime).

Lorenz
System.The Lorenz system is a very well-known third-order chaotic system [28] and can be defined as follows: where , , and  make up the system state; , , and  are the system parameters.Figures 3(a), 3(b), and 3(c) show the time series of the Lorenz system's state variables , , and , respectively.To generate these time series, the fourth-order Runge-Kutta method is used with a step of 0.01 (this step size is used for all numerical experiments related to Lorenz system) and main parameter values  = 10.0,  = 28.0, and  = 2.66.(i) The first fixed fluctuation state for the scale interval 2 4.75 ≤ scale < 2 7 where the values of the scaling exponent are 1.57, 1.437, and 1.437 for the variables , , and , respectively.
(ii) The second fixed fluctuation state for the scale interval 2 7.5 ≤ scale < 2 12 where the corresponding scaling exponent values are approximately 0.75, 0.72, and 0.3 for , , and , respectively.
(iii) The transition between the previous states.
It is well known that the chaotic systems are very sensitive to initial conditions.In other words, regardless of the closeness to the initial conditions, the chaotic systems orbits will flow for different ways.   of displaying a set of curves with respect to the average, the standard deviation is calculated using the following formula: where  is the total number of curves in the set (generated by varying the system initial conditions or by varying the system parameter values),   () is the value of  of the curve  and the scale , and ⟨  ()⟩ = (1/) ∑  =1   () (the average).Figure 4(d) shows that from the scale axis of value 2 5 onward, the standard deviation for  curves from the times series  and  is identical.For not so large scale values, the graphs are near to the average for time series in the case of all the three variables.Approximately, starting from 2 9 , the value of   () begins to grow for the three cases.
To illustrate the response of the DFA algorithm to the parameters of the system presented by (3), the main parameters are varied as follows: In these parameters' variation, the system stability evolves to one of two equilibrium point attractors with two scrolls.curves for the three variation cases.Figure 5(d) shows the standard deviation of the three cases where it is noticeable that the forms of all curves are similar.
A comparison of DFA graphs under different parameter values for times series of variable  and time series of variable  is shown in Figure 6.Thus, it contains the average of the difference between detrended fluctuations of time series  and detrended fluctuation of time series  for the three Lorenz main parameters.
The results for the time series of the variable  are presented in Figure 7.It is shown in Figure 7(a) that the DFA curves that correspond to variations of parameter  are above the DFA curves that correspond to the parameter .Furthermore, they have approximately similar shapes with different vertical positions.Thus, when the parameter  or  varies, it causes a change in the vertical position of the DFA curve.Figure 7(b) includes the average of each set of the DFA curves determined when the parameter value is varied.The average curves which correspond to every Lorenz parameters begin to be closer when the scale value increases.Meanwhile the standard deviations of the time series fluctuations trend to be the same for the three parameters in large scales (Figure 7(c)), and for larger parameter  values, the standard deviation tends to be bigger than the standard deviation values from the other parameters (Figure 7(d)).
To find the relation between the DFA curves, the correlation coefficient () between each pair of DFA curves obtained from time series of an arbitrary variable but at different parameter is defined.The correlation between DFA curves   and   can be defined as where  is the number of data in the scale.The relationship between th curve with the rest of other −1 curves generated at different values of a specific parameter is defined as follows: The correlation coefficient is calculated for DFA curves from time series of the variables  and .For the variable  the results are 0.9976, 0.9910, and 0.9958 when the parameters , , and  are varied, respectively.For the time series  we have 0.9903, 0.9927, and 0.9907 corresponding to the parameters , , and  variations, respectively.The presented numerical results indicate that the DFA algorithm has a constant response and changing the values of the parameter does not affect the stability and the asymptotic behavior of the system.

Rossler
System.The Rossler system can be defined using (8) with the following changes: where , , and  are the system parameters.From the first sight to the DFA algorithm response toward the chaotic time series from Rossler system, it can be seen that the DFA curves of time series  and  are quite similar.On the other hand, the dashed vertical line in Figure 8(d) presents the separation between different behaviors in the time series detected by the DFA algorithm.For scale values smaller than 2 7.25 , the values of  slopes are approximately 1.4, 1.5, and 0.92 for the time series , , and , respectively.As a result, at this time scale the DFA algorithm describes the time series  and  as persistent data and the time series from variable  as a Brownian motion.At scale values greater than 2 8 the values of the  slopes are approximately 0.0 for the three time series which leads to conclude that, for large scales, the DFA algorithm characterizes the Rossler time series as a pure antipersistent processes represented by an alternating values sequence.
Figure 9 presents the study of the DFA algorithm response for Rossler time series initialized with the set −0.8 ≤ (0) ≤ 10 and −10.0 ≤ (0) ≤ 8 at increments of 1 while  variable initial value is fixed to (0) = 3.0.Figures 9(a), 9(b), and 9(c) contain the set of the DFA graphs under these initial conditions for time series , , and , respectively.As seen in Figure 9, the DFA response for the time series of variable  is similar to the DFA response for time series of variable.Besides, their respective standard deviations from each set of DFA curves are plotted in Figure 9(d).As expected, the standard deviation of the DFA curves is bigger when the scale value is bigger likewise.In spite of the DFA standard deviation of the curves from the time series  being bigger in comparison with the standard deviation for DFA curves from the time series  and , the standard deviation of the time series  tends to a small and constant value.

Multistable System.
In the previous subsections, we explore the DFA algorithm response to time series from chaotic system under the case when (⋅) is a nonlinear function.In this section, the time series from a system defined under the case when (⋅) is a piecewise linear function; namely, the multistable system is considered.The multistable system is defined as where  31 ,  32 , and  33 are system parameters, and  : R → R is the piecewise linear function: In [29], it is shown that based on the number of linear pieces in the function (), the orbit of the system ( 14) and ( 15) could flow with the same number of scrolls.
For the numerical integration made with multistable oscillator, the form ( 16) is applied and in most cases, the assigned values to the system (14) parameters are chosen to generate orbit solutions with shapes of three scroll oscillators.Figures 10(a), 10(b), and 10(c) present the time series of the multistable oscillator corresponded to the variables , , and  generated with a step size ℎ = 0.1 and parameter values  31 = 1.5,  32 = 1.0, and  33 = 1.0 (ℎ = 0.1 is used in each numerical integration given in this subsection).In Figure 10(d), the DFA corresponding curves are presented where one can see that the persistency of each time series changes as a function of the scale.On the scale interval [2 4 , 2 6 ), the scaling exponents are 1.75, 1.64, and 1.54 for time series , , and , respectively; that is, the DFA algorithm characterizes them as time series with long-range correlation.Meanwhile on the scale interval (2 7 , 2 12 ), each time series has different characterization.The values of the scaling exponents are 0.8772 (persistent series), 0.2087 (antipersistent), and 0.0144 (antipersistent series oscillating between two values).
Figure 11 shows the results for the DFA response when the initial conditions ((0), (0), (0)) of the multistable numerical integration are varied as −1.0 ≤ (0) ≤ 1.0, −0.6 ≤ (0) ≤ 0.6, and a constant value for the initial condition (0) = 1.0 with the same parameter values mentioned above.The DFA curves for time series of the three variables are plotted in the same figure (Figure 11).In the interval [2 4 , 2 6 ], the curves are almost overlapped; meanwhile in the interval [2 6 , 2 12 ] the set of curves corresponding to time series variable  are located in the upper position, DFA curves for time series  are in the middle position, and DFA curves for time series  are located in the below position.The error average percentage coefficients  = (1/) ∑  =1  2 () for each curves set are 2.94%, 0.58%, and 0.37% which correspond to time series of variables , , and , respectively.It is noticeable that these percentages have very small values which indicate that the forms of the curves are similar among them.
Figure 12 shows numerical examples of  graphs from the time series of the three variables of the multistable oscillator when the parameter values are varied and the initial conditions are fixed to (1.5, 1.0, 1.0).The parameter value intervals are selected as a set of values so that the numerically generated orbits with systems ( 14), (15), and ( 16) are strange attractors.As in the previous oscillators time series analysis,  graphs from the time series generated under various parameter values have similar form.Figure 12(a) corresponds to the case for  graphs of time series  generated when the parameter  31 is varied and the parameters  32 and  33 are fixed to 1.0.It is easy to see that from the values of scale > 2 6 , the graphs begin to separate in two groups.Graphs located in the lower positions are generated by varying the parameter  31 in the range 1.0 <  31 < 1.05 where the generated orbits with these values are strange attractors with disk form like the one shown in Figure 12(d  Furthermore, to describe the similitude and the closeness among DFA curves, the percent relative standard deviation is used and formulated as  A more complete analysis for the time series of the system defined by ( 14), (15), and ( 16) is listed in Table 2 where the average of scaling exponents and the relative standard deviations as a function of the parameter value interval (second column) and the scale interval (third column) are presented.As we can see, the parameter  31 has only two parameter intervals.The interval [1.0, 1.04) corresponds to multistable oscillator graphs with only one scroll like the plot in Figure 12(d).Meanwhile, the interval [1.05, 1.53] corresponds to triple scroll oscillators.As a consequence, the scaling exponents are different from each other.Besides, Table 2 shows that from the small scale ([2 4 , 2 6 ]) to the large scale ([2 7 , 2 12 ]), each time series decreases its scaling exponent value.With respect to the parameter  31 in the interval [1.0, 1.04], the scaling exponent average of the sets of time series for variables  and  becomes zero at large scales because their time series are oscillating signals with amplitude that varies very slowly.For time series , the DFA algorithm characterizes the change from series with long-range correlation (small scales) to antipersistent series (large scales).For time series  from triple scroll oscillator, the DFA algorithm characterizes the change from persistent series (small scales) to series oscillating between two values (large scale).Finally, for the cases with  rel > 20%, this value does not mean that the DFA curve slopes are greatly expanded over the average.The real reason behind this case is that the scaling exponent averages are very close to zero ( → 0).

DFA Algorithm for Fine Spectrum Sensing
As mentioned before, the DFA algorithm is used for fine SS (PU signal detection) in two-stage spectrum sensor.This idea is proposed based on the DFA algorithm behavior analysis presented in Sections 3 and 4. The DFA algorithm performance and response are stable and invariant under different initial conditions and for different chaotic system parameters.The last important property (DFA response robustness against the changes of the initial conditions and parameters) made the use and implementation of the DFA algorithm for fine SS a very practical solution under various types of imperfections that change the initial conditions of the SS process such as spatial correlation between the antenna array elements and noise variance uncertainty.Additionally, one of the main conclusions from the presented analysis in Sections 3 and 4 that the DFA algorithm stability, reliability, and immunity toward the fluctuations in the initial conditions is maintained at the expense of relatively large number of samples.As a result, the DFA is implemented for fine SS and not for coarse SS.A brief summary about the last discussion can be found in Figure 13.These two sensing stages can be performed in serial or parallel modes.For serial mode, if the ED decides that there is PU signal in a certain radio channel, the fine sensing stage is eliminated or aborted.Otherwise, the fine sensing stage is performed on the same radio channel for the final decision.In the case of parallel mode, the two sensing stages are performed in parallel and the final decision about the PU signal presence or absence is made by combining the decisions of the two stages.The block diagram of the proposed spectrum sensor is shown in Figure 14 where the ED is implemented for the coarse or fast SS and the fine SS is performed by the DFA.

System Model.
In the bandwidth of the spectrum to be sensed there are number of radio channels equal to  ch and all these channels are serially sensed.The proposed spectrum sensor is equipped by antenna array with the number of elements equal to  and each antenna array element receives  samples during the sensing time.The ED accumulates the energy of  samples and compares it to the ED threshold (THR ED ) to decide the PU signal presence or absence.The times required for coarse sensing  1 and for fine sensing  2 per each channel are given by where  1 and  2 are the numbers of samples needed for coarse and fine sensing, respectively, and  is the channel bandwidth.Since  2 >  1 then it is axiomatic that  2 >  1 .
In the case if the two sensing stages are performed in parallel mode, the total sensing time for all radio channels is defined as where  ch is the total number of radio channels to be sensed and  fs is the mean number of the reported radio channels for the fine sensing stage.The PU signal detection problem (spectrum sensing) can be presented using the conventional binary hypothesis test: = 1, . . ., ;  = 0, . . .,  − 1, where The DFA algorithm is more suitable for fine spectrum sensing not for coarse spectrum sensing The DFA algorithm can be used for spectrum sensing (PU signal detection) Figure 13: A brief summary about the DFA algorithm studying results and the related application.

Media access control
Receiver front end Antenna array (M) Figure 14: Two-stage spectrum sensing using ED and DFA.
obeying the circularly symmetric complex Gaussian distribution with zero mean and variance equal to  2 ℎ , that is, ℎ  [] ∼ CN(0,  2 ℎ ); and   [] is the discrete-time PU signal.The channel coefficients ℎ  [𝑘] are spatially correlated between each other.The coefficient of spatial correlation  sc between the adjacent antenna array elements (0 ≤  sc ≤ 1) can be given as [30] where  is the wavelength;  is the distance between two adjacent antennas; and Λ is defined as where  is the angular spread.The components of the  ×  antenna array element correlation matrix C can be presented in the following form: The  × 1 signal vector X that collects all the observed signal samples during the sensing time can be defined using the following form: where  denotes a transpose.The data distribution of the matrix X can be expressed as where   is the average energy of the PU signal and I is the  ×  identity matrix.

Fine Spectrum Sensing
Using DFA.In general, the radio channel with low power has a high probability to be unused by the PU.The fine sensing stage in this scheme is performed owing to the fact that the ED threshold definition is susceptible to noise power variations and as a consequence, its detection performance.The DFA algorithm behaves differently if the input data, in this case X, contains noise samples only (the hypothesis  0 ) or contains samples of PU signal plus noise (the hypothesis  1 ).The test statistics based on the DFA algorithm is highly attached to the output parameter or the scaling exponent  given by (6).Thus, the DFA test statistics  DFA (X) can be defined using the following form: Under the hypothesis  0 , the value of the scaling exponent  of the DFA algorithm indicates that the input data corresponds to random or stochastic process like colored noise The final decision about the PU signal presence or absence is obtained by combining the decisions of the two stages using the decision making rule.For example, if the "OR" rule is used, Table 3 presents the two sensing stages final decision making results.
In the case of two-stage spectrum sensing scheme, the probability of false alarm  FA and the probability of detection   are defined with the respect to the two stages test statistics.Under the hypothesis  0 , if  ED (X) > THR ED or if  ̸ = 0.5 when  ED (X) ≤ THR ED , the false alarm occurs, where  ED (X) is the ED test statistics given by The correct PU detection occurs under the hypothesis  1 , if  ED (X) > THR ED or if  ̸ = 0.5 given that  ED ≤ THR ED .Thus, the overall value can be defined as where  ED FA and  ED  are the probability of false alarm and the probability of detection of the first sensing stage using ED, respectively, and  DFA FA and  DFA  are the probabilities of false alarm and detection of the second stage using DFA, respectively.

Fine Spectrum under Noise Power Uncertainty.
The second type of imperfections after the spatial correlation between the adjacent antenna array elements is the noise power or variance uncertainty.The sensitivity to noise power uncertainty (the fluctuations on the noise variance value as a function of time) is one of the most common problems among the spectrum sensors such as the energy detector (ED), matched filter (MF), and even the cyclostationary detector under the low signal-to-noise ratio (SNR) [31,32].The negative effect of noise power uncertainty is quantified by the SNR wall location.If the SNR value is less than the SNR wall, the PU signal detector (spectrum sensor) will fail to achieve the desired detection performance and maintain the immunity against power noise uncertainty regardless of the sensing time length.This problem negatively affects the receiver operation characteristic (ROC).
The SNR wall phenomenon can be presented using the sampling complexity of the spectrum sensor that is defined as the number of samples  as a function of the SNR, the probability of false alarm  FA , and probability of miss  miss :  =  (SNR,  FA ,  miss ) .
In general, any spectrum sensor (or the PU signal detector) must minimize the number of samples  required to achieve the desired detection performance.The lowest SNR satisfying the probability of false alarm  FA and the probability of miss  miss constraints is called the detector sensitivity.Under noise power uncertainty initial conditions, the actual noise power or variance  2 at the spectrum sensor input can be determined only within the limits of a definite range [31]: where  is the uncertainty parameter defined as where  is the parameter used to define the amount of nonprobabilistic uncertainty in the noise power.Clearly, small values of  (least uncertainty case) are preferred.As an example, the sampling complexity of the ED can be defined using the uncertainty parameter  as follows: where is the Gaussian -function.From (32), we can define the ED SNR wall in the following form [31]: With the new initial conditions when there is noise power uncertainty, the data distribution of the matrix X given by ( 25) should be modified to consider the uncertainty parameter  (the actual noise power  2 ): As a result, the sample complexity tends to approach infinity as the SNR decreases to approach the SNR wall: It is important to mention that the SNR at the input of the spectrum sensor is defined using the following form: As seen from ( 30) and (37), the noise power uncertainty leads to SNR value uncertainty at the spectrum sensor input.

Numerical Results.
The main objective of this section is to show how the DFA algorithm behaves differently in accordance with the input data, in this case the data matrix X of the observed received signal samples at the spectrum sensor input.In other words, the DFA response and the value of the scaling exponent (test statistics) differ with a considerable way when X contains data samples from noise only or data samples from PU signal plus noise (see Figure 15).This fact about the DFA algorithm behavior is illustrated in Figure 15 where in Figure 15(a) the data matrix (DFA input) is formed from noise samples only and in Figure 15(b) the data matrix is formed from PU signal plus noise samples.The simulation results are obtained using the parameters given in Table 4.Note that in this section, the value of scale is equal to the number of samples .The DFA performance as a fine spectrum sensor is shown by following the same evaluation approach presented in Sections 3 and 4. As shown in Figure 15, the DFA algorithm responses and behaves differently in accordance with the input data matrix X contents in the presented two cases, namely, noise only (a) and PU signal plus noise (b).As a direct observation about the test statistics presented in (26), the value of the scaling exponent  in the case of noise is 0.5 ((a),  = 0.5) while in the case of noise plus signal this value is not equal to 0.5 ((b),  ̸ = 0.5).It is important to mention that the value of the scaling exponent  changes in Figure 15(b) according to the scale range.For instance, when the Scale ≤ 2 6 the scaling exponent  = 0.3 and when the Scale > 2 6 the scaling exponent  = 1.0.
Under two types of imperfections, namely, spatially correlated antenna array elements and noise power uncertainty, the DFA response is illustrated in Figure 16 when the value of the  is equal to 1 dB (this parameter value is chosen to roughly give a big noise power uncertainly when in the case of ED we have that SNR ED wall = −3 dB).From Figure 16, it is confirmed that the DFA algorithm performance differs according to the input data type of the matrix X.For instance, in the case of noise only as shown in Figure 16(a), the value of the scaling exponent  is equal to 0.5 ( = 0.5) even under uncertainty case.The attention should be paid to the case of PU signal plus noise presented in Figure 16(b) at two types of imperfections, namely, spatially correlated antenna array elements and noise power uncertainty when the scaling exponent has two values ( ̸ = 0.5) according to the number of samples (the scale range).Approximately, when the Scale ≤ 2 9 the scaling exponent  = 0.24 and when the Scale > 2 9 the scaling exponent  = 2.0.Thus, the DFA algorithm is still able to distinguish between the basic binary hypothesis cases ( 0 and  1 ) even under the combined extreme conditions of spatial correlation and noise uncertainty.Comparing the last two cases presented in Figure 15 (spatial correlation case) and Figure 16 (combined spatial correlation and noise uncertainty case), it is clear that the scale and the scaling exponent have different limits and values.
A complete analysis for the DFA employment in fine spectrum sensing (full theoretical analysis with detailed performance comparison with other spectrum sensing algorithms) requires a special dedicated work that can be considered as a future work for the authors.

Conclusions
In this work, a numerical study about the DFA algorithm response and behavior toward the periodical and chaotic signals is discussed.According to the presented results, each signal can have more than one characterization that depends on the time scale under analysis.For example, in the case of short signal intervals (small scales), the DFA algorithm characterizes the signals as periodic or persistent.For periodical signals like sinusoidal, the slope of the detrended fluctuation graphs changes its value to zero when the scale reaches the first period of the signal.For chaotic signals (not periodic), the DFA characterizes them as periodic at small time scales.
The presented results have evidently demonstrated that the DFA response is invariant and robust to initial condition and chaotic system parameters.For several numerical proofs, the results were very similar when the initial conditions of the orbit are varied and the system parameters are changed.As a result, it is possible to characterize the persistency (or antipersistency) of time series from real systems (like physical systems, social system, biological, etc.) without the need of performing and repeating large number of experiments.
A primitive proposal about using the DFA algorithm for fine spectrum sensing along with ED (coarse stage) in twostage topology is presented.The DFA does not require any prior knowledge of the PU signal parameters and relies only on the data matrix of the sampled received signal.The test statistics of the DFA fine spectrum sensing is simple and based on the main algorithm output, namely, the scaling exponent.The numerical outcomes show promising results about futuristic use of DFA algorithm for spectrum sensing in CR systems under practical imperfections like noise power uncertainty and spatially correlated antenna array elements.

Figure 1 (
a) shows the amplitude of the sinusoidal signal as a function of the time scale, and Figure 1(b) contains the respective  fluctuations.

and 2 12
correspond to the time at which the signal flows from positive to negative values (see Figure1(b)).

Figure 3 (
d) presents the DFA response for these time series.In this figure it is observed that the fluctuations for  and  are very similar and that the general behavior of the fluctuation has three states, two fixed states and a transition between them:

Figure 4 :
Figure 4: The DFA analysis of the Lorenz time series under different conditions for the variables: (a) , (b) , and (c)  and the related standard deviation of the DFA graphs in (a) dotted line, (b) circles line, and (c) continuous line.

Figure 5 :
Figure 5: The DFA analysis of the Lorenz time series under different parameters values: (a) , (b) , and (c)  and the related standard deviation as a function of the scale in (a) dotted line, (b) circles line, and (c) continuous line.

Figure 5
shows the DFA curves for time series when the three parameters are varied.Figures 5(a), 5(b), and 5(c) present the

Figure 6 :
Figure 6: The average of the differences between the DFA of time series () and the DFA of time series () from the Lorenz system at different parameters values:  is varied (continuous line),  is varied (dashed line), and  is varied (circles line).

Figure 7 :
Figure 7: The DFA of the variable  (Lorenz system) at different parameter values; in (a) the dashed line and the dotted line correspond to the DFA curves when the parameters  and  are varied, respectively; in (b), (c), and (d) the dotted line corresponds to the parameter , the dashed line with circles corresponds to the parameter , and the dashed line with diamonds corresponds to the parameter .The average over the DFA graphs when all Lorenz parameters are varied is presented in (b).The standard deviation as a function of the scale and the parameters value range are presented in (c) and (d), respectively.

Figure 8 :
Figure 8: The time series of the Rossler system variables  (a),  (b), and  (c), the DFA of the time series in (d) where the dashed line corresponds to the variable , the continuous line to variable , and circles line to the variable .

Figure 9 :
Figure 9: The DFA of Rossler system time series under various conditions: (a) time series of the variable , (b) time series of the variable , and (c) the time series of the variable .The standard deviation of each set of curves is presented in (d) where the dotted line corresponds to (a), the dashed line to (b), and continuous line to (c).

Figure 10 :ScaleFigure 11 :
Figure 10: The multistable system time series: (a) variable , (b) variable , and (c) variable .(d) shows the DFA for the time series in (a) with continuous line, in (b) with dashed line, and in (c) with dotted line.

Figure 12
(b) shows the fluctuation graphs of time series  when the parameter  32 value is varied as 1.0 <  32 < 1.41 and the other parameters are fixed to  31 = 1.5 and  33 = 1.0.Meanwhile Figure 12(c) shows the fluctuation graphs for time series  generated when the parameter  33 varies as 1.0 <  33 < 1.47.Clearly, the DFA curve sets of all figures are adjacent or spliced among them and they change their slope approximately for scale > 2 7 values.

Figure 12 :
Figure 12: The DFA of the time series multistable system generated at different parameters values: (a) variable  31 , (b) variable  32 , and (c) variable  33 .(d) presents the multistable strange attractor which revolves around one fixed point.

Figure 15 : 8 Coefficient 1
Figure 15: The DFA algorithm behavior under two cases: (a) the data matrix X contains samples from noise only; (b) the data matrix X contains samples from signal plus noise under spatial correlation between antenna array elements.

Figure 16 :
Figure 16: The DFA algorithm behavior under two cases: (a) the data matrix X contains samples from noise only; (b) the data matrix X contains samples from signal plus noise under two types of imperfections (spatial correlation and noise uncertainty).
(period  1 = 1) has 216= 65536 samples by cycle, then the period of the signal with  = 16 is equal to  16 = 2 16 /16 = 2 16 /2 4 = 212, and for the signal with  = 8 we have  8 = 216/8 = 2 16 /2 3 = 213.Thus, for sinusoidal signal, the DFA algorithm gives a saturation value of the scaling exponent  = 0 which corresponds to a signal that oscillates between two values.The dents positioned at scale values 211

Table 1 :
Relative standard deviation of the DFA curves calculated with different parameter values. is the number of scales, and  2 is the variance.In Table 1, the relative standard deviations of each set of DFA curves are presented where a set of curves is obtained by analyzing the time series of the same variable at different parameter values.As we can see in the right column of Table 1 ( rel (%)), the DFA curves are not dispersed.

Table 2 :
The scaling exponent average of the DFA curves and their respective relative standard deviation for multistable system time series generated with various parameter values.The behavior column abbreviations: P: persistent; AP: antipersistent; AP2: antipersistent oscillating between two values; LRC: long range correlation.The symbol "+" in the case of LRC and P refers to the set of the DFA curves contain time series that behaves as Brownian motion.

Table 3 :
The final decision making rule of two stages SS using "OR" rule.