A Robust FLOM Based Spectrum Sensing Scheme under Middleton Class A Noise in IoT

Accessibility to remote users in dynamic environment, high spectrum utilization, and no spectrum purchase make Cognitive Radio (CR) a feasible solution of wireless communications in the Internet of Things (IoT). Reliable spectrum sensing becomes the prerequisite for the establishment of communication between IoT-capable objects. Considering the application environment, spectrum sensing not only has to cope with man-made impulsive noises but also needs to overcome noise fluctuations. In this paper, we study the Fractional LowerOrderMoments (FLOM) based spectrum sensingmethod underMiddleton Class A noise and incorporate a Noise Power Estimation (NPE) module into the sensing system to deal with the issue of noise uncertainty. Moreover, the NPE process does not need noise-only samples.The analytical expressions of the probabilities of detection and the probability of false alarm are derived. The impact on sensing performance of the parameters of the NPE module is also analyzed. The theoretical analysis and simulation results show that our proposed sensing method achieves a satisfactory performance at low SNR.


Introduction
The Internet of Things (IoT) has to construct comprehensive connections among variety of objects distributed over an extensive area.So the resources allocation to this large number of objects has to be resolved carefully to maintain a satisfactory Quality-of-Service (QoS) [1].Generally, the frequency spectrum is one of the most important resources in wireless communications, and the problem of spectrum scarcity is getting worse as a result of the large number of applications [2].Therefore, the available spectrum has to be carefully utilized by the IoT to ensure plenty of reliable connections between different objects.Fortunately, Cognitive Radio (CR) allows unlicensed users to utilize licensed bands opportunistically and enable them to reuse the frequency bands that are not heavily occupied by Primary Users (PU).Being able to address the spectrum scarcity issue, CR as a promising solution to exploit the available spectrum for the IoT has been proposed [3][4][5].The ability of spectrum sensing to measure or sense the presence and absence of PU signal is essential because the operation of CR starts with detecting spectrum holes [6].Spectrum Sensing methods proposed for identifying the presence of signal transmissions can be categorized as energy detector (ED) based sensing [7], matched filtering based sensing [8], waveform based sensing [9], cyclostationarity based sensing [10], radio identification based sensing [11], and so forth.Among these methods, ED based approaches are the most commonly used because of their low computational and implementation complexities [12].
Most of the previous studies on spectrum sensing only focused on signals contaminated by Additive White Gaussian Noise (AWGN).However, this assumption fails to model the behavior of certain noise types in IoT applications.Considering the applications of IoT such as Machine to Machine (M2M) networks and smart grids, a key challenge in establishing the IoT is wireless communication in the vicinity of vehicles, machines, or electrical power equipment which often radiates electromagnetic waves from switching power electronics components.In particular, this kind of waves in the form of impulse noise and high power transients disrupt wireless communication [3,13].Middleton Class A noise model is one of the widely investigated statistical distributions that are used to model this kind of man-made interference and the narrow band impulsive noise in different systems [14].Being different from AWGN hypothesis, ED based detector is no longer an optimal detector and it has poor performance.Besides, Generalized Likelihood Ratio Test (GLRT) based detector as the optimal detector has a very complex structure, which will be explained later, when it is used under Middleton Class A noise.Recently, a large number of spectrum sensing approaches under different non-Gaussian noises have been proposed [15][16][17].However, the implementation of these detectors in the IoT remains challenging because multiple antennas were used or the noise uncertainty was not considered.
Fractional Lower Order Moments (FLOM) demonstrated its capability in signal processing under non-Gaussian noise in [18].When FLOM is applied to spectrum sensing, the test statistic has a similar expression as that of ED based sensing.Nevertheless, determination of the threshold also depends on the noise parameters in FLOM based sensing as in the case of ED based sensing, and a small noise uncertainty will cause significant performance loss [19,20].To this end, a robust FLOM based sensing method should be studied as a promising solution of spectrum sensing under Middleton Class A noise, especially at low SNR.As is shown in our previous work [21], uncertainty of noise power is really destructive, while small estimation errors on other parameters of the noise do not have a strong effect on spectrum sensing performance.So in this paper, the focus is put on the spectrum sensing method under Middleton Class A noise with Noise Power Estimation (NPE), making the following contributions.
(i) We study the problem of spectrum sensing under Middleton Class A noise adopting FLOM based detector and derive the analytical expressions of the probability of false alarm  fa and the probability of detection   .Then, we analyze the relationship between the sensing performance enhancement and the noise parameters.
(ii) We propose such an NPE based structure to deal with the issue of noise uncertainty that noise-only samples are not necessary in the estimation process.The performance of the proposed structure is analyzed, which relates the accuracy of the estimator to the estimation duration and the order of the estimator.
The following parts of this paper are organized as follows.The signal and noise models are defined in Section 2. In Section 3, FLOM are introduced to spectrum sensing as a suboptimal detector.The NPE based sensing structure is proposed in Section 4 and the derivation and analysis are presented in the same section.Section 5 includes the simulation and the results, and the conclusion is drawn in Section 6.

Signal and Noise Models
In spectrum sensing, the PU signal to be sensed is considered as a random process (called Bayesian model) in some works; and it is also considered as an unknown deterministic signal (called classical model) in others [22].Lacking in the knowledge of the PU signals, we choose Bayesian model and consider a source  with a zero-mean Gaussian probability density function (pdf) and  is transmitted over a channel impaired by a Middleton Class A noise , whose pdf is where   =

Spectrum Sensing under Middleton Class A Noise
Depending on the idle state and busy state of the PU, with the presence of the noise, detecting the presence of PU is usually considered as the following binary hypothesis testing problem [23]: in which  = 1, 2, 3, . . ., ,  is the number of observed samples; () is the signal observed by sensing receiver with () and () denoting the PU signal and the additive impulsive noise respectively. 0 means that the PU signal is absent and  1 means that the PU signal is present.
According to the Neyman-Pearson (NP) theorem, GLRT can maximize detection probability when the probability of false alarm is fixed.So we attempted to use GLRT as an optimal method first.With the signal and noise models described in Section 2, the globally optimal detector can be expressed as where y is the vector of the received samples.If Λ GO ≥ 0, it means that the PU signal is present.Otherwise, it means that the PU signal is absent.Substituting (3) into (4), we have in which E{⋅} denotes statistical expectation.
Obviously, Λ GO is with respect to the pdf of the PU signal which may not be obtained by unlicensed users.However, CR always operates in low SNR (i.e.,  2  / 2  → 0) circumstance, especially in IoT.By making the low SNR assumption, we can obtain a locally optimal detector from the globally optimal detector.Equation ( 5) can be simplified by using Taylor series [24], where we drop the time dependence  for clarity.
The second term of the right side in ( 6) equals 0 when taking the expectation.Hence, where we assume E{ 2 } =  2  .Then the locally optimal detector under low SNR can be expressed as Here we use log(1 + ) ≈  for  → 0. Obviously, this locally optimal detector only requires the pdf of the noise.
For the case of AWGN, the differential part of ( 8) can be simplified into The locally optimal detector equals the traditional ED.Unfortunately, (()) cannot be easily simplified under Middleton Class A noise hypothesis.From (8), it can be seen that the locally optimal GLRT detector is with respect to the pdf of the noise which contains infinite summation.Moreover, the structure of the detector has to change with the change of the noise parameters which makes the implementation of GLRT detector impossible.Consequently, (()) should be converted into a simpler nonlinear operator so that it can be implemented practically.Under impulsive noise hypothesis, the presence of impulses increases the false alarm during sensing process.To improve the sensing performance, the impact of randomly appearing large amplitudes in the noise should be reduced.Inspired by the capability of FLOM in signal processing under non-Gaussian noise and the expression of energy detector, we use FLOM as a suboptimal detector and the corresponding test statistic is given in where 0 <  < 2.  Through fractional power operation, a nonlinear operation, the large impulse amplitude in the noise can be reduced, while the small values almost remain unchanged.As a result, a good sensing performance can be obtained.In addition, with the similar expression of ED, FLOM based detector can be implemented practically when the parameter  is determined.

FLOM Based Spectrum Sensing with
Noise Power Estimation 4.1.FLOM Based Spectrum Sensing.As for the FLOM based detector, the structure is shown in Figure 1.According to the central limit theorem [25], when  is large enough, the metric Λ SO () in ( 10) can be approximated as a Gaussian random variable, 1,  ) (11) in which The th order central absolute moments of a Gaussian random variable  ∼ N(0,  2 ) can be expressed as where Γ(⋅) is the usual gamma function defined by Assuming that the noise () and the PU signal () are independent of each other, we can obtain that Then it can be easily concluded that The optimal value of  can be determined by the deflection coefficient which is usually defined as (16) to measure the performance for binary hypothesis testing problems [26].
The probability of false alarm  fa and probability of detection   can be given in terms of the  function by where (⋅) is defined by () = (1/ √ 2) ∫ ∞  exp(− 2 /2) and  is the threshold chosen for judging the presence or absence of the PU signal.
The threshold is normally chosen to satisfy a target  fa_DES according to the Neyman-Pearson criterion in practice [27], which requires the noise parameters to be known, as shown in  Substituting (18) into the second equation of ( 17), we can obtain the final probability of detection.
As shown in (19), increasing the number of the samples  can improve the sensing ability but will take longer time.So in practice, compromise should be made between sensing performance and processing time.

FLOM Detector with Noise Power Estimation.
As mentioned above, FLOM based sensing outperforms ED based methods.However, the noise power fluctuation is generally encountered in practice.So the noise power uncertainty should be considered to ensure the implementation of the sensing method.If the estimated noise power is assumed to be in an interval σ2 ∈ [(1/) 2  ,  2  ],  > 1 is a parameter that quantifies the amount of the uncertainty.To achieve a target  fa and   robustly, the following equations need to hold: From ( 21), we have Figure 3 displays how many samples are needed to achieve target  fa = 0.1 and   = 0.9 under different SNR with  =  = 0.1 (note that the vertical axis is a logarithmic scale).It is obvious that the number of samples which needed  will tend to ∞ when SNR approaches a specific value with any  > 1, which is called SNR Wall phenomenon.And the SNR wall will become higher with the uncertainty  becoming higher.It means that the detector cannot robustly detect the signal if the signal-noise power ratio is lower than the specific value.So the performance of an FLOM detector with a fixed threshold also degrades like ED when the noise power uncertainty occurs [19].
To keep a satisfactory sensing performance, the real time value of noise power should be used to determine the threshold in every sensing period.Fortunately, we found that the higher order moment of the samples can be a maximum likelihood (ML) noise power estimator, although it is not able to act as a good sensing statistic.In this paper, an adaptive threshold sensing structure is adopted with a higher order  moment detector as a real time noise power estimator and updating thresholds during each sensing process.From the structure of the NPE based detector shown in Figure 4, it can be seen that the estimator uses the same samples as the sensor and noise-only samples are not necessary in the noise estimation process.Without the noise-only slots, the transmitter can make full use of time slots and further improves the spectrum utilization.
The statistic of the noise power estimator is shown in in which  ≥ 2 is the higher order and  ≥  is the estimation time window.Due to the fact that the expectations of |()|  and |() + ()|  have very small differences at low SNR, this th order estimator can operate an unbiased estimation of the th order absolute moment of () no matter whether the PU signal is present or absent.The th order absolute moment of () is estimated by PW; that is, where  2  = ((/ + )/(1 + )) 2  is substituted.It is worth mentioning that the estimated noise power σ2  in the proposed estimator is not obtained directly when  > 2. But the value we need for calculating the threshold is σ , so estimating σ2 or σ  can both achieve this purpose.The only difference between them is the estimation accuracy which will be discussed in the following section.
Substituting (15) into (18) and replacing    by the estimated σ  in (24), the corresponding threshold can be obtained where Utilizing the threshold γ, the real probability of false alarm and probability of detection can be expressed as in which  0 (⋅) and  1 (⋅) denote the pdf of PW under  0 and  1 , respectively.The previous discussion indicates that PW also obeys Gaussian distribution with a large .So where Substituting ( 28) and ( 29) into (27), the value of  fa_ and  _ can be got.

Simulation and Results
In this section, we give the simulations and analysis of the theoretical results in Section 4. when SNR is as low as −15 dB.However, the FLOM detector shows large improvement and a satisfactory   can be easily achieved when an appropriate value for  is chosen,  = 0.4 as an example.The result also indicates that the performance improvement also depends on .The larger improvement will be obtained when  is getting smaller.In other words, the FLOM detector is more suitable for the circumstances under which the noise is more impulsive.From another perspective, Figure 6 shows several Receiver Operating Characteristic (ROC) curves of ED and FLOM based sensing in both theory and simulation with  =  = 0.1;  = 1000 at SNR = −15 dB.The simulation results are well matched with theory results.Obviously, FLOM based sensing outperforms ED based one and a better performance will be achieved when the value of  is smaller.
The detection ability of the FLOM detector under different SNR is plotted in Figure 7.The similar conclusion is drawn as those in Figures 5 and 6.In addition, take   = 0.9 as a target probability of detection; FLOM based sensing can work with SNR as low as −15 dB when  = 0.4, but ED based sensing will only work when SNR is not lower than −5 dB.

FLOM Detector with Noise Power Estimation.
The performance of the NPE based FLOM sensing approach must be related to the estimation accuracy, so the relationship between  _ ,  fa_ , and the order  of the estimator is given in From Figure 8, it can be seen that  _ and  fa_ are both increasing with different slopes as the value of  increases, but the change of  fa_ is slight.After  ≥ 3,  _ and  fa_ tend to be stable.Moreover, the order of the estimator can be set as  = 2 as a tradeoff to reduce computational complexity.Figure 9 demonstrates the theory and the simulation results of  _ and  fa_ versus SNR of our proposed FLOM detector.The threshold γ is calculated using σ  and the target probability of false alarm  fa_DES = 0.1,  = 1000,  =  = 0.1,  = 0.4, and  = 2.In these curves, ideal FLOM Mobile Information Systems represents the noise power perfectly known.It is obvious that a better sensing performance can be obtained when the estimation time window  is getting larger, while it also takes a longer processing time and the performance improvement will slow down when  is larger than 10.
It is worth mentioning that the number of observed samples has impact on the energy consumption and the sensing time.As an example, the proposed sensing structure works well at SNR = −15dB with the set of parameters ( =  = 0.1) using in the simulation.In practice, the noise power ratio  usually varies within [10 −4 , 0.1].Figures 5 and 7 indicate that the better performance will be obtained when  is getting smaller and SNR is getting larger.So in the case that SNR is not as low as −15 dB or the noise is more impulsive ( ≪ 0.1), the number of the observed samples could be drastically reduced.As a result, the energy consumption and the sensing time will be reduced and the proposed detector can be applied in IoT devices.

Conclusion
Cognitive Radio can be a helpful technology for utilizing and allocating frequency spectrum in the IoT.The utilization of the FLOM can successfully achieve spectrum sensing task for CR under Middleton Class A noise in the IoT.We derive the analytical expressions of the probability of false alarm  fa and the probability of detection   demonstrating the relationship between the performance enhancement and the value of the fractional order .For all circumstances, reasonable choices of  can be made to achieve a satisfactory performance using the deflection coefficient.After an NPE module being added to the sensing system, the problem of noise uncertainty at low SNR has been resolved.Both the power estimation time window  and the order  of the NPE module impact sensing performance.By choosing reasonable  and , simulation shows that our proposed sensing scheme has almost the same performance as that when the noise power is perfectly known.

Figure 2
Figure 2 demonstrates the distribution of Λ()|  0 and Λ()|  1 in ED based sensing and FLOM based sensing at SNR = −5 dB, respectively.Here we choose  =  = 0.1,  = 0.4, and  = 100 as an example.The two distribution curves in FLOM are more clearly separated than those of ED.It intuitively explained why the performance of FLOM based sensing is much better than ED based sensing under Middleton Class A noise.The optimal value of  can be determined by the deflection coefficient which is usually defined as(16) to measure the performance for binary hypothesis testing problems[26].

Figure 4 :
Figure 4: Structure of FLOM detector with noise power estimation module.