Next Article in Journal
Dynamics of Hierarchical Urban Green Space Patches and Implications for Management Policy
Next Article in Special Issue
Two Novel Two-Stage Direction of Arrival Estimation Algorithms for Two-Dimensional Mixed Noncircular and Circular Sources
Previous Article in Journal
Improved Short-Term Clock Prediction Method for Real-Time Positioning
Previous Article in Special Issue
Research and Analysis on the Localization of a 3-D Single Source in Lossy Medium Using Uniform Circular Array
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two-Dimensional DOA Estimation for Coherently Distributed Sources with Symmetric Properties in Crossed Arrays

National Digital System Engineering and Technological Research R&D Center, Zhengzhou 450001, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(6), 1300; https://doi.org/10.3390/s17061300
Submission received: 21 March 2017 / Revised: 1 June 2017 / Accepted: 5 June 2017 / Published: 6 June 2017

Abstract

:
In this paper, a novel algorithm is proposed for the two-dimensional (2D) central direction-of-arrival (DOA) estimation of coherently distributed (CD) sources. Specifically, we focus on a centro-symmetric crossed array consisting of two uniform linear arrays (ULAs). Unlike the conventional low-complexity methods using the one-order Taylor series approximation to obtain the approximate rotational invariance relation, we first prove the symmetric property of angular signal distributed weight vectors of the CD source for an arbitrary centrosymmetric array, and then use this property to establish two generalized rotational invariance relations inside the array manifolds in the two ULAs. Making use of such relations, the central elevation and azimuth DOAs are obtained by employing a polynomial-root-based search-free approach, respectively. Finally, simple parameter matching is accomplished by searching for the minimums of the cost function of the estimated 2D angular parameters. When compared with the existing low-complexity methods, the proposed algorithm can greatly improve estimation accuracy without significant increment in computation complexity. Moreover, it performs independently of the deterministic angular distributed function. Simulation results are presented to illustrate the performance of the proposed algorithm.

1. Introduction

In recent decades, the problem of direction-of-arrival (DOA) estimation, which plays an important role in radar, sonar and wireless communication systems, has attracted a lot of attention. The most commonly considered system model in the DOA finding techniques is the point source model, where the signals are assumed to arrive at the array via a single path [1,2,3,4]. When dealing with a point source, conventional subspace-based algorithms, such as the multiple signal classification (MUSIC) algorithm [5,6] and the estimation of signal parameters via rotational invariance techniques (ESPRIT) algorithm [7,8], have high DOA estimation resolution. However, in many practical applications, the signals will reach the array through many rays reflected or scattered from the vicinity, which causes angular spreading. In these cases, directly applying the MUSIC and ESPRIT algorithms may lead to biased estimations. Therefore, some researchers have considered a more realistic signal model called the spatially distributed source model [9,10]. Depending on the correlation among different rays, distributed sources are classified into two types: coherently and incoherently distributed (CD and ID) sources. In this paper, we only consider the DOA estimation of the CD sources.
Many DOA estimation techniques for CD sources have been published. Since conventional subspace-based methods cannot be directly applied to a distributed source, some researchers studied modifications to the MUSIC algorithm, which gave rise to the distributed signal parameter estimator [11], dispersed parametric estimator [12] and Vec-MUSIC estimator [13]. These three methods were established based on the assumption that the distribution shapes of multiply distributed sources are identical and known. In addition, the computational complexity is high because of the need for a two-dimensional (2D) search. The literature [14] used an unstructured model for the part of covariance matrix, where the 2D search problem was replaced by a successive one-dimensional (1D) search. With even lower computational complexity, a search-free algorithm called the spread root-MUSIC algorithm [15] was proposed to fit a two-ray model of the data. However, low-complexity was obtained in the special case where only one distributed source existed. In [16], the authors considered two identical and closely spaced sub-arrays. When the distance between the two sub-arrays was far shorter than the wavelength, an approximate rotational invariance relation between the two sub-arrays was obtained based on Taylor series approximation, and finally the central DOA of the CD source can be estimated by total least square estimation parameter via rotational invariance techniques (TLS-ESPRIT) using the generalized array manifold (GAM) model. All of the aforementioned works [11,12,13,14,15,16] were designed for 1D DOA and angular spread estimation of distributed sources. However, when the distributed source and the receiving sensor array are not in the same plane, it is reasonable to instead model the source as a 2D distributed source.
Since a 2D CD source is characterized by four parameters: the central azimuth direction, the azimuth angular extension, the central elevation direction and the elevation angular extension, the conventional optimum estimators will be computationally expensive owing to high-dimensional parameters [17]. Consequently, it is very necessary for 2D distributed sources to find some suboptimum algorithms to reduce the computational cost. To date, several low-complexity DOA estimation algorithms for 2D CD sources have been proposed. Specifically, the authors in [18] considered a pair of uniform circular arrays (UCAs). Preliminary estimations of central elevation DOAs were obtained using TLS-ESPRIT. Next, by using the estimated elevation DOAs, a sequential one-dimensional searching (SOS) method was proposed to estimate the central azimuth DOAs. In [19], using two parallel uniform linear arrays (ULAs), a low-complexity algorithm without searching was proposed for CD sources. Similarly, central elevation DOAs are obtained based on the approximate rotational invariance relation between the two ULAs. Instead of SOS, the quadric rotational invariance property (QRIP) of the GAM was used to estimate the central azimuth DOAs. Finally, a parameter matching approach was given to obtain the correct DOA estimation. In [20], the central elevation and azimuth DOAs were both estimated based on TLS-ESPRIT, which used two parallel ULAs, and the parameter matching method was also required. In order to avoid the parameter matching procedure, the literature [21] estimated the central elevation and central azimuth DOAs by applying the singular value decomposition method to the cross-correlation (CC) matrix of the received data in the double parallel ULAs. However, all the algorithms in [18,19,20,21] were all based on the special array geometry composed of two sub-arrays. The approximate rotational invariance relation between the two sub-arrays was obtained by using the one-order Taylor series approximation, which may introduce additional errors and affect the estimation accuracy.
In array processing, the crossed array is a commonly used 2D array geometry [22]. Compared to other 2D arrays such as the UCA and plane rectangular array, the crossed array can provide a larger aperture and hence offer better resolution for a given number of elements. Moreover, the crossed array consists of two intersecting ULAs working independently, thus the computational complexity is only double that of a single dimensional array. Several 2D estimation algorithms based on the crossed array have been proposed [23,24]. However, they are all based on the point source model. To the best of our knowledge, there have been few reports about the DOA estimation for CD sources in a crossed array.
In this paper, we consider a crossed array and divide it into two sub-ULAs. In particular, instead of using the Taylor series approximation, we prove the symmetric property of the angular signal distributed weight (ASDW) vector for an arbitrary centrosymmetric array, and use this property to establish the generalized rotational invariance relations inside the GAMs for the two sub-ULAs. Resorting to such relations, the central elevation and azimuth DOAs are estimated based on polynomial-root-based search-free method, respectively. Then simple parameter matching is accomplished by searching the minimums of the cost function of the estimated 2D angular parameters. The proposed algorithm does not require that the angular distribution functions of the multiple distributed sources are the same and known. In addition, it does not suffer additional errors induced by Taylor series approximation and high computational complexity brought about by spectrum-peak searching.
The rest of this paper is organized as follows: Section 2 presents the data model. In Section 3, we describe the proposed algorithm in detail. Some simulation results which illustrate the validity and performance of the proposed method are given in Section 4. Section 5 concludes the paper.
The following notations will be used throughout this paper. Superscript ( ) , ( ) T , and ( ) H represent the conjugate, transpose and conjugate transpose operations, respectively. The symbol denotes the Schur-Hadamard product; E [ ] stands for the mathematical expectation and det ( ) is the matrix determinant; d i a g [ ] is a diagonal matrix and the values in the brackets are the diagonal elements.

2. Signal Model

Let us consider the plane crossed array presented in Figure 1. The array is centered at the origin of the three-dimensional coordinate system with two ULAs directed along the y-axis and z-axis. The ULAs Y a and Z a are composed of M y and M z omni-directional antenna elements, respectively. The distance between adjacent sensors is d in the two ULAs. We assume that there are D narrowband CD sources impinging on the crossed array. The observation vectors of Y a and Z a at time t are given by [17,18,19,20,21].
y ( t ) = i = 1 D a y ( θ , γ ) s i ( θ , γ , t ; μ i ) d θ d γ + n y ( t ) ,
z ( t ) = i = 1 D a z ( θ , γ ) s i ( θ , γ , t ; μ i ) d θ d γ + n z ( t ) ,
where y ( t ) is the M y × 1 array output vector of the sub-array Y a ; z ( t ) is the M z × 1 array output vector of the sub-array Z a ; s i ( θ , γ , t ; μ i ) is the complex random angular signal density function of the i-th source. The vector μ i = ( θ i , σ θ i , γ i , σ γ i ) determines the central azimuth DOA θ i , the azimuth angular spread σ θ i , the central elevation DOA γ i and the elevation angular spread σ γ i of the i-th sensor; n y ( t ) and n z ( t ) are Gaussian white noise with zero-mean and variance, while σ n 2 ; a y ( θ , γ ) and a z ( θ , γ ) are the two array manifold vectors at direction ( θ , γ ) :
a y ( θ , γ ) = [ e j η 0.5 ( M y 1 ) sin θ sin γ , e j η 0.5 ( M y 3 ) sin θ sin γ , ,   e j η 0.5 ( M y 1 ) sin θ sin γ ] T , a z ( θ , γ ) = [ e j η 0.5 ( M z 1 ) cos γ , e j η 0.5 ( M z 3 ) cos γ , , e j η 0.5 ( M z 1 ) cos γ ] T ,
where η = 2 π d / λ , and λ is the wavelength of the coming signal.
For a 2D CD source, the angular signal density function s i ( θ , γ , t ; μ i ) can be written as:
s i ( θ , γ , t ; μ i ) = s i ( t ) ρ i ( θ , γ ; μ i ) ,
where s i ( t ) is a random variable and ρ i ( θ , γ ; μ i ) is the deterministic angular distribution function.
Define the GAM vectors of distributed source for subarray Y a and Z a as follows:
b y ( μ i ) = a y ( θ , γ ) ρ i ( θ , γ ; μ i ) d θ d γ , b z ( μ i ) = a z ( θ , γ ) ρ i ( θ , γ ; μ i ) d θ d γ .
For small angular extensions, we have the following closed forms [11,20]:
b y ( μ i ) = a y ( θ i , γ i ) g y ( μ i ) , b z ( μ i ) = a z ( θ i , γ i ) g z ( μ i ) ,
where g y ( μ i ) and g z ( μ i ) are the ASDW vectors. The observation vectors in (1) and (2) can be written as:
y ( t ) = B y ( μ ) s ( t ) + n y ( t ) , z ( t ) = B z ( μ ) s ( t ) + n z ( t ) ,
where s ( t ) = [ s 1 ( t ) , s 2 ( t ) , , s D ( t ) ] T is a D × 1 signal vector, and B y ( μ ) and B z ( μ ) are the GAM matrices, which are composed of D GAM vectors:
B y ( μ ) = [ b y ( μ 1 ) , b y ( μ 2 ) , , b y ( μ D ) ] , B z ( μ ) = [ b z ( μ 1 ) , b z ( μ 2 ) , , b z ( μ D ) ] .
The total array output vector is expressed as:
x ( t ) = [ y ( t ) z ( t ) ] .

3. The Proposed Algorithm

This section consists of four parts. Firstly, the symmetric property of an ASDW vector is identified in a centro-symmetric array. Then, by making use of the symmetric property of the ASDW vectors in the two sub-ULAs Z a and Y a , we establish two generalized rotational invariance relations for the GAM vectors. On the premise of such relations, the central elevation and azimuth DOAs are obtained by using a polynomial-root-based search-free approach, respectively. Afterwards, and when multiple CD sources exist, a simple parameter matching approach is addressed. Finally, we provide the algorithm’s realization steps and an analysis of the computational complexity.

3.1. Symmetric Property of an ASDW Vector in a Centro-Symmetric Array

In this part, the symmetric property of an ASDW vector in a centro-symmetric array is derived in detail. Let us consider a centro-symmetric array consisting of M identical antenna elements centered at the coordinate origin, where the m-th sensor is placed at ( x m , y m , z m ) for m = 1 , 2 , M . The array manifold vector in direction ( θ , γ ) is expressed as:
a ( θ , γ ) = [ e j ( 2 π / λ ) ( x 1 cos θ sin γ + y 1 sin θ sin γ + z 1 cos γ ) e j ( 2 π / λ ) ( x 2 cos θ sin γ + y 2 sin θ sin γ + z 2 cos γ ) ...... e j ( 2 π / λ ) ( x M cos θ sin γ + y M sin θ sin γ + z M cos γ ) ] .
If we define θ = θ i + θ ˜ and γ = γ i + γ ˜ , in which θ i , γ i are the central azimuth DOA and the central elevation DOA of the i-th source, and θ ˜ , γ ˜ are the corresponding random angular deviations , the GAM vector can be presented as:
b ( μ i ) = a ( θ , γ ) ρ ( θ , γ ; μ i ) d θ d γ = a ( θ i + θ ˜ , γ i + γ ˜ ) ρ ( θ ˜ , γ ˜ ; μ i ) d θ ˜ d γ ˜ ,
For small angular extensions, the m-th element of b ( μ i ) can be written as [11]:
[ b ( μ i ) ] m [ a ( θ i , γ i ) ] m e j ς m ρ ( θ ˜ , γ ˜ ; μ i ) d θ ˜ d γ ˜ .
Thus, the m-th element of ASDW vector is given by:
[ g ( u ) ] m = e j ς m ρ ( θ ˜ , γ ˜ ; μ i ) d θ ˜ d γ ˜ ,
where:
ς m = ( 2 π / λ ) [ x m ( θ ˜ sin θ i sin γ i + γ ˜ cos θ i cos γ i ) +   y m ( θ ˜ cos θ i sin γ i + γ ˜ sin θ i cos γ i ) + z m ( γ ˜ sin γ i ) ] .
We have ( x m , y m , z m ) = ( x m + M / 2 , y m + M / 2 , z m + M / 2 ) in the centrosymmetric array, thus ς m = ς m + M / 2 . Respecting the fact that ρ ( θ ˜ , γ ˜ ; μ i ) is an even function (see Appendix A), we can obtain the symmetric property of the ASDW vector such as:
[ g ( u ) ] m = [ g ( u ) ] m + M / 2 .
It is obvious that Y a and Z a are centrosymmetric arrays, thus we have:
[ g y ( u ) ] m = [ g y ( u ) ] m + M / 2 , [ g z ( u ) ] m = [ g z ( u ) ] m + M / 2 .

3.2. Derivation

3.2.1. Central Elevation DOA Estimation

For sub-array Z a , and owing to the symmetric property of the ASDW vector in (16), we can establish the following generalized rotational invariance relation of the GAM vector:
Π M Z b z ( μ i ) = Ψ ( γ i ) b z ( μ i ) ,
where Π M Z is the M z × M z exchange matrix with ones on its anti-diagonal and zeros elsewhere. Ψ ( γ i ) is an M z × M z diagonal matrix which is given by:
Ψ ( γ i ) = d i a g [ e j η ( M z 1 ) cos γ i , e j η ( M z 3 ) cos γ i , , e j η ( M z 1 ) cos γ i ] .
If we define the complex variable k = e j η cos γ i , Ψ ( γ i ) can be written as:
Ψ ( k ) = d i a g [ k ( M z 1 ) , k ( M z 3 ) , , k ( M z 1 ) ] .
According to the observation vector z ( t ) in (7), the covariance matrix of z ( t ) is expressed as:
R z = E { z ( t ) z H ( t ) } = B y ( μ ) R s B y H ( μ ) + σ n 2 I M z ,
where R s = E { s ( t ) s H ( t ) } is the signal covariance matrix of the CD sources. The eigenvalue decomposition of R z is given by:
R z = U z Σ z U z H + U n Σ n U n H ,
where U z is the signal subspace matrix, whose columns are the eigenvectors corresponding to the D largest eigenvalues of R z . When R s is of full rank, the subspace spanned by the columns of U z is equal to the subspace spanned by the columns of B z ( μ ) . At this point, there exists a unique non-singular D × D matrix T 1 such that U z = B z ( μ ) T 1 . According to the generalized rotational invariance relation in (17), we can formulate a matrix F z ( k ) :
F z ( k ) = Π M Z U z Ψ ( k ) U z = Π M Z B z ( μ ) T 1 Ψ ( k ) B z ( μ ) T 1 = [ ( Ψ ( k 1 ) Ψ ( k ) ) b z ( μ 1 ) , ( Ψ ( k 2 ) Ψ ( k ) ) b z ( μ 2 ) ,   , ( Ψ ( k D ) Ψ ( k ) ) b z ( μ D ) ] T 1 .
Therefore, when Ψ ( k i ) = Ψ ( k ) for i = 1 , 2 , , D , the i-th column of F z ( k ) is a zero vector. Hence, F z ( k ) is rank deficient and the determinant of F z H ( k ) F z ( k ) is zero. The central elevation DOA estimations γ ^ i ( i = 1 , 2 , D ) can be obtained by finding the highest D peaks of the following spectrum function:
H z ( k ) = 1 / det ( F z H ( k ) F z ( k ) ) .
However, estimator (23) involves computationally intensive spectral-peak searching. In order to reduce the complexity, we derive a polynomial-root-based search-free approach. The denominator of (23) can be written as the following polynomial:
h z ( k ) = det ( F z H ( k ) F z ( k ) ) .
It is obvious that the central elevation DOAs can be obtained by rooting this polynomial. Note that the roots of (24) appear in conjugate reciprocal pairs, as in the conventional root-MUSIC [25]. To find the D central elevation DOAs, we select the D roots k i ( i = 1 , 2 , D ) inside the unit circle that maximize (23). Finally, the central elevation DOA estimates are obtained by:
γ ^ i = a r cos ( λ 2 π d arg ( k i ) ) .

3.2.2. Central Azimuth DOA Estimation

The method is similar to that of the central elevation DOA estimate, thus we simplify the process of deduction. For the centrosymmetric sub-array Y a , we also have the generalized shift invariance relation:
Π M y b y ( μ i ) = Ω ( θ i , φ i ) b y ( μ i ) ,
where Ω ( θ i , φ i ) is a M y × M y diagonal matrix which is given by:
Ω ( θ i , φ i ) = d i a g [ e j η ( M y 1 ) sin θ i sin γ i , e j η ( M y 3 ) sin θ i sin γ i ,   , e j η ( M y 1 ) sin θ i sin γ i ] .
If we define the complex variable l = e j η sin θ i sin γ i , Ω ( θ i , φ i ) can be written as: Ω ( l ) = d i a g [ l ( M y 1 ) , l ( M y 3 ) , , l ( M y 1 ) ] .
Let U y be the signal subspace matrix, whose columns are the eigenvectors corresponding to the D largest eigenvalues of R y = E { y ( t ) y H ( t ) } . Similarly, there exists a unique non-singular D × D matrix T 2 such that U y = B y ( μ ) T 2 . Thus, we introduce a matrix F y ( l ) which is expressed as:
F y ( l ) = Π M y U y Ω ( l ) U y = Π M y B y ( μ ) T 3 Ω ( l ) B y ( μ ) T 2 = [ ( Ω ( l 1 ) Ω ( l ) ) b y ( μ 1 ) , ( Ω ( l 2 ) Ω ( l ) ) b y ( μ 2 ) ,   , ( Ω ( l D ) Ω ( l ) ) b y ( μ D ) ] T 2 .
The central azimuth DOA estimations θ ^ i ( i = 1 , 2 , D ) can be obtained by rooting the following polynomial:
h y ( l ) = det { F y H ( l ) F y ( l ) } .
Similarly, we use the roots inside the unit circle, and select the D roots l i ( i = 1 , 2 , D ) that maximize the spectral function such as:
H y ( l ) = 1 / det ( F y H ( l ) F y ( l ) ) .
The values of sin θ ^ i sin γ ^ i for i = 1 , 2 , , D are obtained as:
sin θ ^ i sin γ ^ i = λ 2 π d arg ( l i ) .

3.2.3. The Parameter Matching Method

For only one CD source, the central elevation and azimuth DOAs can be estimated directly using (25) and (31). However, when multiple CD sources exist, the estimated elevation and azimuth DOAs are required to be matched. To perform the pair-matching procedure, we need to consider the GAM vector of the whole cross array such as:
b x ( μ i ) = [ b y ( μ i ) b z ( μ i ) ] .
Let J be an ( M y + M z ) × ( M y + M z ) selection matrix which is defined as:
J = [ Π M y 0 M y × M Z 0 M z × M y Π M Z ] .
Based on the symmetric property of the ASDW vector in (16), we have the following generalized rotational invariance relation:
J b x ( μ i ) = Φ ( θ i , γ i ) b x ( μ i ) ,
where Φ ( θ i , γ i ) is an ( M y + M z ) × ( M y + M z ) diagonal matrix given by:
Φ ( θ i , γ i ) = d i a g [ e j η ( M y 1 ) sin θ i sin γ i , e j η ( M y 3 ) sin θ i sin γ i , , e j η ( M y 1 ) sin θ i sin γ i , e j η ( M z 1 ) cos γ i , e j η ( M z 3 ) cos γ i , , e j η ( M z 1 ) cos γ i ] .
Let U x be the signal subspace matrix of R x = E { x ( t ) x H ( t ) } . Similarly, there exists a unique non-singular D × D matrix T 3 such that U x = B x ( μ ) T 3 . We can also introduce a matrix F x ( θ , γ ) which is written as:
F x ( θ , γ ) = J U x Φ ( θ , γ ) U x = J B x ( μ ) T 2 Φ ( θ , γ ) B x ( μ ) T 3 = [ ( Φ ( θ 1 , γ 1 ) Φ ( θ , γ ) ) b x ( μ 1 ) , ( Φ ( θ 2 , γ 2 ) Φ ( θ , γ ) ) b x ( μ 2 ) , , ( Φ ( θ D , γ D ) Φ ( θ , γ ) ) b x ( μ D ) ] T 3 .
Similarly, when Φ ( θ i , γ i ) = Φ ( θ , γ ) for i = 1 , 2 , , D , the i-th column of F x ( θ , γ ) is a zero vector. Therefore, the central elevation and azimuth DOA estimations can be matched by minimizing of the following cost function:
f ( θ , γ ) = det ( F x H ( θ , γ ) F x ( θ , γ ) ) .
If we pick γ ^ i from the elevation DOA estimations { γ ^ 1 , γ ^ 2 , γ ^ D } , there will be D pairs of 2D central DOAs for γ ^ i , which is given by { ( θ ^ i 1 , γ ^ i ) , ( θ ^ i 2 , γ ^ i ) , ( θ ^ i D , γ ^ i ) } . We then substitute the DOA estimations into (37) and calculate the function value f ( θ , γ ) . If f ( θ ^ i j , γ ^ i ) for j = 1 , 2 , D is the largest, then ( θ ^ i j , γ ^ i ) is the correct match.

3.2.4. Algorithm Implementation and Complexity Analysis

Now, we can summarize the proposed algorithm as follows:
Step 1:
Calculate the covariance matrix R z = E { z ( t ) z H ( t ) } . Through the eigen-decomposition of R z , obtain the signal subspace matrix U z .
Step 2:
Construct the matrix F z ( k ) in (22), and root the polynomial in (24) to obtain the central elevation DOA estimations γ ^ i for i = 1 , 2 , D . It is noted that the roots are inside a unit circle and maximize (23).
Step 3:
Calculate the covariance matrix R y = E { y ( t ) y H ( t ) } . Through the eigen-decomposition of R y , obtain the signal subspace matrix U y .
Step 4:
Construct the matrix F y ( k ) in (28), and root the polynomial in (29) to obtain sin θ ^ i sin γ ^ i for i = 1 , 2 , D . It is noted that the roots are inside a unit circle and maximize (30).
Step 5:
Compute all the possible 2D DOAs ( θ ^ i j , γ ^ i ) for the elevation DOA estimations γ ^ i . Calculate the function values f ( θ ^ i j , γ ^ i ) for j = 1 , 2 , D in (37). The largest one is the correct match.
Step 6:
Repeat the process in Step 5 to match all the parameters.
Next, when the number of sensor elements M and the number of snapshots L change, we analyze the computational complexity of the proposed algorithm in comparison with the SOS algorithm in [18], the CC algorithm in [21] and Zheng’s algorithm in [20]. The main computational cost of the proposed algorithm is mostly made of four operations: the estimation of the covariance matrix, the eigen-decomposition of the covariance matrix, the polynomials rooting, and the pair-matching procedure. Specifically, the cost involved by the estimation of covariance matrices R x , R y and R z is found to be in O ( 6 M 2 L ) . The eigen-decomposition of the covariance matrices R x , R y and R z needs O ( 10 M 3 ) multiplications. The computational complexity of rooting polynomials h z ( k ) and h y ( k ) is found to be in O ( 2 M 3 ) , and the pair-matching procedure adds O ( D 5 + D 4 M ) multiplications to the proposed algorithm. In above, the computational complexity of the proposed algorithm is O ( 12 M 3 + 6 M 2 L + D 5 + D 4 M ) . Moreover, the main computational complexity of the SOS algorithm, the CC algorithm and Zheng’s algorithm is given in Table 1 (N denotes the number of searching points).
When the number of searching points N is large, it is clear to see that the propose algorithm provides lower computational cost compared to the SOS algorithm. Although the computational complexity of the proposed algorithm is higher than Zheng’s algorithm and the CC algorithm, it is not significant increment since the proposed algorithm does not require any spectrum searching. In addition, and unlike the SOS algorithm, Zheng’s algorithm and the CC algorithm, the proposed algorithm does not use the Taylor series approximation to establish the rotational invariance relation, as this approximation may introduce additional errors.

4. Simulation Results and Performance Analysis

In the following experiments, noise is a complex Gaussian process with zero mean. The number of snapshots is 200. We use the root mean squared error (RMSE) to evaluate the estimation performance, where the RMSE of the central azimuth and elevation DOAs ( RMSE ( θ ) and RMSE ( γ ) ) are defined as:
RMSE ( θ ) = E [ 1 D i = 1 D ( θ ^ i - θ i ) 2 ] ,
RMSE ( γ ) = E [ 1 D i = 1 D ( γ ^ i - γ i ) 2 ] ,
where θ ^ i and θ i are the estimated and true central azimuth DOA of the i-th source, respectively. Additionally, γ ^ i and γ i are the estimated and true central elevation DOA of the i-th source, respectively.
In the following simulations, the signal power of sources is assumed to be the same. In addition, signal-to-noise ratio (SNR) is defined as:
SNR = 10 log σ s 2 σ n 2 ,
where σ s 2 is the signal power of sources, while σ n 2 is the variance of noise.

4.1. Effect of Different Deterministic Angular Distributed Functions

In this part, we examine if the proposed algorithm works properly for different angular distributed functions. The sub-arrays Y a and Z a are both composed of M y = M z = 16 sensors. The distance between adjacent sensors is 0.5 λ . The parameters of two CD sources are μ 1 = ( 20 ° , 3 ° , 20 ° , 5 ° ) and μ 2 = ( 60 ° , 4 ° , 80 ° , 4 ° ) . Their deterministic angular distributed functions are Gaussian and uniform shaped, respectively. The SNR is 15 dB. For 30 independent trials, the central DOA estimations of CD sources are plotted in Figure 2. It can be seen that the proposed algorithm can give the correct DOA estimations for cases where different CD sources have different deterministic angular distributed functions, or unknown deterministic angular distributed functions.

4.2. Performance Comparison

In this part, we compare the estimation accuracy of the proposed algorithm with the SOS algorithm in [18], the CC algorithm in [21] and Zheng’s algorithm in [20] with respect to SNR from 0 dB to 30 dB. The Cramér-Rao lower bound (CRLB) is also used as a benchmark [26]. The sub-arrays Y a and Z a of the crossed array in proposed algorithm are both composed of M y = M z = 16 sensors. The arrays in the SOS algorithm and the CC algorithm are both composed of 32 sensors. Since the number of antenna elements in Zheng’s algorithm must be odd, we set it to 33. In these algorithms, the distance between adjacent sensors in a sub-array is 0.5 λ , the vertical distance between the two sub-arrays is 0.5 λ , and the radius of the UCA is λ / ( 4 sin ( π / 16 ) ) . The parameters of two Gaussian-shaped CD sources are μ 1 = ( 20 ° , 2 ° , 60 ° , 3 ° ) and μ 2 = ( 15 ° , 3 ° , 70 ° , 2 ° ) , respectively. Based on 500 Monte Carlo experiments, the RMSE curves of the central DOA estimations versus SNR are shown in Figure 3. It is clearly indicated that the estimation accuracy of the proposed algorithm is higher than the SOS algorithm, the CC algorithm and Zheng’s algorithm, which arises from the fact that the proposed algorithm does not suffer additional errors brought about by Taylor series approximation.

4.3. Effect of Snapshots

In this part, we illustrate the influence of the number of snapshots on the performance of the proposed algorithm. The number of snapshots varies from 100 to 900. The SNR is fixed to 15 dB and the other parameters are the same as in Section 4.2. Based on 500 Monte Carlo experiments, the RMSE curves for different algorithms are shown in Figure 4, from which we can draw similar conclusion as in Section 4.2.

4.4. Effect of the Central Elevation and Azimuth DOAs

In the last example, we examine the performance of the proposed method versus the central elevation and azimuth DOAs for a Gaussian-shaped CD source with σ θ i = σ γ i = 1 ° . First, let us consider in Figure 5a the influence of the central azimuth DOA on performance. Assume that γ i = 30 ° , SNR = 10 dB and the number of snapshots L = 200 . As can be seen from the Figure 5a, the RMSE of θ ^ i estimated by the proposed method increases dramatically when the central azimuth DOA approaches the lower bound ( θ i = 90 ° ) or the upper bound ( θ i = 90 ° ), but our method can still estimate effectively the central azimuth DOA. Next, the influence of the central elevation DOA on performance is examined in Figure 5b. At this time, we assume that θ i = 30 ° , SNR = 10 dB and L = 200 . Similarly, the RMSE of γ ^ i estimated by the proposed method also increases dramatically when the central elevation DOA approaches the lower bound ( γ i = 0 ° ), and our method has still satisfying estimation performance for approaching the lower bound.

5. Conclusions

In this paper, we have presented a new approach for estimating the 2D central DOA of CD sources using a centrosymmetric crossed array. Instead of using the Taylor series approximation, we derive the symmetric property of an ASDW vector in a centrosymmetric array, based on which the generalized shift invariance relations inside the GAMs are established in the two sub-ULAs. Resorting to such relations, the central elevation and azimuth DOAs are estimated based on the polynomial-root-based search-free method, respectively. After that, the pair-matching method is presented. The proposed algorithm performs independently of the deterministic angular distributed function. Compared to the existing low-complexity algorithms, the proposed algorithm does not suffer additional errors brought by the Taylor series approximation, which allows it to achieve a higher estimation accuracy. Moreover, the proposed algorithm does not suffer high computational complexity brought by spectrum-peak searching.

Acknowledgments

The work was supported by the National Natural Science Foundation of China (grant no. 61401513).

Author Contributions

Zhengliang Dai has contributed to the scientific part of this work. Zhengliang Dai, Weijia Cui and Bin Ba wrote and corrected the manuscript. Daming Wang contributed in the conception and design of the study. Youming Sun has critically reviewed the paper. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The most commonly used the deterministic angular distributed functions for a 2D CD source are presented as follows:
Gaussian shaped:
ρ ( θ ˜ , γ ˜ ; μ i ) = 1 2 π σ θ σ γ e 1 2 ( θ ˜ 2 / σ θ 2 + γ ˜ 2 / σ γ 2 ) ,
Uniform shaped:
ρ ( θ ˜ , γ ˜ ; μ i ) = { 1 12 σ θ σ γ | θ ˜ | 3 σ θ   a n d   | γ ˜ | 3 σ γ 0 otherwise ,
Laplacian shaped:
ρ ( θ ˜ , γ ˜ ; μ i ) = 1 2 σ θ σ γ e ( 2 | θ ˜ | / σ θ + 2 | γ ˜ | / σ γ ) ,
where the small deviation θ ˜ and γ ˜ are defined as θ ˜ = θ θ i and γ ˜ = γ γ i , respectively.

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. Mccloud, M.L.; Scharf, L.L. A new subspace identification algorithm for high-resolution DOA estimation. IEEE Trans. Antennas Propag. 2002, 50, 1382–1390. [Google Scholar] [CrossRef]
  3. Ciuonzo, D.; Romano, G.; Solimene, R. Performance analysis of time-reversal MUSIC. IEEE Trans. Signal Process. 2015, 63, 2650–2662. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Xu, X.; Sheikh, Y.A.; Ye, Z. A rank-reduction based 2-D DOA estimation algorithm for three parallel uniform linear arrays. Signal Process. 2016, 120, 305–310. [Google Scholar] [CrossRef]
  5. Yan, F.; Jin, M.; Qiao, X. Low-complexity DOA estimation based on compressed MUSIC and its performance analysis. IEEE Trans. Signal Process. 2013, 61, 1915–1930. [Google Scholar] [CrossRef]
  6. Vallet, P.; Mestre, X.; Loubaton, P. Performance analysis of an improved MUSIC DOA estimator. IEEE Trans. Signal Process. 2015, 63, 6407–6422. [Google Scholar] [CrossRef]
  7. Zhang, X.; Xu, D. Low-complexity ESPRIT-based DOA estimation for colocated MIMO radar using reduced-dimension transformation. Electron. Lett. 2011, 47, 283–284. [Google Scholar] [CrossRef]
  8. Suleiman, W.; Pesavento, M.; Zoubir, A.M. Performance analysis of the decentralized eigendecomposition and ESPRIT algorithm. IEEE Trans. Signal Process. 2016, 64, 2375–2386. [Google Scholar] [CrossRef]
  9. Besson, O.; Stoica, P. Decoupled estimation of DOA and angular spread for spatially distributed sources. Conference Record of the Thirty-Third Asilomar Conference on Signals, Systems, and Computers, Pacific Grove, CA, USA, 24–27 October 1999; pp. 253–257. [Google Scholar]
  10. Cao, R.; Gao, F.; Zhang, X. An angular parameter estimation method for incoherently distributed sources via generalized shift invariance. IEEE Trans. Signal Process. 2016, 64, 4493–4503. [Google Scholar] [CrossRef]
  11. Valaee, S.; Champagne, B.; Kabal, P. Parametric localization of distributed sources. IEEE Trans. Signal Process. 1995, 43, 2144–2153. [Google Scholar] [CrossRef]
  12. Meng, Y.; Stoica, P. Estimation of direction of arrival of spatially dispersed signals in array processing. IEEE Proc. Sonar Navig. 1996, 143, 1–9. [Google Scholar] [CrossRef]
  13. Wu, Q.; Wong, K.M.; Meng, Y.; Read, W. DOA estimation of point and scattered sources—Vec-Music. Processing of the IEEE Seventh SP Workshop on Statistical Signal and Array Processing, Quebec, QC, Canada, 26–29 June 1994; pp. 365–368. [Google Scholar]
  14. Besson, O.; Stoica, P. Decoupled estimation of DOA and angular spread for a spatially distributed source. IEEE Trans. Signal Process. 2008, 48, 1872–1882. [Google Scholar] [CrossRef]
  15. Bengtsson, M.; Ottersten, B. Low-complexity estimators for distributed sources. IEEE Trans. Signal Process. 2000, 48, 2185–2194. [Google Scholar] [CrossRef]
  16. Shahbazpanahi, S.; Valaee, S.; Bastani, M.H. Distributed source localization using ESPRIT algorithm. IEEE Trans. Signal Process. 2001, 49, 2169–2178. [Google Scholar] [CrossRef]
  17. Boujemâa, H. Extension of COMET algorithm to multiple diffuse source localization in azimuth and elevation. Eur. Trans. Telecommun. 2005, 16, 557–566. [Google Scholar] [CrossRef]
  18. Lee, J.; Song, I.; Kwon, H.; Lee, S.R. Low-complexity estimation of 2D DOA for coherently distributed sources. Signal Process. 2003, 83, 1789–1802. [Google Scholar] [CrossRef]
  19. Guo, X.S.; Wan, Q.; Yang, W.L.; Lei, X. Low-complexity 2D coherently distributed sources decoupled DOAs estimation method. Sci. Chin. Inf. Sci. 2009, 52, 835–842. [Google Scholar] [CrossRef]
  20. Zheng, Z.; Li, G.; Teng, Y. Simplified Estimation of 2D DOA for Coherently Distributed Sources. Wirel. Pers. Commun. 2012, 53, 907–922. [Google Scholar] [CrossRef]
  21. Yang, X.; Li, G.; Zheng, Z.; Zhong, L. Low-complexity 2d central DOA estimation of coherently distributed sources with cross-correlation matrix. Electron. Lett. 2014, 50, 1118–1120. [Google Scholar] [CrossRef]
  22. Liu, S.; Yang, L.; Wu, D.; Huang, J. Two-dimensional DOA estimation using a co-prime symmetric cross array. Prog. Electromagn. Res. C 2014, 54, 67–74. [Google Scholar] [CrossRef]
  23. Shi, Y.; Huang, L.; Qian, C.; So, H. Direction-of-arrival estimation for noncircular sources via structured least squares–based esprit using three-axis crossed array. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 1267–1278. [Google Scholar] [CrossRef]
  24. Wang, X.; Mao, X.; Wang, Y.; Zhang, N.; Li, B. A novel 2-D coherent DOA estimation method based on dimension reduction sparse reconstruction for orthogonal arrays. Sensors 2016, 16, 1496. [Google Scholar] [CrossRef] [PubMed]
  25. Yan, F.; Shen, Y.; Jin, M.; Qiao, X. Computationally efficient direction finding using polynomial rooting with reduced-order and real-valued computations. J. Syst. Eng. Electron. 2016, 27, 739–745. [Google Scholar]
  26. Ghogho, M.; Besson, O.; Swami, A. Estimation of directions of arrival of multiple scattered sources. IEEE Trans. Signal Process. 2001, 49, 2467–2480. [Google Scholar] [CrossRef]
Figure 1. Geometry of the considered crossed array.
Figure 1. Geometry of the considered crossed array.
Sensors 17 01300 g001
Figure 2. The 2D central DOA estimation results of the proposed algorithm (30 trials).
Figure 2. The 2D central DOA estimation results of the proposed algorithm (30 trials).
Sensors 17 01300 g002
Figure 3. (a) RMSE curves of the central azimuth DOA estimations versus SNR; (b) RMSE curves of the central elevation DOA estimations versus SNR.
Figure 3. (a) RMSE curves of the central azimuth DOA estimations versus SNR; (b) RMSE curves of the central elevation DOA estimations versus SNR.
Sensors 17 01300 g003
Figure 4. (a) RMSE curves of the central azimuth DOA estimations versus the number of snapshots; (b) RMSE curves of the central elevation DOA estimations versus the number of snapshots.
Figure 4. (a) RMSE curves of the central azimuth DOA estimations versus the number of snapshots; (b) RMSE curves of the central elevation DOA estimations versus the number of snapshots.
Sensors 17 01300 g004
Figure 5. (a) RMSE curves versus the central azimuth DOA; (b) RMSE curves versus the central elevation DOA.
Figure 5. (a) RMSE curves versus the central azimuth DOA; (b) RMSE curves versus the central elevation DOA.
Sensors 17 01300 g005
Table 1. Comparison of different algorithms in computational complexity.
Table 1. Comparison of different algorithms in computational complexity.
AlgorithmMain Computational Complexity
Proposed algorithm O ( 12 M 3 + 6 M 2 L + D 5 + D 4 M )
SOS algorithm O ( 8 M 3 + 4 M 2 L + N ( D 3 + 3 D 2 M ) )
CC algorithm O ( M 3 + M 2 L + D 3 )
Zheng’s algorithm O ( ( 2 M + 1 ) 3 + ( 2 M + 1 ) 2 L + D 3 )

Share and Cite

MDPI and ACS Style

Dai, Z.; Cui, W.; Ba, B.; Wang, D.; Sun, Y. Two-Dimensional DOA Estimation for Coherently Distributed Sources with Symmetric Properties in Crossed Arrays. Sensors 2017, 17, 1300. https://doi.org/10.3390/s17061300

AMA Style

Dai Z, Cui W, Ba B, Wang D, Sun Y. Two-Dimensional DOA Estimation for Coherently Distributed Sources with Symmetric Properties in Crossed Arrays. Sensors. 2017; 17(6):1300. https://doi.org/10.3390/s17061300

Chicago/Turabian Style

Dai, Zhengliang, Weijia Cui, Bin Ba, Daming Wang, and Youming Sun. 2017. "Two-Dimensional DOA Estimation for Coherently Distributed Sources with Symmetric Properties in Crossed Arrays" Sensors 17, no. 6: 1300. https://doi.org/10.3390/s17061300

APA Style

Dai, Z., Cui, W., Ba, B., Wang, D., & Sun, Y. (2017). Two-Dimensional DOA Estimation for Coherently Distributed Sources with Symmetric Properties in Crossed Arrays. Sensors, 17(6), 1300. https://doi.org/10.3390/s17061300

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