Next Article in Journal
Hybrid Game Trading Mechanism for Virtual Power Plant Based on Main-Side Consortium Blockchains
Next Article in Special Issue
Enhancing Underwater Image Quality Assessment with Influential Perceptual Features
Previous Article in Journal
Explicit Representation of Mechanical Functions for Maintenance Decision Support
Previous Article in Special Issue
THANet: Transferring Human Pose Estimation to Animal Pose Estimation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Enhanced DOA Estimation Method for Coherent Sources via Toeplitz Matrix Reconstruction and Khatri–Rao Subspace

1
School of Information and Electronics, Beijing Institute of Technology, Beijing 100080, China
2
School of Computer Science, Xiamen University, Xiamen 361000, China
*
Author to whom correspondence should be addressed.
Electronics 2023, 12(20), 4268; https://doi.org/10.3390/electronics12204268
Submission received: 8 September 2023 / Revised: 7 October 2023 / Accepted: 10 October 2023 / Published: 16 October 2023
(This article belongs to the Special Issue Advances in Computer Vision and Multimedia Information Processing)

Abstract

:
The Toeplitz matrix reconstruction methods are capable of resolving coherent signals, playing a crucial role in the direction-of-arrival (DOA) estimation of acoustic sources. However, the decoherence processing sacrifices the array aperture and further results in a reduced resolution capability for the number of identifiable sources. To solve this issue, we propose an enhanced method using the Khatri–Rao subspace to resolve more coherent sources than that of the existing Toeplitz matrix reconstruction methods. Firstly, a full set of Toeplitz matrices with full rank is obtained. Then, the virtual array aperture can be obtained using the Khatri–Rao product of the array response, and the degree of freedom provided inherently in the virtual array structure is about twice the size of that of the existing Toeplitz methods. Next, linear processing is further used to achieve complexity reduction without losing the effective degree of freedom. Finally, the DOA estimation for more coherent sources can be achieved by combining it with conventional methods. Numerical simulations verify the superiority of the proposed method.

1. Introduction

The ocean is a vast and largely unexplored frontier, replete with an abundance of resources within its aquatic realm. This circumstance has engendered an escalating demand for the development of underwater sensing and communication systems, based upon advanced multimedia information processing technologies. Central to the functionality of these systems is the process of direction-of-arrival (DOA) estimation, which serves as the linchpin for source localization, beamforming, and target tracking [1,2,3,4,5]. In practice, coherent signal sources are usually encountered owing to multipath propagation or various interference [6,7]. In such cases, the subspace-based methods such as estimation of signal parameters via invariance techniques (ESPRIT) [8] and multiple signal classification (MUSIC) [9] will suffer from performance degradation owing to the rank deficiency in the source covariance matrix. To restore its rank, the methods based on spatial smoothing were developed to solve this problem [10,11,12,13,14,15]. However, because of the reduced array aperture in the overlapped subarrays, a common drawback is that they reduce the ability to identify the source number that can be resolved.
Methods such as maximum likelihood [16] and subspace fitting [17] have been proposed to resolve coherent sources. The coherent DOA estimation problem can be achieved via multi-dimensional parameter optimization without using eigenvalue decomposition. Thus, they are insensitive to signal coherence. However, the optimal solution is obtained via alternating maximization, which further results in computational burden. Meanwhile, compressed sensing methods [18,19,20] have been developed for coherent signal DOA estimation. These methods exploit the spatial sparse characteristics of the source signals and then solve the coherent signal DOA estimation through sparse vector recovery. However, the optimal solution needs to be implemented interactively, which results in high computational complexity. The high-order cumulants methods [21,22] have also been employed to resolve the coherent signals, which produce better estimation performance because of the enhanced noise suppression ability. Meanwhile, this process can produce an extended virtual array that can distinguish more sources. Nonetheless, the high-order cumulants estimation accuracy usually requires more snapshots and further results in increased computational complexity.
Recently, the idea of the so-called Toeplitz matrix reconstruction technique with moderate computational cost was also proposed to obtain the DOA estimation of the coherent sources [23,24,25,26]. The source covariance matrix in the reconstructed Toeplitz matrix corresponding to its row vector of the array output covariance matrix is a diagonal matrix with the restored full rank. Then, the DOAs are obtained by combining with conventional subspace-based methods. In [27], the beamspace-based Toeplitz matrix reconstruction method based on beamspace processing was proposed for its noise suppression ability. Furthermore, it performs better estimation in cases of a low SNR regime. However, similar to spatial smoothing methods, the existing Toeplitz matrix reconstruction methods realize the signal decorrelation at the cost of reducing the array aperture, which further results in resolution capability degradation for resolving more source signals.
To address the aforementioned problems, an enhanced coherent DOA estimation method via Toeplitz matrix reconstruction and Khatri–Rao subspace is designed to resolve more sources. Firstly, a set of Toeplitz matrices associated with each row vector of the array output covariance matrix is constructed, and then the reconstructed matrices are vectored and stacked to achieve the covariance matrix reconstruction of the array output. Furthermore, the virtual array aperture for the proposed method is about twice as big as that of the existing Toeplitz matrix methods, which further provides us with the capability of resolving more sources. Next, linear processing is employed to reduce the dimension, which leads to a lower computational cost in implementation. Meanwhile, it also achieves noise-whitening processing. In the end, the estimated DOA of coherent sources can be obtained via the conventional subspace-based methods. Numerical results also show that the proposed method can resolve more sources compared to the existing methods.
This paper is structured as follows: Section 2 and Section 3 present the signal model and the existing method, and the hypothesis proposed in this paper is further introduced in Section 4. In Section 5, several analysis experiments are then presented. Finally, the paper is concluded in Section 6.
Notation: (.)H represents conjugate transpose, (.)T denotes transpose, (.)* is conjugate. X, x, and x denote the matrix, vector, and scalar, respectively. E[.] is the expectation operation. I M denotes M × M identity matrix.

2. Signal Model

Consider a uniform linear array (ULA) that contains N = 2 M + 1 sensors with half-wavelength ( d = λ / 2 ) as shown in Figure 1, and P narrowband far-field signals with distinct directions θ i , i = 1 , 2 , , P are received by the array. Moreover, the complex envelope for source signals can be denoted as s i ( t ) , i = 1 , 2 , , P at the sensor 0, which is set as the reference.
The array output vector y ( t ) at time t can be written as
y ( t ) = [ y M ( t ) , , y 0 ( t ) , , y M ( t ) ] T       = x ( t ) + n ( t )       = A s ( t ) + n ( t )
where the s ( t ) = [ s 1 ( t ) , , s P ( t ) ] T denotes P × 1 source vector. The ( 2 M + 1 ) × P array steering matrix is expressed as A = [ a ( θ 1 ) , , a ( θ P ) ] with the i-th steering vector a ( θ i ) = [ e j ( 2 π / λ ) M d sin θ i , , 1 , , e j ( 2 π / λ ) M d sin θ i ] T . n ( t ) = [ n M ( t ) , , n M ( t ) ] T represents the ( 2 M + 1 ) × 1 noise vector with Gaussian white characteristics which are assumed to be uncorrelated to signal sources with zero mean and variance σ n 2 .
Based on the assumptions described in (1), the covariance matrix of the array output R can be expressed as
R = E [ y ( t ) y H ( t ) ]       = A E [ s ( t ) s H ( t ) ] A H + E [ n ( t ) n H ( t ) ]       = A R s A H + σ n 2 I ( 2 M + 1 )
where R s = E [ s ( t ) s H ( t ) ] is the P × P source covariance matrix. When performing eigen-decomposition, the rank loss of the matrix R s will result in performance degradation in DOA estimation or even failure [12,13]. Therefore, the key processing for the coherent DOA estimation when combined with the subspace-based methods is to restore the rank of the source covariance matrix [25,26,27].

3. Toeplitz Matrix Reconstruction Method

The Toeplitz matrix reconstruction methods [24,25,26,27] exploit a set of Toeplitz matrices with restored rank associated with row vectors to achieve the matrix reconstruction of the array output which restores the rank of the source covariance matrix. Thus, the DOAs of sources can be further obtained when using the subspace-based methods. In a manner akin to spatial smoothing techniques, the Toeplitz matrix reconstruction process also reduces the array aperture, resulting in a decrease in resolution for more signals.
Using the result in [26], the reconstructed matrix R Y i with full rank corresponding to the i-th row vector of the array output covariance matrix R is given as
R Y i = E R Y ( t ) y i ( t )     = E R X ( t ) x i ( t ) + E R N ( t ) n i ( t )     = R X i + σ n 2 I ˜ ( M + 1 ) , i
R Y ( t ) = y 0 ( t ) y 1 ( t ) y M ( t ) y 1 ( t ) y 0 ( t ) y M 1 ( t ) y M ( t ) y M + 1 ( t ) y 0 ( t )     = R X ( t ) + R N ( t )
R X ( t ) = x 0 ( t ) x 1 ( t ) x N ( t ) x 1 ( t ) x 0 ( t ) x N 1 ( t ) x N ( t ) x N + 1 ( t ) x 0 ( t )     = A ˜ S ( t ) A ˜ H
R N ( t ) = n 0 ( t ) n 1 ( t ) n N ( t ) n 1 ( t ) n 0 ( t ) n N 1 ( t ) n N ( t ) n N + 1 ( t ) n 0 ( t )
where R Y ( t ) is the reconstructed ( M + 1 ) × ( M + 1 ) Toeplitz matrix via y ( t ) . R X ( t ) and R N ( t ) denotes the reconstructed signal matrix and noise matrix with Toeplitz characteristics, respectively. A ˜ = a ˜ ( θ 1 ) , , a ˜ ( θ P ) represents the ( M + 1 ) × P array steering matrix with a ˜ ( θ i ) = 1 , , e j ( 2 π / λ ) M d sin θ i T corresponding to the reconstructed matrix R X ( t ) , and S ( t ) = d i a g s 1 ( t ) , , s P ( t ) is the reconstructed P × P source covariance matrix containing each signal on its diagonal elements.
Furthermore, the reconstructed correlation matrix based on R X ( t ) and x i ( t ) can be given by [24]
R X i = E R X ( t ) x i ( t ) = A ˜ E S ( t ) x i ( t ) A ˜ H = A ˜ S ˜ i A ˜ H
where S ˜ i denotes the correlation matrix between S ( t ) and the obtained signal from the i-th sensor.
In a similar way, the reconstructed matrix R Y i under a noise field can be written as
R Y i = A ˜ E S ( t ) x i ( t ) A ˜ H + σ n 2 I ˜ ( M + 1 ) , i = A ˜ S ˜ i A ˜ H + σ n 2 I ˜ ( M + 1 ) , i
We can conclude from (3) and (8) that the Toeplitz matrix processing named the ESPRIT-like method transforms the source matrix R s described in (1) to a diagonal matrix S ˜ i with restored full rank described in (7) to achieve signal decoherence. Then, the coherent DOAs can be obtained using the conventional subspace-based methods. However, the matrix reconstruction process reduces the effective aperture of the array since the dimension of the reconstructed matrix R Y i is smaller than that of the original matrix R , which further reduces the ability to identify the sources that can be resolved.
To avoid only employing partial information of the covariance matrix of the array output in the ESPRIT-like method, many enhanced methods including the forward and backward partial Toeplitz matrices reconstruction named FB-PTMR [24] and the multiple Toeplitz matrices reconstruction called MTOEP [25] were proposed to employ more information to achieve the matrix reconstruction. Thus, they can produce a better estimation of the reconstructed matrix and then improve the DOA estimation performance. The reconstructed covariance matrix of the array output for the FB-PTMR and MTOEP can be expressed as
R F P P T M R = i = 0 M R Y i R Y i H
R M T O E P = i = M M R Y i R Y i H
We can observe from (9) and (10) that the improved methods employ a set of Toeplitz matrices corresponding to several rows to reconstruct the covariance matrix of the array output through quadratic spatial smoothing. However, similar to the ESPRIT-like method, the reconstructed covariance matrix of the array output for the FB-PTMR and MTOEP methods still have the ( M + 1 ) × ( M + 1 ) characteristics. Thus, this will result in a reduced number of distinguishable signal sources.

4. The Proposed Method

To resolve more source signals, an enhanced DOA estimation method for coherent sources via Toeplitz matrix reconstruction and Khatri–Rao subspace is proposed. We first construct a full set of Toeplitz matrices with restored rank, and then the Khatri–Rao processing is further employed to achieve the aperture extension of the reconstructed matrix. Therefore, it improves the ability to resolve more coherent signals than that of the existing ESPRIT-like method.

4.1. Khatri–Rao Processing Criterion

Firstly, similar to the vectorization process in [28,29], the reconstructed covariance matrix of the array output shown in (8) can be rewritten as
y i vec R Y i = ( A ˜ A ˜ ) d i + vec σ n 2 I ˜ ( M + 1 ) , i
where d i = E [ s 1 ( t ) x i ( t ) ] , , E [ s P ( t ) x i ( t ) ] T is the P × 1 source vector for the reconstructed matrix; y i . vec { } denotes vectorization. represents the Khatri–Rao product
Then, when stacking Y y M , , y 0 , , y M , we obtain
Y y M , , y 0 , , y M     = ( A ˜ A ˜ ) d M , , d 0 , , d M + σ n 2 vec I ˜ ( M + 1 ) , M , , vec I ˜ ( M + 1 ) , 0 , , vec I ˜ ( M + 1 ) , M     = ( A ˜ A ˜ ) ψ + σ n 2 G
where A ˜ A ˜ is the virtual array steering response. ψ denotes the P × ( 2 M + 1 ) signal matrix. G is the ( M + 1 ) 2 × ( 2 M + 1 ) noise matrix. ψ and G have the following expressions
ψ = d M , , d 0 , , d M = E [ s 1 ( t ) x M ( t ) ] E [ s 1 ( t ) x M 1 ( t ) ] E [ s 1 ( t ) x M ( t ) ] E [ s 2 ( t ) x M ( t ) ] E [ s 2 ( t ) x M 1 ( t ) ] E [ s 2 ( t ) x M ( t ) ] E [ s P ( t ) x M ( t ) ] E [ s P ( t ) x M 1 ( t ) ] E [ s P ( t ) x M ( t ) ]
G = vec I ˜ ( M + 1 ) , M , , vec I ˜ ( M + 1 ) , 0 , , vec I ˜ ( M + 1 ) , M = 0 0 0 0 0 0 1 0 0 0 1 0 0 1 1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 1 0 0 0 0 0 0

4.2. The Number of Sources Which Can Be Estimated

It has been proven in [28] that for A ˜ C ( M + 1 ) × P as described in (12), the Khatri–Rao product A ˜ A ˜ is a full-column rank matrix assuming that the source number satisfies P 2 M + 1 . That is
R a n k A ˜ A ˜ = P
where R a n k { } denotes the standard rank.
Furthermore, using the result in [28], the signal matrix ψ C P × ( 2 M + 1 ) can maintain a full rank assuming that the source power distributions over the columns are different and the source number satisfies P 2 M . The i-th row for the matrix ψ in (13) can be written as
ψ i = E [ s i ( t ) x M ( t ) ] , E [ s i ( t ) x M 1 ( t ) ] , , E [ s i ( t ) x M ( t ) ]
Firstly, substituting (1) into (16), we have
ψ i = l = 1 P E [ s i ( t ) s l ( t ) ] e j ( 2 π / λ ) M d sin θ l , , l = 1 P E [ s i ( t ) s l ( t ) ] e j ( 2 π / λ ) M d sin θ l
Note that as defined in (1), because the source DOAs are distinct ( θ i θ j , i j ), the source power distributions within the vector ψ i are different. Meanwhile, by combining the assumption P 2 M as shown in (16), the reconstructed source matrix ψ is of full-column rank [28,29], which has the following expression
R a n k ψ = P
Then, inserting (15) and (18) into (12), we have the following expression [28,29]
R a n k Y = R a n k A ˜ A ˜ = P
When performing the singular value decomposition (SVD), we obtain
Y = [ U s   U n ] Σ s 0 0 0 V s H V n H
where U n C ( M + 1 ) 2 × ( ( M + 1 ) 2 P ) is the left singular matrices associated with zero singular values, and V n C ( 2 M + 1 ) × ( ( M + 1 ) 2 P ) denotes its corresponding right singular matrices. Meanwhile, U s C ( M + 1 ) 2 × P and V s C ( 2 M + 1 ) × P represent the left and right matrices corresponding to their nonzero singular values, respectively. Σ s R P × P is a diagonal matrix that contains nonzero singular values.
According to the array signal processing theory in [1,8,9], we have
R A ˜ A ˜ = R U s = R U n
where R ( ) represents the range space.
Finally, the estimated DOAs can be derived as follows
U n A ˜ A ˜ k = U n a ˜ ( θ k ) a ˜ ( θ k ) = 0
where denotes the Kronecker product.
By employing the conclusions in [28,29], the actual angles θ k , k = 1 , 2 , , P via the criterion in (22) based on Khatri–Rao processing can be estimated when the number of source signals P satisfy
P 2 ( M + 1 ) 2 2 M
From (8)–(10) and (23), we can conclude that the maximum number of sources that can be resolved in our proposed method, given by 2 M , are greater than that of the existing Toeplitz matrix, which is M . This is because of the proposed method based on Khatri–Rao subspace processing, which can extend the array aperture as shown in (12).

4.3. Dimension Reduction

It can be seen in (12) and (22) that the array dimension in our proposed method, given by ( M + 1 ) 2 , is greater than the existing Toeplitz-based methods [23,24,25,26,27], that is ( M + 1 ) , which increased the computational burden when it is combined with subspace-based methods. Furthermore, the noise contribution non-scalar matrix G in (14) needs to be further processed, which also impacts the computational efficiency. For these problems, we can solve them by performing linear processing on the reconstructed matrix Y , which reduces the dimension and further achieves complexity reduction. Meanwhile, the noise is also whitened during this processing. Therefore, the proposed method does not require noise elimination, making it less computationally complex.
Firstly, according to the reference [29], the virtual array response matrix A ˜ A ˜ in (12) can be computed as follows
A ˜ A ˜ = G B
where B = [ b ( θ 1 ) , b ( θ 2 ) , , b ( θ P ) ] denotes the ( 2 M + 1 ) × P dimension-reduced array steering matrix with b ( θ i ) = e j ( 2 π / λ ) M d sin θ i , , 1 , , e j ( 2 π / λ ) M d sin θ i T .
Then, substituting (24) into (12), the reconstructed matrix Y can be rewritten as
Y = G B ψ + σ n 2 G
where σ n 2 denotes the noise energy.
Next, it can be verified from (14) that the column vectors contained in the matrix G with full-column rank are orthogonal. Thus, the diagonal matrix W can be further obtained, which has the following expression [29]
W = G T G = diag 1 , 2 , , M , M + 1 , M , , 2 , 1
By combining (25) and (26), the reconstructed matrix with dimension reduction can be achieved through a linear operation [29]. We can write that
Y ˜ = W 1 G T Y     = W 1 G T ( G B ψ + σ n 2 G )     = W 1 G T G B ψ + σ n 2 W 1 G T G     = B ψ + σ n 2 I ( 2 M + 1 )
According to the results from (12) and (27), when performing linear processing, the reconstructed matrix is transformed from Y C ( M + 1 ) 2 × ( 2 M + 1 ) into Y ˜ C ( 2 M + 1 ) × ( 2 M + 1 ) , which significantly reduces the complexity burden in implementations without losing the effective degree. Meanwhile, noise-whiten processing is also achieved, which further results in better DOA estimation accuracy than that in color-noise cases [30,31,32].
The auto-correlation and spatial smoothing for the dimension-reduced matrix Y ˜ are employed to improve the noise suppression ability [26,33]. We obtain
R ˜ = Y ˜ Y ˜ H
R p r o p o s e d = R ˜ + J ( R ˜ ) J
where J is the exchange matrix.
Finally, the DOAs for coherent sources can be estimated via the total least squares ESPRIT (TLS-ESPRIT) method [34].
The major steps of the proposed method are summarized in Algorithm 1.
Algorithm 1: Pseudocode of the proposed DOA estimation method
Input: Array output vector, y ( t ) .
Output: Estimated signal source DOAs, θ d , d = 1 , , P .
1:Apply R = l = 1 L y ( t ) y H ( t ) instead of (2) to calculate the output covariance matrix.
2:Construct the Toeplitz matrix R Y i , i = M , , M with full rank corresponding to the i-th row vector of the covariance matrix R using (3).
3:Perform the Khatri–Rao processing criterion to obtain the reconstructed matrix with aperture extension Y using (11) and (12).
4:The reconstructed matrix Y ˜ with dimension reduction is obtained using (25)–(27).
5:Obtain the reconstructed covariance R P r o p o s e d with better noise suppression ability using (28) and (29).
6: Perform TLS-ESPRIT (total least squares ESPRIT) to estimate the DOAs of the signal sources.
For the TLS-ESPRIT method, the computational complexity of our proposed method mainly lies in the calculation of Y , dimension reduction, and DOA estimation. Thus, we add these processes to obtain the computational burden. 2 M + 1 denotes the number of elements in the array, L represents the snapshot number, X denotes the smoothing number of FOSS/FBSS methods, and P represents the number of source signals. The main computational complexities of these methods are shown in Table 1.
Compared with the existing Toeplitz and spatial smoothing methods, our proposed method exploits the Khatri–Rao processing to achieve the aperture extension of the reconstructed matrix, which increases the computational burden. However, it can provide the capability to identify more sources. Meanwhile, the DOA estimation performance is also improved.
In addition, the proposed method in this paper mainly emphasizes that the array aperture extension can be achieved by combing with Khatri–Rao subspace processing, and further provides better DOA estimation accuracy in cases of low SNR when compared to the existing Toeplitz matrix reconstruction methods. Moreover, the aperture extension essentially provides us with the ability to resolve more source signals than that of the existing Toeplitz matrix reconstruction methods. For simplicity of analysis, we select the ESPRIT-like [23] method as an estimator to reflect the advantages of Khatri–Rao subspace processing. Thus, when combined with another high resolution-method [26,27,30], it will also provide better performance on estimation and resolution. The DOA estimation performance is also affected by element spacing, array patterns, and mutual coupling [1,2,7]. However, this paper mainly focuses on the advantages of Khatri–Rao subspace processing. Thus, DOA estimation under these conditions will be the topic of our future research.

5. Simulation Results

Numerical simulations are presented to illustrate the superiority of our proposed method compared with the spatial smoothing methods including FBSS [11], FOSS [10], and the Toeplitz reconstruction methods such as MTOEP [25], ESPRIT-like [23], and FB-PTMR [24] based on the basis of the probability of resolution (POR) and the root-mean-square error (RMSE). The RMSE and POR of the DOA estimates are defined as
R M S E = 1 W P i = 1 P w = 1 W θ ˜ i , w θ 2
P O R = 1 W P i = 1 P N i × 100 %
where W represents the Monte Carlo trial number, which is set to 2000 in our simulations. θ ˜ i , w is the i-th estimated value of θ i in the w-th trial. P denotes the source number assumed to be known or accurately estimated [35,36,37]. The resolution is considered successful when the estimated value satisfies | θ ˜ i , w θ | 2 , and N i counts the number of successful resolutions for the i-th signal source.

5.1. Estimated Source Number Verfication

Firstly, we verify the maximum source number that can be resolved for our proposed method. As shown in Figure 2, the case where ( N , P ) = ( 5 , 4 ) is considered. A group of two coherent sources located at [ 40 , 20 ] and another group contains two coherent sources received from [ 10 , 45 ] . The spacing between the adjacent elements is half-wavelength ( d = λ / 2 ), and the number of snapshots is 200. The noise is zero-mean white Gaussian with power σ n 2 and the SNR varies from −10 dB to 12 dB according to S N R = 10 log ( σ n 2 ) [38].
The existing spatial smoothing and Toeplitz matrix reconstruction methods can not resolve the coherent sources under this condition. However, our proposed method is still able to identify the sources because our proposed method exploits the Khatri–Rao subspace processing that provides an extended virtual array, which is consistent with the analysis in Section 4. Thus, our proposed method provides the ability to resolve the angle estimation problem for more signals.

5.2. Estimation Performance versus SNR

We then consider an overdetermined case where ( N , P ) = ( 7 , 4 ) and our proposed method provides the ability to provide better estimation accuracy than other estimators in different SNRs. The true DOAs are { θ 1 , θ 2 , θ 3 , θ 4 } = { 20 , 10 , 12 , 24 } , where { 20 , 10 } is a group of coherent signals and { 12 , 24 } is another group. The number for subarrays and snapshots is set to 3 and 100, respectively. The other simulation conditions are similar to Example 5.1. Moreover, the Cramér–Rao bound (CRB) [39,40] is also presented.
As observed in the simulation results of Figure 3, our proposed method can still achieve coherent DOA estimation when compared with the existing methods based on Toeplitz matrix reconstruction. This is because the array dimension for the Toeplitz methods is less than the number of sources, while our proposed method, based on Khatri–Rao subspace processing, employs a virtual array structure to extend the physical array aperture, which results in resolving more signals. Meanwhile, the Khatri–Rao processing also eliminates the color noise. This aligns with the theoretical analysis presented in Section 4.
Figure 3 also shows that when S N R 6   dB , our proposed method provides better estimation performance on accuracy and resolution. Meanwhile, the DOA estimation performance is closer to CRB than other methods. Furthermore, the proposed method can produce a slightly interior RMSE performance at higher SNRs than the FBSS method, but the resolution probability in our method obviously outperforms the FBSS method. Particularly when S N R = 0   dB , our proposed algorithm improves estimation accuracy ( R M S E 7 ) and resolution probability ( P O R 74 % ) than that of other methods ( R M S E 10 , P O R 27 % ).

5.3. Estimation Performance versus Snapshots

The DOA estimation performance concerning the snapshot numbers is also investigated. The simulation conditions are similar to those in Example 5.2, except that the snapshot number varies from 32 to 256 and the SNR is fixed at 0   dB .
Figure 4 illustrates that our proposed method provides good estimation accuracy over the entire snapshot regime. When compared with the other methods, our proposed method can realize the array aperture extension using the Khatri–Rao subspace processing. Thus, our proposed method yields better estimation accuracy and resolution under the same snapshot number. Figure 4 also shows that the DOA estimation performance for our proposed method is significantly improved faster than other methods as an increase in the number of snapshots. This is especially true for a low number of snapshots such as 96; our proposed method can provide better estimation accuracy ( R M S E 7 ) than other methods ( R M S 11 ). Moreover, the resolution of probability in our proposed method is about 74 % , while those of the others are below 26 % .

5.4. Estimation Performance versus Angular Separation

We tested the estimation performance versus the angular separation. Two groups of coherent signals are also considered, which are set to { 30 , 30 + Δ θ } and { 0 , 0 + Δ θ } , respectively. The angular separation Δ θ varies from 2 to 18 , and the snapshot number is set to 200. The other conditions are the same as those in Example 5.1.
From the simulation results as shown in Figure 5, we can conclude that our proposed method provides better estimation performance than that of other estimators across the entire angular separation regime. Compared to the FOSS and FBSS methods, our proposed method achieves the extension of the array physical aperture when combined with the Khatri–Rao subspace processing. Therefore, this will result in better estimation accuracy. Furthermore, Figure 5 also shows that when the angular separation is enlarged, the RMSEs of the spatial smoothing methods gradually approach and are still lower than our proposed method. Meanwhile, the resolution probability of our proposed method is also significantly improved. In particular, when the angular separation Δ θ = 18 , the estimation accuracy of the FBSS is similar to our proposed method, while the resolution probability ( P O R 90 % ) in our scheme is higher than that of the FBSS method ( P O R 53 % ).

6. Conclusions

An enhanced ESPRIT-like method that combines Khatri–Rao product processing and Toeplitz reconstruction was proposed to solve the DOA estimation for coherent sources. The proposed method firstly constructs the Toeplitz matrix with full rank to achieve decorrelation. Then, the Khatri–Rao characteristic contained in the array response is inherently further employed to achieve the array aperture extension, which further provides the capability to identify more sources that can be resolved by the other methods. The numerical results confirm the validity of our proposed method. Furthermore, our analysis also demonstrates that the proposed method can produce better estimates of the DOA in cases of low snapshots, low SNR regime, and closely spaced sources.

Author Contributions

Conceptualization, B.Q.; Methodology, B.Q., X.L., D.D. and R.H.; Software, B.Q.; Validation, B.Q.; Formal analysis, X.L, D.D. and R.H.; Investigation, R.H.; Resources, B.Q. and X.L.; Data curation, D.D. and R.H.; Writing—original draft, B.Q. and R.H.; Writing—review& editing, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Krim, H.; Viberg, M. Two decades of array signal processing research: The parametric approach. IEEE Signal Process. Mag. 1996, 13, 67–94. [Google Scholar] [CrossRef]
  2. Godara, L.C. Application of antenna arrays to mobile communications. II. Beam-forming and direction-of-arrival considerations. Proc. IEEE 1997, 85, 1195–1245. [Google Scholar] [CrossRef]
  3. Kim, S.; Oh, D.; Lee, J. Joint DFT-ESPRIT Estimation for TOA and DOA in Vehicle FMCW Radars. IEEE Antennas Wirel. Propag. Lett. 2015, 14, 1710–1713. [Google Scholar] [CrossRef]
  4. Lee, H.; Ahn, J.; Kim, Y.; Chung, J. Direction-of-Arrival Estimation of Far-Field Sources Under Near-Field Interferences in Passive Sonar Array. IEEE Access 2021, 9, 28413–28420. [Google Scholar] [CrossRef]
  5. Amin, M.G.; Wang, X.; Zhang, Y.D.; Ahmad, F.; Aboutanios, E. Sparse Arrays and Sampling for Interference Mitigation and DOA Estimation in GNSS. Proc. IEEE 2016, 104, 1302–1317. [Google Scholar] [CrossRef]
  6. Yan, H.; Fan, H. On source association of DOA estimation under multipath propagation. IEEE Signal Process. Lett. 2005, 12, 717–720. [Google Scholar]
  7. Xie, W.; Wen, F.; Liu, J.; Wan, Q. Source Association, DOA, and Fading Coefficients Estimation for Multipath Signals. IEEE Trans. Signal Process. 2017, 65, 2773–2786. [Google Scholar] [CrossRef]
  8. Roy, R.; Kailath, T. ESPRIT-estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef]
  9. Schmidt, R. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
  10. Shan, T.J.; Wax, M.; Kailath, T. On spatial smoothing for direction-of-arrival estimation of coherent signals. IEEE Trans. Acoust. Speech Signal Process. 1985, 33, 806–811. [Google Scholar] [CrossRef]
  11. Pillai, S.; Kwon, B. Forward/backward spatial smoothing techniques for coherent signal identification. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 8–15. [Google Scholar] [CrossRef]
  12. Choi, Y.-H. On conditions for the rank restoration in forward/backward spatial smoothing. IEEE Trans. Signal Process. 2002, 50, 2900–2901. [Google Scholar] [CrossRef]
  13. 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]
  14. Du, W.; Kirlin, R. Improved spatial smoothing techniques for DOA estimation of coherent signals. IEEE Trans. Signal Process. 1991, 39, 1208–1210. [Google Scholar] [CrossRef]
  15. Pan, J.; Sun, M.; Wang, Y.; Zhang, X. An Enhanced Spatial Smoothing Technique with ESPRIT Algorithm for Direction of Arrival Estimation in Coherent Scenarios. IEEE Trans. Signal Process. 2020, 68, 3635–3643. [Google Scholar] [CrossRef]
  16. Vincent, F.; Besson, O.; Chaumette, E. Approximate Unconditional Maximum Likelihood Direction of Arrival Estimation for Two Closely Spaced Targets. IEEE Signal Process. Lett. 2015, 22, 86–89. [Google Scholar] [CrossRef]
  17. Yao, B.; Zhang, W.; Wu, Q. Weighted Subspace Fitting for Two-Dimension DOA Estimation in Massive MIMO Systems. IEEE Access 2017, 5, 14020–14027. [Google Scholar] [CrossRef]
  18. Das, A. Deterministic and Bayesian Sparse Signal Processing Algorithms for Coherent Multipath Directions-of-Arrival (DOAs) Estimation. IEEE J. Ocean. Eng. 2019, 44, 1150–1164. [Google Scholar] [CrossRef]
  19. Uehashi, S.; Ogawa, Y.; Nishimura, T.; Ohgane, T. Prediction of Time-Varying Multi-User MIMO Channels Based on DOA Estimation Using Compressed Sensing. IEEE Trans. Veh. Technol. 2019, 68, 565–577. [Google Scholar] [CrossRef]
  20. Das, A.; Hodgkiss, W.S.; Gerstoft, P. Coherent Multipath Direction-of-Arrival Resolution Using Compressed Sensing. IEEE J. Ocean. Eng. 2017, 42, 494–505. [Google Scholar] [CrossRef]
  21. Chen, H.; Hou, C.-P.; Wang, Q.; Huang, L.; Yan, W.-Q. Cumulants-Based Toeplitz Matrices Reconstruction Method for 2-D Coherent DOA Estimation. IEEE Sens. J. 2014, 14, 2824–2832. [Google Scholar] [CrossRef]
  22. Shi, H.; Li, Z.; Liu, D.; Chen, H. Efficient method of two-dimensional DOA estimation for coherent signals. EURASIP J. Wirel. Commun. Netw. 2017, 53, 1–10. [Google Scholar] [CrossRef]
  23. Han, F.-M.; Zhang, X.-D. An ESPRIT-like algorithm for coherent DOA estimation. IEEE Antennas Wirel. Propag. Lett. 2005, 4, 443–446. [Google Scholar]
  24. Zhang, W.; Han, Y.; Jin, M.; Li, X.-S. An Improved ESPRIT-Like Algorithm for Coherent Signals DOA Estimation. IEEE Commun. Lett. 2020, 24, 339–343. [Google Scholar] [CrossRef]
  25. Zhang, W.; Han, Y.; Jin, M.; Qiao, X. Multiple-Toeplitz Matrices Reconstruction Algorithm for DOA Estimation of Coherent Signals. IEEE Access 2019, 7, 49504–49512. [Google Scholar] [CrossRef]
  26. Qi, B.B.; Li, W. An Improved Multiple-Toeplitz Matrices Reconstruction Algorithm for DOA Estimation of Coherent Signals. Radioengineering 2021, 30, 532–539. [Google Scholar] [CrossRef]
  27. Qi, B. DOA estimation of the coherent signals using beamspace matrix reconstruction. Signal Process. 2021, 191, 108349. [Google Scholar] [CrossRef]
  28. Ma, W.; Hsieh, T.; Chi, C. DOA estimation of quasi-stationary signals via Khatri-Rao subspace. In Proceedings of the IEEE International Conference on Acoustics, Speech and Signal Processing, Taipei, Taiwan, 19–24 April 2009; pp. 2165–2168. [Google Scholar]
  29. Ma, W.-K.; Hsieh, T.-H.; Chi, C.-Y. DOA Estimation of Quasi-Stationary Signals With Less Sensors Than Sources and Unknown Spatial Noise Covariance: A Khatri–Rao Subspace Approach. IEEE Trans. Signal Process. 2010, 58, 2168–2180. [Google Scholar] [CrossRef]
  30. Fang, Y.; Zhu, S.; Gao, Y.; Zeng, C. DOA Estimation for Coherent Signals With Improved Sparse Representation in the Presence of Unknown Spatially Correlated Gaussian Noise. IEEE Trans. Veh. Technol. 2020, 69, 10059–10069. [Google Scholar] [CrossRef]
  31. Viberg, M.; Stoica, P.; Ottersten, B. Maximum likelihood array processing in spatially correlated noise fields using parameterized signals. IEEE Trans. Signal Process. 1997, 45, 996–1004. [Google Scholar] [CrossRef]
  32. Wang, Y.; Trinkle, M.; Ng, B.W.H. Two-stage DOA estimation of independent and coherent signals in spatially coloured noise. Signal Process. 2016, 128, 350–359. [Google Scholar] [CrossRef]
  33. Mei, D.; Shouhong, Z.; Xiangdong, W.; Huanying, Z. A High Resolution Spatial Smoothing Algorithm. In Proceedings of the 2007 International Symposium on Microwave, Antenna, Propagation and EMC Technologies for Wireless Communications, Hangzhou, China, 16–17 August 2007; pp. 1031–1034. [Google Scholar]
  34. Bachl, R. The forward-backward averaging technique applied to TLS-ESPRIT processing. IEEE Trans. Signal Process. 1995, 43, 2691–2699. [Google Scholar] [CrossRef]
  35. Qi, B.; Zhang, H.; Zhang, X. Time-frequency DOA estimation of chirp signals based on multi-subarray. Digit. Signal Process. 2021, 113, 103031. [Google Scholar] [CrossRef]
  36. Huang, L.; Long, T.; Mao, E.; So, H.C. MMSE-Based MDL Method for Accurate Source Number Estimation. IEEE Signal Process. Lett. 2009, 16, 798–801. [Google Scholar] [CrossRef]
  37. Chang, C.-I.; Du, Q. Estimation of Number of Spectrally Distinct Signal Sources in Hyperspectral Imagery. IEEE Trans. Geosci. Remote Sens. 2004, 42, 608–619. [Google Scholar] [CrossRef]
  38. Belouchrani, A.; Amin, M.G. Time-frequency MUSIC. IEEE Signal Process. Lett. 1999, 6, 109–110. [Google Scholar] [CrossRef]
  39. Stoica, P.; Nehorai, A. MUSIC, maximum likelihood, and Cramer-Rao bound. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 720–741. [Google Scholar] [CrossRef]
  40. Stoica, P.; Nehorai, A. Performance comparison of subspace rotation and MUSIC methods for direction estimation. IEEE Trans. Signal Process. 1991, 39, 446–453. [Google Scholar] [CrossRef]
Figure 1. Antenna array model.
Figure 1. Antenna array model.
Electronics 12 04268 g001
Figure 2. DOA estimation performance in ( N , P ) = ( 5 , 4 ) case. (a) RMSE versus SNR. (b) POR versus SNR.
Figure 2. DOA estimation performance in ( N , P ) = ( 5 , 4 ) case. (a) RMSE versus SNR. (b) POR versus SNR.
Electronics 12 04268 g002
Figure 3. DOA estimation performance versus SNR. (a) RMSE versus SNR. (b) POR versus SNR.
Figure 3. DOA estimation performance versus SNR. (a) RMSE versus SNR. (b) POR versus SNR.
Electronics 12 04268 g003
Figure 4. DOA estimation performance versus snapshot number. (a) RMSE versus snapshot number. (b) POR versus snapshot number.
Figure 4. DOA estimation performance versus snapshot number. (a) RMSE versus snapshot number. (b) POR versus snapshot number.
Electronics 12 04268 g004
Figure 5. DOA estimation performance versus angular separation. (a) RMSE versus angular separation. (b) POR versus angular separation.
Figure 5. DOA estimation performance versus angular separation. (a) RMSE versus angular separation. (b) POR versus angular separation.
Electronics 12 04268 g005
Table 1. Computational complexity analysis.
Table 1. Computational complexity analysis.
AlgorithmComputational Complexity
FOSS/FBSS O ( 2 M + 1 ) 2 L + ( 2 M X + 2 ) 3 + 3 ( 2 M X + 1 ) P 2 + 2 P 3
FB-PTMR O ( 2 M + 1 ) 2 L + ( M + 1 ) 4 + ( M + 1 ) 3 + 3 M P 2 + 2 P 3
MTOEP O ( 2 M + 1 ) 2 L + ( M + 1 ) 3 ( 2 M + 1 ) + ( M + 1 ) 3 + 3 M P 2 + 2 P 3
ESPRIT-like O ( 2 M + 1 ) 2 L + ( M + 1 ) 3 + 3 M P 2 + 2 P 3
The proposed method O ( 2 M + 1 ) 2 L + ( 2 M + 1 ) 2 ( M + 1 ) 2 + ( 2 M + 1 ) 3 + 3 ( 2 M + 1 ) P 2 + 2 P 3
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

Qi, B.; Liu, X.; Dou, D.; Zhang, Y.; Hu, R. An Enhanced DOA Estimation Method for Coherent Sources via Toeplitz Matrix Reconstruction and Khatri–Rao Subspace. Electronics 2023, 12, 4268. https://doi.org/10.3390/electronics12204268

AMA Style

Qi B, Liu X, Dou D, Zhang Y, Hu R. An Enhanced DOA Estimation Method for Coherent Sources via Toeplitz Matrix Reconstruction and Khatri–Rao Subspace. Electronics. 2023; 12(20):4268. https://doi.org/10.3390/electronics12204268

Chicago/Turabian Style

Qi, Bingbing, Xiaogang Liu, Daowei Dou, Yan Zhang, and Runze Hu. 2023. "An Enhanced DOA Estimation Method for Coherent Sources via Toeplitz Matrix Reconstruction and Khatri–Rao Subspace" Electronics 12, no. 20: 4268. https://doi.org/10.3390/electronics12204268

APA Style

Qi, B., Liu, X., Dou, D., Zhang, Y., & Hu, R. (2023). An Enhanced DOA Estimation Method for Coherent Sources via Toeplitz Matrix Reconstruction and Khatri–Rao Subspace. Electronics, 12(20), 4268. https://doi.org/10.3390/electronics12204268

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