A Spectrum Sensing Method Based on Empirical Mode Decomposition and K-Means Clustering Algorithm
To solve the problems of poor performance of traditional spectrum sensing method under low signal-to-noise ratio, a new spectrum sensing method based on Empirical Mode Decomposition algorithm and K-means clustering algorithm is proposed. Firstly, the Empirical Mode Decomposition algorithm and the wavelet threshold algorithm are used to remove the noise components in the spectrum sensing signal, and K-means clustering algorithm is used to determine whether the primary user exists. The method can remove the redundant components such as noise in the nonstationary or nonlinear sampling signal in the real environment and does not need to know the prior information such as signal, channel, and noise, so it can well handle the complicated sensing signal in real environment. This method can reduce the impact of noise on the spectrum sensing system and thus can improve the sensing performance of the system. In the experimental part, the difference between maximum and minimum eigenvalues and the difference between the maximum eigenvalue and the average energy in the random matrix are selected as signal features. Experiments also show that the proposed method is better than the traditional spectrum sensing methods.
Spectrum sensing technology is the key technology of cognitive radio. The main purpose is to detect spectrum holes and improve spectrum utilization. Conventional spectrum sensing methods include energy detection, matched filter, and the cyclostationary feature detection [1–4]. The perceived signals received in the actual environment include noise, which can affect the sensing performance of traditional spectrum sensing methods [5, 6].
Spectrum sensing based on random matrix is a cooperative spectrum sensing method. This method does not require prior knowledge of sampling signals and can reduce the interference of various fading factors in the environment. Therefore, the method of random matrix theory is gradually applied to spectrum sensing. Literature  has extracted a perception method based on the Difference between maximum and minimum eigenvalue (DMM). This method uses DMM as the test statistic and then makes a decision by calculating the threshold. However, the algorithm is susceptible to low signal-to-noise ratio (SNR) and few collaborative users, which will greatly affect the system performance. The literature  proposed the Maximum to Minimum Eigenvalue (MME) spectrum sensing method. The method uses the MME eigenvalue as the test statistic to reduce the influence of noise uncertainty, but only when the number of samples is high does it have good sensing performance. Literature  proposed a spectrum sensing method based on the Difference between Maximum Eigenvalue and the Average Energy (DMEAE). This method uses DMEAE eigenvalue as the test statistic and derive the corresponding threshold to make a decision. Since the statistical covariance of the received signal and noise are usually different, they can be used to distinguish between the presence of the primary user signal and the presence of noise only. However, the traditional method of random matrix theory will have the problem of inaccurate estimation of the threshold. At the same time, in the complex real environment, the SNR of the primary user signal is low and susceptible to the multipath fading and time diffusion of the wireless channel, which will greatly affect the system performance.
In order to solve the problem of inaccurate threshold derivation in traditional stochastic matrix theory, the literature  proposed spectrum sensing method based on signal eigenvalue and clustering algorithm. This method uses a K-means clustering algorithm and signal eigenvalues to train a classifier and use the classifier to determine whether the primary user signal exists. Spectrum sensing can also be considered as a problem of two classifications (the existence of a primary user or the absence of a primary user). Because the clustering algorithm can deal with the two classification problems well, the spectrum sensing method based on machine learning has gradually become a focus of people’s research. Literature  proposed a spectrum sensing method based on K-means clustering. This method takes the energy of the signal as a feature and then divides this feature into a channel available class or a channel unavailable class through a K-means clustering algorithm. Literature  proposed a spectrum sensing method combining support vector machine and MME. The literature  analyzes the spectrum sensing performance under different clustering algorithms. Compared with the traditional spectrum sensing method, the spectrum sensing method based on machine learning is more adaptive and does not need to know the prior knowledge of the sensing environment. However, taking into account the actual environment, the signal often contains complex uncertain noises which will affect the effectiveness of the clustering algorithm, making the spectrum sensing system have poor clustering performance at low SNR and it will affect the sensing performance of the entire system.
In a cognitive radio network, noise in the actual environment will affect the sensing performance of the system. In order to improve the sensing performance of the system, it is necessary to remove redundant information such as noise in the spectrum signal. However, most signals received by the system are nonstationary and nonlinear. In order to deal with nonstationary and nonlinear signals, the traditional methods use short-time Fourier transform, Wigner-Ville distribution, and wavelet transform. However, these methods are essentially Fourier transforms, which are subject to indefinite principles. Therefore, there are certain defects in the traditional signal processing methods. Literature  proposed an Empirical Mode Decomposition (EMD) algorithm that decomposes complex signals into a collection of Intrinsic Mode Functions (IMF). Literature  proposed a spectrum sensing method based on Hilbert-Huang Transform (HHT). First, the spectrum signal is subjected to EMD decomposition to obtain the IMF component. Calculate the Hilbert spectrum and total HHT spectrum for each IMF component and then find the total marginal spectrum. Compare the total marginal spectrum with the threshold to determine whether the primary user exists.
In order to make the system have good adaptability in the actual environment, reduce the impact of noise and accurately determine whether the main user exists. Based on previous studies, we propose a spectrum sensing method (KEMDSS) based on EMD algorithm and K-means clustering algorithm. The KEMDSS method will uses EMD algorithm and wavelet threshold algorithm to process the received spectrum signal, removing the redundant information in the signal, and then use the K-means clustering algorithm to determine whether the primary user exists. To verify the effectiveness of the KEMDSS method, this article uses DMM and DMEAE as signal features. The experimental results also show that KEMDSS method can improve the sensing performance of the system no matter which feature is used.
2. System Model
The KEMDSS method is mainly composed of a signal processing part and a clustering algorithm part. In order to be able to deal with complex spectrum signals in the actual environment, we posed a signal processing method (WEMD) which combines the EMD algorithm and the wavelet threshold algorithm. According to the multiscale filtering characteristics of EMD, the spectrum signal can be decomposed into several high-frequency and low-frequency IMF components. Considering signal noise and other interference factors are mainly concentrated in the high-frequency band, the WEMD method uses wavelet algorithm to remove the noise in high-frequency IMF components. The purpose is to reduce the impact of noise and other factors on the system. The clustering algorithm mainly uses the DMEAE and DMM signal features, collects a certain amount of training sets and test sets, and then uses the K-means clustering algorithm to train the classifiers for judgment, thereby determining whether the primary user exists. The entire KEMDSS flowchart is as shown in Figure 1.
2.1. Cognitive Radio Network Model
KEMDSS is a cooperative spectrum sensing method. The cooperative spectrum sensing system can reduce the influence of uncertain factors such as multipath fading in the actual environment and thus can improve the sensing performance of the entire system. Assume that there are secondary users (SU) in the cognitive radio network and each SU has sampling points. We assume that indicates that the primary user (PU) signal does not exist and indicates that the primary user PU exists. So the model of the received signal under these two assumptions can be expressed as
where represents the PU signal and represents a Gaussian white noise signal with a mean of 0 and a variance of . Assume that and , respectively, indicate the state of the PU. The availability A of the channel can be defined as
where A = 0 indicates that the channel is unavailable and A = 1 indicates that the channel is available. The false alarm probability () and detection probability () of the system can be defined as
2.2. WEMD Signal Processing Algorithm
2.2.1. EMD Algorithm
This paper introduces the EMD algorithm. Compared to traditional signal processing methods, EMD algorithm can handle nonstationary and nonlinear signals well. The EMD algorithm decomposes the original signal into several IMF components based on the local time feature scale of the signal . The specific steps of the EMD decomposition process are as follows.
Step 1. Find all local maxima and local minima in the spectrum signal .
Step 2. Find the maximum envelope and the minimum envelope , and calculate the average :
Step 3. Subtract from the original signal to get the first component :Determine whether satisfies the two conditions of the IMF. If Step 4 is satisfied, perform Steps 1 and 2 for to obtain the mean envelope , and calculate :In order until meets the IMF condition, then is the first IMF component.
Step 4. Subtract from the spectrum signal to obtain residual ; that is,Treat as the original signal and repeat Steps 1~4 to get . And so on until the residual is a monotonous function or a constant. Therefore, the signal can be expressed as
2.2.2. Wavelet Threshold Algorithm
Considering wavelet coefficients of signal and noise are different, so the wavelet threshold algorithm can remove noise components in the signal based on this characteristic. The specific steps of the wavelet threshold denoise algorithm are as follows [16, 17].
Step 1. Wavelet transform the signal to obtain the wavelet coefficient .
Step 2. Threshold the wavelet coefficients to obtain the estimated coefficient .
Step 3. Use to perform wavelet reconstruction and obtain the denoised signal.
This paper will use the soft threshold function, which is expressed as follows:
where is the VisuShrink threshold proposed by Donoho and Johnstone , which satisfy
where is the standard deviation of noise and is the number of sample points of the signal.
2.2.3. WEMD Signal Processing Algorithm
After EMD decomposition, the spectrum signal can be represented by adding several IMF components and residual components. In an actual environment, noise usually occurs at high frequencies and useful signals are generally concentrated at low frequencies. Therefore, some high-band IMF can be removed, and then the remaining IMF can be reconstructed to obtain denoised signals, but this method will lose some useful signals . In order to solve this problem, WEMD signal processing algorithm is proposed in this paper. After the signal is decomposed by EMD, wavelet threshold method is used to denoise high-band IMF to obtain DIMF, then DIMF and the original low-band IMF are superimposed to obtain the reconstructed signal .
The literature  proves that this algorithm can remove the noise component of the original signal well, and it will not lose the useful information in the original signal. Therefore, the use of WEMD signal processing method can reduce the interference of noise on the spectrum sensing system and improve the system sensing performance. Figure 2 shows the WEMD signal processing flow chart.
WEMD signal processing algorithm steps are as follows.
Step 1. The spectrum signal is decomposed by EMD to obtain IMF components.
Step 2. Find the high-frequency components in the IMF and obtain the critical point .
The literature [20, 21] uses the continuous mean square error criterion to find the critical value :Among them, N represents the sampling point number of the signal . The critical point represents the global minimum point of all energy and is also the demarcation point between the low-frequency band and the high-frequency band.
Step 3. Use the wavelet soft threshold algorithm to denoise the high-frequency components to obtain the denoised component .
Step 4. The DIMF is superimposed with the original low-frequency IMF to obtain the reconstructed signal .
2.3. Clustering Algorithm Model
After the spectrum signal is processed by the WEMD algorithm, we will use the K-means clustering algorithm for further processing. The traditional spectrum sensing method generally deduces the corresponding threshold to judge whether the main user exists, but the problem of inaccurate threshold derivation occurs. Spectrum sensing can also be considered as a problem of two classifications (PU presence or PU absence). KEMDSS method will use the K-means clustering algorithm to determine whether the primary user exists. Tradition spectrum sensing methods for machine learning often use energy features to characterize signals. Then the clustering algorithm is used to divide the feature values into channel usable classes and channel unavailable classes. In this paper, we use the eigenvalues of the covariance matrix in the random matrix to describe the signal characteristics.
Assuming that represents the reconstructed sample vector of the i-th SU user, the sensing matrix in a period of time can be expressed as
Therefore, the covariance matrix of the signal can be expressed as . Traditional random matrix spectrum sensing method usually uses DMM eigenvalues and DMEAE eigenvalues as the statistics and then deduced the corresponding thresholds to judge. However, the traditional method has relatively high requirements for the derivation of the threshold, and it cannot adaptively make corresponding changes according to different sensing environments, which will adversely affect the sensing performance of the entire system. Therefore, this paper introduces the method of K-means clustering algorithm on the traditional method of random matrix. The status of the primary user is determined by training the classifier.
In the sensing system, the signal eigenvalues of SU nodes are first collected and use these eigenvalues to construct a training set:
where is a column vector, represents the random matrix eigenvalues, and represents the number of eigenvalues. After this training set is obtained, the K-means clustering algorithm is used to classify these features and train a classifier. Considering the actual situation of spectrum sensing, it is generally divided into two categories. One class represents that the primary user does not exist, and the other represents the primary user exists. After the training process is completed, the corresponding cluster center is obtained. We use formula (17) to make the decision . Where is used to control and of the system, represents the number of clusters. In this paper, the value of is 2.
If the data of the test set satisfies (17), it is determined as A = 1; otherwise it is determined as A = 0.
The specific steps of the K-means clustering algorithm are as follows.
Step 1. Firstly we collect sufficient number of training feature matrix T.
Step 2. Randomly initialize the cluster center .
Step 3. Calculate the distance from each feature to each cluster center , find the smallest and reclassify.
Step 4. Calculate the mean of all points in each cluster and take this as a new cluster center.
Step 5. Calculate . If converges, the algorithm stops; otherwise returns to Step 3.
Step 6. After successful training the classifier, we import the corresponding test set and make a corresponding decision on channel availability A by formula (17).
3. Experimental Part
In order to simulate the actual environment, this experiment uses the FM signal, chirp signal, and mixed signal as experimental simulation signals. The carrier frequency is 702 KHz and the sampling rate is 4 MHz. To ensure the accuracy of the experiment, we obtained 1000 signal features (DMEAE, DMM), 500 of them were used as training sets and 500 were used as test sets.
3.1. Clustering Algorithm Classification Effect
In order to verify the effectiveness of the method, we first simulated the mixed signal. In the experimental part, we adopted two characteristic values of DMEAE and DMM. The clustering results are shown in Figures 4–7.(1)DMEAE clustering effect graph under SNR=-16(2)SNR=-16 DMM clustering effect chart
Figures 4 and 5 show the clustering of DMEAE features before and after the KEMDSS method under SNR=-16, respectively. Figures 6 and 7, respectively, show that in the conditions of SNR=-16, the DMM features are processed before and after the KEMDSS method. It can be seen from the above that after the KEMDSS method is processed, the clustering effect of the two features is better than the original. This is because after the KEMDSS method is processed, the difference between the sampled signal covariance matrix and the noise covariance matrix is more obvious, so the eigenvalue of the covariance matrix has a better clustering effect under the K-means clustering algorithm.
3.2. ROC Curve under Constant False Alarm Condition
In order to reflect the superiority of the KEMDSS method, we used the KEMDSS method and the traditional EMD (TEMD) method to process the spectrum signal. The TEMD method decomposes the spectrum signal and discards some high-band signals. However, the KEMDSS method uses the wavelet threshold method to further process the high-frequency signals. The test results of the two methods are as follows.
(1) ROC Curve of DMEAE Characteristics. Figures 8 and 9 show the ROC curves of the DMEAE features after KEMDSS and TEMD processing at different SNR, respectively. From this we can see that the detection performance of the spectrum signal after KEMDSS method or TEMD method is better than the original one, and the best effect can be achieved by using KEMDSS method. Tables 2 and 3 show the detection probability of each algorithm when the false alarm probability is constant at different SNR. From Table 2, we can see that, under the condition of SNR=-16 and , the using KEMDSS and TEMD as the algorithm has improved the detection performance by 69% and 41%, respectively, compared with the traditional algorithm. This is because both KEMDSS and TEMD can remove the noise in the spectrum signal. Compared with the TEMD method, KEMDSS can remove the redundant components such as noise in the spectrum signal without losing some useful signals and thus can obtain better detection results. Table 3 shows the comparison of different algorithms in SNR = -18.
(2) ROC Curves of DMM Features. Figures 10 and 11 show the ROC curves of DMM features after KEMDSS and TEMD processing at different SNR, respectively. From the figure, it can still be seen that the detection performance of the system has been improved after KEMDSS and TEMD removal of noise components. We can find that KEMDSS method is more effective than the TEMD method. Tables 4 and 5 show the detection probability of each algorithm when the false alarm probability is constant under different SNR. From Table 4, we can calculate that, under the condition of SNR = -16 and , using KEMDSS and TEMD as an algorithm is 88% and 62% higher than the conventional algorithm, respectively.
From the experimental results in Section 3.2 we can know that regardless of the characteristics used, the detection performance of the system is improved after the spectrum signal is processed by KEMDSS and TEMD. And the KEMDSS algorithm will have a best result. The sampled signal is processed by the KEMDSS algorithm, which minimizes the impact of noise uncertainties on the system and thus achieves good sensing performance.
3.3. ROC Curve for Different Cooperative Users
Figures 12 and 13 show the ROC curve of the KEMDSS algorithm under different numbers of users. From this we can see that the sensing performance of the system is better as the number of collaborative users increases. This is because cooperative perception can reduce the interference of various factors such as multipath fading and shadows in the propagation environment. Therefore, the greater the number of collaborating users, the stronger the anti-interference ability of the system and the better the sensing performance.
3.4. ROC Curve under FM Signal and Chirp Signal
In order to verify the practicality of the KEMDSS method, different spectrum signals (FM signals and chirp signals) were used for experiments. The specific experimental results are as follows:(1)The simulation results under the FM signal(2)Simulation results under the chirp signal
Figures 14–17 show the ROC curves of the DMM and DMEAE features under the FM and chirp signals, respectively. It can be seen from the figure that, compared with the traditional spectrum sensing method, the KEMDSS method can achieve better detection performance and it can be applied to various spectrum signals in the actual environment. This is because EMD plus wavelet threshold algorithm can remove the noise information in the spectrum signal. It can reduce the impact of noise and other pieces of redundant information on the sensing system, thereby improving the system sensing performance.
This paper proposes a spectrum sensing method based on EMD algorithm and K-means algorithm. Compared with traditional spectrum sensing methods, KEMDSS can well solve the problem that traditional spectrum sensing methods are affected by noise and the threshold is difficult to determine. The wireless spectrum signal is first processed by the WEMD algorithm, then uses the K-means clustering algorithm to train the classifier, and finally uses the classifier to judge whether the primary user exists. In the experimental part, we analyzed the detection results of two random matrix features DMEAE and DMM under the KEMDSS algorithm. From the experimental results, it can be seen that KEMDSS spectrum sensing method has better sensing performance than the conventional spectrum sensing method.
The data used to support the findings of this study are currently under embargo while the research findings are commercialized. Requests for data, 12 months after publication of this article, will be considered by the corresponding author.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This work was supported in part by Special Funds from the Central Finance to Support the Development of Local Universities under no. 400170044, the State Key Laboratory of Management and Control for Complex Systems, Institute of Automation, Chinese Academy of Sciences, under Grant no. 20180106, the Science and Technology Program of Guangdong Province under Grant no. 2016B090918031, the Degree and Graduate Education Reform Project of Guangdong Province under Grant no. 2016JGXM_MS_26, the Foundation of Key Laboratory of Machine Intelligence and Advanced Computing of the Ministry of Education under Grant no. MSC-201706A, and the Higher Education Quality Projects of Guangdong Province and Guangdong University of Technology.
T. Chen, J. Liu, L. Xiao, and L. Huang, “Anti-jamming transmissions with learning in heterogenous cognitive radio networks,” in Proceedings of the 2015 IEEE Wireless Communications and Networking Conference Workshops, WCNCW 2015, pp. 293–298, March 2015.View at: Google Scholar
X. Y. Wang and Y. G. Lu, “DMM based spectrum sensing method for cognitive radio systems,” Journal of Electronics & Information Technology, vol. 32, no. 11, pp. 2571–2575, 2010.View at: Google Scholar
N. Liu, S. H. Shi, and B. Yang, “Spectrum sensing method based on ME-S-ED,” Measurement & Control Technology, 2016 (Chinese).View at: Google Scholar
Y. Zhang, P. Wan, S. Zhang, Y. Wang, and N. Li, “A spectrum sensing method based on signal feature and clustering algorithm in cognitive wireless multimedia sensor networks,” Advances in Multimedia, vol. 2017, no. 4, pp. 1–10, 2017.View at: Google Scholar
Y. Wang, Research on Spectrum Sensing and Spectrum Allocation in Cognitive Radio Networks [Doctoral, thesis] , Post Doctoral Research Report of South China University of Technology, 2016 (Chinese).
W. Zhang, F. Yu, and H.-M. Guo, “Improved adaptive wavelet threshold for image denoising,” in Proceedings of the 2009 Chinese Control and Decision Conference, CCDC 2009, pp. 5958–5963, June 2009.View at: Google Scholar
Q. Han, Z. Zhang, D. Wang, L. Liu, and X. Hou, “An OTDR event analysis algorithm based on EMD-based denoising and wavelet transform,” in Proceedings of the 12th IEEE International Conference on Electronic Measurement and Instruments, ICEMI 2015, pp. 1034–1038, July 2015.View at: Google Scholar
Y. Yan and Z. Cui, “Noise and zero excursion elimination of electrostatic detection signals based on EMD and wavelet transform,” in Proceedings of the 2009 2nd International Congress on Image and Signal Processing, CISP'09, October 2009.View at: Google Scholar
W. Sun F, Y. Peng H, and U. X. Jian-Hua, “A de-noising method for laser ultrasonic signal based on EMD,” Journal of Shandong University, 2008 (Chinese).View at: Google Scholar