Robust group compressive sensing for DOA estimation with partially distorted observations

EURASIP Journal on Advances in Signal Processing, Dec 2016

In this paper, we propose a robust direction-of-arrival (DOA) estimation algorithm based on group sparse reconstruction algorithm utilizing signals observed at multiple frequencies. The group sparse reconstruction scheme for DOA estimation is solved through the complex multitask Bayesian compressive sensing algorithm by exploiting the group sparse property of the received multi-frequency signals. Then, we propose a robust reconstruction algorithm in the presence of distorted signals. In particular, we consider a problem where the observed data in some frequencies are distorted due to, e.g., interference contamination. In this case, the residual error will follow the impulsive Gaussian mixture distribution instead of the Gaussian distribution due to the fact that some of the estimation errors significantly depart from the mean value of the estimation error distribution. Thus, the minimum least square restriction used in the conventional sparse reconstruction algorithm may lead to a failed reconstruction result. By exploiting the maximum correntropy criterion which is inherently insensitive to the impulsive noise, a weighting vector is derived to automatically mitigate the effect of the distorted narrowband signals, and a robust group compressive sensing approach is developed to achieve reliable DOA estimation. The robustness and effectiveness of the proposed algorithm are verified using simulation results.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://link.springer.com/content/pdf/10.1186%2Fs13634-016-0425-0.pdf

Robust group compressive sensing for DOA estimation with partially distorted observations

Wang et al. EURASIP Journal on Advances in Signal Processing Robust group compressive sensing for DOA estimation with partially distorted observations Ben Wang 0 1 Yimin D. Zhang 0 Wei Wang 1 0 Department of Electrical and Computer Engineering, Temple University , 19122 Philadelphia, PA , USA 1 College of Automation, Harbin Engineering University , Harbin, 150001 Heilongjiang , China In this paper, we propose a robust direction-of-arrival (DOA) estimation algorithm based on group sparse reconstruction algorithm utilizing signals observed at multiple frequencies. The group sparse reconstruction scheme for DOA estimation is solved through the complex multitask Bayesian compressive sensing algorithm by exploiting the group sparse property of the received multi-frequency signals. Then, we propose a robust reconstruction algorithm in the presence of distorted signals. In particular, we consider a problem where the observed data in some frequencies are distorted due to, e.g., interference contamination. In this case, the residual error will follow the impulsive Gaussian mixture distribution instead of the Gaussian distribution due to the fact that some of the estimation errors significantly depart from the mean value of the estimation error distribution. Thus, the minimum least square restriction used in the conventional sparse reconstruction algorithm may lead to a failed reconstruction result. By exploiting the maximum correntropy criterion which is inherently insensitive to the impulsive noise, a weighting vector is derived to automatically mitigate the effect of the distorted narrowband signals, and a robust group compressive sensing approach is developed to achieve reliable DOA estimation. The robustness and effectiveness of the proposed algorithm are verified using simulation results. Compressive sensing; DOA estimation; Group sparsity; Robust group compressive sensing; Maximum correntropy criterion 1 Introduction Direction-of-arrival (DOA) estimation is an important array signal processing technique which finds broad applications in radar, sonar, navigation and wireless communications. In the past decades, a number of high-resolution DOA estimation algorithms, such as Capon estimator [1], multiple signal classification (MUSIC) [2], estimation of parameters via rotation invariance techniques (ESPRIT) [3], and maximum likehood estimator (ML) [4], have been proposed. Recently, compressive sensing and sparse reconstruction methods become attractive in DOA estimation by exploiting the spatial sparsity of the signal arrivals (e.g., [5–8]). Various reconstruction algorithms are available in the literature to obtain desired sparse solutions, including least absolute shrinkage and selection operator (Lasso) [9], orthogonal matching pursuit (OMP) [10], and Bayesian compressive sensing (BCS) [11]. When multiple observations are available due to the use of, e.g., multiple polarizations and/or multiple frequencies, group sparse reconstruction algorithms allow exploitation of multiple observations for better determination of the sparse signal support and thus achieve improved sparse reconstruction performance as compared with conventional methods that do not account for such group sparse property. A number of group sparse reconstruction algorithms, such as group LASSO (gLasso) [12], block orthogonal matching pursuit (BOMP) [13], and multitask Bayesian compressive sensing (MT-BCS) [14] algorithms, have been developed for this purpose. The BCS algorithm based on the relevance vector machine (RVM) [15] has constituted a family of algorithms to recover sparse signals, and it is applied in this paper to solve the sparse reconstruction problem because it © The Author(s). 2016 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. can achieve superior performance and is less sensitive to the coherence of the dictionary entries. To handle complex values involved in the underlying DOA estimation problem, a complex value can be decomposed into the real and imaginary components [16]. In this paper, we exploit the complex multitask compressive sensing (CMTCS) algorithm [17] which achieves an improved performance by exploiting the sparsity pattern shared by the real and imaginary components of the complex-valued observations. In [6], DOA estimation of wideband signals is examined by using a coprime array. In this approach, the wideband signals are divided into multiple subbands. The problem is considered in the context of group sparsity, i.e., the positions of the sparse DOA entries corresponding to the respective DOAs are shared by all the frequency subbands. That is, signal observations obtained from the virtual sensors as a result of difference coarray are fused in the context of group sparse reconstruction. In this paper, we consider a more challenging situation in which a subset of the subbands is distorted due to, e.g., calibration error, filter misfunction, and/or narrowband interference contamination. A direction finding method with sensor gain and phase uncertainties is proposed in [18] which can simultaneously estimate the DOAs of signals and array parameters. However, this method requires that the array perturbation to be small and, due to the joint iteration between DOA and array parameters, it suffers from suboptimal convergence. In [19], a DOA estimation method is developed to utilize partially calibrated array by minimizing a certain cost function. Both methods in [18, 19] have a high computational complexity due to the requirement of spatial search and iterations. By modeling the imperfections of the partially calibrated array, an ESPRIT-based method is proposed in [20] that avoids the spatial search and iterations, thus leading to a low computational complexity. When a subset of the subbands are distorted, the manifold matrix is distorted in these subbands. In this case, the abovementioned DOA estimation methods [18–20] that are developed under the conventional group sparse reconstruction framework and exploit data information across all the subbands will lead to performance degradation. Because the manifold matrix is distorted at the subbands, the residual errors follow an impulsive Gaussian mixture distribution, i.e., some of the error points are far from the mean value of the error distribution which can be viewed as outliers. The effect of such outliers will be amplified in the process of DOA estimation due to the minimum mean square error (MMSE) criterion used in the conventional sparse reconstruction scheme. Similarly, the group sparsity property will be distorted and, as a result, the performance of group sparsity-based DOA estimation will be degraded. In order to achieve robust reconstruction performance, the outliers caused by the subband distortion must be suppressed. As described in [21, 22], the correntropy can be treated as a generalized correlation function characterized by the kernel bandwidth, which controls the observation window and provides an effective mechanism to eliminate the detrimental effect of outliers. Intrinsically different from the threshold-based methods, the correntropy can be applied to nonlinear and non-Gaussian signal processing. Based on the content of correntropy, a new cost function, referred to as maximum correntopy criterion (MCC), is proposed in [22] that achieves robust results for observations corrupted with non-Gaussian noise and other types of outlier observations. The MCC has been successfully applied in various applications, such as radar localization [23] and ellipse fitting problem [24]. Motivated by this fact, the objective of this paper is to develop a novel robust group compressive sensing and sparse reconstruction algorithm that can be applied for DOA estimation when the amplitude and phase information is destroyed for a subset of the subbands. The proposed approach uses the MCC, instead of the conventionally used least square measure, to fuse the data observed at different subbands. By applying a weighting vector that is automatically optimized in the proposed method, the observations at the distorted subbands are adaptively suppressed, leading to a robust sparse reconstruction and DOA estimation performance in the presence of distorted subband signals. The reminder of this paper is organized as following. In Section 2, we present the signal model based on multi-frequency signals. In Section 3, we first introduce a group sparse representation scheme for DOA estimation utilizing multi-frequency signal model and then briefly review the CMT-BCS algorithm, which is used for group sparsity-based DOA estimation. In Section 4, the effect of narrowband interference to the group sparse reconstruction algorithm is examined, and the MCC content is introduced. The proposed robust group sparse reconstruction algorithm for robust DOA estimation utilizing multi-frequency signal is then provided. The performance of the proposed algorithm is evaluated through computer simulations in Section 5. Finally, the conclusion of this paper is drawn in Section 6. Notations: We use lower-case (upper-case) bold characters to denote vectors (matrices). In particular, IN denotes the N × N identity matrix, 1N denotes a N × 1 all-one vector, and 0m×n denotes an all-zero matrix or, when n = 1, an all-zero vector. (·)∗ denotes complex conjugate, and (·)T and (·)H respectively denote the transpose and conjugate transpose of a matrix or vector. vec(·) denotes the vectorization operator that turns a matrix into a vector by concatenating all the columns, and diag{x} denotes a diagonal matrix with the elements of x constituting the diagonal entries. · 2 denotes the 2-norm of a vector, whereas · 1 denotes the 1-norm. E[ ·] denotes the statistical expectation operator. ⊗ denotes the Kronecker product, ◦ presents the Hadamard product, and denotes Khatri-Rao product. x is the smallest integer greater than or equal to x. 2 Signal model Assume there are K far-field sources impinging on an M-element uniform linear array (ULA) with an interelement spacing of d. The signals are assumed to occupy a frequency band consisting of L consecutive or disconnected subbands respectively centering at frequencies f1, f2, · · · , fL. In this section, no signal distortions are considered. The effect of signal distortion of a subset of frequency subbands is discussed in Section 4. At frequency fl, l = 1, · · · , L, the array received signal at the sampling time t can be expressed as k=1 where θk is the DOA of the kth source and s[(lk])(t) is the corresponding signal coefficient at frequency fl. Different source signals s[1l](t), · · · , s[Kl](t) are assumed to be uncorrelated. In addition, for the lth frequency, a[l](θk) is the steering vector corresponding to the kth signal, expressed as a[l](θk) = 1, e−j 2λπld sin θk , · · · , e−j 2λπld (M−1) sin θk where λl = c/fl is the wavelength corresponding to fl and c is the propagation velocity. Furthermore, n˜[l](t) denotes the additive white Gaussian noise vector. After down converting the array received signal vector x˜ [l](t) to baseband, followed by low-pass filtering, the baseband signal corresponding to the lth subband can be represented as k=1 x[l](t) = = A[l]s[l](t) + n[l](t), with l = 1, · · · , L, where A[l] = a[l](θ1), · · · , a[l](θK ) , T s[l](t) = s[(l1])(t), · · · , s[(lK] )(t) , and n[l](t) is the zeromean noise term with variance matrix σ[2l]IM. The covariance matrix of the received signal corresponding to the lth subband is expressed as Rx[l] = E x[l](t)x[Hl](t) = A[l](θ )Rs[l]A[Hl](θ ) + σ[2l]IM, where Rs[l] = E[ s[l](t)s[Hl](t)] = diag{σ12[l], · · · , σK2[l]} is the source covariance matrix with σk2[l] denoting the kth source power at the lth subband. In practice, the covariance matrix in Eq. (4) can be estimated utilizing the collected T samples, i.e., 1 T t=1 3 DOA estimation based on group sparse reconstruction In this section, the group sparse property of the received signals across different subbands is exploited to estimate the DOAs through a group sparsity-based solution. A Bayesian sparse learning method, termed CMT-BCS [17], is utilized to effectively solve the sparse reconstruction problem. 3.1 Group sparse signal representation In order to implement the CS-based DOA estimation, we vectorize the covariance matrix Rˆ x[l] to an M2 × 1 vector, which is expressed as y[l] = vec Rˆx[l] = Aˆ[l](θ )b[l] + σ[2l]i, , and i = vec(IM). = = a[∗l](θk) ⊗ a[l](θk) is the generalized steering vector , b[l] = T Denote the discretized search grid of the DOA as θ˜ = [ θ˜1, θ˜2, · · · , θ˜Q]. We assume that the number of the grid points Q is much greater than the number of sources, i.e., Q K . Then, the manifold dictionary matrix corresponding to the signal at the lth subband can be expressed as A˜[l] = Aˆ [l] θ˜1 , · · · , Aˆ[l](θ˜Q) , where Aˆ [l](θ˜q) = a[∗l](θ˜q) ⊗ a[l](θ˜q). Then, the sparse representation of y[l] in Eq. (6) can be expressed as where vector b˜[l] defines the spatial power spectrum over different grid points. For different subband l, l = 1, · · · , L, the elements of b˜[l] take different values but all L vectors have a common support corresponding to the K sources which sparsely locate in the spatial domain. Thus, the underlying problem can be regarded as a group sparse reconstruction problem of locating the non-zero entries of b˜[l], l = 1, · · · , L. By stacking vectors y[1], · · · , y[L] into a single vector y = [ y[T1], · · · , y[TL]]T , we have the following sparse representation vector y = A˜ b˜ + z, s.t. y − A˜ b˜ b˜ =[ b˜[1], · · · , b˜[L]], and z = [ σ[21]iT , · · · , σ[2L]iT ]T . Note that b˜ has only LK non-zero entries corresponding to the elements of σk2[l], where k = 1, · · · , K and l = 1, · · · , L. The respective positions of non-zero entries in b˜ represent the targets’ DOAs. The elements of b˜ satisfy the group sparse property, i.e., the non-zero entries corresponding to the different frequency combinations share the same support. Thus, the DOA estimation can be achieved by solving the following 1-norm minimization problem where ξ(·) is an operation that obtains the 2-norm of the L-element entries corresponding to each spatial grid points and ε is a user-specific tolerance factor. The CSbased DOA estimation methods have been successfully used for multi-frequency and wideband DOA estimation problems [6, 25]. As we described earlier, a number of group sparse reconstruction algorithms can solve the above 1-norm minimization problem. In this paper, the CMT-BCS algorithm [17], which effectively handles complex group sparse problems in the context of MT-BCS, is utilized due to its superior performance and robustness to dictionary coherence. 3.2 CMT–BCS algorithm The entries of sparse vector b˜[l] in Eq. (8) which characterizes the lth task are assumed to be drawn from the following zero-mean Gaussian distributions: q where b˜[l] is a 2 × 1 vector which is consisted by the realpart coefficient b[rl],q and imagery-part coefficient bi[l],q of the qth element of b˜[l]. αq is the variance of Gaussian probability distribution function (pdf ). It should be note that the parameter α =[ α1, · · · , αQ]T is shared by all the L q groups. Thus, the qth vector b[l] tends to be zero with probability 1 across the L groups when αq is set to zero [14]. To deal with the complex-valued variables in the context of MT-BCS which was originally developed to handle realvalued problems [14], the real and imaginary components of a complex variable are treated as two separate variables in [16] without utilizing the fact that their nonzero entries usually appear at the same positions. The CMT-BCS algorithm treats the real and imaginary components as group sparse and αq is shared and jointly estimated for both real and imagery components [17]. To promote the sparsity of b˜[l], a Gamma prior is placed on αq−1, i.e., (a)−1bax−(a−1)e− bx , and (·) where Gamma(x−1|a, b) = denotes the Gamma function. It has been demonstrated in [15] that a proper selection of the hyper-parameters a and b encourages a sparse representation for the coefficients in b˜[l]. We set a = b = 0 as a default choice which avoids a subjective choice of a and b and leads to simplifications of computation. Considering that the covariance matrix is obtained on basis of received data samples in Eq. (5), a Gaussian prior N (0, β0I2) is also placed on the additive noise. Similarly, β0−1 is placed on the Gamma prior with hyper-parameters c and d, i.e., β0−1 ∼ Gamma(β0−1|c, d). We also let c = d = 0 as a default choice. RI R T I T T R Define b˘[l] = b˘[l] , b˘[l] with b˘ [l] = b[Rl],1, . . . , b[Rl],Q T and b˘[Il] = b[Il],1, . . . , b[Il],Q T respectively denoting the real and imaginary parts. The poste RI rior probability density function for b˘[l] can be evaluated analytically based on the Bayes’ rule as y[RlI] = Re y[l] , Im y[l] T T [l] = β0−1 [Tl] [l] + F−1 −1 , −Im A˜ [l] F = diag α1, · · · , αQ, α1, · · · , αQ . It is noted that, the real and imaginary parts share the same α to ensure their group sparsity. Once the parameters α and β0 are evaluated, the mean and variance of each signal powers corresponding to different spatial grid points can be derived. The logarithm of marginal likehood function is expressed as {α, β0} = arg max L (α, β0) , α,β0 2 [l]μ[l] 2 l=1 1 L = − 2 l=1 + y[RlI] T C[−l]1y[RlI] and C[l] = β0I2M2 + [l]F [Tl]. β0 and α can be evaluated by maximizing Eq. (20) via the expectation maximization (EM) algorithm. α and β0 can be expressed as [17]: C[∗l] −1 RI [l]q denotes the qth column of [l] and b˜[l]q = T is a vector consisted of the real and imaginary part of the corresponding qth entry of b˜[l]. 4 Robust DOA estimation The CMT-BCS sparse reconstruction algorithm described in Section 3 provides effective DOA estimation only when there is no distortions. When the observed signals at some frequencies are contaminated with distortions to yield incorrect amplitude and phase information, these methods will result in erroneous DOA estimation. In this section, we propose an MCC-based robust group compressive sensing method for reliable DOA estimation in such scenarios. In the following, we first analyze the effect of narrowband distortions on the group sparsity-based DOA estimation and then introduce the concept of the MCC [21, 22]. Finally, the proposed robust group compressive sensing method for DOA estimation is presented. 4.1 Effect of narrowband distortion As we discussed earlier, narrowband distortions may result from calibration error and filter misfunction. In addition, local interferer scattering that does not correspond to a clear array manifold would yield a similar effect (Fig. 1). Let L be a subset of the frequency subbands which are contaminated by narrowband distortions or interference. Denote the cardinality of L as |L| = Lo Fig. 1 Array signals with narrowband interference scattering where Lo < L. At subband lo ∈ L, the received signal vector becomes k=1 where a[olo](θk) is the distorted steering vector contaminated by the loth subband interference, which can be expressed by the Hadamard product of steering vector and a N × 1 distortion vector v[lo] with unknown amplitude o and phase, i.e., a[lo](θk) = a[lo](θk) ◦ v[lo] . After converting x˜ [lo](t) to baseband, the corresponding received signal vector can be represented as x[lo](t) = o where A[lo] = A[lo] ◦ V [lo] is the distorted manifold matrix and V = v[lo]1TK is the distortion matrix. The corresponding vectorized covariance matrix of (24) can be written as o o y[lo] = vec Rˆ x[lo] = Aˆ[lo]b[lo] + σ[2l]i, where Rˆ x[lo] = 1/T A[olo] o A[lo], respectively. With the existence of distortion matrix V [lo], the group sparse property in (9) will be destroyed. In addition, under the sparse representation scheme, the distorted manifold matrix leads the residual errors, elo = y[lo] − A˜ [lo] b˜[lo], following the impulsive Gaussian mixture distribution. Note that both the ξ(·) operation and the constraint operation in (10) use l2-norm operations which will generally amplify the effect of outliers so that the estimates DOAs become biased or erroneous unless the distortions is properly suppressed or compensated. 4.2 Maximum correntropy criterion The correntropy between two arbitrary random variables X and Y is defined by [22] (X − Y )2 is the Gaussian kernel function with user-defined kernel size δ that controls the observation window. By adjusting the kernel size, the impulsive Gaussian noise can be effectively eliminated [22]. In practice, the correntropy is calculated by a finite number of samples of xn and yn for n = 1, . . . , N , i.e., 1 N N n=1 y(n) = f (x(n), p) + z(n) In order to derive the MCC, we first give a simple regression model for two random variables y(n) and x(n) as an example: with n ∈[ 1, · · · , N ], where p is a parameter of function f, and z(n) is a noise process. If the probability density function of noise variable z(n) follows the Gaussian distribution, the optimal solution of the above regression problem in the minimum mean square error (MMSE) sense is found by n=1 J(p) = f (x(n), p) − y(n) 2. However, when the noise process z(n) follows a nonGaussian distribution with large outliers, the solution obtained under the above MMSE criterion is no longer optimal. In this case, we consider the optimal solution under the following MCC [22]: J (p) = exp − n=1 f (x(n), p) − y(n) 2 It is shown in [22] that the MCC outperforms MMSE in the impulsive Gaussian mixture noise case since correntropy is inherently insensitive to outliers and the MCC attains the same efficiency compared with MMSE under a Gaussian noise process. 4.3 Proposed robust DOA estimation technique In this subsection, we apply the MCC to develop a robust group sparse reconstruction method, which is used to provide reliable DOA estimation in the presence of distortion in a subset of frequency subbands. Consider that the loth subband, where lo ∈ L, is distorted with the corresponding manifold matrix A˜[lo]. Define yo = where yom is the mth element of yo and A˜m is the mth row of dictionary matrix A˜. Similar as ε in (10), ε1 is a userspecific tolerance parameter. However, the inequality constraint in Eq. (32) is both non-linear and nonconvex, and thus, it is difficult to solve this nonconvex optimization problem. In order to achieve an optimal solution of (32), we introduce the following proposition [26]. Proposition 1 For κ(x) = exp − |xδ2|2 , there exists a conjugate function ϕ such that and for a fixed x, the maximum value is reached at p = −κ(x). By exploiting Proposition 1, we can solve the above nonconvex optimization problem in an iterative method. Using the result presented in Proposition 1, there exists a conjugate function ϕ(ω) for the Gaussian kernel function κ yom − A˜m b˜ ; thus, the constraint function in (32) can be modified as m=1 where ω = ω1, · · · , ωLM2 is a weighting row vector to be optimized and ε2 is a user-specific parameter. Comparing (34) and (33), it is clear that ωm in (34) plays the role of parameter p in (33). Based on Proposition 1, inequality (34) achieves its supreme when ωm is set as . If the block diagonal dictionary matrix designed in Eq. (9) is utilized, the entries of residual error vector yo − A˜ b˜ follow impulsive Gaussian mixture distribution rather than zero-mean Gaussian distribution. Thus, the solution computed from the 1-norm minimization problem in (10) is no longer optimal. In order to solve this problem, the constraint function in (10) needs to be modified. By utilizing the MCC, the constrained 1-norm minimization problem is modified as [23] m=1 |yom − A˜ m b˜|2 m=1 where l = Mm2 . Note that the received signal is divided into L groups corresponding to L subband signals, and it is necessary to automatically adjust the kernel size δ[l] based on the individual subband signals rather than the whole received signals. According to Proposition 1, inequality (34) and the constraint function in (32) are equal when the value of ω is specified [26]. In addition, the conjugate function ϕ(ωm) is a function of parameter ωm. Because ωm keeps constant in one iteration, ϕ(ωm) also keeps constant in each iteration. In this case, the sparse vector b˜ still can be determined when the constraint function in (34) arrives its maximum value with ignoring the value of ϕ(ωm). Hence, the inequality (34) can be reformulated as where ε3 is a user-specific parameter as similar before. For description convenience, the vector ω is divided into L groups, i.e., ω = ω[1], · · · , ω[L] , where ω[l] = ωM2(l−1)+1, · · · , ωM2l corresponds to the lth subband. It is desired that, when there is no distortion in the lth subband, the corresponding ω[l] is a 1 × M2 vector consisting of elements close to 1. Otherwise, for lo ∈ L, ω[l0] should be set to a 1 × M2 vector with small value elements. Thus, the weighting vector ω is adaptively optimized to mitigate the effect of the narrowband distortions. After the weighting vector w is determined, the optimization problem in (32) can be solved through the following optimization problem s.t. − m=1 where ε4 is a user-specific tolerance parameter. The exploitation of the MCC effectively suppresses the effect of the narrowband interference, which leads to a robust DOA estimation based on group sparse reconstruction scheme. In summary, the proposed robust DOA estimation algorithm performs the following iterative process: 1. Initialization: Set the initial kernel size as δ[l] = ∞, and b˜ = 0LM2×1. 2. Compute the initial values of ωn based on Eq. (35). 3. Solve the modified 1-norm minimization problem (37) utilizing group sparse reconstruction algorithm. As shown in Section 3, the CMT-BCS algorithm is applied in this step. 4. Update the initial kernel size δ[l] following the Silverman’s rule [22, 23]: δ[l] = 1.06 × min{δE, R /1.34 } × (M2)−0.4, (38) where δE and R stand for the standard deviation and error interquartile range of (yom − A˜m b˜), respectively, where m ∈[ (l − 1)M2 + 1, lM2] [22]. Repeat steps 2 to 4 until a user-specific number of iteration is reached or the iterations are converged. 4.4 Computational complexity In this subsection, we analyze the computational complexity of the proposed method in terms of the number of multiplications required in each iteration depicted in the previous subsection. The computational complexity to calculate the values of ωm in step 2 requires O{LM2(LQ + 1)} multiplications. In step 3, the numbers of multiplications required for the computation of parameters α and β0 are O{M2LQ(1 + 2M4)} and O L(4Q2M2 + 4Q2M4 + M2 + 2QM2) , respectively. In addition, O{LM2} multiplications are required to update the kernel size. 5 Simulation results In this section, simulation results are presented to demonstrate the performance of the proposed algorithm. In the following simulations, we consider a 7-element ULA. Three far-field targets (K = 3) impinge on the array from directions of −20◦, 5◦, and 25◦, respectively. Assume that each signal occupies five frequency subbands (L = 5). The ratio between the wavelength λl and the interelement spacing d for the five subbbands is respectively 1.2, 1.6, 2.0, 2.4, and 2.8. Assume that the noise power is equal at the five subbands, and the input signal-to-noise ratio (SNR) is 5 dB. The number of snapshot is T = 200, and the dictionary matrix consists of steering vectors corresponding to a uniformly sampled grid between −60◦ to 60◦ with a step size of 0.1◦. We assume that the fifth subband is interfered with an omni-directional narrowband signal with an input interference-to-noise ratio (INR) of 10 dB. The iteration number we used is 10. We first provide the sparse spatial spectrums obtained from the two Bayesian compressive sensing algorithms when there is no narrowband interference in Fig. 2. The true DOAs are marked with red vertical lines. From Fig. 2a, b, it is clear that the CMT-BCS algorithm [17] achieves an outstanding performance compared with the MT-BCS algorithm [11]. Although the DOA estimates obtained from the MT-BCS altorithm are approximately accurate, there exists a spurious peak close to the true -40 -20 0 20 40 Angle (°) (a) Estimated from the MT-BCS algorithm m ru0.8 t c e p laS0.6 it a p dS0.4 e z il a rm0.2 o N m ru0.8 t c e p laS0.6 it a p dS0.4 e z li a rm0.2 o N m ru0.8 t c e p laS0.6 it a p dS0.4 e z li a rm0.2 o N m ru0.8 t c e p laS0.6 it a p dS0.4 e z il a rm0.2 o N -40 -20 0 20 40 Angle (°) (b) Estimated from the CMT-BCS algorithm Fig. 2 Estimated spatial spectrum using the conventional method implemented by the a MT-BCS and b CMT-BCS methods, where no narrowband distortions are present Fig. 3 Estimated spatial spectrum in the presence of narrowband interference using the conventional method implemented by both a MT-BCS and b CMT-BCS algorithms DOA peaks. Correspondingly, the CMT-BCS algorithm achieves an accurate spatial spectrum without fake peaks. In Fig. 3, we provide the estimated spatial spectra using the two BCS algorithms which are not robust to the presence of narrowband signal interference. It is easy to see that both of the algorithms fail to achieve an accurate DOA estimation. From Fig. 3a, which provides the estimated spatial spectrum obtained from the MT-BCS algorithm, we can see that there exists many spurious peaks which have higher peak values than the true DOA peaks. In Fig. 3b, the spatial spectrum estimated from the CMT-BCS algorithm only captures one signal arrival with a bias, whereas it fails to estimate the DOAs of the other two signals. On the other hand, two fake spurious peaks are generated instead. The estimated spatial spectrum obtained from the proposed algorithm is given in Fig. 4. The two sparse reconstruction algorithms are respectively utilized in step 3 of the proposed algorithm. From the two plots in Fig. 4, we can see that the narrowband interference signal is eliminate efficiently and accurate DOA estimations are obtained. However, there exist spurious components in the estimated spatial spectrum when the MT-BCS algorithm is applied. On the other hannd, no such spurious components are in the spatial spectrum obtained from the CMT-BCS algorithm. Define the root mean square error (RMSE) as follows: 1 K k=1 where θˆk is the estimated DOA of the kth source. The RMSE versus SNR is provided in Fig. 5 where the results are averaged over 100 independent trails using the proposed robust DOA estimation algorithm and the input SNR varies from −5 to 25 dB. It is clearly shown that the proposed method achieves a robust estimation and the reconstruction performance is significantly improved by using the CMT-BCS reconstruction algorithm. In the next example, we compare the performance of the proposed method under different proportions of distorted subbands in Fig. 6. We consider that three cases of subband distortion, i.e., 1 subband (20%), 2 subbands (40%), and 3 subbands (60%), are randomly distorted. While the RMSE performance generally degrades as the number of distorted subbands increases, the proposed m ru0.8 t c e p laS0.6 it a p dS0.4 e z li a rm0.2 o N m ru0.8 t c e p laS0.6 it a p dS0.4 e z li a rm0.2 o N ()°0.8 E S M0.6 R 0.4 Fig. 5 Comparison of the RMSE performance versus the input SNR Fig. 4 Estimated spatial spectrum in the presence of narrowband interference using the proposed method implemented by both a MT-BCS and b CMT-BCS algorithms method still provides robust DOA estimation in all the three considered cases. 6 Conclusions In this paper, a robust group compressive sensing algorithm was developed to achieve reliable DOA estimation in the presence of narrowband signal distortion due to interference and calibration errors. We first construct the Proposed method using CMT-BCS Proposed method using MT-BCS group sparse signal model utilizing the received multifrequency signals. Then, we review the complex mutitask Bayesian compressive sensing (CMT-BCS) technique to solve the complex value group sparse reconstruction problem. Consider narrowband signal distortions that affect the manifold matrix in some of the frequency bands. In this case, the residual error between the observations and estimated signals follows the impulsive Gaussian mixture distribution which is characterized as outliers. Thus, it leads to performance degradation of DOA estimation when conventional group sparse reconstruction methods are used. In order to achieve robust DOA estimation, a robust group sparse reconstruction method is developed based on the maximum correntopy criterion. In this approach, an adaptive weight is applied to each subband observation to mitigate such outliers caused by the narrowband distortion. The reliable DOA estimation capability is verified in such challenging situations where conventional group compressive sensing methods fail. Acknowledgments B. Wang was supported by the China Scholarship Council for his stay at the Temple University. The work of Y. D. Zhang was supported in part by the National Science Foundation under grant AST-1547420. The work of W. Wang was supported by the National Natural Science Foundation (61571148), China Postdoctoral Special Funding (2015T80328), and China Postdoctoral Science Foundation Grant (2014M550182). Competing interests The authors declare that they have no competing interests. 1. J Capon, High-solution frequency-wavenumber spectrum analysis . Proc. IEEE . 57 ( 8 ), 1408 - 1418 ( 1969 ) 2. RO Schmidt, Multiple emitter location and signal parameter estimation . IEEE Trans. Antennas Propagat . 37 (12), 276 - 280 ( 1986 ) 3. R Roy, A Paulraj, T Kailath, ESPRIT-a subspae rotation approach to estimation of parameters of Cisoids in noise . IEEE Trans. Acoust. Speech Signal Proc . 34 ( 5 ), 1340 - 1342 ( 1986 ) 4. I Ziskind, M Wax, Maximum likelihood localization of multiple sources by alternate projection . IEEE Trans. Acoust. Speech Signal Process . 36 (10), 1553 - 1560 ( 1988 ) 5. YD Zhang, MG Amin, B Himed, in Proc. IEEE ICASSP . Sparsity-based DOA estimation using co-prime arrays , (Vancouver, Canada, 2013 ) 6. Q Shen, W Liu, W Cui, S Wu, YD Zhang, MG Amin, Low-complexity wideband direction-of-arrival estimation based on co-prime arrays , IEEE/ACM Trans. Audio Speech Language Process . 23 ( 9 ), 1445 - 1456 ( 2015 ) 7. S Qin, YD Zhang, MG Amin, Generalized coprime array configurations for direction-of-arrival estimation . IEEE Trans.Signal Process . 63 ( 6 ), 1377 - 1390 ( 2015 ) 8. A Koochakzadeh , P Pal, Sparse source localization using perturbed arrays via bi-affine modeling . Digital Signal Process . in press. doi:10.1016/j.dsp. 2016 .06.004 9. R Tibshirani, Regression shrinkage and selection via the Lasso . JR Statist.Soc . 58 ( 1 ), 267 - 288 ( 1996 ) 10. JA Tropp, AC Gilbert, Signal recovery from random measurements via orthogonal matching pursuit . IEEE Trans. Inf. Theory . 53 (12), 4655 - 4666 ( 2007 ) 11. S Ji, Y Xue, L Carin, Bayesian compressive sensing . IEEE Trans. Signal Process . 56 ( 6 ), 2346 - 2356 ( 2008 ) 12. M Yuan, Y Lin, Model selection and estimation in regression with grouped variables . J. Royal Stat. Soc. Ser. B . 68 ( 1 ), 49 - 67 ( 2007 ) 13. Y Eldar , P Kuppinger, H Bölcskei, Block-sparse signal: uncertainty relations and efficient recovery . IEEE Trans. Signal Process . 58 ( 6 ), 3042 - 3054 ( 2010 ) 14. S Ji, D Dunson, L Carin, Multitask compressive sensing . IEEE Trans. Signal Process . 57 ( 1 ), 92 - 106 ( 2009 ) 15. ME Tipping, Sparse Bayesian shrinkage and selection learning and the relevance vector machine . J. Royal Statist. Soc. Ser. B . 68 ( 9 ), 211 - 244 ( 2001 ) 16. M Carlin, P Roccca, G Oliveri, F Viani, A Massa, Directions-of-arrival estimation through Bayesian compressive sensing strategies . IEEE Trans. Antennas Propagat . 61 ( 7 ), 3828 - 3838 ( 2013 ) 17. Q Wu, YD Zhang, MG Amin, B Himed, in Proc. IEEE ICASSP. Complex multitask Bayesian compressive sensing, (Florence, Italy , 2014 ), pp. 3375 - 3379 18. A Weiss, B Freidlander, Eigenstructure methods for direction finding with sensor gain and phase uncertainties . Circuits Syst. Signal Process . 9 ( 3 ), 271 - 300 ( 1990 ) 19. A Weiss, B Friedlander, DOA and steering vector estimation using a partially calibrated array . IEEE Trans. Aerospace Electron. Syst . 32 ( 3 ), 1047 - 1057 ( 1996 ) 20. B Liao, S Chan, Direction finding with partly calibrated uniform linear arrays . IEEE Trans. Antennas Propagat . 60 ( 2 ), 922 - 929 ( 2012 ) 21. I Santamaría, PP Pokharel, JC Principe, Generalized correlation function: definition, properties and application to blind equalization . IEEE Trans. Signal Process . 54 ( 6 ), 2187 - 2197 ( 2006 ) 22. W Liu, P Pokharel , JC Principe, Correntropy: properties and application in non-Gaussian signal processing . IEEE Trans. Signal Process . 55 (11), 4634 - 4643 ( 2007 ) 23. J Liang, D Wang, L Su, B Chen, H Chen , HC So, Robust MIMO radar target localization via nonconvex optimization . Signal Process . 122 , 33 - 38 ( 2016 ) 24. J Liang, Y Wang, X Zeng, Robust ellipse fitting via half-quadratic and semidefinite relaxation optimization . IEEE Trans. Image Process . 24 (11), 4276 - 4286 ( 2015 ) 25. YD Zhang, MG Amin, F Ahmad, B Himed, in Proc. IEEE Int. Workshop on Comput. Adv. Multi-Sens. Adaptive Process . DOA estimation using a sparse uniform linear array with two CW signals of co-prime frequencies, ( Saint Martin , 2013 ), pp. 404 - 407 26. X Yuan, B Hu, in Proc. Int. Conf. Machine Learning. Robust feature extraction via information theoretic learning , (Montreal, Canada, 2009 ), pp. 1193 - 1200


This is a preview of a remote PDF: http://link.springer.com/content/pdf/10.1186%2Fs13634-016-0425-0.pdf

Robust group compressive sensing for DOA estimation with partially distorted observations, EURASIP Journal on Advances in Signal Processing, 2016, DOI: 10.1186/s13634-016-0425-0