Industrial Informatics: Applications of Mobile and Wireless Emerging Technologies in Industry 4.0
View this Special IssueResearch Article  Open Access
MedianDifference Correntropy for DOA under the Impulsive Noise Environment
Abstract
The source localization using direction of arrival (DOA) of target is an important research in the field of Internet of Things (IoTs). However, correntropy suffers the performance degradation for direction of arrival when the two signals contain the similar impulsive noise, which cannot be detected by the difference between two signals. This paper proposes a new correntropy, called the mediandifference correntropy, which combines the generalized correntropy and the median difference. The median difference is defined as the deviation between the sampling value and the median of the signal, and it intuitively reflects the abnormality of impulsive noise. Then, the median difference is combined with the generalized correntropy to form a new weighting factor that can effectively suppress the amplitude level of impulsive noise. To improve the robustness of the algorithm, an adaptive kernel size is also integrated into the weighting factor to obtain the optimal local feature. The influence of adaptive kernel sizes on the proposed algorithm is simulated, and the comparison between three typical directionofarrival estimation algorithms is conducted. The results show that the accuracy of the mediandifference correntropy is significantly superior to the correntropybased correlation and the phased fractional lowerorder moment for a wide range of alphastable distribution noise environments.
1. Introduction
The source localization using direction of arrival (DOA) of target is an important research in the field of Internet of Things (IoTs). The directionofarrival (DOA) approaches based on acoustics have many applications including radar, sonar, seismic exploration, navigation, and sound source tracking [1–3]. DOA estimation is usually regarded as a problem of signal matching, and the performance is significantly influenced by noise. A majority of existing DOA estimations are based on the concept that noise follows Gaussian distribution [4, 5]. Since the Gaussian process has secondorder and higherorder statistics, the traditional DOA algorithms can easily evaluate the signal characteristics according to secondorder statistics [6]. The multiple signal classification (MUSIC) algorithm [7, 8] and estimation method of signal parameters via rotational invariance techniques (ESPRIT) are the basic subspace algorithms which have good performance [9, 10]. MUSIC is the representative of the noise subspace algorithm, and ESPRIT is the representative of the signal subspace algorithm.
Because of atmospheric noise, electromagnetic interference, sea clutter, car ignitions, and office equipment, the signal is corrupted by the extremely impulsive noise that exhibits irregularity in time domain. In addition, the probability density functions of impulsive noise decay with heavy tails and do not follow a common Gaussian distribution [11]. Therefore, alphastable distribution is usually used to define impulsive noise [12]. The conventional covariance matrix is calculated from the secondorder statistics of the signal, which may be infinite when the data are corrupted by the extremely impulsive noise [13, 14]. In addition, the conventional DOA algorithms cannot be decomposed into the signal subspace and the noise subspace with covariance matrix. Thus, it has become increasingly important to study DOA under the impulsive noise environment.
The relevant statistical algorithms mainly focus on the existence of statistics. The fractional lowerorder statistics (FLOS) algorithms [15] exhibit a more desirable performance than the secondorder statistics for alphastable distribution [16]. When the signal contains impulsive noise, the signal bears the fractional lowerorder statistics. The FLOS algorithms employ the minimum dispersion (MD) criterion to suppress impulsive noise, such as the robust covariation in ROCMUSIC [17], the fractional lowerorder statistics (FLOM) based MUSIC, and the phased fractional lowerorder moment (PFLOM). The FLOM algorithmbased MUSIC obtains a finite covariance by suppressing one of two crosscorrelation signals when the characteristic exponent of alphastable distribution ranges from 1 to 2 [18]. The PFLOM algorithm gets the accurate DOA estimation with circular symmetrical signals, embedded in the additive impulsive noise [19]. However, investigators demonstrate that the performance of FLOM and PFLOM algorithms depends on the relationship between the parameter of fractional lowerorder moment and the characteristic exponent of alphastable distribution. If the characteristic exponent is unknown, the performance of FLOM and PFLOM algorithms seriously decreases.
The correntropy criterion [20] is a relatively simple method that can measure the local similarity between two signals [21, 22]. Because it has the properties of Mestimation, the correntropy has been widely used in the impulsive noise environment [23]. Zhang et al. [24] investigated a narrowband model based on the generalized correntropy which is called the correntropybased correlation (CRCO) in impulsive noise environment. The CRCO algorithm imposes a correntropy operator on the covariance matrix to depress impulsive noise. The generalized correntropy is suitably for dealing with the template matching between the received signal and the template signal [25]. Since DOA estimation is a matching problem between two signals, the generalized correntropy is incapable of measuring the difference between the outliers and fails to suppress impulsive noise. Thus, when data contain the similar impulsive noise, the crosscorrelation of the generalized correntropy is infinite.
In order to solve the problem that correntropy cannot distinguish the similar impulsive noise, we propose a mediandifference correntropy (MDCO) algorithm. The MDCO depending on the inner product of vectors measures the similarity of multidimensional properties from input space. A weighting factor of a median difference is defined and evaluates the similarity between the sample value and the median of the signal. The median difference intuitively reflects the abnormality of impulsive noise to guarantee that the autocorrelation is finite. Then, MDCO derives the weighting factor of the generalized correntropy from the correntropy criterion which suppresses the larger impulsive noise. Hence, the weighted covariance function of signal employs the generalized correntropy instead of secondorder statistics. These two criteria map the signal from the lowdimensional space to an infinitedimensional reproducing kernel Hilbert space (RKHS) with impulsive noise, thus including higherorder signal statistics. The adaptive kernel size is applied to the weighting factor for MDCO which can obtain the optimal local feature. At convergence, MDCO is unbiased and it is applicable to achieve CRLB under various parameters. The main work is summarized as follows:(1)We propose a mediandifference correntropy (MDCO) algorithm, which can effectively combine the generalized correntropy and the median difference to suppress impulsive noise(2)To improve the robustness of MDCO, we also introduce a novel adaptive kernel size into the weighting factor of the generalized correntropy and the median difference
The rest of this paper is organized as follows: The problem model is defined with DOA in Section 2. In Section 3, we present the mediandifference correntropy (MDCO) for DOA estimation under impulsive noise environment. The performance evaluation is presented in Section 4. Finally, conclusions are drawn in Section 5.
2. Problem Formulation
It is considered that a uniform linear array of M isotropic acoustic sensors receives the farfield signal generated by narrowband sources. Then, Figure 1 illustrates the linear array model.
The first array sensor is a reference sensor, and the received signal of the mth sensor can be expressed aswhere is the pth signal at time t, is the noise of the mth sensor, denotes the time delay of the mth acoustics sensor relative to the reference sensor for the pth signal, and f is the center frequency of signal.
The received signal of acoustics can be expressed aswherein which the superscript represents transpose, is the vector of the signal received by the acoustic array sensors, is the vector of the acoustic source, is the vector of impulsive noise that follows alphastable distribution, and is the array manifolds of the array sensors. There is an assumption that the number M of the array elements is greater than the number P of the acoustic sources. is the steering vector that can be expressed as
Impulsive noise in (2) is usually defined as symmetric alphastable distribution. However, the probability density function of symmetric alphastable distribution cannot be expressed by a general expression [24]. Therefore, it is generally introduced by its characteristic function which can be expressed aswhere α is the characteristic exponent of symmetric alphastable distribution whose range is . Moreover, the smaller the α becomes, the more impulsive the nonGaussian noise will be. γ is the dispersion.
Usually, the received data contain very little impulsive noise and can also be represented aswhere denotes the noise without outliers, denotes outliers with sparse characteristics, , and . Meanwhile, denotes the received data without outliers. We can assume that the noise follows Gaussian distribution. Then, the noise follows symmetric alphastable distribution. The purpose of our algorithm is to eliminate outliers .
Many DOA estimation algorithms usually use the covariance to represent the secondorder statistics of the signal. The conventional covariance matrix can be given aswhere the superscript represents transpose, the notation represents the expectation operation, and and are the covariance matrix of signal and noise, respectively. If the noise contains impulsive noise, the signal and noise cannot be completely orthogonal. The conventional DOA algorithms that describe signal with impulsive noise cannot be decomposed into the signal subspace and the noise subspace by the covariance matrix. Therefore, the covariance matrix (7) can also be represented aswhere denotes the covariance term generated by outliers and the covariance is dominated by impulsive noise when some elements of covariance is much larger than the corresponding elements of . Because natural noise and manmade noise often include outliers, some elements of matrix may have a large value which gives rise to false DOA estimation. From Figure 1, we can see that the received data of the acoustic sensors suffer impulsive noise with large values . For example, because data and which are located on the different sensors are both impulsive noise at the same time, the crosscorrelation with a relatively large value is infinite:where Y and Z are two sampling points of total snapshots and represents the crosscorrelation of data and at time . Assuming that only one acoustic source impinges on the array, the signal is a delayed signal of with time when the signal is received by the array sensor M at time . Furthermore, array sensors m and M both contain impulsive noise at time . At this point, the noise covariance is dominant so that the noise subspace would spread to the signal subspace, causing the characteristics of the signal to be covered. In addition, the autocorrelation of data , , and is also nonexistent.
The goal is to search for an efficient strategy of suppressing impulsive noise that makes it possible towhere indicates the operator of suppressing impulsive noise. In principle, as long as provides a small contribution, we can accurately estimate DOA from the covariance.
3. MedianDifference Correntropy (MDCO)
In order to solve impulsive noise, the structure of this section is as follows: firstly, the mediandifference correntropy (MDCO) is proposed. Next, we summarize some properties for MDCO. Finally, the application of MDCO for DOA estimation is designed.
Normally, if the signal follows Gaussian distribution, there is a probability that results in , where is the collection of Gaussian distribution signal, is the mean value, and is the standard deviation ( and have no relationship to the similar parameters behind). If , the is a zeromean Gaussian process with variance . Therefore, there is a probability of . At this point, we can simply set a threshold ε to measure the similarity between the signals and set .
Definition 1. For two variables Y and Z of total snapshots , where , if , we say that the variable Y is similar to Z.
As shown in Table 1, the similarity of the variables Y and Z is divided into four cases. If , , the variables Y and Z are considered normal values. But, if , where , we consider the variables Y and Z to be the outliers containing impulsive noise. One of our goals is to suppress the three abnormal cases that contain impulsive noise.

3.1. MedianDifference Correntropy (MDCO)
In this section, we present a mediandifference correntropy (MDCO) algorithm for DOA with alphastable distribution. The correntropy is a new method for nonlinear and local optimal measurement of two random variables Y and Z. It is defined as follows:where is a translation function of the shiftinvariant Mercer kernel [26] and is the joint probability density function of Y and Z. In this paper, the Gaussian density function is used as kernel function. The Gaussian kernel creates a RKHS [26] with universal approximating capability. The Gaussian kernel is numerically stable and usually gets reasonable results. In this paper, we use a Gaussian kernel to suppress impulsive noise.
Assuming that the random variables Y and Z follow symmetric alphastable distribution whose characteristic exponent is , the mediandifference correntropy (MDCO) can be considered aswhere σ is the kernel size of Gaussian kernel. The variable B is the average of of the preprocessing data which can approximately represent the median of the received data.
Because impulsive noise contained in the signal is sparse, the variable B can effectively weigh the amplitude of the signal. The variables and reflect the median difference to which the received signal deviates from the median B. Therefore, the variables and are called the weighting factors of the median difference. The correntropy measures the similarity between the signals and is called the weighting factor of the generalized correntropy. The median difference evaluates the deviation between the sampling value and the median of the signal, and it intuitively reflects the abnormality of impulsive noise. Then, the median difference is combined with the generalized correntropy to form a new weighting factor, which can effectively suppress the amplitude level of impulsive noise. Then, the median difference and the generalized correntropy are combined with the traditional covariance, and a novel generalized weighted covariance is obtained.
To simplify (12), the median differences and can be combined. The variable reflects the deviation degree of . The relaxation of (12) can be expressed aswhere and are the absolute value of the signal and and are the kernel sizes that control the scale of the metric.
The shorttime energy method is a common method to detect impulsive noise in time domain [27]. We use a shorttime average energy method to preprocess the received data. However, the data processed by the shorttime average energy method does not directly be used as the input of MDCO. The shorttime average energy indicates that the energy of the signal is average in a short segment. The shorttime average energy in the segment can be expressed as :where is the sampling sequence of the original signal. By the formula above, we can obtain the average energy with I data points.
If the signal energy is much larger than the shorttime average energy, the preprocessing signal can be represented as
Assuming that is the sum of in total snapshots. Therefore, in this study, B can be expressed aswhere N is the total sampling size of snapshots. The signal with impulsive noise does not participate in calculating the median of the signal. Therefore, for a majority of data, the median B is less than due to preprocessing.
The key to the mediandifference correntropy is that the kernel size can work in the confidence interval to eliminate impulsive noise. The kernel size is usually related to the dispersion coefficient γ by considering the local feature of the signal. We take a modified Sigmoid function to make the kernel size adaptive [24, 28–30]. Then, kernel sizes and can be expressed aswhere and indicate the scale of the modified Sigmoid function and and control the monotonicity of the modified Sigmoid function [28]. The shrinkage direction and rate are determined by and , respectively.
The adaptive mediandifference correntropy provides a new metric criterion for impulsive noise. And the Mestimators of the autocorrelation and crosscorrelation of random variables Y and Z are in existence associated with the generalized correntropy and median difference.
In order to prove the effectiveness of MDCO in theory, five crucial properties of MDCO are listed as follows:
Property 1. (C1). When the noise follows temporally stationary zeromean white Gaussian processes, MDCO is reduced to traditional secondorder statistics, and the performance of MDCO is equal to traditional signal subspace algorithms (see Appendix A for this proof).
Property 2 (C2). If the variable Y is an outlier rather than Z, the MDCO can eliminate impulsive noise. of MDCO approaches zero; that is to say, is close to zero at this point (see Appendix B for this proof).
Property 3. (C3, C4). When the variables Y and Z both contain impulsive noise, and both approach zero (see Appendix C for this proof).
Property 4. When the variables Y and Z both contain impulsive noise, the MDCO is more effective than the CRCO for noise suppression (see Appendix D for this proof).
Property 5. Assuming that the random variables Y and Z follow the symmetric alphastable distribution, is bounded and the MDCO has the generalized correntropy statistics (see Appendix E for this proof).
Because of the abovementioned properties, the generalized weighted covariance has finite autocorrelation and crosscorrelation at all moments. When the sampling data contain impulsive noise, of MDCO approaches zero; that is to say, is close to zero at this point. If two signals contain the different degrees of impulsive noise at the same time, MDCO can induce a very small weighted covariance by means of the median difference operator and then obtain a convergent covariance matrix. As shown in Table 2, the performances of the MDCO is listed for different noise cases. In summary, the MDCO has reliably good performance in all cases.

4. Performance Simulations
In this experiment, we compare the MDCOMUSIC to FLOMMUSIC, PFLOMMUSIC, and CRCOMUSIC and test the performance of the proposed MDCOMUSIC under different parameter conditions. In our experiment, we evaluate the performance of the algorithm from five aspects, namely, various kernel sizes of MDCOMUSIC, GSNR, snapshots of sampling, different characteristic exponents of alphastable distribution, and angular separation of two direction angles.
The resolution probability can well define the performance of the four algorithms. We use a popular resolution criterion to measure the spatial spectrum which can be given aswhere and are the two independent direction angles, is the midangle between them, and is the spatial spectrum. The two direction angles are resolvable if the result on the left of the equation is smaller than that on the right; otherwise, the two direction angles are not resolvable.
RMSE is the deviation criterion between the observation and the true value and can be given bywhere L represents the total number of MC run and and are estimations of and in the lth MC run, respectively.
Kozick and Sadler has come to a closedform expression of CRLB for the impulsive noise [31], which can be expressed aswhere is a diagonal matrix of the received signal, the subscript of is omitted from , is the differential of the array manifolds , is the real part, and the coefficient can be expressed as , in which is the modulus of the impulsive noise and is probability density function of x. Note that with and with . For simplicity, the coefficient for can be approximated by firstorder linear interpolation with and [32].
4.1. Experimental Setup
The linear array is set to omnidirectional acoustic sensors which are placed in a half of the wavelength at the center frequency of the signal. We consider the case that the number of narrowband sources is two which follow temporally stationary zeromean white Gaussian processes. The central frequency of signal is set to . In order to distinguish the signal, the power of each signal is different. The symmetric alphastable distribution is used to model impulsive noise with the characteristic exponent and dispersion . The GSNR is set to . The sampling frequency of each array element is . All experiments are carried out through two directions of arrival associated with and . The snapshots are . Assume that the parameter p of the fractional lowerorder statistics is equal to 1.1 in FLOMMUSIC and PFLOMMUSIC algorithms. In every experiment, 500 Monte Carlo (MC) runs are performed. Furthermore, we compare the resolution probability and the rootmeansquare error (RMSE) between the different algorithms.
4.2. Kernel Parameters of MDCO
In this experiment, we test the performance of MDCOMUSIC with respect to the different kernel parameters of the weighting factors for the separate characteristic exponents and . In order to eliminate impulsive noise, it is necessary to satisfy the requirement that the large local weighting factor of the median difference and generalized correntropy associate with the small kernel sizes and ; that is, and must be less than zero. Because the weighting factor of the median difference and generalized correntropy are independent parameters, the performance can be analyzed separately. Figures 2 and 3 show that the appropriate and for the median difference can effectively suppress impulsive noise and evidently decrease RMSE. From Figure 2, we can see that would get the best result with three different values of . The convergence of is not affected by . Empirically, it is revealed that impulsive noise can be effectively alleviated in Figure 3, when .
Figures 4 and 5 illustrate the robustness of MDCOMUSIC in terms of the appropriate and for the generalized correntropy. MDCOMUSIC gains the optimal performance when the kernel parameter varies from 9 to 12. Figure 5 shows that RMSE gains a more significant decrease if . Usually, it is a common method that the optimal kernel parameters are obtained by experiments. Therefore, we set , , , and in other experiments.
4.3. Performance of Different Algorithms
4.3.1. Performance Comparison versus GSNR
Figure 6 analyzes the performance of MDCOMUSIC, FLOMMUSIC, PFLOMMUSIC, and CRCOMUSIC under the average GSNR from −10 dB to 8 dB. We can derive the conclusion that MDCOMUSIC outperforms FLOMMUSIC, PFLOMMUSIC, and CRCOMUSIC algorithms both in resolution probability and in RMSE. The MDCOMUSIC algorithm yields particularly robust estimation to impulsive noise and can accurately estimate the DOA when the GSNR is 2 dB. The RMSE of MDCOMUSIC is much smaller than that of FLOMMUSIC and PFLOMMUSIC. It is clear from Figure 6 that that the performance of MDCO coincides with the CRLB at high SNR.
(a)
(b)
4.3.2. Performance Comparison versus Characteristic Exponent of AlphaStable Distribution
In this simulation, we aim at evaluating the performance of the characteristic exponent on the basis of probability of resolution and RMSE. The characteristic exponent of the symmetric alphastable distribution ranges from 1 to 2. Figure 7 shows that the performance of MDCOMUSIC outperforms the FLOMMUSIC, PFLOMMUSIC, and CRCOMUSIC both in resolution probability and in RMSE. Moreover, the smaller the α becomes, the worse the performance of FLOMMUSIC and PFLOMMUSIC will be. Another observation is that it is beneficial to employ MDCOMUSIC rather than PFLOMMUSIC if the source signal is not circularly symmetrical. Furthermore, when α approaches 1, FLOMMUSIC and PFLOMMUSIC gain a large RMSE if the parameters of the fractional lowerorder statistics cannot automatically adapt. The result shows that the accuracy of the PFLOMMUSIC algorithm is influenced by the relationship of the FLOM parameter and the characteristic exponent α. Because of the contribution of two weighting factors, MDCOMUSIC is not affected by the variety of characteristic exponent and can effectively estimate the DOA in any case.
(a)
(b)
4.3.3. Performance Comparison versus Snapshots of Sampling
In this simulation, the performance of MDCOMUSIC, FLOMMUSIC, PFLOMMUSIC, and CRCOMUSIC about snapshots is tested under symmetric alphastable distribution. Considering that the snapshots start at 50 and end at about 1000. Figure 8 illustrates that MDCOMUSIC gains a more significant enhancement in resolution probability with the increasing number of snapshots and a more evident decrease in RMSE compared to the other three algorithms. Consequently, MDCOMUSIC can work effectively for a small number of snapshots.
(a)
(b)
4.3.4. Performance Comparison versus Angular Separation
Figure 9 depicts the performance of four algorithms when the direction angle of the second signal is gradually away from that of the first angle. The first direction of arrival is . As expected, the MDCOMUSIC can effectively estimate the two directions of arrival with angle separation. Furthermore, MDCOMUSIC requires a smaller angle separation threshold than FLOMMUSIC, PFLOMMUSIC, and CRCOMUSIC associated with a fixed probability of resolution, and then MDCOMUSIC gains a lower RMSE than other algorithms associated with a fixed angular separation threshold.
(a)
(b)
Through the RMSE curve in the above experiment, MDCO can effectively converge to a stable trend under different parameters, which ensure the boundedness and improve the robustness of the algorithm.
5. Conclusion
The issue of DOA in the presence of the abnormal similarity is solved through the adaptive mediandifference correntropy. The MDCO combines the weighting factor of the median difference and generalized correntropy to suppress impulsive noise. The median difference evaluates the deviation between the sampling value and the median of the signal, and it intuitively reflects the abnormality of impulsive noise. The generalized correntropy measures the similarity between the two signals. These two weighting factors map the signal to an infinitedimensional space with impulsive noise, thus including much statistical information. By controlling the adaptive kernel size, MDCO can effectively deal with DOA. We also prove that the MDCO satisfies some robust properties and applies the MDCO to DOA estimation combined with MUSIC. Experimental results illustrate that the MDCOMUSIC is more robust than the FLOMMUSIC, PFLOMMUSIC, and CRCOMUSIC at a much lower GSNR and in the extremely impulsive noise environment.
Appendix
A. Proof of Property 1
In (12), the weighting factors can be represented as
When the noise does not contain impulsive noise, it can be obtained:
The constraint of the weight factor C can be translated into
At this point, a small enough positive number ξ can always be found so that the following formula is established:
Thus, it is easy to arrive at . In other words, the authors can getwhich means that the MDCO is reduced to the traditional secondorder statistics algorithm when the noise follows Gaussian distribution.
B. Proof of Property 2
If the variable Y is an outlier rather than Z, that is, and , then , , and . Therefore, the weighting factor of generalized correntropy and median difference is very small which can effectively suppress outlier Y. Let be expressed aswhere results in a small value that is close to zero.
C. Proof of Property 3
If the variables Y and Z are signals with impulsive noise, the authors can obtain and . If , the variable Y is similar to Z and ; otherwise, . Let us discuss the two cases together, and can be expressed as
Although the variables Y and Z involve impulsive noise, the median difference of and works well on it.
D. Proof of Property 4
To prove Property 4, it is only necessary to prove that when the sampling data contains impulsive noise. According to Property 3, The authors only need to prove .
For the MDCO and CRCO algorithms, the authors have
For the generalized covariance of the two variables Y and Z at a given moment, (D.1) can be rewritten as
When the variables Y and Z contain impulsive noise, the authors can get and . Then,
In summary, the MDCO algorithm performs better than the CRCO algorithm under impulsive noise environments, and results in .
E. Proof of Property 5
Under the premise of Property 4, there exists which makes
Meanwhile, is bounded with [24] and can be expressed as
Combining (E.3) with (E.2), the authors can easily get
Transforming (E.4), the authors havewhich means that is bounded.
Data Availability
All the data generated or analyzed in this study are available from the corresponding author on reasonable request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The paper was sponsored by National Key R&D Program of China (no. 2017YFB0702300) and National Natural Science Foundation of China (no. 61971031).
References
 J. Zheng and M. Kaveh, “Sparse spatial spectral estimation: a covariance fitting algorithm, performance and regularization,” IEEE Transactions on Signal Processing, vol. 61, no. 11, pp. 2767–2777, 2013. View at: Publisher Site  Google Scholar
 Q. Li, B. Liao, L. Huang, C. Guo, G. Liao, and S. Zhu, “A robust STAP method for airborne radar with array steering vector mismatch,” Signal Processing, vol. 128, pp. 198–203, 2016. View at: Publisher Site  Google Scholar
 Z. Xiahou and X. Zhang, “Adaptive localization in wireless sensor network through Bayesian compressive sensing,” International Journal of Distributed Sensor Networks, vol. 11, no. 8, 2015. View at: Publisher Site  Google Scholar
 F.G. Yan, Y. Shen, and M. Jin, “Fast DOA estimation based on a split subspace decomposition on the array covariance matrix,” Signal Processing, vol. 115, pp. 1–8, 2015. View at: Publisher Site  Google Scholar
 F. Ma and X. Zhang, “Wideband DOA estimation based on focusing signal subspace,” Signal, Image and Video Processing, vol. 13, no. 4, pp. 675–682, 2019. View at: Publisher Site  Google Scholar
 M. R. Anbiyaei, W. Liu, and D. C. McLernon, “Performance improvement for wideband DOA estimation with white noise reduction based on uniform linear arrays,” in Proceedings of the 2016 IEEE Sensor Array and Multichannel Signal Processing Workshop (SAM), pp. 1–5, IEEE, Rio de Janeiro, Brazil, July 2016. View at: Publisher Site  Google Scholar
 R. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986. View at: Publisher Site  Google Scholar
 C. Qian, L. Huang, and H. C. So, “Improved unitary rootMUSIC for DOA estimation based on pseudonoise resampling,” IEEE Signal Processing Letters, vol. 21, no. 2, pp. 140–144, 2014. View at: Publisher Site  Google Scholar
 R. Roy and T. Kailath, “ESPRITestimation of signal parameters via rotational invariance techniques,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 7, pp. 984–995, 1989. View at: Publisher Site  Google Scholar
 A. Hu, T. Lv, H. Gao, Z. Zhang, and S. Yang, “An ESPRITbased approach for 2D localization of incoherently distributed sources in massive MIMO systems,” IEEE Journal of Selected Topics in Signal Processing, vol. 8, no. 5, pp. 996–1011, 2014. View at: Publisher Site  Google Scholar
 L. Lu and H. Zhao, “Active impulsive noise control using maximum correntropy with adaptive kernel size,” Mechanical Systems and Signal Processing, vol. 87, pp. 180–191, 2017. View at: Publisher Site  Google Scholar
 W.J. Zeng, H. C. So, and L. Huang, “pMUSIC: robust directionofarrival estimator for impulsive noise environments,” IEEE Transactions on Signal Processing, vol. 61, no. 17, pp. 4296–4308, 2013. View at: Publisher Site  Google Scholar
 P. Wang, T.S. Qiu, F.Q. Ren, and A.M. Song, “A robust DOA estimator based on the correntropy in alphastable noise environments,” Digital Signal Processing, vol. 60, pp. 242–251, 2017. View at: Publisher Site  Google Scholar
 J.F. Zhang and T.S. Qiu, “A robust correntropy based subspace tracking algorithm in impulsive noise environments,” Digital Signal Processing, vol. 62, pp. 168–175, 2017. View at: Publisher Site  Google Scholar
 Z. Chen, X. Geng, and F. Yin, “A harmonic suppression method based on fractional lower order statistics for power system,” IEEE Transactions on Industrial Electronics, vol. 63, no. 6, pp. 3745–3755, 2016. View at: Publisher Site  Google Scholar
 X. Jiang, W.J. Zeng, A. Yasotharan, H. C. So, and T. Kirubarajan, “Minimum dispersion beamforming for nonGaussian signals,” IEEE Transactions on Signal Processing, vol. 62, no. 7, pp. 1879–1893, 2014. View at: Publisher Site  Google Scholar
 P. Tsakalides and C. L. Nikias, “The robust covariationbased MUSIC (ROCMUSIC) algorithm for bearing estimation in impulsive noise environments,” IEEE Transactions on Signal Processing, vol. 44, no. 7, pp. 1623–1633, 1996. View at: Publisher Site  Google Scholar
 T.H. Liu and J. M. Mendel, “A subspacebased direction finding algorithm using fractional lower order statistics,” IEEE Transactions on Signal Processing, , vol. 49, no. 8, pp. 1605–1613, 2001. View at: Publisher Site  Google Scholar
 H. Belkacemi and S. Marcos, “Robust subspacebased algorithms for joint angle/doppler estimation in nonGaussian clutter,” Signal Processing, vol. 87, no. 7, pp. 1547–1558, 2007. View at: Publisher Site  Google Scholar
 I. Santamaria, P. P. Pokharel, and J. C. Principe, “Generalized correlation function: definition, properties, and application to blind equalization,” IEEE Transactions on Signal Processing, vol. 54, no. 6, pp. 2187–2197, 2006. View at: Publisher Site  Google Scholar
 B. Chen and J. C. Principe, “Maximum correntropy estimation is a smoothed MAP estimation,” IEEE Signal Processing Letters, vol. 19, no. 8, pp. 491–494, 2012. View at: Publisher Site  Google Scholar
 W. Ma, B. Chen, J. Duan, and H. Zhao, “Diffusion maximum correntropy criterion algorithms for robust distributed estimation,” Digital Signal Processing, vol. 58, pp. 10–19, 2016. View at: Publisher Site  Google Scholar
 Z. Jinfeng and Q. Tianshuang, “A novel covariation based noncircular sources direction finding method under impulsive noise environments,” Signal Processing, vol. 98, pp. 252–262, 2014. View at: Publisher Site  Google Scholar
 J. Zhang, T. Qiu, A. Song, and H. Tang, “A novel correntropy based DOA estimation algorithm in impulsive noise environments,” Signal Processing, vol. 104, pp. 346–357, 2014. View at: Publisher Site  Google Scholar
 R. Patel, M. P. Janawadkar, S. Sengottuvel, K. Gireesan, and T. S. Radhakrishnan, “Effective extraction of visual eventrelated pattern by combining template matching with ensemble empirical mode decomposition,” IEEE Sensors Journal, vol. 17, no. 7, pp. 2146–2153, 2017. View at: Publisher Site  Google Scholar
 W. Liu, P. P. Pokharel, and J. C. Principe, “Correntropy: properties and applications in nonGaussian signal processing,” IEEE Transactions on Signal Processing, vol. 55, no. 11, pp. 5286–5298, 2007. View at: Publisher Site  Google Scholar
 N. Erdol, C. Castelluccia, and A. Zilouchian, “Recovery of missing speech packets using the shorttime energy and zerocrossing measurements,” IEEE Transactions on Speech and Audio Processing, vol. 1, no. 3, pp. 295–303, 1993. View at: Publisher Site  Google Scholar
 T. Dai, W. Lu, W. Wang, J. Wang, and S.T. Xia, “Entropybased bilateral filtering with a new range kernel,” Signal Processing, vol. 137, pp. 223–234, 2017. View at: Publisher Site  Google Scholar
 M. Zhang and B. K. Gunturk, “Multiresolution bilateral filtering for image denoising,” IEEE Transactions on Image Processing, vol. 17, no. 12, pp. 2324–2333, 2008. View at: Publisher Site  Google Scholar
 I. Frosioa, K. Egiazarian, and K. Pulli, “Machine learning for adaptive bilateral filtering,” in Proceedings of the Image Processing: Algorithms and Systems XIII, San Francisco, CA, USA, February 2015. View at: Publisher Site  Google Scholar
 R. J. Kozick and B. M. Sadler, “Maximumlikelihood array processing in nonGaussian noise with Gaussian mixtures,” IEEE Transactions on Signal Processing, , vol. 48, no. 12, pp. 3520–3535, 2000. View at: Publisher Site  Google Scholar
 B. D. Rao and K. V. S. Hari, “Performance analysis of rootMUSIC,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 37, no. 12, pp. 1939–1949, 1989. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Fuqiang Ma et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.