Next Article in Journal
Open-Set Recognition Model for SAR Target Based on Capsule Network with the KLD
Next Article in Special Issue
Four-Dimensional Parameter Estimation for Mixed Far-Field and Near-Field Target Localization Using Bistatic MIMO Arrays and Higher-Order Singular Value Decomposition
Previous Article in Journal
Accurate Deformation Retrieval of the 2023 Turkey–Syria Earthquakes Using Multi-Track InSAR Data and a Spatio-Temporal Correlation Analysis with the ICA Method
Previous Article in Special Issue
Joint Antenna Scheduling and Power Allocation for Multi-Target Tracking under Range Deception Jamming in Distributed MIMO Radar System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Direction-of-Arrival Estimation in the Presence of Outliers and Noise Nonuniformity

1
State Key Laboratory of Mechanics and Control for Aerospace Structures, Nanjing University of Aeronautics and Astronautics, No. 29 Yudao Street, Nanjing 210016, China
2
Shanghai Aircraft Design and Research Institute, State Key Laboratory of Airliner Integration Technology and Flight Simulation, Shanghai 201210, China
3
College of Electronics and Information Engineering, Shenzhen University, Shenzhen 518060, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(17), 3140; https://doi.org/10.3390/rs16173140
Submission received: 10 July 2024 / Revised: 16 August 2024 / Accepted: 19 August 2024 / Published: 26 August 2024
(This article belongs to the Special Issue Array and Signal Processing for Radar)

Abstract

:
In direction-of-arrival (DOA) estimation with sensor arrays, the background noise is usually modeled to be uncorrelated uniform white noise, such that the related algorithms can be greatly simplified by making use of the property of the noise covariance matrix being a diagonal matrix with identical diagonal entries. However, this model can be easily violated by the nonuniformity of sensor noise and the presence of outliers that may arise from unexpected impulsive noise. To tackle this problem, we first introduce an exploratory factor analysis (EFA) model for DOA estimation in nonuniform noise. Then, to deal with the outliers, a generalized extreme Studentized deviate (ESD) test is applied for outlier detection and trimming. Based on the trimmed data matrix, a modified EFA model, which belongs to weighted least-squares (WLS) fitting problems, is presented. Furthermore, a monotonic convergent iterative reweighted least-squares (IRLS) algorithm, called the iterative majorization approach, is introduced to solve the WLS problem. Simulation results show that the proposed algorithm offers improved robustness against nonuniform noise and observation outliers over traditional algorithms.

1. Introduction

Usually, in the problem of direction-of-arrival (DOA) estimation, the background noise is assumed to be uniform white noise, i.e., all sensor noises constitute a zero-mean Gaussian process with the covariance matrix σ 2 I , where σ 2 is the noise variance and I is the identity matrix [1]. with such an assumption, the problem of DOA estimation can be much simplified. For instance, in subspace-based approaches such as multiple signal classification (MUSIC) and the estimation of signal parameters via rotational invariant techniques (ESPRIT), the signal subspace and noise subspace can be simply separated by the eigenvalue decomposition of the array covariance matrix. In those maximum-likelihood (ML) algorithms [2,3,4], one can concentrate the resultant log-likelihood (LL) function with respect to both signal waveform and noise nuisance parameters. Therefore, the dimension of unknown parameter space and the associated computational burden can be significantly reduced.
However, in practical applications, the noise may be time-varying because of temperature drift, and in some system implementations, the sensor noise variances are not identical to each other due to the difference between the sensor locations and the associated noise environment. As a result, the uniform white noise assumption may be violated, and in general, the classical methods aforementioned cannot provide satisfactory performance in these situations. As discussed in [1], in some applications, such as sparse array systems, the sensor noises are spatially uncorrelated. However, the noise variances of the sensors may be different from each other. This is probably caused by the nonuniformity of sensor noise or the imperfection of array calibration. In this case, the noise covariance matrix can still be modeled as a diagonal matrix, whereas its diagonal elements are no longer identical. This kind of noise is named nonuniform white noise, and it is one of the main aspects that we shall study in this paper.
Many efforts have been devoted to the problem of DOA estimation in the presence of nonuniform noise [1,2,3,4,5,6,7,8]. A series of ML-based algorithms were proposed in [1,2,3,4,5,6], and deterministic/stochastic ML DOA estimators were derived [1,6]. For implementation, the authors gave an iterative procedure, including a stepwise concentration of the LL function with respect to the signal and noise nuisance parameters. For some specific cases, e.g., where the number of sensors is no smaller than three times the number of sources, various methods were proposed for noise covariance matrix estimation, which can be utilized to prewhiten the observations or removed from the array covariance matrix [7,8]. These methods have lower computational complexity compared with the ML DOA estimators. In particular, Liao et al. have made significant contributions on this topic by proposing a series of algorithms for DOA estimation in nonuniform noise, including the iterative ML subspace estimation (IMLSE) and iterative least-squares subspace estimation (ILSSE) algorithms [9], the partly calibrated array with nonuniform noise [10], eigendecomposition- and rank minimization-based approaches [11], the spatial smoothing-based method [12], and the matrix completion-based method [13]. More recently, a sparse reconstruction-based approach was proposed in [14] and a low-rank matrix recovery-based method was reported in [15].
Besides the nonuniformity of sensor noise, the other aspect that may cause significant performance degradation of the conventional Gaussian assumption-based array signal processing methods is the data outliers. Moreover, it is known that in practical systems, one of the key reasons for the presence of data outliers is the existence of unexpected impulsive noise. Due to the importance of this problem, a large number of robust methods against impulsive noise have been developed. In [16], an expectation maximization (EM) algorithm is proposed to estimate the source locations, signal waveforms, and noise distribution parameters. In [17], the Shapiro–Wilk goodness-of-fit W test is utilized to suppress impulsive noise by trimming the outliers prior to covariance estimation, so that the impact of impulsive noise on conventional DOA estimation algorithms can be minimized. It is shown in [18] that as a powerful robust statistical technique, the M-estimator [19,20,21] can be incorporated to the traditional projection approximation subspace tracking (PAST) algorithm [22]. Therefore, the impulse-corrupted data (outliers) can be detected and prevented from corrupting the subspace estimate, and hence, the subspace-based DOA estimators can be properly performed.
In this paper, we develop a more practical and general DOA estimator, which takes both the nonuniform noise and data outliers into account. More precisely, the background noise is assumed to be nonuniform white noise with an unknown covariance matrix. Furthermore, the collected observation data are corrupted by outliers. In order to handle these two aspects simultaneously in the problem of direction finding, a novel robust exploratory factor analysis (EFA)-based DOA estimator is introduced. Firstly, an approach to DOA estimation in the presence of nonuniform noise is derived based on the EFA model [23,24,25,26,27]. Next, with the help of generalized extreme Studentized deviate (ESD) test [28], a modified EFA model is proposed to combat the hostile effect of outliers. We show that the modified EFA model can be deemed a weighted least-squares (WLS) fitting problem. In order to solve this problem, a monotonic convergent iterative reweighted least-squares (IRLS) algorithm [29] is employed. After the subspace is robustly estimated through the above procedure, conventional algorithms such as MUSIC for DOA estimation can be applied directly. Simulation results demonstrate that the proposed method performs well in the presence of nonuniform noise and offers improved robustness against outliers over traditional algorithms.
The remainder of this paper is organized as follows. The EFA-based subspace and DOA estimation algorithms are first introduced in Section 2. The proposed methods for DOA estimation in the presence of nonuniform noise and observation outliers are presented in Section 3 and Section 4. Numerical examples are conducted in Section 5 to evaluate the performance of the proposed method. Finally, Section 6 concludes the paper.

2. Signal Model

Consider an array of M sensors receiving L uncorrelated narrowband source signals. The observation vector can be written as
x ( i ) = As ( i ) + n ( i ) ,
where A = [ a ( θ 1 , ϕ 1 ) , , a ( θ L , ϕ L ) ] denotes the steering matrix with a ( θ , ϕ ) being the steering vector, and s ( i ) = [ s 1 ( i ) , , s L ( i ) ] T and n ( i ) = [ n 1 ( i ) , , n M ( i ) ] T are the signal waveform vector and additive noise measurement vector, respectively. More specifically, the steering vector can be written as
a ( θ , ϕ ) = [ e j w τ 1 ( θ , ϕ ) , , e j w τ M ( θ , ϕ ) ] T ,
where w = 2 π f , f is the carrier frequency, and τ m ( θ , ϕ ) = a T p m a T p m c c , with c being the propagation velocity. e and p m , respectively, are the unite vector and m th sensor position:
e = sin θ cos φ sin θ sin φ cos φ , p m = p x m p y m p z m .
The source signals are assumed to be temporally uncorrelated zero-mean Gaussian processes with E { s ( i ) } = 0 and E { s ( i ) s H ( i ) } = P , where E { · } denotes statistical expectation. Moreover, the covariance matrix P is diagonal, and its diagonal entries represent the signal powers. For the nonuniform white noise considered in this paper, we have
E { n ( i ) n H ( i ) } = Q = d i a g { σ 1 2 , σ 2 2 , , σ M 2 } ,
where σ m 2 denotes the noise power in the m th sensor. It can be seen that if σ 1 2 = σ 2 2 = = σ M 2 , the above model is reduced to that of uniform white noise. For the convenience of the following derivation, we rewrite the noise vector as follows:
n ( i ) = Φ u ( i ) ,
where Φ = Q 1 / 2 = d i a g { σ 1 , σ 2 , , σ M } is a diagonal matrix representing the standard deviation of each sensor noise, and u ( t ) is a spatially and temporally uncorrelated standard normal distribution such that E { u ( i ) } = 0 and E { u ( i ) u H ( i ) } = I M , where I M is an M × M identity matrix. According to (1) and (5), we have
x ( i ) = As ( i ) + Φ u ( i ) .
Assuming that N snapshots are collected, the observation data matrix can be written compactly as
X = S A H + N = S A H + U Φ H ,
where X = [ x ( 1 ) , , x ( N ) ] H C N × M , S = [ s ( 1 ) , , s ( N ) ] H C N × L , N = [ n ( 1 ) , , n ( N ) ] H C N × M , and U = [ u ( 1 ) , , u ( N ) ] H C N × M . Given a sufficiently large N, we have
1 N S H S = P , 1 N N N H = Q , and 1 N U U H = I N .
Making use of the above identities, Equation (7) can be further reformulated as
X = S o A o H + U o Φ o H ,
where S o = N 1 / 2 S P 1 / 2 , A o = N 1 / 2 A P 1 / 2 , U o = N 1 / 2 U , and Φ o = N 1 / 2 Φ . Moreover, we have
S o H S o = I L , U o H U o = I M .
It is assumed that the source signals and noise are uncorrelated, so that we have U H S = 0 M × L , and hence
U o H S o = 0 M × L .
In this paper, we focus on the problem of DOA estimation based on the subspace, which spans the same space as the steering vector matrix A . Once such a subspace estimate is available, the conventional subspace-based DOA estimation algorithms can be applied. Moreover, from the identity A o = N 1 / 2 A P 1 / 2 , it is known that the A o spans the same space as A .

3. EFA-Based DOA Estimation

In this section, we shall focus on estimating A o by taking advantage of the EFA model. More precisely, according to (9)–(11), the EFA model can be described as
X = S o A o H + U o Φ o H , s . t . S o H S o = I L , U o H U o = I M , U o H S o = 0 M × L , and Φ o is diagonal .
Consequently, based on the LS goodness-of-fit criterion, the corresponding problem for subspace (or say A o ) estimation is given by
min X S o A o H U o Φ o H F 2 s . t . S o H S o = I L , U o H U o = I M , U o H S o = 0 M × L Φ o is diagonal .
In order to transform the complex-valued problem (13) into a real-valued one, the following proposition will be applied:
Proposition 1.
Let B be a complex-valued matrix and the corresponding real-valued matrix B ¯ be defined as
B ¯ = Re { B } Im { B } Im { B } Re { B } ;
then, we have
B F 2 = 1 2 B ¯ F 2 .
Further, given another complex-valued matrix C , and letting G = B C H , we have
G ¯ = B ¯ C ¯ T
where C ¯ and G ¯ are real-valued matrices defined according to (14).
Proof. 
See Appendix A. □
Following the above proposition, the objective function in (13) can be rewritten as
X S o A o H U o Φ o H F 2 = 1 2 X ¯ S ¯ o A ¯ o T U ¯ o Φ ¯ o T F 2 ,
where S ¯ o , A ¯ o , U ¯ o , and Φ ¯ o are similarly defined as (14). Moreover, the constraints in (13) can be equivalently represented as
S o H S o = I L S ¯ o T S ¯ o = I 2 L ,
U o H U o = I M U ¯ o T U ¯ o = I 2 M ,
U o H S o = 0 M × L U ¯ o T S ¯ o = 0 2 M × 2 L .
Consequently, the problem in (13) can be equivalently transformed to the following real-value optimization problem:
min X ¯ S ¯ o A ¯ o T U ¯ o Φ ¯ o T F 2 s . t . S ¯ o T S ¯ o = I 2 L , U ¯ o T U ¯ o = I 2 M , U ¯ o T S ¯ o = 0 2 M × 2 L Φ ¯ o is diagonal .
Now, an iterative method will be introduced to solve the above problem. Firstly, we define two matrices E ¯ and F ¯ , respectively, as
E ¯ = [ S ¯ o U ¯ o ] R 2 N × ( 2 L + 2 M ) ,
F ¯ = [ A ¯ o Φ ¯ o ] R 2 M × ( 2 L + 2 M ) .
Then, the problem in (19) can be rewritten as
min X ¯ E ¯ F ¯ T F 2 s . t . E ¯ T E ¯ = I 2 ( M + L ) .
It can be noted that for a given F ¯ , the problem in (21) is an orthogonal Procrustes problem, and its solution is given by [30]
E ¯ = M N T ,
where M and N are obtained from the economy SVD of the matrix X ¯ F ¯ as
X ¯ F ¯ = SVD MD N T ,
where D is a diagonal matrix composed of the singular values.
According to (20a), once E ¯ has been obtained, S ¯ o and U ¯ o can be respectively determined from the first 2 L and the last 2 M columns of E ¯ , i.e., S ¯ o = E ¯ : , 1 : 2 L , U ¯ o = E ¯ : , 2 L + 1 : 2 L + 2 M . In addition, taking advantage of the properties of S ¯ o and U ¯ o as shown in (18a), (18b) and (18c), one obtains
S ¯ o T X ¯ = S ¯ o T ( S ¯ o A ¯ o T + U ¯ o Φ ¯ o T ) = A ¯ o T ,
U ¯ o T X ¯ = U ¯ o T ( S ¯ o A ¯ o T + U ¯ o Φ ¯ o T ) = Φ ¯ o T .
This indicates that A ¯ o and Φ ¯ o can be respectively updated as
A ¯ o = X ¯ T S ¯ o , Φ ¯ o = X ¯ T U ¯ o .
Meanwhile, the matrix F ¯ = [ A ¯ o Φ ¯ o ] can be further updated. It can be seen that E ¯ and F ¯ are estimated iteratively, and the procedure should be stopped after a certain convergence criterion is met, such as the difference between two consecutive f L S s (where f L S = X ¯ E ¯ F ¯ T F 2 ) being less than the prescribed threshold or the maximal iterations being reached. Finally, A o can be extracted from A ¯ o as
A o = A ¯ o ( 1 : M , 1 : L ) + j A ¯ o ( 1 : M , L + 1 : 2 L ) ,
and it is then applied to the conventional high-resolution direction-finding algorithms. In this paper, the MUSIC algorithm is employed, and hence, the DOAs can be estimated from the following spectrum:
G ( θ , ϕ ) = 1 | a H ( θ , ϕ ) ( I M A ˜ o A ˜ o H ) a ( θ , ϕ ) | ,
where A ˜ o is an orthonormal matrix of A o . In summary, the main steps of the proposed EFA-based method for DOA estimation are listed in Algorithm 1.
Algorithm 1 Proposed EFA-based Method for DOA Estimation in Nonuniform Noise
1:
Initialization: Set S ¯ o , U ¯ o , A ¯ o = X ¯ T S ¯ o , E ¯ = [ S ¯ o U ¯ o ] , Φ ¯ o = d i a g ( X ¯ T U ¯ o ) , and F ¯ = [ A ¯ o Φ ¯ o ] , f = X ¯ F 2 , f L S = X ¯ E ¯ F ¯ T F 2
2:
while | f f | > ε or max iteration number unreached do
3:
    X ¯ F ¯ = S V D MD N T (economy SVD)
4:
    E ¯ = M N T
5:
    S ¯ o = E ¯ : , 1 : 2 L
6:
    U ¯ o = E ¯ : , 2 L + 1 : 2 L + 2 M
7:
    A ¯ o = X ¯ T S ¯ o , Φ ¯ o = X ¯ H U ¯ o
8:
    E ¯ = [ S ¯ o U ¯ o ] , F ¯ = [ A ¯ o Φ ¯ o ]
9:
    f = f , f = X ¯ E ¯ F ¯ T F 2
10:
end while
11:
A o = A ¯ o ( 1 : M , 1 : L ) + j A ¯ o ( 1 : M , L + 1 : 2 L )
12:
Orthonormalize A o to obtain A ˜ o
13:
Spectrum G ( θ , ϕ ) = 1 | a H ( θ , ϕ ) ( I M A ˜ o A ˜ o H ) a ( θ , ϕ ) |
14:
Output: DOA estimates.

4. Robust EFA-Based DOA Estimation against Outliers

In the above section, we show that the subspace—and hence, DOAs—can be estimated in cases where the background noise is nonuniform white. However, it is worth noting that the above approach does not take the data outliers into account. Unfortunately, the standard LS problem (21) is not robust against—and is unstable to—the data outliers. Outlying data usually result in an effect so strong in the minimization that the parameters and estimates are distorted. In this paper, to combat the hostile effect of data outliers, a simple robust statistical technique is used to detect the outliers. Then, the detected outliers with extremely large values are treated as missing data, which are then replaced by some mild values. Moreover, relatively small weights are attached to these data points in order to indicate their importance, and hence, a WLS problem is formulated. The resulting approach is thus named the robust EFA-based method, which extends the EFA-based method to the scenarios with both nonuniform noise and data outliers.

4.1. Outlier Detection

Many methods have been developed in the past several decades to detect data outliers; the interested reader is referred to [28,31,32,33,34] and related references therein. In this paper, a method called the generalized ESD test is employed [28]. The main reason for choosing this approach is that it is unnecessary to exactly specify the number of outliers. Instead, only an upper bound on the suspected number of outliers is required, and hence, it is the recommended test when the exact number of outliers is not known. As a result, the generalized ESD test is performed both on the real and imaginary parts of the samples of each array channel.
Let x m r ( i ) denote the real part of the i th observation of the m th sensor, and
x m = [ x m r ( 1 ) , x m r ( 2 ) , , x m r ( N ) ] .
To detect whether there are outliers in the sample set x m , the hypothesis is defined as
H 0 : No outliers in x m , H 1 : Up to K outliers in x m ,
where K denotes the upper bound of the number of the outliers, and it is assumed to be known. The statistics R 1 , R 2 ,…, R K , i.e., extreme Studentized deviates, are computed according to the following rule:
R 1 = max | x m r ( 1 ) x ¯ m r | , , | x m r ( N ) x ¯ m r | s ,
where x ¯ m r and s are the sample mean and sample standard deviation, respectively, of x m , i.e.,
x ¯ m r = 1 N i = 1 N x m r ( i ) ,
s = 1 N 1 i = 1 N x m r ( i ) x ¯ m r 2 .
Then, the i th sample x m r ( i ) associated with R 1 is removed from the sample set x m . This yields to a reduced set x m with N 1 samples, and R 2 is calculated similarly as R 1 from x m . The above process is repeated until R k is computed. Next, critical values of the test are determined as
λ k = ( N k ) t p , N k 1 ( N k 1 + t p , N k 1 2 ) ( N k + 1 ) ,
where k = 1 , , K , t p , v represents the percentile of a t distribution with v degrees of freedom, and p is given by
p = 1 α 2 ( N k + 1 ) ,
where α is the prescribed significance level, which is typically selected to be 0.05.
with the statistics R 1 , R 2 , , R K and critical values λ 1 , λ 2 , , λ K , one can then verify the outliers in the sample set x m . If R k λ k for k { 1 , , K } ; then, one can declare that there are no outliers in x m , i.e., H 0 is chosen. Otherwise, if R k > λ k , then the sample in x m associated with R i is flagged as an outlier, and hence, H 1 is chosen. It is worth noting that we do not focus on outlier detection. Instead, further processing of these outliers should be performed to guarantee that the future estimation problem can be carried out properly. Thus, one needs to consider how to deal with outliers. Much attention has been dedicated to the problem of what to do with identified outliers [31]. In this paper, a common method is employed, i.e., the outlying data points are removed and then replaced by some mild values, such as the median of the data set. Therefore, the observation data matrix X ¯ is replaced by a trimmed data matrix X ¯ , where the outlying data points have been preprocessed.

4.2. Robust EFA-Based Subspace Estimation

In this section, the problem of robust subspace estimation using the trimmed data matrix X ¯ is presented. Firstly, we reconsider the objective function in the standard LS problem (21), and obtain
f L S = X ¯ E ¯ F ¯ T F 2 = i = 1 N j = 1 M | e i j | 2 ,
where e i j denotes the ( i , j ) th entry of the residual matrix X ¯ E ¯ F ¯ T . However, when there are outliers in the data matrix X ¯ , the corresponding residuals will tend to be abnormal. In this case, minimizing f L S will lead to an unsatisfactory solution. with the modified data matrix X ¯ , the corresponding standard LS problem is readily given by
f L S = X ¯ E ¯ F ¯ T F 2 = i = 1 N j = 1 M | e i j | 2 ,
if the trimmed data points in X ¯ are treated equally as others. On the other hand, a more natural way is that relatively small weights are attached to those trimmed data points to indicate the importance of them. Consequently, a WLS problem is formulated as
f W L S = ( X ¯ E ¯ F ¯ T ) W F 2 = i = 1 N j = 1 M | w i j e i j | 2 ,
where W denotes the weight matrix, with w i j being a non-negative weight value associated with the residual e i j ; and ∘ denotes the Hadamard product (matrix element-wise product). Generally, the weight value w i j can be given as
w i j = γ , X ¯ i , j is a trimmed value , 1 , otherwise ,
where 0 < γ < 1 . Therefore, according to (21), the resultant problem is
min ( X ¯ E ¯ F ¯ T ) W F 2 s . t . E ¯ T E ¯ = I 2 ( M + L ) .
Similar to the normal case as described in Section 3, once the above problem is solved, the matrix A o standing for the subspace can be obtained.

4.3. Solution of the WLS Problem

In this subsection, a monotonic convergent IRLS algorithm, i.e., the so-called iterative majorization approach [29], is utilized to solve the optimization problem in (38). In this approach, one can minimize a majoring function instead of minimizing the original complicated objective function directly. Let Z be the parameter matrix lying in a certain domain Ω and Z c be the current estimate of Z (called the supporting point): the majorizing function m a j ( Z | Z c ) and the original objective function f ( Z ) must satisfy
m a j ( Z | Z c ) f ( Z ) , Z Ω ,
m a j ( Z c | Z c ) = f ( Z c ) .
Hence, given the supporting point Z c , the corresponding problem for updating the estimate can be formulated as
Z u = arg min m a j ( Z | Z c ) .
According to (39a), (39b) and (40), it can be found that f ( Z u ) m a j ( Z u | Z c ) m a j ( Z c | Z c ) = f ( Z c ) . This implies that by iteratively minimizing the majorizing function m a j ( Z | Z c ) , a sequence of monotonically decreasing loss function values of f ( Z ) can be obtained. For the function bounded below, such an iterative procedure terminates when no significant gains are obtained.
In the considered problem, we aim to minimize the function f W L S ( E ¯ , F ¯ ) . Thus, if the corresponding majorizing function m a j ( E ¯ , F ¯ | E ¯ c , F ¯ c ) is available, the iterative majorization approach can be used. As shown in [29], a majorizing function of f W L S ( E ¯ , F ¯ ) is given by
m a j ( E ¯ , F ¯ | E ¯ c , F ¯ c ) = w max 2 X ¯ ˜ E ¯ F ¯ T F 2 + constant ,
where w max is the maximum element of the weight matrix W ; in this paper, we have w max = 1 , as shown in (37), and X ¯ ˜ is given by
X ¯ ˜ = E ¯ c F ¯ c T + w max 2 W W X ¯ E ¯ c F ¯ c T .
It can be seen from (41) that the minimization of m a j ( E ¯ , F ¯ | E ¯ c , F ¯ c ) is equivalent to minimizing X ¯ ˜ E ¯ F ¯ T F 2 . As a result, given the current estimate ( E ¯ c , F ¯ c ) , in each iteration, the updated estimate ( E ¯ u , F ¯ u ) can be obtained from the following problem:
min X ¯ ˜ E ¯ F ¯ T F 2 s . t . E ¯ T E ¯ = I 2 ( M + L ) ,
which is similar to the original one (21), except that the data matrix X ¯ is replaced by X ¯ ˜ . Again, the iterative algorithm can be employed to solve the above problem. In brief, given a certain estimate F ¯ , The solution of (43) is given by
E ¯ = M ˜ N ˜ T ,
where M ˜ and N ˜ are obtained from the economy SVD of the matrix X ¯ ˜ F ¯ as
X ¯ ˜ F ¯ = S V D M ˜ D ˜ N ˜ T ,
where D ˜ is a diagonal matrix composed of the singular values.
Similarly to the EFA algorithm, the above iterative procedure stops when the difference between two consecutive f W L S ’s is negligible or the maximal iterations are reached. Therefore, A o can be estimated from F ¯ according to (20a), (20b) and (26), and the DOA can be estimated from the MUSIC spectrum as (27). Finally, the proposed robust EFA algorithm for DOA estimation against outliers is summarized in Algorithm 2.

4.4. Complexities of the Proposed Methods

In the EFA-based approach for subspace estimation, the main complexity comes from the SVD of the matrix X ¯ F ¯ with dimensions 2 N × ( 2 M + 2 L ) in each iteration. Since N is commonly larger than M + L , the complexity is O ( N ( M + L ) 2 ) . Similarly, the complexity of the robust EFA-based approach is O ( N ( M + L ) 2 ) , due to the SVD of the matrix X ¯ ˜ F ¯ with dimensions 2 N × ( 2 M + 2 L ) in (45). The complexity is comparable to the conventional ILSSE algorithm with complexity O ( M 3 ) , Moreover, it is computationally more efficient than the nonuniform ML methods, which require a multidimensional parameter search. Note that to determine the DOAs based on the subspace, a spectrum search is required if the MUSIC algorithm is used.
Algorithm 2 Proposed Robust EFA-based DOA Estimation Against Outliers
1:
Use the generalized ESD test for outlier detection.
2:
Initialization:  E ¯ c , F ¯ c , f = X ¯ E ¯ c F ¯ c T F 2 , f W L S = ( X ¯ E ¯ c F ¯ c T ) W F 2
3:
while  | f f W L S | > ε or max iteration number unreached do
4:
   Calculate X ¯ ˜ = E ¯ c F ¯ c T + w max 2 W W X ¯ E ¯ c F ¯ c T
5:
   Solve the problem min X ¯ ˜ E ¯ F ¯ T F 2 s . t . E ¯ T E ¯ = I 2 ( M + L )
6:
   Obtain the solution E ¯ u = [ S ¯ o u U ¯ o u ]
7:
   Update A ¯ o u = X ¯ ˜ T S ¯ o u , Φ ¯ o u = d i a g ( X ¯ ˜ T U ¯ o u )
8:
   Update F ¯ u = [ A ¯ o u Φ ¯ o u ]
9:
    F ¯ c = F ¯ u and E ¯ c = E ¯ u
10:
    f = f W L S , f W L S = ( X ¯ E ¯ c F ¯ c T ) W F 2
11:
end while
12:
A o = A ¯ o ( 1 : M , 1 : L ) + j A ¯ o ( 1 : M , L + 1 : 2 L )
13:
Orthonormalize A o to obtain A ˜ o
14:
Spectrum G ( θ , ϕ ) = 1 | a H ( θ , ϕ ) ( I M A ˜ o A ˜ o H ) a ( θ , ϕ ) |
15:
Output: DOA estimates.

5. Simulation Results

To evaluate the performance of the proposed method, in this section, we consider the problem of 2D DOA estimation using a 4 × 3 rectangular planar array, and hence, M = 12 . The sensors are spaced by half-wavelength. For illustrative purposes, we assume that the source signal impinges on θ [ 90 , 90 ] , φ [ 0 , 90 ] . More precisely, three uncorrelated narrowband signals with identical power impinge on the array from the far field. The DOAs of them are assumed to be ( 25 , 65 ) , ( 50 , 65 ) , and ( 30 , 30 ) , respectively. Unless otherwise specified, the covariance matrix of the nonuniform noise is assumed to be
Q = d i a g { 0 . 5 , 30 , 2 , 2 . 5 , 1 , 8 , 1 . 5 , 10 , 0 . 8 , 1 , 0 . 5 , 1 } .
Accordingly, the signal-to-noise ratio (SNR) is defined as follows [9]: SNR = σ s 2 M m = 1 M 1 σ m 2 , where σ s 2 denotes the signal power. For the special case of uniform white noise, we set σ 1 2 = σ 2 2 = σ M 2 = σ n o i s e 2 , so that SNR = σ s 2 σ s 2 σ n o i s e 2 σ n o i s e 2 , which is the traditional definition.
Except for the last example in Section 5.2, the number of snapshots is 100. Moreover, if outliers are considered, we assume that they occur at the ( 55 , 1 ) th , ( 75 , 1 ) th , ( 60 , 10 ) th , and ( 100 , 8 ) th entries of the data matrix X . In the simulations, random values produced from normal distribution with large variance are added to those entries to simulate outliers.

5.1. Uniform White Noise without Outliers

As a special case of nonuniform white noise, uniform white noise is considered in this example. Firstly, we assume that no outliers exist. The SNR is 10 dB. Figure 1 shows the resultant spectra and their contour plots of the MUSIC- and EFA-based approaches. It should be mentioned that in all simulations, the MUSIC algorithm is carried out by obtaining the subspace through the eigenvalue decomposition of the covariance matrix estimate, R ^ = N 1 X H X . From Figure 1, we notice that in such an ideal case, i.e., in uniform white noise and in the absence of data outliers, the two methods perform very well; the elevation as well as azimuth angles of the sources can be correctly estimated.

5.2. Uniform White Noise with Outliers

In this subsection, we assume that the data are corrupted by outliers. Figure 2a,b show the spectra and their contour plots of MUSIC and EFA in such a situation. It can be found that the data outliers can significantly deteriorate the performance of these nonrobust methods. Specifically, the spectra do not form proper peaks at the true DOAs. On the contrary, it can be seen from Figure 3a that the hostile effect of data outliers can be effectively suppressed by the proposed robust EFA method and hence, the DOAs can be still correctly estimated from the spectra. In the robust EFA method, the weight for the trimmed data point is set to γ = 0.8 . For comparison, the performance of the conventional robust PAST algorithm [18], where the forgetting factor is set to 0.95, is also tested and the corresponding result is shown in Figure 3b. We notice that the robust PAST algorithm also provide satisfactory performance in this case.
To further examine the performance of the robust EFA algorithm, the root-mean-square error (RMSE) of DOA estimation is compared with those of other methods (MUSIC and robust PAST) in Figure 4. The results are obtained from 100 independent runs. The performance of EFA without outliers is also shown as benchmark. It can be seen that the nonrobust MUSIC algorithm performs poorly in this scenario. Moreover, the proposed robust EFA algorithm generally outperforms the robust PAST algorithm. One reason is that the robust PAST algorithm does not update the subspace when the data are corrupted by outliers. In other words, the corrupted data are completely excluded for subspace estimation, and the number of effective snapshots in the robust PAST algorithm is less than that in the proposed method. In addition, it is found that the robust EFA algorithm is capable of achieving similar performance as the ideal case without outliers.
It is also worth noting the robust PAST algorithm requires several snapshots to converge [18]. Hence, compared with the robust EFA algorithm, it may be further affected when the sample number is relatively small, or outliers exist in the first few snapshots. For illustration, we reduce the number of snapshots to 50, and the outliers occur at the ( 11 , 1 ) th, ( 12 , 3 ) th, ( 20 , 8 ) th, and ( 30 , 6 ) th entries of the data matrix. The resulting spectra of the robust PAST and robust EFA algorithms are shown in Figure 5. It is seen that the robust PAST algorithm does not perform well in this case, whereas the robust EFA algorithm still offer good performance.

5.3. Nonuniform White Noise without Outliers

In this subsection, the noise is assumed to be nonuniform white. Firstly, we assume that there are no outliers in the data matrix X . The SNR is 10 dB according to (46). Figure 6 shows the resultant contour plots of the spectra of various methods including EFA, robust EFA, ILSSE [9], reduced covariance matrix (RCM) [11], robust PAST [18], and MUSIC. Obviously, we can notice from Figure 6e,f that two peaks of the MUSIC spectrum and robust PAST spectrum are merged together. This is because in the nonuniform white noise environment, the subspace cannot be well estimated by these approaches. On the contrary, the EFA-based approaches (EFA and robust EFA), ILSSE, and RCM take the noise nonuniformity into account, so that the subspaces and DOAs can be properly determined.
In Figure 7, the RMSEs of DOA estimation of different methods are compared. Since outliers are not considered in this example, the robust EFA and robust PAST are excluded, and we focus on comparing the performance of the proposed EFA-based method with the RCM and ILSSE algorithms. It is observed that when only nonuniform noise exists, the DOA estimation accuracy of these methods is similar, and the performance gap among these methods is small.
In order to examine the influence of the degree of the noise nonuniformity on DOA estimation, we change the maximum noise power, i.e., σ 2 2 in (46), from 25 to 75, and keep σ min 2 = 0.5 . Therefore, the worst noise power ratio (WNPR), σ max 2 / σ min 2 is varied from 50 to 150. For each WNPR, the SNR remains to be 10dB. It is seen from Figure 8 that the performance of all the tested methods degrades with the increase of WNPR. In general, the proposed EFA-based method achieves better performance than the other two methods.

5.4. Nonuniform White Noise with Outliers

Following the previous setup, both noise nonuniformity and data outliers are considered. The resultant contour plots of the spectra of the proposed robust EFA algorithm and various conventional robust and nonrobust methods are shown in Figure 9. Obviously, it can be seen that the proposed robust EFA method can clearly identify the three source signals. As expected, the nonrobust methods, e.g., MUSIC, EFA, EFA, ILSSE, and RCM, are greatly affected by the outliers and they cannot provide satisfactory performance in this case. Furthermore, though the conventional robust PAST algorithm performs well in uniform noise with outliers, it does not offer satisfactory performance in nonuniform white noise. On the contrary, the issues of nonuniformity and data outliers can be simultaneously handled by the robust EFA method.
To further show the ability of the robust EFA against outliers, the RMSEs of azimuth and elevation estimation using this method are compared to various methods, including ILSSE, RCM, robust PAST, and the benchmark method (EFA without outliers) in Figure 10. We notice that the hostile effects of outliers can be successfully suppressed in all SNR levels tested, and the performance is quite closed to the case without outliers. Moreover, although the robust PAST algorithm is robust against outliers, it cannot suppress the nonuniform noise effectively, and hence, it is outperformed by the proposed method. For the ILSSE and RCM algorithms, although they can perform well in nonuniform noise, they cannot provide robustness against outliers and perform even worse than the robust PAST algorithm.

6. Conclusions

In this paper, we investigate the problem of DOA estimation in the presence of nonuniform noise as well as data outliers. Differently to many traditional DOA estimation algorithms, which use a uniform white noise mode, it is assumed that the noise is nonuniform white noise, and the noise covariance matrix is diagonal but the diagonal elements are not identical. In order to estimate the subspace, and hence, DOAs, in this case, a new algorithm based on the EDA model is developed. Next, the data outliers are also considered. A simple method, i.e., the ESD test, is employed for outlier detection. Based on this detection, a trimmed data matrix can be obtained. According to the trimmed data matrix, a WLS problem is formulated to determine the subspace. Moreover, an iterative majorization approach, which is monotonic convergent, is introduced to solve the WLS problem. Simulation results show that the proposed robust DOA estimator outperforms traditional algorithms in the presence of nonuniform noise and outliers.

Author Contributions

Conceptualization, B.G. and B.L.; methodology, B.G. and B.L.; software, B.L. and X.S.; validation, B.G., X.S., Z.L. and B.L.; formal analysis, B.L., X.S. and B.L.; investigation, X.S., Z.L. and B.L.; resources, X.S. and Z.L.; writing—original draft preparation, B.G. and B.L.; writing—review and editing, B.G., X.S., Z.L. and B.L.;supervision, X.S. and Z.L.; project administration, X.S., Z.L. and B.L.; funding acquisition, X.S., Z.L. and B.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded in part by the National Natural Science Foundation of China under Grants 62171292 and 62101340, and in part by the Guangdong Basic and Applied Basic Research Foundation under Grant 2022A1515010188.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

Author Zhengqiang Li was employed by the company Shanghai Aircraft Design and Research Institute. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Appendix A

According to the definition of the real-valued matrix B ¯ in (14), it can be readily calculated that
B ¯ F 2 = 2 Re { B } F 2 + 2 Im { B } F 2 = 2 Re { B } + j Im { B } F 2 = 2 B F 2 .
Therefore, we have B F 2 = 1 2 B ¯ F 2 . This completes the proof of (15).
Since G = B C H , we have
Re { G } = Re { B C H } = Re { ( Re { B } + j Im { B } ) ( Re { C } T j Im { C } T ) } = Re { B } Re { C } T + Im { B } Im { C } T
and
Im { G } = Im { B C H } = Im { ( Re { B } + j Im { B } ) ( Re { C } T j Im { C } T ) } = Re { B } Im { C } T + Im { B } Re { C } T .
As a result, the real-valued matrix M ¯ is given by
G ¯ = Re { B C H } Im { B C H } Im { B C H } Re { B C H } = Re { B } Re { C } T + Im { B } Im { C } T Re { B } Im { C } T + Im { B } Re { C } T Re { B } Im { C } T Im { B } Re { C } T Re { B } Re { C } T + Im { B } Im { C } T = Re { B } Im { B } Im { B } Re { B } × Re { C } T Im { C } T Im { C } T Re { C } T = B ¯ C ¯ T .
This completes the proof of (16).

References

  1. Pesavento, M.; Gershman, A.B. Maximum-likelihood direction-of-arrival estimation in the presence of unknown nonuniform noise. IEEE Trans. Signal Process. 2001, 49, 1310–1324. [Google Scholar] [CrossRef]
  2. Böhme, J.F. Estimation of spectral parameters of correlated signals in wavefields. Signal Process. 1986, 11, 329–337. [Google Scholar] [CrossRef]
  3. Jaffer, A.G. Maximum likelihood direction finding of stochastic sources: A separable solution. In Proceedings of the ICASSP-88, International Conference on Acoustics, Speech, and Signal Processing, New York, NY, USA, 11–14 April 1988; pp. 2893–2894. [Google Scholar]
  4. Stoica, P.; Nehorai, A. On the concentrated stochastic likelihood function in array signal processing. Circuits Syst. Signal Process. 1995, 14, 669–674. [Google Scholar] [CrossRef]
  5. Stoica, P.; Viberg, M.; Wong, K.M.; Wu, Q. Maximum-likelihood bearing estimation with partly calibrated arrays in spatially correlated noise fields. IEEE Trans. Signal Process. 1996, 44, 888–899. [Google Scholar] [CrossRef]
  6. Chen, C.E.; Lorenzelli, F.; Hudson, R.E.; Yao, K. Stochastic maximum-likelihood DOA estimation in the presence of unknown nonuniform noise. IEEE Trans. Signal Process. 2008, 56, 3038–3044. [Google Scholar] [CrossRef]
  7. Wu, Y.; Hou, C.; Liao, G.; Guo, Q. Direction-of-arrival estimation in the presence of unknown nonuniform noise fields. IEEE J. Ocean. Eng. 2006, 31, 504–510. [Google Scholar] [CrossRef]
  8. Liao, B.; Liao, G.; Wen, J. A method for DOA estimation in the presence of unknown nonuniform noise. J. Electromagn. Waves Appl. 2008, 22, 2113–2123. [Google Scholar] [CrossRef]
  9. Liao, B.; Chan, S.C.; Huang, L.; Guo, C. Iterative methods for subspace and DOA estimation in nonuniform noise. IEEE Trans. Signal Process. 2016, 64, 3008–3020. [Google Scholar] [CrossRef]
  10. Liao, B.; Wen, J.; Huang, L.; Guo, C.; Chan, S.C. Direction finding with partly calibrated uniform linear arrays in nonuniform noise. IEEE Sens. J. 2016, 16, 4882–4890. [Google Scholar] [CrossRef]
  11. Liao, B.; Huang, L.; Guo, C.; So, H.C. New approaches to direction-of-arrival estimation with sensor arrays in unknown nonuniform noise. IEEE Sens. J. 2016, 16, 8982–8989. [Google Scholar] [CrossRef]
  12. Wen, J.; Liao, B.; Guo, C. Spatial smoothing based methods for direction-of-arrival estimation of coherent signals in nonuniform noise. Digit. Signal Process. 2017, 67, 116–122. [Google Scholar] [CrossRef]
  13. Liao, B.; Guo, C.; Huang, L.; Wen, J. Matrix completion based direction-of-arrival estimation in nonuniform noise. In Proceedings of the 2016 IEEE International Conference on Digital Signal Processing (DSP), Beijing, China, 16–18 October 2016; pp. 66–69. [Google Scholar]
  14. Fang, Y.; Zhu, S.; Zeng, C.; Gao, Y.; Li, S. DOA Estimations with Limited Snapshots Based on Improved Rank-One Correlation Model in Unknown Nonuniform Noise. IEEE Trans. Veh. Technol. 2021, 70, 10308–10319. [Google Scholar] [CrossRef]
  15. Xu, H.; Jin, M.; Guo, Q.; Jiang, T.; Tian, Y. Direction of Arrival Estimation with Gain-Phase Uncertainties in Unknown Nonuniform Noise. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 9686–9696. [Google Scholar] [CrossRef]
  16. Kozick, R.J.; Sadler, B.M. Maximum-likelihood array processing in non-Gaussian noise with Gaussian mixtures. IEEE Trans. Signal Process. 2000, 48, 3520–3535. [Google Scholar]
  17. Lim, C.H.; See, S.C.M.; Zoubir, A.M.; Ng, B.P. Robust adaptive trimming for high-resolution direction finding. IEEE Signal Process. Lett. 2009, 16, 580–583. [Google Scholar]
  18. Chan, S.C.; Wen, Y.; Ho, K.L. A robust past algorithm for subspace tracking in impulsive noise. IEEE Trans. Signal Process. 2005, 54, 105–116. [Google Scholar] [CrossRef]
  19. Huber, P.J. Robust Statistics; John Wiley & Sons: Hoboken, NJ, USA, 2004; Volume 523. [Google Scholar]
  20. Zou, Y.; Chan, S.; Ng, T. A recursive least M-estimate (RLM) adaptive filter for robust filtering in impulse noise. IEEE Signal Process. Lett. 2000, 7, 324–326. [Google Scholar] [CrossRef]
  21. Zou, Y.; Chan, S.C.; Ng, T.S. Fast least mean M-estimate algorithms for robust adaptive filtering in impulse noise. In Proceedings of the 2000 10th European Signal Processing Conference, Tampere, Finland, 4–8 September 2000; pp. 1–4. [Google Scholar]
  22. Yang, B. Projection approximation subspace tracking. IEEE Trans. Signal Process. 1995, 43, 95–107. [Google Scholar] [CrossRef]
  23. Mardia, K.; Kent, J.; Bibby, J. Multivariate Analysis (Probability and Mathematical Statistics); Academic Press: Cambridge, MA, USA, 1980. [Google Scholar]
  24. McDonald, R.P. The simultaneous estimation of factor loadings and scores. Br. J. Math. Stat. Psychol. 1979, 32, 212–228. [Google Scholar] [CrossRef]
  25. Knott, M.; Bartholomew, D.J. Latent Variable Models and Factor Analysis; Edward Arnold: New York, NY, USA, 1999. [Google Scholar]
  26. Trendafilov, N.T.; Unkel, S. Exploratory factor analysis of data matrices with more variables than observations. J. Comput. Graph. Stat. 2011, 20, 874–891. [Google Scholar] [CrossRef]
  27. Unkel, S.; Trendafilov, N.T. Simultaneous parameter estimation in exploratory factor analysis: An expository review. Int. Stat. Rev. 2010, 78, 363–382. [Google Scholar] [CrossRef]
  28. Rosner, B. Percentage points for a generalized ESD many-outlier procedure. Technometrics 1983, 25, 165–172. [Google Scholar] [CrossRef]
  29. Kiers, H.A. Weighted least squares fitting using ordinary least squares algorithms. Psychometrika 1997, 62, 251–266. [Google Scholar] [CrossRef]
  30. Gower, J.C.; Dijksterhuis, G.B. Procrustes Problems; Oxford University Press: Oxford, UK, 2004. [Google Scholar]
  31. Osborne, J.W.; Overbay, A. The power of outliers (and why researchers should always check for them). Pract. Assess. Res. Eval. 2019, 9, 6. [Google Scholar]
  32. Barnett, V.; Lewis, T. Outliers in Statistical Data; Wiley: New York, NY, USA, 1994. [Google Scholar]
  33. Hawkins, D.M. Identification of Outliers; Springer: Berlin/Heidelberg, Germany, 1980; Volume 11. [Google Scholar]
  34. Hodge, V.; Austin, J. A survey of outlier detection methodologies. Artif. Intell. Rev. 2004, 22, 85–126. [Google Scholar] [CrossRef]
Figure 1. Spectra and their contour plots of different methods in uniform white noise without outliers. Crosses denote the true DOAs. (a) left column: MUSIC; (b) right column: EFA.
Figure 1. Spectra and their contour plots of different methods in uniform white noise without outliers. Crosses denote the true DOAs. (a) left column: MUSIC; (b) right column: EFA.
Remotesensing 16 03140 g001
Figure 2. Spectra and their contour plots of different methods in uniform white noise with outliers. Crosses denote the true DOAs. (a) left column: MUSIC; (b) right column: EFA.
Figure 2. Spectra and their contour plots of different methods in uniform white noise with outliers. Crosses denote the true DOAs. (a) left column: MUSIC; (b) right column: EFA.
Remotesensing 16 03140 g002
Figure 3. Spectra and their contour plots of different methods in uniform white noise with outliers. Crosses denote the true DOAs. (a) left column: Robust EFA; (b) right column: Robust PAST.
Figure 3. Spectra and their contour plots of different methods in uniform white noise with outliers. Crosses denote the true DOAs. (a) left column: Robust EFA; (b) right column: Robust PAST.
Remotesensing 16 03140 g003
Figure 4. RMSE of angle estimation in the uniform white noise with outliers (a) Azimuth angle; (b) Elevation angle.
Figure 4. RMSE of angle estimation in the uniform white noise with outliers (a) Azimuth angle; (b) Elevation angle.
Remotesensing 16 03140 g004
Figure 5. Comparison of the spectra of robust EFA and robust PAST algorithms in the uniform white noise with outliers. The number of snapshots is 50.
Figure 5. Comparison of the spectra of robust EFA and robust PAST algorithms in the uniform white noise with outliers. The number of snapshots is 50.
Remotesensing 16 03140 g005
Figure 6. Contour plots of the spectra of various methods in nonuniform white noise without outliers.
Figure 6. Contour plots of the spectra of various methods in nonuniform white noise without outliers.
Remotesensing 16 03140 g006
Figure 7. Comparison of the RMSEs versus SNR in the nonuniform white noise without outliers.
Figure 7. Comparison of the RMSEs versus SNR in the nonuniform white noise without outliers.
Remotesensing 16 03140 g007
Figure 8. Comparison of the RMSEs versus WNPR in the nonuniform white noise without outliers.
Figure 8. Comparison of the RMSEs versus WNPR in the nonuniform white noise without outliers.
Remotesensing 16 03140 g008
Figure 9. Contour plots of the spectra of various methods in nonuniform white noise with outliers.
Figure 9. Contour plots of the spectra of various methods in nonuniform white noise with outliers.
Remotesensing 16 03140 g009
Figure 10. Comparison of the RMSEs of angle estimation versus SNR in nonuniform white noise with outliers: (a) azimuth angle; (b) elevation angle.
Figure 10. Comparison of the RMSEs of angle estimation versus SNR in nonuniform white noise with outliers: (a) azimuth angle; (b) elevation angle.
Remotesensing 16 03140 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, B.; Shen, X.; Li, Z.; Liao, B. Robust Direction-of-Arrival Estimation in the Presence of Outliers and Noise Nonuniformity. Remote Sens. 2024, 16, 3140. https://doi.org/10.3390/rs16173140

AMA Style

Gao B, Shen X, Li Z, Liao B. Robust Direction-of-Arrival Estimation in the Presence of Outliers and Noise Nonuniformity. Remote Sensing. 2024; 16(17):3140. https://doi.org/10.3390/rs16173140

Chicago/Turabian Style

Gao, Bin, Xing Shen, Zhengqiang Li, and Bin Liao. 2024. "Robust Direction-of-Arrival Estimation in the Presence of Outliers and Noise Nonuniformity" Remote Sensing 16, no. 17: 3140. https://doi.org/10.3390/rs16173140

APA Style

Gao, B., Shen, X., Li, Z., & Liao, B. (2024). Robust Direction-of-Arrival Estimation in the Presence of Outliers and Noise Nonuniformity. Remote Sensing, 16(17), 3140. https://doi.org/10.3390/rs16173140

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

Article Metrics

Back to TopTop