Next Article in Journal
A Comparative Study of Post-Quantum Cryptographic Algorithm Implementations for Secure and Efficient Energy Systems Monitoring
Next Article in Special Issue
An Innovative Design of Isoflux Scanning Digital Phased Array Based on Completely Shared Subarray Architecture for Geostationary Satellites
Previous Article in Journal
Multimodal Interaction Evaluation and Optimization Design of Smart Landscape Devices for Elderly People in Seaside Parks
Previous Article in Special Issue
A Low-RCS 2D Multi-Layer Van Atta Array at X-Band
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Multi-Antenna Spectrum Sensing Method Based on CEEMDAN Decomposition Combined with Wavelet Packet Analysis

1
School of Electrical & Information Engineering, Lanzhou University of Technology, Lanzhou 730050, China
2
Department of Computer Science and Computer Engineering, Universite de Technologie Belfort-Montbeliard, 90010 Belfort, France
*
Authors to whom correspondence should be addressed.
Electronics 2023, 12(18), 3823; https://doi.org/10.3390/electronics12183823
Submission received: 6 August 2023 / Revised: 6 September 2023 / Accepted: 7 September 2023 / Published: 9 September 2023
(This article belongs to the Special Issue Applications of Array Antenna in Modern Wireless Systems)

Abstract

:
In many practical communication environments, the presence of uncertain and hard-to-estimate noise poses significant challenges to cognitive radio spectrum sensing systems, especially when the noise distribution deviates from the Gaussian distribution. This paper introduces a cutting-edge multi-antenna spectrum sensing methodology that synergistically integrates complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN), wavelet packet analysis, and differential entropy. Signal feature extraction commences by employing CEEMDAN decomposition and wavelet packet analysis to denoise signals collected by secondary antenna users. Subsequently, the differential entropy of the preprocessed signal observations serves as the feature vector for spectrum sensing. The spectrum sensing module utilizes the SVM classification algorithm for training, while incorporating elite opposition-based learning and the sparrow search algorithm with genetic variation to determine optimal kernel function parameters. Following successful training, a decision function is derived, which can obviate the need for threshold derivation present in conventional spectrum sensing methods. Experimental validation of the proposed methodology is conducted and comprehensively analyzed, conclusively demonstrating its remarkable efficacy in enhancing spectrum sensing performance.

1. Introduction

In recent years, static spectrum allocation strategies struggled to meet the growing demand for spectrum resources. Within dedicated frequency bands, substantial portions of spectrum resources remain unused, both in terms of time and frequency domains [1]. While static allocation strategies effectively manage interference among communication systems operating in separate frequency bands, they exhibit inherent inflexibility. When primary users either refrain from utilizing their allocated spectrum for extended periods or use it intermittently, these idle frequency bands become inaccessible to other radio users. To address this contradiction, Dr. Mitola proposed the concept of cognitive radio (CR) [2]. As a new wireless communication technology, CR became an ideal solution to improve spectrum utilization, and spectrum sensing plays a crucial role in CR.
Several common spectrum sensing methods were proposed in references [3,4,5]. The energy detection algorithm compares the received signal energy with a threshold to determine the presence of a signal [3]. The matched filter detection correlates the received signal with a known transmitted signal to detect the existence of a primary user signal [4]. The cyclostationary feature detection adjusts the number of samples involved in real time based on changes in signal-to-noise ratio (SNR), analyzing the corresponding spectral correlation characteristics to determine the presence of a primary user [5]. Moreover, these detection methods have their limitations. The energy detection method is susceptible to interference from noise uncertainty and faces challenges in determining the decision threshold. Additionally, these detection algorithms are primarily evaluated under the assumption of additive Gaussian white noise for computational convenience.
In complex electromagnetic environments, the influence of interference and noise often follows a generalized Gaussian distribution (GGD), as observed in communication systems such as ultra-wide band (UWB) receivers [6]. At low SNR, the accuracy of spectrum sensing algorithms is diminished due to the uncertainty arising from noise. Consequently, this uncertainty affects the estimation of the signal by spectrum sensing algorithms. To address this issue, several traditional filtering denoising methods were proposed, including mean filtering, adaptive filtering, and wavelet analysis. However, these methods come with certain limitations. For instance, the effectiveness of adaptive filtering in noise suppression relies on the power spectral differences between the signal and the noise, making it less effective when the power spectra of the signal and noise are similar [7]. Empirical mode decomposition (EMD) decomposes signals into intrinsic mode functions (IMFs) without requiring the specification of wavelet basis functions. However, it struggles to handle non-stationary signals effectively [8]. Ensemble empirical mode decomposition (EEMD) was designed to address the mode mixing problem in EMD decomposition by adding white noise to the original signal for analysis [9]. Although the EEMD method solves the mixed mode problem, it is essentially a decomposition method that reduces the purity of the original signal, which will lead to the distortion of the signal. Furthermore, the residual noise introduced during decomposition contributes to increased signal reconstruction errors. To mitigate this limitation, Torres et al. introduced the complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) denoising method [10]. In the CEEMDAN method, adaptive white noise is added to each EMD stage, and its specific residual is calculated to obtain each modal component. Compared with the EEMD method, it overcomes the mode aliasing and effectively reduces the reconstruction error [11]. Although the CEEMDAN algorithm eliminates the mode aliasing and the reconstruction error is small, it will cause the loss of useful signal and incomplete echo signal in the high-frequency component.
To address the spectrum sensing problem, researchers employed machine learning classification algorithms to convert it into a binary classification problem, i.e., whether the primary user is using the licensed spectrum or not. Sarikhani et al. proposed a deep reinforcement learning method for spectrum sensing classification based on a collaborative spectrum sensing model, achieving promising results [12]. Thilina et al. adopted signal energy as a feature vector and applied the machine learning methods such as K-means, neural networks, and support vector machines (SVM), achieving effective spectrum sensing [13]. Tan et al. proposed a fusion algorithm based on random forests, incorporating sensor velocity, signal energy, and average eigenvalues of the covariance matrix as feature parameters, resulting in enhanced detection performance [14]. Wang et al. proposed a cooperative spectrum sensing method based on empirical mode decomposition and information geometry. They employed a K-medoids clustering algorithm for training in the spectrum sensing module, resulting in improved detection probability. However, under low SNR conditions, the detection probability significantly decreases [15]. Saravanan et al. used the differential entropy estimates in the received observations as feature vectors for the SVM algorithm spectrum sensing model training, but the spectrum sensing performance degrades rapidly as the shape parameter β increases [16]. Sheng et al. proposed a method to extract the trace of the covariance matrix and the quadratic covariance matrix as the combination of feature vectors for spectrum sensing training samples. They employed a genetic algorithm to optimize SVM parameters, achieving good sensing results, although the genetic algorithm optimization method is considered outdated [17].
In recent years, many scholars used heuristic algorithms to optimize the relevant parameters of machine learning algorithms, such as the sparrow search algorithm (SSA) [18]. Gao et al. used SSA to optimize the parameters of variational mode decomposition (VMD) to solve the problem that it is difficult to manually set the parameters to achieve the optimal VMD decomposition [19]. Li et al. used the sparrow search algorithm to optimize the parameters of the SVM and established a nonlinear relationship model between the melting process and technological indicators, which improved the model accuracy [20]. Qu et al. proposed a hybrid fault diagnosis method based on wavelet packet energy spectrum and SSA-SVM. The sparrow search algorithm was used to optimize the penalty factor and kernel parameters globally to improve the classification accuracy of the SVM [21]. However, the original sparrow search algorithm used in the above-related literature has the problems of low convergence accuracy, poor optimization efficiency, single search direction, and falling into local optimum.
Based on existing research, this paper proposes a multi-antenna spectrum sensing method named CEEMDAN-DE, which combines complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) decomposition and differential entropy (DE). Firstly, the signal collected by the secondary user is decomposed by CEEMDAN to obtain multiple IMF components. Then, the IMF components with high noise are screened out by calculating the correlation coefficient. Taking advantage of the fact that wavelet packet decomposition can provide more detailed signal frequency characteristics [22], the IMF component containing noise is denoised by wavelet packet analysis. Finally, the pure signal is obtained by component reconstruction. This joint noise removal method preserves the integrity and particularity of the original signal while eliminating the noise, and avoids the distortion problem and the loss of useful signals in high-frequency components caused by CEEMDAN. Finally, considering the shortcomings of the original sparrow search algorithm, the sparrow search algorithm combining elite opposition-based learning and genetic variation is established to determine the kernel function parameters of the SVM algorithm. This optimized SVM algorithm is employed for spectrum sensing. The experimental section includes the validation of the proposed method using receiver operating characteristic (ROC) curves, demonstrating its effectiveness.
The rest of this paper is organized as follows: Section 2 describes the multi-antenna spectrum sensing model and the CEEMDAN feature extraction method. Section 3 first introduces the optimization process of the sparrow search algorithm based on elite opposition-based learning and genetic variation, and then introduces the complete SVM spectrum sensing model. In Section 4, we substantiate the efficacy of the ISSA algorithm via the heuristic algorithm detection function. Then, the effectiveness of the CEEMDAN-wavelet packet analysis denoising method is verified under different SNRs. Next, the detection probability of the CEEMDAN-DE method under different SVM kernel function models is verified and compared. Under different SNRs, the superiority of the proposed method is proved by comparing the ROC curves of different feature extraction methods. Section 5 discusses the main work and contributions of this paper.

2. Multi-Antenna Spectrum Sensing in the Complex Communication Environment

In cognitive radio systems, users can be classified into two categories based on the usage priority of the licensed spectrum: primary users (PUs) and secondary users (SUs). PUs have absolute priority to use the licensed spectrum, while SUs continuously detect the presence of PUs through spectrum sensing, which ensures that SUs can access the PUs spectrum for communication tasks without causing interference with the normal communication of PUs. Therefore, an efficient and fast spectrum sensing technology is a prerequisite for building cognitive radio systems.
In comparison to single–antenna spectrum sensing, multi-antenna spectrum sensing offers numerous advantages [23]. It not only enhances detection performance by incorporating spatial information, but also reliably perceives primary user activities even in the lack of signal and noise variance information.

2.1. Spectrum Sensing Problem Model

The spectrum sensing problem is classified as a signal detection problem, the core of which is to determine whether the primary user signal exists. This can be effectively described using a binary hypothesis test:
H 0 :   Primary user signal does not exist H 1 :   Primary user signal exists   .
Consider a cognitive radio system consisting of a primary user with a single antenna and secondary users with M   ( M > 1 ) receiving antennas. The secondary users utilize their receiving antennas to detect whether the primary user is transmitting a signal. Under a binary hypothesis test, the received signals from the antennas can be represented as follows:
x m n = w m n ,     H 0 s n + w m n , H 1     n = 1,2 , , N .
In the scenario considered, we have two hypotheses: H 0 and H 1 . When H 0 is true, the channel solely consists of samples of noise denoted as w m n , which follows a generalized Gaussian distribution in a complex environment. In this case, there is no presence of a primary user signal s n , and the secondary user (SU) can utilize the channel. Conversely, when H 1 holds true, the channel encompasses both noise and the primary user’s signal, rendering the channel unavailable. N represents the number of observed samples.
For measuring the performance of a sensing system, common reference metrics include false alarm probability   P f   and detection probability P d . The false alarm probability P f represents the event probability that the fusion center misjudges the presence of the primary user when the primary user does not exist.
P f = P r o b H 1 H 0
The detection probability P d represents the probability that the primary user exists and is correctly detected.
P d = P r o b H 1 H 1

2.2. Feature Extraction Based on CEEMDAN

During the process of feature extraction, the observations collected by M   ( M > 1 ) receiving antennas are initially subjected to CEEMDAN combined with wavelet packet analysis for denoising, reducing noise, and redundant information in the acquired signals, thereby enhancing the overall performance of the sensing system. Subsequently, the estimated differential entropy (DE) values based on the observed data are utilized as feature vectors for signal feature extraction.

2.2.1. CEEMDAN Decomposition Principle

The EEMD and complete ensemble empirical mode decomposition (CEEMD) decomposition algorithms both address the mode mixing problem in EMD decomposition by adding Gaussian white noise or paired positive–negative Gaussian white noise to the signal to be decomposed [24]. However, neither of these algorithms isolate the residual noise, resulting in the added white noise signal being able to shift from high frequencies to low frequencies. As a result, a portion of the white noise signal remains in the decomposed intrinsic mode components, which can affect subsequent signal analysis and processing.
CEEMDAN decomposition addresses the above-mentioned issues in two aspects:
  • Adding the IMF component with auxiliary noise after EMD decomposition, rather than adding Gaussian white noise directly to the original signal;
  • Both EEMD and CEEMD methods employ a strategy of averaging the acquired mode components following empirical mode decomposition to tackle the mode mixing problem. However, CEEMDAN introduces a unique approach. It performs global averaging on the first-order IMF to obtain the final first-order IMF and then iterates this process on the residual component. This innovative approach effectively addresses the issue of noise transfer from high frequencies to low frequencies.
Let E i represent the i th IMF obtained after EMD. The i th IMF obtained through CEEMDAN is denoted as C i n ¯ , where v j represents Gaussian white noise following a standard normal distribution, ε denotes the standard deviation of the noise, and x m n represents the original signal received by the antenna. The steps of CEEMDAN decomposition are as follows:
A new signal x m n + 1 q ε v i n is obtained by adding plus–minus pairs of Gaussian white noise to the original signal x m n . The first order intrinsic mode component C 1 is obtained by EMD decomposition of the new signal:
x m n + 1 q ε v i n = C 1 j n + r j , q = 1   o r   2
Summing and averaging the resulting L modal components, the following is obtained.
C 1 n ¯ = 1 L j = 1 L C 1 j n
Calculate the residual error after removing first modal component:
r 1 n = x m n C 1 n ¯ .
Adding pairs of positive and negative Gaussian white noise to r 1 n is a new signal; with a new signal as the carrier for the EMD decomposition, we can obtain the phase 1 modal component D 1 j , thus the available modal component phase 2:
C 2 n ¯ = 1 L j = 1 L D 1 j n .
Calculate to remove the second residual modal component:
r 2 n = r 1 n C 2 n ¯ .
Repeat the above steps until the obtained residual signal is a monotone function and cannot be decomposed further, and the algorithm ends. At this time, if the number of eigenmode components obtained is K , then the original signal x m n is decomposed into:
x m n = k = 1 K C K n ¯ + r K n .

2.2.2. Wavelet Packet Analysis of Noise Reduction Principle

Wavelet packet analysis decomposes a noisy signal into multiple sub-band signals, allowing for better separation of different frequency components. By applying noise reduction techniques to these sub-bands, it becomes possible to selectively remove noise while preserving useful information in the signal. Compared to wavelet denoising, wavelet packet analysis offers higher frequency resolution, enabling further elimination of noise residues in the high-frequency portion.
Define the subspace U j n to be the closure space of the function U n t and U j 2 n to be the closure space of the function U 2 n t , such that U n t is satisfied,
U 2 n + 1 t = 2 K Ζ g k U n 2 t k   U 2 n t = 2 K Ζ q k U n 2 t k  
where q k and g k are orthogonal filter coefficients and g k = 1 k q 1 k . When u 0 = t , the sequence constructed by Equation   11   is called an orthogonal wavelet packet.
If g j n U j n , then g j n t can be expressed as:
g j n = l d l j , n U n 2 j t 1   .
The result of wavelet packet decomposition can be obtained:
d l j , 2 n = k a k 2 l d k j + 1 , n     d l j , 2 n + 1 = k s k 2 l d k j + 1 , n   .
Wavelet packet reconstruction for the inverse process of wavelet packet decomposition algorithm, by d l j , 2 n and d l j , 2 n + 1 construct d l j + 1 , n .
d l j + 1 , n = k q l 2 k d l j , 2 n + g l 2 k d l j , 2 n + 1  

2.2.3. CEEMDAN Decomposes Joint Wavelet Packet Analysis

As shown in Figure 1, the initial step involves utilizing the CEEMDAN algorithm to decompose the original signal into multiple intrinsic mode components. By comparing the correlation coefficient between each mode component and the original signal, the intrinsic mode components containing noise are screened out, and the spectral distribution characteristics and variance contribution rate of the intrinsic mode components are checked. Finally, the modal components containing noise are selected by wavelet packet analysis for threshold denoising, which avoids the problem that IMF components containing more noise are directly discarded, resulting in the loss of effective information. The signal x m ~ n is obtained by reconstructing the denoised modal component and the unprocessed modal component.

2.2.4. Feature Extraction Method Based on GGD Distribution

Generalized Gaussian noise finds relevance in various applications, including image processing, signal processing, and communication systems, among others. In these domains, signals are influenced by diverse factors, including weather-related interference, resulting in signals comprising different types and varying intensities of noise. Consequently, the noise distribution deviates from the conventional Gaussian distribution, assuming a more intricate non-Gaussian form.
If the M observed signal samples from the secondary users (SUs) are independently and identically distributed, they follow a generalized gaussian distribution (GGD). In this distribution, the parameter α   α > 0 controls the spread of the density function, while the parameter β   0 < β 2 determines the shape [25]. The probability density function of the GGD is expressed as shown in Equation (15):
f X x = 1 2 α Γ 1 β e x p x β α , x R .
Differential entropy (DE) is an extension of information entropy that measures the instantaneous entropy of a signal. It provides a measure of the speed of change and uncertainty in a signal. For a continuous random variable x with a probability density function   p ( x ) , the DE is defined as:
h X =   p x   ln p x   d x .
The DE feature of samples generated from GGD noise is given by its maximum likelihood estimation [25] and can be expressed as:
T C E E M D A N D E = 1 β log β 2 Γ 1 β + 1 β log β N i = 1 N x m ~ n X ^ β
where 0 < β 2 controls the “shape” of the generalized Gaussian distribution (GGD), which determines the rate of decay. X ^ represents the sample mean of the received observations and can be expressed as:
X ^ = 1 N n = 1 N x m ~ n .
Let X be a logarithmic random variable. According to the concept of geometric moments in [26], we define the geometric moment of X as follows:
S 0 = S 0 X = e E log X .
Suppose that x m ~ 1 , x m ~ 2 , , x m ~ n is an independent sample sequence from the geometric power distribution S 0 , then the statistical result can be expressed as:
T C E E M D A N G P = e x p 1 N n = 1 N log x m ~ n = n = 1 N x m ~ n 1 / N .
In multi-antenna spectrum sensing, the signal samples are received by the m th secondary user’s antenna and denoted as x m n   n = 1,2 , , N . Based on the feature extraction method, we compute the differential entropy and geometric power of x m ~ n . The results are denoted as T C E E M D A N D E , m and T C E E M D A N G P , m   m = 1 , 2 , , M . Consequently, two sets of feature vectors are obtained T C E E M D A N D E , 1 , , T C E E M D A N D E , M and T C E E M D A N G P , 1 , , T C E E M D A N G P , M . These feature vectors capture relevant information about the signal characteristics and are obtained through the CEEMDAN-DE and CEEMDAN-GP methods, respectively. They provide a comprehensive representation of the observed signals from different antennas in the multi-antenna spectrum sensing scenario. The feature extraction process based on CEEMDAN is shown in Figure 2.

3. Multi-Antenna Spectrum Sensing Based on ISSA-SVM

3.1. Improvement Sparrow Search Algorithm

Addressing the shortcomings of the original sparrow search algorithm (SSA) [15], we propose an improved sparrow search algorithm (ISSA) based on the integration of elite opposition-based learning and genetic variation. To overcome SSA limitations, including low convergence accuracy and poor optimization efficiency, we introduce the elite opposition-based learning mechanism. Furthermore, we address the issue of SSA’s single search direction by incorporating a multi-directional learning strategy. The genetic variation strategy enhances the algorithm’s capability to escape local optima by introducing genetic variations to the local optimal individual.
Elite opposition solution definition: Let x i t be a solution of the t iteration, and the reverse solution is x ¯ i t . When f x i t f x ¯ i t , x i t is called the elite individual of the t iteration; when f x i t f x ¯ i t , x ¯ i t is the elite individual of the t iteration, and x i j is the value of the ordinary individual x i on the j th dimension, then its reverse solution is
x ¯ i j t = m l i j t + u i j t x i j t
where m is a random number between 0 and 1, called the elite opposition coefficient. l i j t and u i j t are the minimum and maximum values of x ¯ i t on the j th dimension, and l i j , u i j are constructed intervals by elite groups.
In the traditional SSA, followers only learn from a single explorer in each iteration, which increases the likelihood of getting trapped in local optima and weakens the algorithm’s optimization ability. To address this issue, a multi-directional learning strategy is applied to the followers, offering relatively weaker individuals among the sparrows more opportunities to explore different regions [27]. With the improvement from the multi-directional learning strategy, the update formula for follower positions becomes as follows:
X i , j t + 1 = Q · e x p X w o r s t X i , j t i 2 ,   i > n 2   ω z · X z , j t + ω o · X o , j t + ω b · X b , j t ω z + ω o + ω b , otherwise .  
The equation can be expressed as follows, where X i j represents the position information of the i th sparrow in the j th dimension ( j = 1,2 , , i t e r m a x ), t denotes the current iteration count, Q is a random number following a normal distribution, X w o r s t represents the global worst position, and ω z , ω o , and ω b represent the weights assigned to the randomly selected sparrow individuals z ,   o , and b .
The iterative process of SSA aims to approach the optimal individual gradually. However, if the current individual becomes the local optimum, SSA is easily trapped in local optima. To overcome this issue, the mutation operator from the genetic algorithm is introduced into ISSA. This operator mutates the optimal individual within the sparrow population, enhancing the algorithm’s global search capability. Additionally, to avoid prolonged search times and emphasize the effectiveness of the mutation operator, an adaptive mutation probability is utilized in ISSA [27]. The mutation probability P m increases with the number of iterations. This adaptive approach helps individual sparrows evade local optima during the later iterations and improves the algorithm’s convergence accuracy. In ISSA, set s = 0.5   d = 1 ,   e = 0.03 ,   a n d   θ = 2 .
P m = e + s d i t e r m a x t i t e r m a x θ
The pseudo-code of the ISSA algorithm is shown in Algorithm 1.
Algorithm 1: ISSA pseudo-code.
Input: sparrow population G , the number of producers PD, proportion of scaredy birds SD, warning value ST, the maximum iterations i t e r m a x , w i   l b X i , u b X i
Output: sweet spot   X b e s t , optimal value   f g
1: Initialize the population using the elite opposition-based learning strategy
2: Calculate the fitness value to find the current best individual and the worst individual
3: While  t < i t e r m a x
4: R = r a n d 1
5: for  i = 1 : P D G
6: Update the discoverer position X i , j t X i , j t + 1
7: end for
8: for  i = P D G + 1 : G
9: Update the follower position according to Equation (22)
10: end for
11: for  i = 1 : S D
12: Update scout position
13: end for
14: Calculate the mutation probability P m , generate a random number p = r a n d   ( 1 )
15: if  p < P m
16: Mutation is performed according to X ˇ i   X i + w i obtain the new sparrow position
17: else
18: Keep the previous sparrow position X i
19: end if
20: Calculate the fitness value of the new position and the original position, and compare it
21: if f n e w < f o l d
22: Preserve current position
23: end if
24: t = t + 1
25: end while

3.2. Spectrum Sensing Based on ISSA-SVM

Since support vector machines (SVM) have certain advantages in nonlinear data processing [28], we use the SVM algorithm for spectrum sensing. The SVM classification problem is to find a hyperplane to maximize the classification interval on the feature space. The SVM mathematical model is as follows:
argmin w , b , ξ   1 2 w 2 + C i = 1 n ξ i s . t .   y i ω T x i + b 1 ξ i , ξ i 0 , i = 1 , , n
where w is the normal vector of the hyperplane, C is the penalty parameter, ξ i is the relaxation variable, x i is the eigenvector, y i is the result label, b is the threshold, and n is the number of samples.
Considering that spectral sensing data are often linear and indivisible, it is necessary to introduce kernel function to improve generalization ability. Therefore, the Gaussian kernel function with parameter σ is selected in this paper, and the expression is as follows:
K x 1 , x 2 = e x p σ x 1 x 2 2 .
When the classifier is successfully trained, the decision function is obtained according to the SVM algorithm as follows:
f x = s g n i = 1 D α i * y i K x 1 , x 2 + b * .
where s g n · is a symbolic function as shown in (27) and α i * is a Lagrange factor. After feature extraction of the classified signal sample, it is substituted into the decision function. If f x = 1 , it is considered that the primary user signal does not exist. If f x = 1 , the primary user signal is considered to exist.
s g n x = 1 ,   x < 0 0 ,   x = 0   1 ,   x > 0
The complete spectrum sensing algorithm is shown in Algorithm 2.
Algorithm 2: Spectrum sensing algorithm based on CEEMDAN-DE.
Input: Signal sequence x m n , number of observed samples n , number of data sets L
Output: P f and P d
     // CEEMDAN decomposition
1: i m f , its  ceemdan ( x m n )
2: for  i   = 1 : 14
3: Calculate the correlation coefficient of the i m f and the original signal
     C C ( i ) c o r r ( i m f ( i , : ) , x m n , t y p e , P e a r s o n )
4: end for
     // wavelet denoise and reconstruct
5: for  i   = 1 : 14 // the i m f component is selected for wavelet packet decomposition
6: x m n _ i m f   =   i m f [ i ]
7: w p t w p d e c ( x m n _ i m f , 3 , d b 8 ) // perform wavelet packet decomposition for each i m f [ i ] , using Daubechies wavelet ‘db8’ and decompose 3 layers
8: n o d e s = g e t ( w p t , t n )
9: N _ c f s = l e n g t h ( n o d e s )
10: f o r   j = 1 : N _ c f s
11: y s o f t ( : , j ) w t h r e s h ( r e x 3 ( : , j ) , s , t h r ) // soft threshold processing
12: end for
13: x m n _ i m f ( : , i ) = s u m ( y s o f t , 2 )
14: end for
15: x m ~ n = s u m ( x m n _ i m f , 2 ) + s u m ( i m f ( : , 6 : 8 ) , 2 ) //acquired reconstruction signal
     // feature extraction
16: for  i = 1 : L
17: A sample of n observations is randomly selected in x m ~ n each time
18: For x m ~ n , d a t a s e t i = T C E E M D A N D E , 1 , T C E E M D A N D E , 2 is obtained is obtained from Equation (17)
19: end
20: Divide d a t a s e t i into training sets and test sets
21: The training set is used for the SVM model parameter training, and ISSA is used to determine kernel function parameters C , σ
22: The decision function (26) is obtained after successful training
23: Input the test set according to the decision function 26
24: P f and P d are calculated according to Equations (3) and (4)

4. Simulation Results and Performance Analysis

The simulation environment for CEEMDEN-DE feature extraction is created using MATLAB2021a, and ROC curve performance analysis and ISSA algorithm verification are conducted using Python 3.8.9. The processor is an Intel(R) Core (TM) i7-8550U CPU @ 1.80 GHz processor, and the RAM is a 12 GB 64-bit computer.

4.1. ISSA Algorithm Optimization Ability Test

To assess the performance of the improved sparrow search algorithm (ISSA) and compare it with the sparrow search algorithm (SSA), we conducted performance testing using six typical benchmark functions shown in Table 1. The maximum number of iterations of the experiment i t e r m a x = 1000 , and the parameter settings [29] of the comparative experiment are shown in Table 2. The convergence curves are shown in Figure 3, Figure 4 and Figure 5. Through the analysis of the convergence curves for functions F 1 to F 6 , it can be observed that the proposed ISSA demonstrates faster convergence speed and better global search ability compared to SSA when solving complex problems.

4.2. Experimental Verification of Spectrum Sensing Method

Then, in this section, the simulated PU signal is obtained by superposition of sine and cosine functions. The simulation signal is represented as follows:
y 1 t = 0.35 × c o s 0.735 π × 50 t   y 2 t = 0.25 × s i n 2 π × 50 t × 1 + 1.5 s i n 0.55 π × 50 t y t = y 1 t + y 2 t + w t  
where t 0 ,   2 , the signal sampling interval t = 2 × 10 4   s ; and w t is the noise signal following the generalized Gaussian distribution.
The parameters related to the CEEMDAN method are set as follows: the standard deviation of the added Gaussian white noise N s t d = 0.2 , number of times noise is added N R = 100 , and the maximum number of iterations M a x I t e r = 500 . Simulation signal y t is decomposed by the CEEMDAN method under different signal-to-noise ratios (SNRs).
The simulation signals of the primary users are undergoing CEEMDAN decomposition under various SNRs, yielding a series of intrinsic mode function (IMF) components ranging from high to low frequencies, as demonstrated in Table 3. The correlation coefficients among these IMF components are being calculated using the “corr” function in MATLAB. As an example, for S N R = 12   d B , IMF1–IMF6 components with correlation coefficients exceeding 0.2 are categorized as high-frequency components with significant noise content. Subsequently, these components undergo denoising through wavelet packet analysis. Conversely, IMF7–IMF14 components with correlation coefficients less than 0.2 are identified as low-noise components and are preserved.
Subsequently, the selected high-frequency noisy components are denoised using the wavelet packet analysis method. The denoising process utilizes the db8 wavelet basis function, a soft thresholding function, and a decomposition level of 3. After denoising, the high-frequency IMF components, now free of noise, are combined with the initially low-noise low-frequency IMF components. This reconstruction step finalizes the CEEMDAN-wavelet packet thresholding composite denoising of the primary user signal y t .
In order to more fully reflect the difference of the primary user signal y t before and after noise reduction, two indicators of continuous mean square error and signal-to-noise ratio are introduced to evaluate the signal after noise reduction.
  • Root mean squared error (RMSE): in the signal noise reduction metric, the root mean square error is defined as the expected value of the squared difference between the un-denoised signal and the denoised signal recalculated as shown in Equation (29).
R M S E = 1 N i = 1 N y i t y i * t 2   i = 1,2 , n
  • SNR: The signal-to-noise ratio is defined as shown in Equation (30).
    S N R = 10 × l g i = 1 N y i t 2 i = 1 N y i t y i * t 2 d B
    where y i * t   is the denoised signal and n is the signal length.
To comprehensively demonstrate the superiority of the CEEMDAN-wavelet packet analysis denoising method, we conducted a comparative analysis, contrasting the performance of the proposed method with that of the EMD-wavelet packet analysis denoising and conventional EEMD denoising methods across various SNRs. The denoising efficacy for both methods is evaluated using the RMSE and SNR.
In the wavelet packet threshold denoising experiment, db8 wavelet basis function and soft threshold function are used again, and the decomposition level is still set to 3.
From the comparative results in Table 4 and the evaluation criteria for denoising performance, it is evident that the proposed CEEMDAN-wavelet packet analysis-combined denoising method outperforms both the EMD-wavelet packet analysis-combined denoising method and the EEMD denoising method for the primary user signal. Additionally, the RMSE values for the CEEMDAN-wavelet packet analysis method are consistently lower than those of the other two methods. Moreover, the SNR achieved by the CEEMDAN-wavelet packet analysis-combined denoising method is notably higher than that achieved by the other two methods.
In the training of the spectrum sensing module using the SVM classification algorithm, as presented in literature [13] and literature [16], only the “linear” and “polynomial” kernel functions are considered. Regrettably, the more appropriate “RBF” kernel function, which is well suited for addressing nonlinear classification problems, is overlooked. These methods are denoted as linear-SVM, poly-SVM, and RBF-SVM, respectively. It is noteworthy that the “RBF” kernel function introduces an element of uncertainty due to the random generation of penalty factors C and kernel function parameters σ during each training session. Therefore, the improved sparrow search algorithm (ISSA) is used to optimize the RBF kernel parameters   C and σ , which are referred to as ISSA-SVM. To comprehensively showcase the superiority of the CEEMDAN-DE-ISSA-SVM algorithm in spectrum sensing classification, it is compared with two methods from the literature.
During the simulation, the parameters are configured as follows: β = 1.8 and α = 5.2 , where α controlled the SNR of the generated signals. Each training sample consisted of N = 100 sampling points from the secondary user’s received signal. The chosen signal feature extraction methods encompassed CEEMDAN-DE and DE [16]. Based on the size of the generated data set, the average training time of the machine-learning-based multi-antenna spectrum sensing methods is presented in Table 5. Among the data samples, normalization is performed and 50% is selected as the training set for the SVM model and the remaining 50% as the test set for the SVM model. The relevant parameter settings when the ISSA algorithm optimizes the SVM kernel function are shown in Table 6.
As can be seen from Table 5, among the three different kernel functions, the training time of the CEEMDAN-DE-SVM model under various data sets is shorter than that of the DE-SVM model in the literature. This observation validates that the proposed denoising method in this paper aids SVM in swiftly identifying an appropriate hyperplane for data partitioning. Moreover, it is observable that the “linear” and “polynomial” kernel functions result in low detection probabilities P d when applied to small and non-denoised data sets. In contrast, ISSA-SVM maintains normal detection probabilities even in such cases.
Among the three kernel functions of the CEEMDAN-DE-SVM model, CEEMDAN-DE-poly-SVM has the longest training time and CEEMDAN-DE-Linear-SVM has the shortest training time, but the detection probability is low. The CEEMDAN-DE-ISSA-SVM algorithm has the highest classification accuracy and detection probability. Its increased computational complexity is offset by its superior classification accuracy and detection probability. Then, under the same α and N conditions, we analyzed and compared the impact of β   change on the detection probability.
We then proceeded to extract 2000 signal features using the CEEMDAN-DE feature extraction method and conducted a comparative analysis of the detection probabilities between the ISSA-SVM and RBF-SVM algorithms. This analysis is performed under various β conditions while keeping α and N constant. As shown in Table 7, the detection probability of the RBF-SVM algorithm is significantly lower than that of the ISSA-SVM under different β . It can be found that the suitability of C and σ values have a great impact on the detection probability performance of the SVM spectrum sensing model.
All the subsequent experiments are implemented based on the ISSA-SVM, and the relevant parameters of the ISSA algorithm are shown in Table 6. To ensure the accuracy of the experiment, we obtained 2000 signal features based on the CEEMDAN-DE feature extraction method and normalized them, among which, 1000 are used as the training set and 1000 are used as the test set. A test set is used to verify the spectrum sensing performance of the proposed method.
ROC curves are used to compare the performance of different feature extraction methods. Among them, DE and GP represent the feature extraction methods mentioned in references [16,26]. EMD-DE and EMD-GP denote the feature extraction method using the EMD-joint wavelet packet analysis denoising method, and ED is the energy detection algorithm.
Figure 6 shows the ROC curves of different feature extraction methods when M = 2 ,   β = 0.8 ,   α = 1 , and the number of sampling points N = 25 for each training sample of the secondary user’s received signal. As indicated in Table 8, it can be observed that the detection probability based on the CEEDMAN-DE algorithm, at P f = 0.1 , the CEEMDAN-DE method improved by 355% compared to the original DE method, and by 95% compared to the GP method. It can be observed that, under the same N , the multi-antenna spectrum sensing algorithm based on SVM with the CEEMDAN-DE feature extraction method proposed in this paper achieves the best performance.
Figure 7 illustrates the classification results obtained through feature extraction using CEEMDAN-DE. The red “+” symbols represent the feature vectors corresponding to the presence of primary user signals, while the blue “o” symbols represent feature vectors corresponding to the absence of primary user signals. The black curve represents the classification hyperplane derived after training the SVM. It is evident that the classification hyperplane successfully separates the samples into two distinct categories, with only a small number of overlapping data points.
Figure 8 shows the ROC curves for different feature extraction methods when M = 2 ,   β = 1.7 ,   α = 3.2 ,   S N R = 14   d B , and N = 100 . As indicated in Table 9, it can be observed that the detection probability based on the CEEDMAN-DE algorithm, at P f = 0.1 , exhibits a significant improvement compared to the performance of the DE method mentioned in the literature, with an increase of 588%. Furthermore, compared to the EMD-DE method, the CEEDMAN-DE algorithm demonstrates a performance improvement of 17%. It can be observed that the spectrum sensing algorithms based on DE and GP mentioned in the literature significantly decrease the spectrum sensing performance when 1 < β < 2 , but the method proposed in this paper can still maintain a high spectrum sensing performance.
In Table 10, we compare the classification performance of the ISSA-SVM spectrum sensing algorithm with K-means, decision tree, and random forest spectrum sensing algorithms when S N R = 16   d B . The metrics are area under the curve (AUC) and classification accuracy rate. Where AUC is the area under the ROC curve, which is used to measure the performance of a classifier. A higher AUC value indicates better performance of the classifier, while a lower AUC value indicates poorer performance.
We can see from Table 10 that the ISSA-SVM spectrum sensing algorithm can maintain the highest classification performance regardless of 0 < β < 1 or 1 < β < 2 . Other algorithms, such as the K-means algorithm, even fail in classification performance.
Figure 9 presents the ROC curves for the SVM-based multi-antenna spectrum sensing method under generalized Gaussian distribution noise, with varying numbers of sampling points N . As shown in the figure, under the conditions of M = 2 ,   β = 0.8 ,   a n d   α = 1 , with an increase in the number of sampling points N for the secondary user, the spectral sensing performance of all algorithms exhibits noticeable improvement. However, the performance of these algorithms remains lower than that of the CEEDMAN-DE method proposed in this study. This validates the effectiveness of the proposed method in scenarios with a small number of sampling points.

5. Conclusions

This paper proposes a CEEDMAN-DE method to effectively improve spectrum sensing performance in complex electromagnetic environments. In terms of feature extraction, the proposed method first utilizes CEMMDAN decomposition combined with wavelet packet analysis to denoise the signals collected by the secondary user (SU). This process enables finer removal of noise components and avoids the loss of primary user information. Then, the denoised signals are subjected to differential entropy feature extraction to obtain features for spectrum sensing. Finally, in the spectrum sensing module, the SVM classification algorithm is employed for training, combined with elite opposition-based learning and a sparrow search algorithm with genetic variation, to determine the optimal kernel function parameters. In the experiments, the performance of the CEEDMAN-DE algorithm is validated and analyzed through ROC curve analysis. When 0 < β < 1 , the proposed method ensures accurate spectrum sensing even with a small value of N, outperforming the spectrum sensing methods in the literature. Moreover, when 1 < β 2 the performance of the spectrum sensing methods in the literature deteriorates, the CEEDMAN-DE algorithm maintains accurate spectrum sensing classification. In our future work, we plan to explore a spectrum sensing algorithm that combines CEEMDAN-DE feature extraction with deep learning techniques. Additionally, we aim to investigate methods for reducing algorithmic complexity to achieve more efficient spectrum sensing.

Author Contributions

Y.H. researched the literatures, provided the mathematical models, designed the algorithm, completed numerical simulations, took charge of the original draft preparation, and edited the manuscript; S.L. conceived the study concepts, improved the systematic research and analysis methodology, and supervised the completion of the refinement of the paper; J.G., S.Y. and Q.Y. gave valuable suggestions for revision, checked formula deducing and English grammar. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the National Natural Science Foundation of China (61663024), in part by the Hongliu First Class Discipline Development Project of Lanzhou University of Technology (25–225305).

Data Availability Statement

The data supporting this article are from previously reported studies and datasets, which have been cited.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Joykutty, A.M.; Baranidharan, B. Cognitive Radio Networks: Recent Advances in Spectrum Sensing Techniques and Security. In Proceedings of the 2020 International Conference on Smart Electronics and Communication, Trichy, India, 10–12 September 2020. [Google Scholar] [CrossRef]
  2. Mitola, J. The Software Radio Architecture. IEEE Commun. Mag. 1995, 33, 26–38. [Google Scholar] [CrossRef]
  3. Bhavana, B.; Sabat, S.L.; Namburu, S.; Panigrahi, T. Energy Detector for Spectrum Sensing Using Robust Statistics in Non-Gaussian Noise Environment. In Proceedings of the 2023 15th International Conference on Communication Systems & Networks, Bangalore, India, 3–8 January 2023. [Google Scholar] [CrossRef]
  4. Kumar, R.V.R. Performance Analysis of a Novel Matched Filter Detector with Reduced out-of-Band Response for Narrowband Signals. IEEE Trans. Signal Process. 2006, 54, 2691–2703. [Google Scholar] [CrossRef]
  5. Sepidband, P.; Entesari, K. A CMOS Spectrum Sensor Based on Quasi-Cyclostationary Feature Detection for Cognitive Radios. IEEE Trans. Microw. Theory Tech. 2015, 63, 4098–4109. [Google Scholar] [CrossRef]
  6. Ahmed, Q.Z.; Park, K.H.; Alouini, M.S. Ultrawide Bandwidth Receiver Based on a Multivariate Generalized Gaussian Distribution. IEEE Trans. Wirel. Commun. 2014, 14, 1800–1810. [Google Scholar] [CrossRef]
  7. Habibi, Z.; Zayyani, H. Markovian Adaptive Filtering Algorithm for Block-Sparse System Identification. IEEE Trans. Circuits Syst. II Express Briefs 2021, 68, 3032–3036. [Google Scholar] [CrossRef]
  8. Wang, Y.H.; Cheng, S.H. Boundary Effects for EMD-Based Algorithms. IEEE Signal Process. Lett. 2022, 29, 1032–1036. [Google Scholar] [CrossRef]
  9. Xue, W.; Dai, X.; Zhu, J.; Luo, Y.; Yang, Y. A Noise Suppression Method of Ground Penetrating Radar Based on EEMD and Permutation Entropy. IEEE Geosci. Remote Sens. Lett. 2019, 16, 1625–1629. [Google Scholar] [CrossRef]
  10. Torres, M.E.; Colominas, M.A.; Schlotthauer, G.; Flandrin, P. A Complete Ensemble Empirical Mode Decomposition with Adaptive Noise. In Proceedings of the 2011 IEEE International Conference on Acoustics, Speech and Signal Processing, Prague, Czech Republic, 22–27 May 2011. [Google Scholar] [CrossRef]
  11. He, C.; Niu, P.; Yang, R.; Wang, C.; Li, Z.; Li, H. Incipient Rolling Element Bearing Weak Fault Feature Extraction Based on Adaptive Second-order Stochastic Resonance Incorporated by Mode Decomposition. Measurement 2019, 145, 687–701. [Google Scholar] [CrossRef]
  12. Sarikhani, R.; Keynia, F. Cooperative Spectrum Sensing Meets Machine Learning: Deep Reinforcement Learning Approach. IEEE Commun. Lett. 2020, 24, 1459–1462. [Google Scholar] [CrossRef]
  13. Thilina, K.M.; Choi, K.W.; Saquib, N.; Hossain, E. Machine Learning Techniques for Cooperative Spectrum Sensing in Cognitive Radio Networks. IEEE J. Sel. Areas Commun. 2013, 32, 2209–2221. [Google Scholar] [CrossRef]
  14. Tan, C.; Chen, J.; Chen, S.; Li, C.; Liu, H.; Zheng, M. Combination Spectrum Sensing Algorithm for Wireless Sensor Network Based on Random Forest. In Proceedings of the 2022 10th International Conference on Intelligent Computing and Wireless Optical Communications, Chongqing, China, 10–12 June 2022. [Google Scholar] [CrossRef]
  15. Wang, Y.; Zhang, S.; Zhang, Y.; Wan, P.; Li, J.; Li, N. A Cooperative Spectrum Sensing Method Based on Empirical Mode Decomposition and Information Geometry in Complex Electromagnetic Environment. Complexity 2019, 2019, 15. [Google Scholar] [CrossRef]
  16. Saravanan, P.; Chandra, S.S.; Upadhye, A.; Gurugopinath, S. A Supervised Learning Approach for Differential Entropy Feature-based Spectrum Sensing. In Proceedings of the 2021 Sixth International Conference on Wireless Communications, Signal Processing and Networking, Chennai, India, 25–27 March 2021. [Google Scholar] [CrossRef]
  17. Sheng, J.; Liu, F.; Zhang, Z.; Huang, C. TCVQ-SVM Algorithm for Narrowband Spectrum Sensing. Phys. Commun. 2022, 50, 101502. [Google Scholar] [CrossRef]
  18. Xue, J.; Shen, B. A Novel Swarm Intelligence Optimization Approach: Sparrow Search Algorithm. Syst. Sci. Control Eng. 2020, 8, 22–34. [Google Scholar] [CrossRef]
  19. Gao, X.; Guo, W.; Mei, C.; Sha, J.; Guo, Y.; Sun, H. Short-term Wind Power Forecasting Based on SSA-VMD-LSTM. Energy Rep. 2023, 9, 335–344. [Google Scholar] [CrossRef]
  20. Li, X.; Yang, C.; Zhang, X.; Li, G. Prediction of Matte Grade in The Oxygen-rich Top Blown Smelting Based on WD-SSA-SVM Algorithm. In Proceedings of the 2023 IEEE 12th Data Driven Control and Learning Systems Conference, Xiangtan, China, 12–14 May 2023. [Google Scholar] [CrossRef]
  21. Qu, J.; Ma, B.; Ma, X.; Wang, M. Hybrid Fault Diagnosis Method based on Wavelet Packet Energy Spectrum and SSA-SVM. Int. J. Adv. Comput. Sci. Appl. 2022, 13, 52–60. [Google Scholar] [CrossRef]
  22. Li, L.; Wang, F.; Shang, F.; Jia, Y.; Zhao, C.; Kong, D. Energy Spectrum Analysis of Blast Waves Based on an Improved Hilbert–Huang Transform. Shock Waves 2017, 27, 487–494. [Google Scholar] [CrossRef]
  23. Liu, M.; Zhang, X.; Chen, Y.; Tan, H. Multi-antenna Intelligent Spectrum Sensing in The Presence of Non-Gaussian Interference. Digit. Signal Process. 2023, 140, 104135. [Google Scholar] [CrossRef]
  24. Yeh, J.R.; Shieh, J.S.; Huang, N.E. Complementary Ensemble Empirical Mode Decomposition: A Novel Noise Enhanced Data Analysis Method. Adv. Adapt. Data Anal. 2010, 2, 135–156. [Google Scholar] [CrossRef]
  25. Gurugopinath, S.; Muralishankar, R.; Shankar, H.N. Differential Entropy-Driven Spectrum Sensing Under Generalized Gaussian Noise. IEEE Commun. Lett. 2016, 20, 1321–1324. [Google Scholar] [CrossRef]
  26. Gonzalez, J.G.; Paredes, J.L.; Arce, G.R. Zero-Order Statistics: A Mathematical Framework for the Processing and Characterization of Very Impulsive Signals. IEEE Trans. Signal Process. 2006, 54, 3839–3851. [Google Scholar] [CrossRef]
  27. Chai, Y.; Sun, X.; Ren, S. Chaotic Sparrow Search Algorithm Based on Multi-Directional Learning. Comput. Eng. Appl. 2023, 59, 81–91. [Google Scholar] [CrossRef]
  28. Akram, Z.; Fernández, M.; Cernadas, E.; Barro, S. Fast Support Vector Classification for Large-Scale Problems. IEEE Trans. Pattern Anal. Mach. Intell. 2021, 44, 6184–6195. [Google Scholar] [CrossRef] [PubMed]
  29. Li, Z.; Guo, J.; Gao, X.; Yang, X.; He, Y.L. A Multi-strategy Improved Sparrow Search Algorithm of Large-scale Refrigeration System: Optimal Loading Distribution of Chillers. Appl. Energy 2023, 349, 121623. [Google Scholar] [CrossRef]
Figure 1. CEEMDAN-wavelet packet analysis joint noise reduction.
Figure 1. CEEMDAN-wavelet packet analysis joint noise reduction.
Electronics 12 03823 g001
Figure 2. Feature extraction process based on CEEMDAN.
Figure 2. Feature extraction process based on CEEMDAN.
Electronics 12 03823 g002
Figure 3. (a) The F1 test function optimizes the contrast curve; (b) the F2 test function optimizes the contrast curve.
Figure 3. (a) The F1 test function optimizes the contrast curve; (b) the F2 test function optimizes the contrast curve.
Electronics 12 03823 g003
Figure 4. (a) The F3 test function optimizes the contrast curve; (b) the F4 test function optimizes the contrast curve.
Figure 4. (a) The F3 test function optimizes the contrast curve; (b) the F4 test function optimizes the contrast curve.
Electronics 12 03823 g004
Figure 5. (a) The F5 test function optimizes the contrast curve; (b) the F6 test function optimizes the contrast curve.
Figure 5. (a) The F5 test function optimizes the contrast curve; (b) the F6 test function optimizes the contrast curve.
Electronics 12 03823 g005
Figure 6. ROC curves of different feature extraction methods at SNR = −9 dB.
Figure 6. ROC curves of different feature extraction methods at SNR = −9 dB.
Electronics 12 03823 g006
Figure 7. Classification effect diagram.
Figure 7. Classification effect diagram.
Electronics 12 03823 g007
Figure 8. ROC curves of different feature extraction methods at SNR = −14 dB.
Figure 8. ROC curves of different feature extraction methods at SNR = −14 dB.
Electronics 12 03823 g008
Figure 9. ROC curves for different number of sampling points N.
Figure 9. ROC curves for different number of sampling points N.
Electronics 12 03823 g009
Table 1. Test functions.
Table 1. Test functions.
FunctionsFunction NameDimension DomainOptimal Value
F 1 Sphere30 [ 100,100 ] 0
F 2 Schwefel’s Problem 2.2230 [ 10,10 ] 0
F 3 Ackley30 [ 32,32 ] 0
F 4 Schwefel’s Problem 2.2130 [ 100,100 ] 0
F 5 Griewank30 [ 600,600 ] 0
F 6 Rastrigin30 [ 5.12,5.12 ] 0
Table 2. Parameters setting of comparison algorithm.
Table 2. Parameters setting of comparison algorithm.
AlgorithmGPDSDST
SSA300.20.10.8
ISSA300.20.10.8
Table 3. Correlation coefficient value corresponding to each order component of the primary user signal with different SNRs.
Table 3. Correlation coefficient value corresponding to each order component of the primary user signal with different SNRs.
SNR = −9 dBSNR = −12 dB
IMF10.60940.6271
IMF20.42710.4541
IMF30.34650.3659
IMF40.27660.2907
IMF50.27530.2271
IMF60.25010.2023
IMF70.22030.1886
IMF80.22910.1881
IMF90.15990.1121
IMF100.01940.0375
IMF110.03660.0146
IMF120.03420.0161
IMF130.01980.0143
IMF140.01640.0151
Table 4. Performance indicators of different noise reduction methods.
Table 4. Performance indicators of different noise reduction methods.
Initial SNREvaluation IndexEMD-Wavelet Packet Analysis-Combined DenoisingEEMDCEEMDAN-Wavelet Packet Analysis-Combined Denoising
SNR = −9 dB RMSE0.23830.28070.1968
SNR/dB3.5322.10954.9743
SNR = −12 dBRMSE0.31540.30390.2644
SNR/dB1.09871.42892.6297
Table 5. Comparison of SVM kernel function models under different feature extraction methods when SNR = −16 dB.
Table 5. Comparison of SVM kernel function models under different feature extraction methods when SNR = −16 dB.
Kernel Function Model CategoryExperimental Data Set SizeTraining Time (s) P d P f = 0.1 A c c u r a c y   R a t e
CEEMDAN-DE-linear-SVM20000.5623370.9340.90
30000.4998870.9520.91
CEEMDAN-DE-poly-SVM20003.8218290.9360.931
30003.4019490.9550.937
CEEMDAN-DE-ISSA-SVM20000.8592160.940.93225
30000.7963240.960.9405
DE-linear-SVM20001.0310030.090.7545
30001.3278540.1250.68925
DE -poly-SVM200017.3393490.0060.7545
300015.2062000.020.68925
DE-ISSA-SVM20001.5774600.790.83575
30002.1522440.6960.769
Table 6. ISSA related parameter settings.
Table 6. ISSA related parameter settings.
G i t e r m a x Warning Value PDSD
30300.60.70.2
Table 7. Comparison of detection probability between ISSA-SVM and RBF-SVM under different β .
Table 7. Comparison of detection probability between ISSA-SVM and RBF-SVM under different β .
Model Category β P d P f = 0.1 Penalty   Factors   ( C ) Kernel   Function   Parameters   ( σ )
ISSA-SVM0.80.9635.207797484.23560914
10.9690.014.08559084
1.50.9390.871509950.08230669
20.8930.106259310.90453521
RBF-SVM0.80.96111.010.5030644
10.87111.010.5069423
1.50.86911.010.4894858
20.73411.010.4988999
Table 8. AUC and P d values of each feature extraction method when S N R = 9   d B .
Table 8. AUC and P d values of each feature extraction method when S N R = 9   d B .
AlgorithmsAUC P d P f = 0.1 Penalty   Factors   ( C ) Kernel   Function   Parameters   ( σ )
CEEMDAN-DE0.9474610.9790.050190980.35761508
CEEMDAN-GP0.8300640.8010.285401710.88132311
EMD-DE0.8343990.7910.949712811.63474901
EMD-GP0.8627920.8096.963368710.67814947
DE0.6662860.2156.945636246.78160817
GP0.7737800.5014.982294435.49424544
Table 9. AUC and P d values of each feature extraction method when S N R = 14   d B .
Table 9. AUC and P d values of each feature extraction method when S N R = 14   d B .
AlgorithmsAUC P d P f = 0.1 Penalty   Factors   ( C ) Kernel   Function   Parameters   ( σ )
CEEMDAN-DE0.9626330.98395.0093994310
CEEMDAN-GP0.8826640.8748.2303817919.5211308
EMD-DE0.8716920.8390.292435046.28011474
EMD-GP0.6795000.5072.593375231.24635901
DE0.5405750.1435.929159855.31179191
GP0.5076790.14359.831491239.78577177
Table 10. Classification performance comparison of spectrum sensing algorithms when SNR = −16 dB.
Table 10. Classification performance comparison of spectrum sensing algorithms when SNR = −16 dB.
Spectrum Sensing Algorithm β AUCAccuracy Rate
ISSA-SVM0.80.9118710.93225
1.50.9188630.90075
Random forest0.80.9071700.924
1.50.8149990.789
K-means0.80.6004180.495
1.50.6778110.7595
Decision tree0.80.889720.907
1.50.7992650.77025
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, S.; Han, Y.; Gaber, J.; Yang, S.; Yang, Q. A Multi-Antenna Spectrum Sensing Method Based on CEEMDAN Decomposition Combined with Wavelet Packet Analysis. Electronics 2023, 12, 3823. https://doi.org/10.3390/electronics12183823

AMA Style

Li S, Han Y, Gaber J, Yang S, Yang Q. A Multi-Antenna Spectrum Sensing Method Based on CEEMDAN Decomposition Combined with Wavelet Packet Analysis. Electronics. 2023; 12(18):3823. https://doi.org/10.3390/electronics12183823

Chicago/Turabian Style

Li, Suoping, Yuzhou Han, Jaafar Gaber, Sa Yang, and Qian Yang. 2023. "A Multi-Antenna Spectrum Sensing Method Based on CEEMDAN Decomposition Combined with Wavelet Packet Analysis" Electronics 12, no. 18: 3823. https://doi.org/10.3390/electronics12183823

APA Style

Li, S., Han, Y., Gaber, J., Yang, S., & Yang, Q. (2023). A Multi-Antenna Spectrum Sensing Method Based on CEEMDAN Decomposition Combined with Wavelet Packet Analysis. Electronics, 12(18), 3823. https://doi.org/10.3390/electronics12183823

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop