A PARALIND Decomposition-Based Coherent Two-Dimensional Direction of Arrival Estimation Algorithm for Acoustic Vector-Sensor Arrays

In this paper, we combine the acoustic vector-sensor array parameter estimation problem with the parallel profiles with linear dependencies (PARALIND) model, which was originally applied to biology and chemistry. Exploiting the PARALIND decomposition approach, we propose a blind coherent two-dimensional direction of arrival (2D-DOA) estimation algorithm for arbitrarily spaced acoustic vector-sensor arrays subject to unknown locations. The proposed algorithm works well to achieve automatically paired azimuth and elevation angles for coherent and incoherent angle estimation of acoustic vector-sensor arrays, as well as the paired correlated matrix of the sources. Our algorithm, in contrast with conventional coherent angle estimation algorithms such as the forward backward spatial smoothing (FBSS) estimation of signal parameters via rotational invariance technique (ESPRIT) algorithm, not only has much better angle estimation performance, even for closely-spaced sources, but is also available for arbitrary arrays. Simulation results verify the effectiveness of our algorithm.


Introduction
Compared with traditional acoustic pressure sensor arrays, acoustic vector sensors can measure the acoustic pressure as well as all three orthogonal components of the acoustic particle velocity at a single point in space, which offers certain significant advantages in collecting more acoustic information and enhancing the system performance [1][2][3][4][5][6]. Since the measurement model of an acoustic vector sensor array was developed in [2], researchers mainly turned to the study of incoming signal direction of arrival (DOA) estimation and have proposed many DOA estimation algorithms, which include the Capon technique [4], estimation of signal parameters via rotational invariance technique (ESPRIT) [7][8][9], root multiple signal classification (MUSIC) [10], self-initiating MUSIC [11], hypercomplex [12], propagator method (PM) [13], trilinear decomposition or parallel factor (PARAFAC) [14], as well as others [15][16][17][18][19][20][21]. The subspace-based methods, such as ESPRIT [7][8][9] and MUSIC [10,11], require eigenvalue decomposition (EVD) of the cross correlation matrix and singular value decomposition (SVD) of the received data to obtain the signal subspace or noise subspace, which implies fairly high computational complexity, while the propagator method (PM) is considerably less demanding because the PM does not require any EVD of the cross correlation matrix and SVD of the received data [22]. However, only in high-snapshots situation and signal to noise ratio (SNR), can the PM algorithm provide a better estimation performance. In most current algorithms for DOA estimation, some precise a priori knowledge, including the sensor locations, gain, phase response and mutual coupling of the receiver array is needed, but in realistic situations this can seldom be pre-known. For example, Capon [4] and the MUSIC algorithm can be used for arbitrary arrays, however they need to pre-know the array geometry. Besides, peak searching is also required for Capon and MUSIC algorithms, which renders a heavier computational cost. EVD of the cross spectral matrix or SVD of the received data to obtain the signal subspace is needed for the ESPRIT algorithm [7], which has been used for two dimensional (2D) DOA estimation for arbitrarily spaced arrays at unknown locations based on the acoustic vector-sensor property. However, a problem existing in [7] is that the ESPRIT algorithm needs an extra pair matching which increases the computational load, and usually fails to work in lower SNR when it is used for 2D-DOA estimation. Reference [14] proposes a trilinear decomposition-based 2D-DOA estimation algorithm for acoustic vector sensor arrays, which provides DOA estimation for the arbitrarily spaced sensor arrays and doesn't require knowledge of sensor locations and extra pair matching.
The reflected signal of the same sources through different propagation paths will produce multipath signals. Therefore, it is of significant importance to study the coherent angle estimation problem. The angle estimation algorithms listed above are all proposed for incoherent sources. When it comes to coherent sources, the coherency of sources will result in serious degradation or invalidity of the above algorithms. Some conventional coherent angle estimation algorithms, including forward backward spatial smoothing (FBSS) [23,24], only work for uniform arrays.
This paper combines the acoustic vector-sensor array parameter estimation problem with the so-called PARAllel profiles with LINear Dependencies (PARALIND) model, and proposes a blind coherent 2D-DOA estimation algorithm for arbitrarily spaced acoustic vector-sensor arrays subject to unknown locations by exploiting the PARALIND decomposition approach. Our algorithm can provide coherent and incoherent two-dimensional angle estimation for arbitrary arrays, and it automatically archives paired azimuth and elevation angles, and the paired correlated matrix of the sources can also be acquired. Compared with conventional coherent angle estimation algorithms such as the FBSS-ESPRIT algorithm which only works for uniform arrays, our algorithm has much better angle estimation performance. Furthermore, our algorithm performs considerably well for angle estimation of closely spaced sources. We also derive the Cramér-Rao bound (CRB) of angle estimation for arbitrarily spaced acoustic vector-sensor arrays. Simulation results verify the effectiveness of the proposed algorithm.
The trilinear decomposition, also known as PARAFAC analysis [25,26], has been naturally related to angle estimation for arbitrarily spaced acoustic vector-sensor arrays at unknown locations [14]. However, the PARAFAC angle estimation solution is usually non-unique when the coherent sources exist. The PARALIND model [27,28] is a generalization of PARAFAC suitable for solving problems with linear dependent factors where PARAFAC analysis will fail to provide meaningful results. Our work links the coherent angle estimation problem to the PARALIND model, and proposes a PARALIND decomposition-based coherent angle estimation algorithm for arbitrary arrays, which can be regarded as an extension of the work presented in [14]. Reference [27] proposed the PARALIND model for application in biology and chemistry; the present paper expands this model to the acoustic vector-sensor array signal processing problem to estimate coherent DOA and automatically achieve paired two-dimensional angle estimations, which is an innovation.
Although the ESPRIT algorithm [7] and our algorithm can be used for DOA estimation for arbitrarily spaced acoustic vector-sensor arrays, there are some differences between our paper and Reference [7]. The latter proposed an ESPRIT algorithm for DOA estimation for arbitrarily spaced three-dimensional arrays of vector hydrophones, but it fails to work well for coherent sources. Our work exploits the PARALIND decomposition approach to estimate 2D-DOA in arbitrarily spaced acoustic vector-sensor arrays with unknown locations, and our algorithm is suitable for coherent sources.
The present paper is structured as follows: Section 2 develops the data model for arbitrarily spaced acoustic vector-sensor arrays at unknown locations; Section 3 establishes our PARALIND decomposition-based coherent 2D-DOA angle estimation algorithm in addition with the identifiability issues and complexity analysis; In Section 4, simulation results are presented to verify effectiveness of the proposed algorithm, while the final conclusions are made in Section 5.

Data Model
We assume that a total of K narrowband plane waves impinge on an array equipped with M acoustic vector sensors, which are all located at arbitrary unknown three dimensional positions as shown in Figure 1. We define the location of the mth vector sensor as r m = (x m , y m , z m ).
We also assume the signals in the far-field. K sources, including K 1 incoherent sources and K-K 1 coherent sources, are considered. Assume that the noise is additive white Gaussian, which is independent of the sources. The kth signal arrives from direction (φ k , ϕ k ), where φ k and ϕ k respectively stand for the azimuth angle and the elevation angle. Let θ k = [φ k , ϕ k ] T as the 2D-DOA of the kth source.
According to [2], the output at the acoustic vector sensors array can be expressed as: T is a column vector consisting of amplitudes and phases of the K 1 incoherent sources. Γ is the correlated matrix with K × K. n(t) is the received additive white Gaussian noise vector.
Therefore, the output of J snapshots can be given by: X can be compactly expressed as: where S = [s(t 1 ), s(t 2 ),…, s(t J )], N = [n(t 1 ), n(t 2 ),…, n(t J )] is an 4M × J matrix composed of J snapshots of received additive white Gaussian noise.

PARALIND Decomposition
Define T = Y X % , according to the signal model in Equation (4), Y % can be expressed as: According to [27], least fitting for the signal model in Equation (5) amounts to: In a no-noise case, according to Equation (5), we have: According to [7], stacking these vectors leads to: Equation (9) can be expressed compactly as: According to Equation (10), we can obtain: For attaining vec(Γ T ), then Γ can be easily updated via transforming the column vector to its original column-wise matrix.
According to Equation (6), the least square (LS) update for S T is given by: According to Equation (7), we have: A has full column rank, so we can acquire A * via: Similarly, we have: Extracting the diagonal elements of the matrices in two sides of the equation, we get: Then we get: The matrix H can be straightforwardly obtained via diag −1 (D n (H)), n = 1,2,3,4. According to Equations (11), (12), (14), and (17), we can show PARALIND algorithm applied in the data model established in this paper in detail as follows: According to Equation (11), the update for ( ) T vec Γ is given by: where Y % is the noisy signal. Ĥ , Â , and ˆT S are the previously obtained estimates of H, A, and S T , respectively. According to Equation (12), the LS update for S T is obtained via: where Ĥ , Â , and Γ are the previously obtained estimates of H, A, and Γ, respectively. According to Equation (14), the update for A * is shown as: where n Y % denotes the noisy signal. Ĥ , Ŝ , and Γ are the previously obtained estimates of H, S, and Γ, respectively. According to Equation (17), the update for D n (H) is: where Â , Ŝ , and Γ are the previously obtained estimates of A, S, and Γ, respectively. Finally, the update of H , noted as Ĥ , can be straightforwardly obtained via , where e ji is the (j, i) element of the matrix E. According to Equations (18), (19), (20) and (21), the matrices A, H, Γ, and S are updated until the SSR ≤ 10 −8 , finally we obtain Â , Ĥ , Ŝ , and Γ .

Uniqueness of PARALIND Decomposition
According to [27], we derive the uniqueness of PARALIND decomposition in an acoustic vector-sensor array. The signal matrix in Equation (7) can be transformed to another equivalent matrix via column and row exchanging, which can be expressed as: The two slices Z i and Z j (i ≠ j) in Equation (22) are represented as: where . Then we form the following matrix: H being full row rank assures that ( ) span span The matrix U can be denoted as: where T is a nonsingular matrix. Construct auto and cross correlation matrices as follows: where G = T H H H H, R 1 and R 2 are full rank. According to Equation (26), we obtain: Λ and G consist of the eigenvalues and the corresponding eigenvectors of the matrix R 2 R 1 −1 . Λ is unique, and G is recovered with the scale ambiguity and permutation ambiguity.
Notably, scale ambiguity and permutation ambiguity are inherent to the separation problem. However, the scale ambiguity can be resolved easily by normalization, while the existence of permutation ambiguity is not considered for angle estimation.

Two-Dimensional Angle Estimation
By imposing PARALIND decomposition for the received data matrix, we get the estimate of matrix  (28) where ˆk h is the normalized kth column vector of the estimated Ĥ . Then the elevation angle estimate of kth source can be easily obtained via: 1ŝ in ( (4)) Finally, the azimuth angle estimate of kth source can be attained via: Obviously, the azimuth angle and elevation angle are automatically paired. We also obtain the estimate of the correlated matrix Γ from the PARALIND decomposition, because Γ is a matrix whose elements are 0 or 1, so we can easily decide it by a decision function to obtain accurate estimated Γ .

The Procedures of the Proposed Algorithm
Till now, we have achieved the proposal for the proposed coherent 2D-DOA estimation algorithm for the acoustic vector-sensor array. We show major steps of the proposed algorithm as follows: Step 1: Obtain Y % from the received data matrix X, and then initialize for the matrices H, A, Γ, and S.
Remark A: Our algorithm obtains the estimate of correlated matrix Γ from PARALIND decomposition, even for partly coherent sources, the correlated matrix Γ can be acquired as well.

Remark B:
The PARALIND decomposition brings the same permutation ambiguity for the estimated Ĥ , Â , and Γ , so the elevation and azimuth angles and the correlated matrix are automatically paired.
Remark C: If the number of sources K is unknown, it can be estimated by performing SVD for a slice Z i in Equation (22) and finding the number of largest singular values.

CRB and Advantages of the Proposed Algorithm
According to [29], we can derive the CRB of coherent angle estimation for the acoustic vector-sensor array with unknown locations as: where 2 σ stands for the covariance of the noise and J denotes the number of snapshots: , , , , , Γs .
The advantages of the proposed algorithm can be presented as follows, which can be verified by the simulation results in Section 4: (1) The proposed algorithm is effective for coherent and incoherent two-dimensional angle estimation.

Simulation Results
In most of the following simulations, we assume that there are 3 sources in which only source 1 and source 3 are coherent, namely the correlated matrix is 1 0 The sources are located at angles (φ 1 ,ϕ 1 ) = (−15°, 10°), (φ 2 ,ϕ 2 ) = (25°,20°), and (φ 3 ,ϕ 3 ) = (35°,30°), respectively. M, J, and K denote the number of receive sensors, snapshots, and sources, respectively. We present 1000 Monte Carlo simulations to assess the angle estimation performance of the proposed algorithm. Define root mean squared error (RMSE) as follows: where φ k and ϕ k denote the perfect azimuth and elevation angle of kth source, respectively. ,k l φ and ,k l ϕ are the estimates of φ k and ϕ k d in the lth Monte Carlo trail.
In order to present the angle estimation performance comparison of the proposed algorithm and FBSS-ESPRIT algorithm, we assume that the acoustic vector-sensor array is a uniform linear array (ULA) in Figures 2 and 3. Figure 2 depicts the two-dimensional angle estimation with M = 12, J = 100, K = 3, and SNR = 15 dB. It illustrates that our algorithm is effective for paired two-dimensional angle estimation using ULA. Figure 3 presents the angle estimation performance comparison of the proposed algorithm, FBSS-ESPRIT algorithm, and CRB with M = 12, J = 100, and K = 3. It is shown that our algorithm possesses much better angle estimation performance than the FBSS-ESPRIT algorithm in the ULA situation.  The following Figures 4-8 are for the case of arbitrarily spaced acoustic vector-sensor arrays subject to unknown locations, and we assume that the column full rank receive direction matrix is generated randomly. Figure 4 presents the two-dimensional angle estimation of the proposed algorithm with M = 12, J = 100, K = 3, and SNR = 5 dB, and Figure 5 depicts the angle estimation with M = 12, J = 100, K = 3, and SNR = 15 dB. Figures 4 and 5 illustrate that our algorithm is effective for paired two-dimensional angle estimation using arbitrary acoustic vector-sensor arrays. Figure 6 presents the angle estimation performance comparison of the proposed algorithm and CRB with M = 12, J = 100, and K = 3, while Figure 7 depicts the angle estimation performance comparison with M = 10, J = 50, and K = 3.     . It illustrates that the estimation performance becomes better in collaboration with J increasing and the proposed algorithm is effective in small snapshots. Figure 9 shows angle estimation performance    We assume two coherent closely spaced sources located at angles (φ 1 ,ϕ 1 ) = (0°,30°), (φ 2 ,ϕ 2 ) = (2°,28°). Figure 10 displays angle estimation of closely spaced sources exploiting the proposed algorithm with M = 12, J = 100, K = 2, and SNR = 15 dB. It implies that our angle algorithm has considerably good performance for angle estimation of closely spaced sources. For non-coherent sources, the proposed algorithm becomes the PARAFAC algorithm [14]. In non-coherent sources, we compared our algorithm against the ESPRIT algorithm [7] and CRB. Figure 11 presents angle estimation performance comparisons with K = 3, M = 8 and L = 150. It is indicated our algorithm has better DOA estimation performance than the ESPRIT algorithm [7]. Figure 11. Angle estimation performance comparison for non-coherent sources.

Conclusions
In the present paper we have expanded the PARALIND model, originally applied to biology and chemistry, to acoustic vector-sensor array signal processing. The PARALIND model is used for coherent 2D-DOA estimation for arbitrarily spaced acoustic vector-sensor arrays subject to unknown locations. We have derived the PARALIND decomposition and its uniqueness in acoustic vector-sensor arrays, which is convenient for simulation, performance analysis and further study. Our algorithm obtains automatically paired 2D-DOA estimation and the correlated matrix of the sources. Compared with the FBSS-ESPRIT algorithm which only works for uniform arrays, our algorithm has much better angle estimation performance. Furthermore, our algorithm has considerably good performance for angle estimation of closely spaced sources.