Abstract

In this article, a difference-coarray-based direction of arrival (DOA) method is introduced, which utilizes the uniform linear array (ULA) in a novel fashion to address the problem of DOA estimation for coherent signals. Inspired by the coarray-based estimators employed in cases of sparse arrays, we convert the sample covariance matrix of the observed signals into the difference coarray domain and process the signals using a spatial smoothing technique. The proposed method exhibits good accuracy and robustness in both the uncorrelated and coherent cases. Numerical simulations verify that the ULA difference coarray- (UDC-) based method can achieve good DOA estimation accuracy even when the SNR is very low. In addition, the UDC-based method is insensitive to the number of snapshots. Under extremely challenging conditions, the proposed UDC-ES-DOA method is preferred because of its outstanding robustness, while the UDC-MUSIC method is suitable for most moderate cases of lower complexity. Due to its demonstrated advantages, the proposed method is a promising and competitive solution for practical DOA estimation, especially for low-SNR or snapshot-limited applications.

1. Introduction

Antenna-array-based direction of arrival (DOA) estimation algorithms have been widely used in radar, navigation, measurement, and control systems. Various theories and techniques have been developed for array signal processing related to DOA estimation [1]. Uniform linear arrays (ULAs) are often considered in both theoretical research and engineering applications due to their symmetrical and simple array geometries. In this context, subspace-based methods are widely used, including multiple signal classification (MUSIC) [2] and the estimation of signal parameters via rotation invariance techniques (ESPRIT) [3]. In some practical scenarios, such as indoor environments, the signals may be highly correlated (coherent) due to multipath propagation. In this case, the signal covariance matrix is no longer nonsingular, and as a result, the performance of the abovementioned subspace-based DOA algorithms will deteriorate dramatically.

To address this problem, the most common method is to pretreat the observed signals using a spatial smoothing (SS) scheme to recover the reduced rank of the covariance matrix [4]. To cope with K coherent source signals received by a ULA, the forward spatial smoothing (FSS) method requires at least 2K sensors, while the forward-backward spatial smoothing (FBSS) method requires 3K/2 sensors [5]. Compared to the uncorrelated case, SS approaches realize the decorrelation of coherent signals at the cost of the number of DOFs, and these approaches often lead to worse DOA estimation due to the reduction of the array aperture [6]. Although SS-based approaches are simple and easy to implement, their accuracy and robustness are not satisfactory since the estimation results are sensitive to the SNR and the number of snapshots.

In [7], the Improved MUSIC algorithm was introduced, which enables accurate DOA estimation by reconstructing the covariance matrix of the received signals. In [8], another algorithm, named ES-DOA, was proposed based on an investigation of Improved MUSIC. The ES-DOA algorithm can achieve high and robust performance by making use of both the signal subspace and noise subspace characteristics. The main drawback of both Improved MUSIC and ES-DOA algorithms is that they are valid only if there is just one coherent signal. An ESPRIT-like algorithm is proposed to solve the coherent problem by eliminating the rank loss of the spatial covariance matrix with a reconstructed Toeplitz matrix [911]. In [12], an MSE-based regularization method is proposed to detect the number of coherent signals and estimate the DOAs. Based on the KR subspace approach, Masaaki et al. proposed the Kronecker-MUSIC method [13], which can cope with N − 1 coherent source signals by means of N sensors. However, the source signals are restricted to being quasistationary, thus limiting the widespread application of this method. In addition, there are some methods that are developed based on the spatial differencing technique to address the DOA estimation problem of the mixed uncorrelated and coherent sources in multipath environment [1416]. Compressive sensing- (CS-) based approaches, such as the l1-SVD algorithm [17], the JLZA-DOA algorithm [18], and the SBL-DOA algorithm [19], were introduced to identify the DOAs. The CS-based methods are very attractive in that they can cope with the coherent signals with different topologies of the antenna arrays and identify the number of sources automatically [2022]. Nevertheless, there are basis mismatch problem and heavy computational complexity drawbacks for the CS-based DOA methods.

Recently, approaches based on sparse arrays, such as coprime arrays [23] and nested arrays [24], have been suggested. It is possible to achieve O (N2) DOFs using O (N) sensors based on a nested array, which is formed by combining two or more ULAs with successively increasing spacing. A number of follow-up methods have been developed to further increase the number of DOFs by modifying the array configurations [2528]. The difference coarray concept is the foundation of these sparse-array-based approaches [29]. However, for the existing sparse array estimators, SS schemes cannot be utilized to cope with coherent source signals [30].

Inspired by the coarray-based estimators employed in cases of sparse arrays, we convert the sample covariance matrix of the signals observed by a ULA into the difference coarray domain and pretreat the signals using an SS method. Since the difference coarray of a ULA is another virtual ULA, we can obtain almost twice the original number of DOFs. Furthermore, we utilize the invariance of the ULA shift to handle coherent source signals because the original and virtual arrays both have ULA structures.

The major contributions of this paper include the following:A novel ULA difference-coarray-based strategy is presented to address the general coherent signal DOA estimate issue. This approach can achieve good accuracy and robustness of DOA estimation in both the uncorrelated and coherent cases.An enhanced method named UDC-ES-DOA is discussed to solve the DOA estimate problem under extremely challenging conditions. This method can achieve high accuracy with limited snapshots and demonstrates significant superiority in low-SNR scenarios. It is favorable for practical applications.

In addition, numerical simulations are conducted to evaluate the validity of the proposed method. The results show the superiority of the proposed method to the conventional approaches; meanwhile, it can achieve similar performance as the CS-based methods with significantly less computational complexity.

The remainder of the paper is organized as follows. Section 2 introduces the signal model and the conventional methods used in this research. The proposed method is introduced in Section 3. The details of the numerical simulations are discussed in Section 4. Concluding remarks are given in Section 5.

Notations: vectors and matrices are denoted by and lower-case letters (k) and boldface capital letters (k), respectively. Diag{•} is a diagonal matrix operator, and vec (•) denotes an operator that returns the vectorization result of a matrix. (_)T, (_)H, and (_) stand for transpose, conjugate transpose, and complex conjugate, respectively. ⊙ depicts the KR product, and ⊗ denotes the Kronecker product.

2. Preliminaries

To facilitate the presentation of the proposed method, we shall briefly present some preliminaries regarding the signal model, the FBSS method, difference-coarray-based DOA estimation, and the ES-DOA algorithm.

2.1. Signal Model

Suppose that there are K (K < N) narrowband plane waves impinging on an N-element ULA with the source matrix . Hence, the observed signal is modeled aswhere is the array manifold matrix with Vandermonde-structured distinct columns (N > K) and is of rank K, whereas n(t)∈ℂN represents a white Gaussian noise vector and is not correlated with the sources. The steering vector of the array can be expressed aswhere θi ∈ (−π/2, π/2) denotes the DOA of source i, d is the interval of the array, and λ is the signal wavelength. If the signal and noise components are stationary, zero-mean, uncorrelated random processes, then the covariance matrix of the received signals can be defined aswhere P is the source autocorrelation matrix of S(t), denotes the noise power, and I denotes the identity matrix. The source covariance matrix R is a full-rank matrix as long as the incident signals are uncorrelated. The real value of the sample covariance matrix R is not able to be obtained in practice; instead, the maximum-likelihood estimate of the sample covariance matrix is obtained from M snapshots:where M is the number of snapshots and i is an index representing different snapshots.

2.2. Forward-Backward Spatial Smoothing

As mentioned above, if the received signals are coherent, then conventional subspace-based algorithms cannot yield accurate DOA estimates because the source autocorrelation matrix P is not of full rank. By utilizing the invariance of the ULA shift, a spatial smoothing scheme can be used to address this problem.

The basic principle of an SS algorithm is to divide the original ULA into L overlapping uniform subarrays and introduce phase shifts between these subarrays. There are N − L + 1 sensors in each subarray. The FBSS covariance matrix is given bywhere J denotes the exchange matrix, with ones on its antidiagonal and zeros elsewhere.

It has been proven that when the number of sensors in the ULA is more than 3K/2, the DOA estimates for all coherent signals are available by using the FBSS scheme [5]. However, the SS process significantly reduces the effective array aperture, which often results in a deterioration of the resolution performance.

2.3. Difference-Coarray-Based DOA Method

For sparse linear arrays, such as nested arrays and coprime arrays, the number of DOFs can be significantly increased by exploiting the coarray structure. The difference set of a coarray is based on the KR subspace of the physical array manifold. The KR product has the following important property [31].

Property 1. Let A∈ℂn×k, B∈ℂn×k, and d∈ℂk. In addition, let D = Diag (d). Then,Under the assumption that all the incident signals are uncorrelated, from (3) and (7), we can obtain the vectorized as follows:where and , with ei being a column vector consisting of 1 in the i-th position and zeros elsewhere. The new array manifold has dimensions of N2 × D.
A coarray with a filled ULA structure that is centered at the origin can be obtained after the sorting and replacement of repeated rows. For a ULA, AA is similar to a Vandermonde matrix with (2N − 1) distinct rows. The observation vector of the virtual ULA is given bywhere W is the coarray selection matrix and denotes the (2 + 1) × K steering matrix of the (2 + 1)-element virtual ULA.
Since the difference coarray is a filled ULA, intuitively, it can be applied for DOA estimation. However, the equivalent source signal vector p consists of the powers of the actual source signals and behaves like the coherent sources in (6).
SS processing can be employed to solve the rank reconstruction problem, as demonstrated in [24]. As in the SS algorithm, the virtual array is divided into overlapping subarrays, each containing sensors. The i-th subarray has sensors located at and is represented bywhere A1i is an  × K matrix consisting of the ( − i + 1)th to (2 − i)th rows of and ei is a vector consisting entirely of zeros except for 1 in the i-th position. From (5) and (10), we can obtain the smoothed covariance matrix as follows:where .
As seen from (11), Rss can be treated as the covariance matrix of an N-element ULA whose array manifold is represented by A11. Thus, many subspace-based DOA estimation methods, such as MUSIC, can be applied.

2.4. ES-DOA Algorithm

Compared to the classical MUSIC algorithm, the ES-DOA algorithm not only utilizes the noise subspace characteristics for DOA estimation but also makes use of the signal subspace characteristics.

The source covariance matrix is reconstructed into a Toeplitz matrix as follows:where J denotes the N × N exchange matrix with ones on its antidiagonal and zeros elsewhere. The reconstructed matrix Rx can be decomposed aswhere the subscripts S and U represent the signal subspace part and the noise subspace part, respectively.

The generalized inverse matrix of the signal subspace part can be calculated as

The spatial frequency spectrum function of ES-DOA is defined as

The ES-DOA algorithm can achieve high resolution and demonstrates a strong antinoise capability when dealing with coherent array signals.

3. Proposed UDC-DOA Method

In this section, we will introduce the UDC-DOA method for addressing coherent signals received by a ULA.

Suppose that there are K coherent signals impinging on an N-element ULA at different azimuths {θ1, θ2,…, θK}. Since coherent signals can be regarded as phase-delayed amplitude-weighted replicas of each other, the observed signal can be expressed aswhere represents the complex attenuation vector of the coherent source signals with respect to the first signal s1(t). Suppose that ; then, using (3)–(5), the covariance matrix reduces to

It is clear that the signal covariance matrix Rco is no longer nonsingular in this case. From (4) and (16), we can rewrite the covariance matrix as

Inspired by the coarray-based estimators employed in cases of sparse arrays, we attempt to convert the sample covariance matrix of the signals observed by the ULA into the difference coarray domain. Since the product of is no longer diagonal in general in this case, we cannot vectorize Rco by using Property 1 directly.

For convenience in obtaining the equivalent observed signal in the difference coarray mode for the ULA, we introduce the difference set D as follows:

The difference value d ranges from 1 − N to N − 1, and the weight function of the difference value can be expressed as

The difference coarray of the ULA can be constructed via the following steps:(1)The covariance matrix Rco is rotated clockwise by 90°. As shown in Figure 1, each column corresponds to a distinct difference value, ranging from (N − 1) to (1 − N) when moving from left to right. Therefore, the elements in each column share the same difference value.(2)By calculating the mean value of each column, the new observation value can be obtained as follows:(3)It is clear that the resulting difference coarray is equal to another longer, virtual ULA with (2N − 1) elements. Following the process described above, the observation vector z of the virtual ULA can be expressed aswhere the matrix represents the new array manifold of the virtual ULA, q is the equivalent source signal vector, and represents the transition function of the matrix Rco that is mentioned above.

Then, we can obtain the covariance matrix of the new observed signals as follows:

Since the SS method can be used as a technique for building up the rank of an observation matrix for any array whose difference coarray is a filled ULA [24], this method can also naturally be applied to a ULA.

By applying the SS method, the virtual array is divided into N overlapping subarrays, each containing N sensors. According to (3) and (22), the covariance matrix of the i-th subarray can be expressed aswhere B1 is the submatrix consisting of the first to N-th rows of B. ϕi−1 denotes the i-th power of the K×K diagonal matrix

From (5), (6), and (23), the spatially smoothed covariance matrix is defined as the mean of the forward-backward subarray covariance matriceswhere represents the selection matrix for the i-th subarray.

Now, the DOA spectrum can be estimated using a general subspace-based algorithm such as MUSIC. In this paper, we prefer to use the ES-DOA algorithm to calculate the DOA spectrum because of its antinoise capability and robustness.

Notably, the original purpose of using SS method is to be able to exploit the DOFs of the difference coarray instead of decorrelating the coherent sources. In addition, the sparse arrays, such as nested arrays and coprime arrays, can commonly handle only uncorrelated signals. However, it is exciting and interesting that this restriction to uncorrelated signals no longer exists when the difference coarray is converted from a ULA.

As mentioned in Section 2, the ES-DOA algorithm can increase the rank of the source covariance matrix only by one in the coherent case, and it cannot be applied in the case of a sparse array. However, by virtue of the characteristics of ULAs, the ES-DOA algorithm can enhance the resolution and robustness performance in cases involving UDC signals.

Thus, our proposed method demonstrates improved accuracy and robustness in both the uncorrelated and coherent cases.

4. Numerical Examples

In this section, we employ numerical simulations to evaluate the validity of the proposed UDC-DOA method and compare its performance with that of conventional approaches. We assume that there are 11 far-field signals impinging on a ULA composed of 23 sensors. The locations of the signals range from −50° to +50°, in increments of 10°. The spatial noise is zero-mean uniform white complex Gaussian noise. We simulate the DOA spectra under different SNR conditions and with different numbers of snapshots in both the uncorrelated case and coherent cases.

4.1. Uncorrelated Case

First, we evaluate the performance of the proposed method in the uncorrelated case. Suppose that all the incident signals are completely uncorrelated. The corresponding DOA spectra estimated using the classic MUSIC algorithm, the FBSS-MUSIC algorithm, and the UDC-based MUSIC (UDC-MUSIC) algorithm are compared under different SNR conditions in Figure 2 and under different snapshot conditions in Figure 3.

4.1.1. Different SNR Conditions

As seen in Figure 2(a), with 200 signal snapshots, the spectra peaks obtained with all three algorithms are in good agreement with the true DOAs under good noise conditions, that is, with an SNR of 10 dB. As the SNR deteriorates to 0 dB, small but obvious peak deviations occur in the spectrum of the FBSS-MUSIC method, as shown in Figure 2(b). In this case, the spectral peaks obtained using the proposed method become sharper than those of the classic MUSIC algorithm. As shown in Figure 2(c), as the SNR continues to decline to −10 dB, the spectrum of the FBSS-MUSIC algorithm significantly deteriorates, with the peaks deviating considerably from the true DOAs. Meanwhile, for the classic MUSIC algorithm, the peaks become blurred and inconspicuous while the spectral peaks found with the proposed algorithm remain accurate and sharp.

The results showed in Figure 2 indicate that the proposed method can achieve better resolution and robustness compared to the traditional algorithm, especially under severe noise interference conditions.

4.1.2. Different Number of Signal Snapshots

In Figure 3, the SNR condition is set to 5 dB, and the DOA spectra in (a), (b), and (c) are calculated with different numbers of signal snapshots: 100, 20, and 5, respectively. By comparing the results, we find that the proposed method and FBSS-MUSIC both demonstrate similar resolution performance that is superior to that of the classic MUSIC algorithm when the number of signal snapshots is sufficient. However, the FBSS-MUSIC algorithm is sensitive to the number of snapshots, causing the peaks to severely deviate from the true values as the number of snapshots decreases. Moreover, although there are no obvious false peaks in the DOA spectrum of the classic MUSIC algorithm, the spectrum shows significant deterioration in resolution and even loses some peaks when the number of snapshots is very limited. In contrast, the proposed method shows reliable performance even with only a few signal snapshots. As shown in Figure 3(c), even with only five snapshots, the DOA peaks yielded by the proposed method remain accurate and sharp.

The results shown in Figure 3 indicate that the proposed method can maintain a high resolution and robustness even when there are only a few signal snapshots, while the traditional methods cannot function properly under these conditions.

4.2. Coherent Case

Second, we evaluate the performance of the proposed method in the coherent case. Suppose that all the incident signals are coherent, while the other simulation conditions are similar to the uncorrelated case. Since the classic MUSIC algorithm does not work in this case, we present only the performance comparison between the FBSS-MUSIC algorithm and the UDC-MUSIC algorithm.

4.2.1. Different SNR Conditions

As shown in Figure 4, the FBSS-MUSIC algorithm achieves better performance than the proposed algorithm under good signal conditions, such as SNR = 10 dB. However, the performance of the FBSS-MUSIC algorithm significantly deteriorates with the degradation of the SNR. When the SNR reaches −10 dB, the peaks in the spectrum obtained with the FBSS-MUSIC algorithm severely deviate from the true DOAs, while the results of the proposed UDC-MUSIC algorithm are still acceptably accurate.

4.2.2. Different Signal Snapshots

As shown in Figure 5(a), with 100 signal snapshots, both the UDC-MUSIC method and the FBSS-MUSIC method show good accuracy performance, and the peaks obtained with the latter are even sharper. However, similar to the uncorrelated case, the DOA peaks in the FBSS-MUSIC spectrum rapidly deviate from the true values as the number of snapshots decreases, while the proposed method shows a superior capability to cope with very few snapshots.

4.3. UDC-Based ES-DOA Algorithm

In the previous simulations, we have used the classic MUSIC algorithm to address the signals from the UDC. In this section, we evaluate the DOA estimation performance of the UDC-based ES-DOA (UDC-ES-DOA) method in both the uncorrelated and coherent cases.

As mentioned in Section 2, the ES-DOA algorithm can achieve a high resolution and has a strong antinoise capability. Thus, the following simulations correspond to challenging conditions, with the SNR set to 10 dB and only 50 snapshots.

Figure 6 shows the results for the uncorrelated case, while Figure 7 shows the results for the coherent case. In both figures, the UDC-ES-DOA method is denoted by UDC-ES. It is clear that the DOA spectrum obtained with UDC-ES-DOA has the highest resolution in both cases. Although the UDC-MUSIC algorithm shows much better performance than the classic MUSIC and FBSS-MUSIC algorithms, we can still observe some blurred and inconspicuous peaks in its spectra. Meanwhile, the UDC-ES-DOA method identifies the DOAs of all signals with little deviation.

4.4. RMSE Evaluation via Monte Carlo Simulations

To verify the advantages of our proposed method, without loss of generality, we report Monte Carlo simulations carried out to evaluate the average root-mean-square error (RMSE) performance for DOA estimation and compare the results with those of the other algorithms. The ULA topology and the incident signals remain the same as in the previous simulations.

We evaluate the RMSEs of the proposed method in both the uncorrelated and coherent cases. In both cases, the results are based on 200 independent Monte Carlo trials implemented to obtain the RMSE curves under different SNR and snapshot conditions. In the uncorrelated case, we compare the RMSE performances of the proposed UDC-ES-DOA and UDC-MUSIC algorithms, the FBSS-MUSIC algorithm, and the classic MUSIC algorithm. In the coherent case, we compare the RMSE performances of the proposed UDC-ES-DOA and UDC-MUSIC algorithms, the FBSS-MUSIC algorithm, and the l1-SVD algorithm. The l1-SVD algorithm is a kind of CS-based DOA approach which utilizes the l1 norm minimization techniques by assuming the incoming DOAs specified on the discrete frequency grid [5]. In our simulation, the grid ranges from −90° to 90° with the width of one degree.

Figure 8 depicts the RMSE curves with respect to the SNR when the number of snapshots is 200. From Figure 8(a), we find that the UDC-ES-DOA and UDC-MUSIC methods are obviously superior and that the UDC-ES-DOA method is significantly more robust when the SNR is below −5 dB. In contrast, the classic MUSIC algorithm is too sensitive to noise interference. The results for the coherent case are shown in Figure 8(b). Similar to the uncorrelated case, the UDC-ES-DOA and UDC-MUSIC methods perform much better than the FBSS method within the low-SNR range. Although the l1-SVD method is extremely reliable when the SNR is below −10 dB, it is surpassed by the proposed UDC-ES-DOA and UDC-MUSIC methods at the crossing points that occur at −7 dB and −5 dB, respectively.

Figure 9 depicts the RMSE curves with respect to the number of snapshots when the SNR is 5 dB. From Figure 9(a), we find that the UDC-ES-DOA and UDC-MUSIC methods demonstrate significant superiority when the snapshots are limited, and the UDC-ES-DOA method performs more impressively with only a few snapshots while the classic MUSIC algorithm cannot estimate correctly. It can be seen that the above three algorithms have similar accuracy performances when the number of snapshots is more than 30. In contrast, the RMSE curve of the FBSS-MUSIC algorithm indicates that it needs many more snapshots than the other three algorithms to achieve the same accuracy performance. The results for the coherent case are shown in Figure 9(b). Similar to the uncorrelated case, the UDC-ES-DOA and UDC-MUSIC methods perform much better than the FBSS method. Although the l1-SVD method is extremely reliable when the SNR is below −10 dB, it is surpassed by the proposed UDC-ES-DOA and UDC-MUSIC methods at the crossing points of −7 dB and −5 dB, respectively.

4.5. Computational Complexity

We evaluate the computational complexity of the proposed method and compare the results with those of the other methods. For a quantitative comparison, the CPU running time (RT) of different algorithms is computed with the same data in both the uncorrelated and coherent cases. The CPU type is Intel Core i7-4790 @ 3.6 GHz and we implemented the simulation in the MATLAB environment. We fix K = 11, N = 23, and SNR = 5 dB. And the simulation parameters and running time results are tabulated in Table 1.

It can be seen from the results that the computational complexity of the proposed method is slightly higher than the classic MUSIC and the FBSS method, while the CS-based algorithm takes several hundred times longer.

4.6. Remarks

From the simulation results above, it can be concluded that the proposed UDC-based methods are effective and impressive under both uncorrelated and coherent conditions. Compared to the traditional methods, the greatest advantage of the UDC-based methods is their outstanding robustness in cases of noise interference and a limited number of snapshots. Meanwhile, in the coherent source cases, the UDC-based methods can achieve similar DOA performance as the CS-based methods with significantly less computational complexity. According to the above evaluation and comparison based on Monte Carlo trials, it is fair to state that the UDC-based approach constitutes a promising and attractive method that is especially suitable for low-SNR or snapshot-limited practical applications.

5. Conclusions

In this paper, we have introduced a difference-coarray-based (UDC-based) DOA method to address the problem of DOA estimation for coherent signals received by a ULA. Inspired by the coarray-based estimators employed in cases of sparse arrays, we convert the sample covariance matrix of the signals observed by the ULA into the difference coarray domain and process the signals with an SS method. The proposed method exhibits good accuracy and robustness in both the uncorrelated and coherent cases. Numerical simulations verify that the UDC-based method can achieve significantly improved DOA estimation accuracy even when the SNR is very low. In addition, the UDC-based method is insensitive to the number of snapshots and thus can achieve high accuracy with a limited number of snapshots. Under extremely challenging conditions, the UDC-ES-DOA method is preferred because of its robustness, while the UDC-MUSIC method is suitable for most moderate cases of lower complexity. Due to its demonstrated advantages, the UDC-based method is a promising and competitive approach for practical DOA estimation, especially for low-SNR or snapshot-limited applications.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Science Foundation of China under Grant no. 61801497.