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
, where
is the noise variance and
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
where
denotes the steering matrix with
being the steering vector, and
and
are the signal waveform vector and additive noise measurement vector, respectively. More specifically, the steering vector can be written as
where
,
f is the carrier frequency, and
, with
c being the propagation velocity.
and
, respectively, are the unite vector and
sensor position:
The source signals are assumed to be temporally uncorrelated zero-mean Gaussian processes with
and
, where
denotes statistical expectation. Moreover, the covariance matrix
is diagonal, and its diagonal entries represent the signal powers. For the nonuniform white noise considered in this paper, we have
where
denotes the noise power in the
sensor. It can be seen that if
, 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:
where
is a diagonal matrix representing the standard deviation of each sensor noise, and
is a spatially and temporally uncorrelated standard normal distribution such that
and
, where
is an
identity matrix. According to (
1) and (
5), we have
Assuming that
N snapshots are collected, the observation data matrix can be written compactly as
where
,
,
, and
. Given a sufficiently large
N, we have
Making use of the above identities, Equation (
7) can be further reformulated as
where
,
,
, and
. Moreover, we have
It is assumed that the source signals and noise are uncorrelated, so that we have
, and hence
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 . Once such a subspace estimate is available, the conventional subspace-based DOA estimation algorithms can be applied. Moreover, from the identity , it is known that the spans the same space as .
3. EFA-Based DOA Estimation
In this section, we shall focus on estimating
by taking advantage of the EFA model. More precisely, according to (
9)–(
11), the EFA model can be described as
Consequently, based on the LS goodness-of-fit criterion, the corresponding problem for subspace (or say
) estimation is given by
In order to transform the complex-valued problem (
13) into a real-valued one, the following proposition will be applied:
Proposition 1. Let be a complex-valued matrix and the corresponding real-valued matrix be defined asthen, we have Further, given another complex-valued matrix , and letting , we havewhere and are real-valued matrices defined according to (14). Following the above proposition, the objective function in (
13) can be rewritten as
where
,
,
, and
are similarly defined as (
14). Moreover, the constraints in (
13) can be equivalently represented as
Consequently, the problem in (
13) can be equivalently transformed to the following real-value optimization problem:
Now, an iterative method will be introduced to solve the above problem. Firstly, we define two matrices
and
, respectively, as
Then, the problem in (
19) can be rewritten as
It can be noted that for a given
, the problem in (
21) is an orthogonal Procrustes problem, and its solution is given by [
30]
where
and
are obtained from the economy SVD of the matrix
as
where
is a diagonal matrix composed of the singular values.
According to (
20a), once
has been obtained,
and
can be respectively determined from the first
and the last
columns of
, i.e.,
,
. In addition, taking advantage of the properties of
and
as shown in (
18a), (
18b) and (
18c), one obtains
This indicates that
and
can be respectively updated as
Meanwhile, the matrix
can be further updated. It can be seen that
and
are estimated iteratively, and the procedure should be stopped after a certain convergence criterion is met, such as the difference between two consecutive
s (where
) being less than the prescribed threshold or the maximal iterations being reached. Finally,
can be extracted from
as
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:
where
is an orthonormal matrix of
. 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 , , , , , and , , - 2:
while or max iteration number unreached do - 3:
(economy SVD) - 4:
- 5:
- 6:
- 7:
- 8:
- 9:
- 10:
end while - 11:
- 12:
Orthonormalize to obtain - 13:
Spectrum - 14:
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
rectangular planar array, and hence,
. The sensors are spaced by half-wavelength. For illustrative purposes, we assume that the source signal impinges on
,
. 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
,
, and
, respectively. Unless otherwise specified, the covariance matrix of the nonuniform noise is assumed to be
Accordingly, the signal-to-noise ratio (SNR) is defined as follows [
9]:
where
denotes the signal power. For the special case of uniform white noise, we set
, so that
, 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
,
,
, and
entries of the data matrix
. 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,
. 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
. 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
th,
th,
th, and
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
. 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.,
in (
46), from 25 to 75, and keep
. Therefore, the worst noise power ratio (WNPR),
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.