Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
The Electrostatically Formed Nanowire: A Novel Platform for Gas-Sensing Applications
Previous Article in Journal
A Wireless Sensor System for Real-Time Monitoring and Fault Detection of Motor Arrays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computationally Efficient 2D DOA Estimation with Uniform Rectangular Array in Low-Grazing Angle

1
Air and Missile Defense College, Air Force Engineering University, Xi’an 710051, China
2
Electronic information college, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
3
College of Information Science and Engineering, Shandong Agricultural University, Tai’an 271018, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(3), 470; https://doi.org/10.3390/s17030470
Submission received: 13 December 2016 / Revised: 7 February 2017 / Accepted: 22 February 2017 / Published: 26 February 2017
(This article belongs to the Section Physical Sensors)

Abstract

:
In this paper, we propose a computationally efficient spatial differencing matrix set (SDMS) method for two-dimensional direction of arrival (2D DOA) estimation with uniform rectangular arrays (URAs) in a low-grazing angle (LGA) condition. By rearranging the auto-correlation and cross-correlation matrices in turn among different subarrays, the SDMS method can estimate the two parameters independently with one-dimensional (1D) subspace-based estimation techniques, where we only perform difference for auto-correlation matrices and the cross-correlation matrices are kept completely. Then, the pair-matching of two parameters is achieved by extracting the diagonal elements of URA. Thus, the proposed method can decrease the computational complexity, suppress the effect of additive noise and also have little information loss. Simulation results show that, in LGA, compared to other methods, the proposed methods can achieve performance improvement in the white or colored noise conditions.

1. Introduction

Two-dimensional (2D) direction of arrival (DOA) estimation including azimuth and elevation angles with different array geometry has been widely applied in wireless communications, radar and sonar signal processing [1]. Recently, the 2D DOA estimation with uniform rectangular arrays (URAs) has attracted widespread concern [2,3,4,5]. Various algorithms have been developed for improving the estimation performance, such as the multiple signal classification (MUSIC) [6], estimation of signal parameters via rotational invariance techniques (ESPRIT) [7] and the matrix pencil (MP) method [8]. However, in a low-grazing angle (LGA) condition [9,10,11,12,13], the coherency between direct and reflected signals of each target can cause the sample covariance matrix to be rank-deficient. To address this problem, both spatial smoothing [14,15,16,17] and spatial differencing [18,19,20] techniques are developed for 2D DOA estimation with URAs.
For spatial smoothing techniques, Yeh et al. [15] developed the spatial smoothing 2D MUSIC algorithm by using the covariance matrices of overlapping rectangular subarrays and Chen et al. [16] presented an analysis of a special smoothing scheme extended in conjunction with the eigenstructure technique. To reduce the computational complexity, a partial spectral search based method (PSS) is proposed for limiting the searching region into a small sector [21]. Then, a tree structure one-dimensional (1D) algorithm [22] was developed based on a repeated use of the 1D MUSIC algorithm. However, it requires a large number of eigenvalue decompositions (EVDs) and does not perform well under a low signal-to-noise ratio (SNR). For spatial differencing techniques, Liu et al. [18] constructed a new spatial differencing matrix to suppress the white noise by using the difference between the first subarray and the spatial backward subarray. By exploiting the difference between the neighboring forward subarrays and backward subarrays, the method in [19] can suppress the colored noise. Additionally, an asymmetric spatial difference smoothing method [20] was used to reduce noise for coherent sources location, especially when the number of targets is odd.
However, due to the 2D peaks search, 2D or 1D EVD operations, all the aforementioned methods suffer from great computational complexity, especially for a large size of subarrays. In other words, these methods can provide better performance at the cost of great computation. Therefore, we concentrate on reducing the computational complexity caused by the 2D operations, while maintaining a high estimation performance. In this paper, we propose a spatial differencing matrix set (SDMS) method for 2D DOA estimation with URA in LGA. Employing the overlapping column or row subarrays along the x- or y-direction, we build the SDMSs by rearranging the auto-correlation and cross-correlation matrices in turn among different subarrays. In addition, to suppress the data loss, we only perform difference for the auto-correlations and the cross-correlations are kept completely. Then, the two parameters are estimated independently by using the 1D subspace-based estimation technique, the pair-matching of which is achieved by extracting the diagonal elements of the URA. Simulation results verify the effectiveness of the proposed method.
The advantages of the proposed method can be given as follows:
  • The methods in [3,4,15,16,17,21] involve the 2D EVD or 2D peak search, while the proposed method can estimate the parameters with 1D subspace-based estimation techniques.
  • The method in [22] can only use the auto-correlations of different subarrays, while the proposed method can use more information including auto-correlations and cross-correlations.
  • The spatial differencing techniques in [18,19] perform a difference operation on the whole subarrays, while the proposed method is only for the auto-correlations and the cross-correlations are kept completely. Thus, the SDMS method has little data loss.
The rest of this paper is given as follows. In Section 2, the basic received signal model of the URA in LGA is developed. Then, we derive the SDMS method using a 1D subspace-based estimation technique and achieve pair-matching by extracting the diagonal elements of the URA in Section 3. Simulation results are given in Section 4, and we conclude the whole paper in Section 5.
In this paper, operators ( · ) T , ( · ) * and ( · ) H represent transpose, conjugation, and conjugate transpose, respectively. I N denotes an N × N identity matrix and J M denotes an M × M exchange matrix with ones on its anti-diagonal and zeros elsewhere. ⊕ and ⊗ represent Hadamard product and Kronecker product, respectively; d i a g ( · ) and b l k d i a g ( · ) denote the diagonal matrix or the block diagonal matrix operator. E [ · ] and v e c ( · ) denote expectation and vectorization, respectively.

2. System Model

As described in Figure 1, we regard the multipath effect as ideal specular reflection, where both curved earth and atmosphere refraction are not considered. We also consider K narrowband far-field signals s k ( t ) ( k = 1 , 2 , , K ) impinging on a URA with MN well calibrated and identically polarized sensors parallel to the xoy plane. Here, both x- and y- directions of the URA are separated by half a wavelength, the height of which is set as h. Since the received signals include two paths, i.e., direct path, reflected path, the output can be given as [11,12]
X ( t ) = k = 1 K a x ( α k , θ d k ) a y T ( α k , θ d k ) s k ( t ) + k = 1 K a x ( α k , θ r k ) a y T ( α k , θ r k ) β k s k ( t ) + Z ( t ) ,
where θ d k and θ r k are the direct and reflected elevation angle for the kth target ( θ d k - θ r k = θ k ), respectively, α k is the azimuth angle and β k is the multipath reflection coefficient from the kth signal to receiver array; Let β k = exp [ j ( π 2 π Δ R k 2 π Δ R k λ λ ) ] , Δ R k 2 h sin θ k for simplicity, Δ R k is the difference value between direct path and reflected path in LGA; a x ( α k , θ k ) = a x ( u k ) = [ 1 , e j π u k , , e j π ( M 1 ) u k ] T , a y ( α k , θ k ) = a y ( v k ) = [ 1 , e j π v k , , e j π ( N 1 ) v k ] T , u k = sin θ k cos α k , and v k = sin θ k sin α k . The elements of Z ( t ) are temporally and spatially complex white Gaussian noises with zero-mean and variance σ 2 . Then, by vectorizing the matrix X ( t ) , a composite data vector can be constructed as
x ( t ) = v e c ( X ( t ) ) = ( A x A y ) s ( t ) + z ( t ) ,
where A x = [ a x ( u 1 ) , a x ( u 1 ) , , a x ( u K ) , a x ( u K ) ] T , A y = [ a y ( v 1 ) , a y ( v 1 ) , , a y ( v K ) , a y ( v K ) ] T , s ( t ) = [ s 1 ( t ) , β 1 s 1 ( t ) , , s K ( t ) , β K s K ( t ) ] 2 K × 1 T , z ( t ) = v e c ( Z ( t ) ) .
By forming Φ y = d i a g [ e j π v 1 , e j π v 1 , , e j π v K , e j π v K ] , the model in (2) can be rewritten as
x ( t ) = x 1 x 2 x N = A x A x Φ y A x Φ y N 1 s ( t ) + z ( t ) .
From (2) and (3), the incident signals can be divided into K parts and the signals in each part are correlated. With L snapshots ( t = 1 , 2 , , L ), the sample covariance matrix can be calculated as
R 0 = E x ( t ) x H ( t ) = 1 L t = 1 L x ( t ) x H ( t ) = A R s A H + σ 2 I M N ,
where A = A x A y , R s = E [ s ( t ) s H ( t ) ] denotes the correlation matrix of coherent signals. Thus, due to rank-deficiency of the covariance matrix R 0 , the classic methods for 2D DOA estimation will lose efficiency [6,7,8].

3. 2D DOA Estimation in LGA

In this section, we proposed a SDMS method that uses the differencing matrices among different spatial smoothing subarrays for 2D DOA estimation in LGA, where the parameters u k and v k are estimated independently.

3.1. 1D Estimation of Parameter u k

Here, as described in Figure 2a, we divide the URA with M × N sensors into P overlapping forward rectangular subarrays of size Q × N along the x-direction. Each rectangular subarray has N column subarrays with Q sensors, such as the shaded areas for the first rectangular subarray. Then, the nth column of pth rectangular subarray can be given as z p n ( t ) = X p y n ( t ) = X p G n y ( t ) , where X p = [ 0 Q × ( p 1 ) I Q 0 Q × ( P p ) ] , P = M-Q + 1, p = 1 , 2 , , P , n = 1 , 2 , , N , and G n = [ 0 M × ( n 1 ) M I M 0 M × ( N n ) M ] . As a result, using the matrix pencil of auto-correlation and cross-correlation matrices between different column subarrays, we can build the new matrix as
R x p = E [ z p 1 z p 1 H ] , , E [ z p 1 z p N H ] , N m a t r i c e s o f s i z e Q × N E [ z p 2 z p 2 H ] , , E [ z p 2 z p N H ] , N 1 m a t r i c e s o f s i z e Q × N , E [ z p N z p N H ] 1 ,
In (5), because of the calculation of covariance matrices in turn, the data information can be used more effectively. In addition, the cross-correlation matrix E [ z p d z p n H ( t ) ] ( d n ) can also restrain the effect of additive noise. However, the auto-correlation matrix E [ z n p z n p H ( t ) ] has the noise covariance matrix σ 2 I Q , which will certainly decrease the performance. Then, in order to suppress the effect of noise, we first build an initial matrix that has the same noise matrix as R x p , and we have
C x = { E [ z 11 z 11 H ] , 0 , 0 N m a t r i c e s o f s i z e Q × N , E [ z 11 z 11 H ] , 0 , , 0 N 1 m a t r i c e s o f s i z e Q × N , , E [ z 11 z 11 H ] 1 } .
Combining (5) and (6), the forward SDMS for the x-direction (SDMS-x) can be defined as
D x p = C x J Q R x p * F Q = A x Q D ¯ x p d i a g A x Q H , , A x Q H ,
where F Q = b l k d i a g [ J Q , J Q , , J Q ] , the number of J Q in F Q is N ( N + 1 ) / 2 , D ¯ x p = [ R s Θ Q p + 1 R s * Θ 1 + p Q , , Θ Q p + 1 R s * Φ N 1 Θ 1 + p Q , , R s Θ Q p + 1 Φ 1 N R s * Φ N 1 Θ 1 + p Q ] , A x Q is the submatrix of the array response matrix A x consisting of the first Q rows, X p A x = A x Q Θ p 1 , and J Q A x Q * = A x Q Θ Q , Θ = d i a g [ e j π u 1 , e j π u 1 , e j π u K , , e j π u k ] . Following the forward backward (FB) technique, the new FB SDMS-x can be given as
D x = 1 2 P p = 1 P D x p + J Q D x p * F Q .
Then, based on the definition in (8), we can prove that the new FB SDMS-x has the following property.
Theorem 1.
Assume that there are 2K narrowband coherent signals impinging on the URA ( M × N sensors). As described in Figure 2a, the URA is divided along the x-direction and the number of sensors in each column subarray is Q. Then, we define the new FB SDMS-x D x as in (8). If Q 2 K , the rank of D x is equal to the number of the signals, namely, r a n k ( D x ) = 2 K.
Proof. 
See the Appendix A. ☐
Under the Theorem 1, we divide A x Q into two submatrices as A x Q = [ A x Q 1 T , A x Q 2 T ] T , where A x Q 1 and A x Q 2 consist of the first 2K rows and the last Q 2 K rows, respectively. Since the matrix A x Q is a Vandermonde matrix with full rank, we can get a 2 K × ( Q 2 K ) linear operator Γ x for A x Q 2 = Γ x H A x Q 1 . Then, the matrix D x can be divided into two submatrices as
D x = D x 1 D x 2 2 K Q 2 K ,
where D x Q 2 = Γ x H D x Q 1 . In addition, the operator Γ x can be calculated as
Γ x = A x Q 1 H A x Q 2 H = D x 1 D x 1 H 1 D x 1 D x 2 H .
Therefore, by constructing the matrix Ω x = [ Γ x H , I Q 2 K ] H , we can get Ω x H A x Q = 0 ( Q 2 K ) × 2 K , which can be used to estimate the parameter u k . By letting a u k = [ 1 , e j π u k , , e j π ( Q 1 ) u k ] T , the parameter u k can be estimated by minimizing the following cost function
f u k = a H u k Π x a u k ,
where Π x = Ω x Ω x H Ω x 1 Ω x H .

3.2. 1D Estimation of Parameter v k

Similarly, as described in Figure 2b, we divide the URA into F overlapping forward rectangular subarrays along the y-direction. Each rectangular subarray has M row subarrays with Q sensors. Then, the mth row of fth rectangular subarray can be set as w f m ( t ) = X f G m T e y ( t ) , where T e is the row permutation matrix, X f = [ 0 Q × ( f 1 ) I Q 0 Q × ( F f ) ] , f = 1 , 2 , , F , m = 1 , 2 , , M , F=N-Q+1, and G m = [ 0 N × ( m 1 ) N I N 0 N × ( M m ) N ] . Then, we can build the new FB SDMS for the y-direction (SDMS-y) as
D y = 1 2 F f = 1 F D y f + J Q D y f * F Q ,
where
D y f = C y J Q R y f * F Q = A y Q D ¯ y f d i a g A y Q H , , A y Q H ,
R y f = E [ w f 1 w f 1 H ] , , E [ w f 1 z f M H ] , M m a t r i c e s o f s i z e Q × M E [ w f 2 w f 2 H ] , , E [ w f 2 w f M H ] , M 1 m a t r i c e s o f s i z e Q × M , E [ w f M w f M H ] 1 ,
C y = { E [ w 11 w 11 H ] , 0 , 0 M m a t r i c e s o f s i z e Q × M , E [ w 11 w 11 H ] , 0 , , 0 M 1 m a t r i c e s o f s i z e Q × M , , E [ w 11 w 11 H ] 1 } ,
and A y Q is the submatrix of the array response matrix A y consisting of the first Q rows, X f A y = A y Q Φ f 1 , J Q A y Q * = A y Q Θ Q , D y f = [ R s Φ Q f + 1 R s * Φ 1 + f Q , Φ Q f + 1 R s * Θ M 1 Φ 1 + f Q , , R s Φ Q f + 1 Θ 1 M R s * Θ M 1 Φ 1 + f Q ] . Likewise, the new FB SDMS-y in (12) has the following property.
Theorem 2.
Assume that there are 2K narrowband coherent signals impinging on the URA ( M × N sensors). As described in Figure 2b, the number of sensors in each row subarray is Q, and the new FB SDMS-y is defined in (12). If Q 2 K , then the rank of D y is equal to the number of the signals, namely, r a n k ( D y ) = 2 K .
Based on the Theorem 2, we can divide the matrix D y as
D y = D y 1 D y 2 2 K Q 2 K ,
where D y Q 2 = Γ y H D y Q 1 and Γ y is the linear operator. By letting a v k = [ 1 , e j π v k , , e j π ( Q 1 ) v k ] T , the parameter v k can be estimated by minimizing the following cost function
f v k = a H v k Π y a v k ,
where Π y = Ω y Ω y H Ω y 1 Ω y H , Ω y = [ Γ y H , I Q 2 K ] H , Γ y = A y Q 1 H A y Q 2 H = D y 1 D y 1 H 1 D y 1 D y 2 H .

3.3. Pair-Matching of Parameters u k and v k

Since the estimated parameters u k and v k are calculated independently, the pair-matching is very important for multiple targets. In the case of M<N, by extracting the diagonal elements of URA, we can write
r d i a g ( t ) = A v k A u k · s ( t ) + N 0 ,
where A v k = A x , A u k is the submatrix of the array response matrix A y consisting of the first M rows, and N 0 is the subarray of N . Since both A v k and A u k are Vandermode matrices, A v k A u k is also a Vandermode matrix. Using the FBSS method, we can divide r d i a g ( t ) into P overlapping forward subarrays with Q sensors. Then, the smoothing matrix can be given as
D d i a g = 1 P p = 1 P R d i a g p + J Q ( R d i a g p ) * J Q ,
where R d i a g p = X p R d i a g X p H and R d i a g = E [ r d i a g ( t ) r d i a g H ( t ) ] . Hence, we can get the estimated parameters by minimizing the following cost function
f v k , u k = a H v k , u k Π d i a g a v k , u k ,
where a v k , u k = [ 1 , e j π v k e j π u k , , e j π ( Q 1 ) v k e j π ( Q 1 ) u k ] T , Π d i a g = Ω d Ω d H Ω d 1 Ω d H , Ω d = [ Γ d H , I Q 2 K ] H , Γ d = A d Q 1 H A d Q 2 H = D d 1 D d 1 H 1 D d 1 D d 2 H , and D d 1 and D d 2 consist of the first 2K rows and the last Q-2K rows of D d i a g , respectively. From (20), the parameters u k and v k can be matched by repeating the following minimization for i = 1 , , 2 K and j = 1 , , i 1
i , k i = arg min k i f v i , u k i , s u b j e c t t o k i k j ,
and the constraint condition k i k j can avoid the different u k paired with the same v k . Then, the azimuth angle and elevation angle can be estimated as
α k = tan 1 u k u k v k v k , θ k = sin 1 v k 2 + u k 2 .

3.4. Implementation of the Proposed Method

From (11), (17) and (22), the azimuth and elevation angles can be estimated by using the spatial differencing method, where the computational burdensome EVD is avoided and the effect of additive noise is also suppressed. Then, with the finite array data y ( t ) t = 1 L , the proposed method can be implemented as follows:
  • Calculate the estimated sample covariance matrix R ^ in (4) as
    R ^ 0 = 1 1 L L t = 1 L y ( t ) y H ( t ) ,
                   L M 2 N 2 flops.
  • Form the FB SDMS-x D ^ x in (8) and the FB SDMS-y D ^ y in (12) as
    D ^ x = 1 2 P p = 1 P D ^ x p + J Q D ^ x p * F Q , D ^ y = 1 2 F f = 1 F D ^ y f + J Q D ^ y f * F Q ,
    where D ^ x p and D ^ y f can be calculated by using the covariance matrix R ^ 0
                     2 Q 3 flops.
  • Estimate the orthogonal projectors Π ^ x in Section 3.1 and Π ^ y in Section 3.2.
             2 N ( N + 1 ) Q 2 K 2 + 2 M ( M + 1 ) Q 2 K 2 + 2 [ 3 Q 2 ( Q 2 K ) + O ( 2 K ) 3 ] flops.
  • Estimate the parameters u k and v k by finding the phases of the p zeros of the polynomial p u ( z ) and p v ( z ) using (11) and (17), where p u ( z ) = Δ z Q 1 a H z Π ^ x a z and p v ( z ) = Δ z Q 1 a H z Π ^ y a z , z = Δ e j π u k or z = Δ e j π v k
             2 [ ( Q 1 ) 2 + O ( Q 1 ) 3 ] flops.
  • Perform the pair-matching of the parameters v k and u k by using (18)–(21) and estimate the azimuth and elevation angles by (22)
             2 Q 2 ( 2 K ) 2 + 3 Q 2 ( Q 2 K ) + O ( 2 K ) 3 flops.
As shown above, the cost of each step is roughly indicated in terms of the number of MATLAB flops. Furthermore, in the case of Q 2 K , the computational complexity of the proposed method mainly includes the calculation of covariance matrix and SDMSs, the estimation of parameters and the pair-matching, which is about L M 2 N 2 + 8 N 2 Q 2 + 8 M 2 Q 2 + 12 Q 3 + 10 Q 2 .
Remark 1.
The proposed method can estimate the parameters independently, and the null space is calculated by using the linear propagator based on the partition of the array response matrix. The polynomial roots can be obtained by the Linsey–Fox root finding algorithm, which is less than the MATLAB function roots. However, the 2D FBSS-MUSIC method in [15,16], the 2D FBSS based DOA matrix (FBSS-DOAM) method in [23], the conventional spatial differencing (CSD) method in [18,19], and the tree structure one-dimensional (1D) based (TSOD) algorithm in [22] all involve the EVD to obtain the signal subspace or noise subspace. Furthermore, as shown in Table 1, FBSS-MUSIC and CSD both need 2D spectrum peak searching. TSOD needs ( 2 P + 8 K + 1 ) 1D EVD, while FBSS-DOAM needs two times of 1D spectrum peak searching. It is easily seen that the proposed method is computationally more efficient than other recently developed methods.
Remark 2.
As described in (7) and (12), the proposed method can form the row and column FB SDMSs by using the auto-correlation and cross-correlation matrices among different subarrays, where we only perform the difference for the auto-correlations, and the cross-correlations are kept completely. In addition, we also employ the FB technique to improve the estimation performance. However, the TSOD method can only use the auto-correlations of different subarrays and CSD performs the difference on the whole subarrays. Thus, SDMS can achieve performance improvement over the methods in [18,19,22].

3.5. Cramér–Rao Bounds (CRB)

As described in Section 2, according to [24], the CRB can be obtained as
C R B = σ 2 2 L Re D H Π A D R ^ s T 1 ,
where
D = a 1 θ 1 , , a K θ K , a 1 ϕ 1 , , a K ϕ K
and R ^ s = R s R s R s R s , Π A = I M N A ( A H A ) 1 A H , a k is the kth column of A , k = 1, ⋯, 2K.

4. Simulation Results

We now evaluate the estimation performance of SDMS by using some numerical experiments. In LGA, the heights of the URA is set as h = 20 m and the number of sensors is M = N = 20. The colored noise is of a second-order AR model with coefficients a = [ 1 , 0.7 , 0.6 ] [25] and the targets are located at α = [ 10 , 20 , 40 ] , θ = [ 20 , 35 , 50 ] . The wavelength of received signals is set as 1 m and the estimation performance is examined over 300 Monte Carlo runs.
Experiment 1: Effectiveness of proposed method. In this experiment, we mainly examine the effectiveness of SDMS in the presence of white noise and colored noise conditions, where the number of sensors in each subarray is Q = 16 and the total number of snapshots is chosen to be L = 500. Here, the SNR is set as 10 dB. Figure 3 shows the estimation results of the proposed method with 100 Monte Carlo runs. As expected, all the 2D DOAs can be estimated effectively and accurately for the white noise and colored noise conditions.
Experiment 2: Performance versus SNR. In Figure 4 and Figure 5, we evaluate the performance in terms of SNR in the white noise and colored noise conditions, where we assume L = 500 and Q = 16. Three methods are performed for comparison, including the CSD method in [18], the TSOD method in [22] and the proposed method. Moreover, the CRB is also provided. It can be observed that, due to the use of more information and the application of the difference operation, the performance of the proposed method is better than those of methods in [18,22]. In the white noise condition, the performance of CSD is better than TSOD for relatively low SNR, whereas it is the opposite with high SNR. It illustrates that the difference operation can reduce the effect of white noise in the low SNR condition, while the data loss will also decrease the performance in a high SNR condition. In the colored noise condition, since the non-diagonal elements of the noise covariance matrix have significant value, CSD is always superior to TSOD. Comparing SDMS and CSD, SDMS can achieve better performance by using more data information. To sum up, SDMS can achieve great performance improvement in the low SNR condition by using more information and performing the difference operation.
Experiment 3: Performance versus the number of snapshots. In this experiment, we evaluate the estimation performance with different methods (described in Experiment 2) in terms of the number of snapshots, where we assume the SNR is 2 dB and Q = 16. Figure 6 and Figure 7 show the RMSE versus the number of snapshots in the white noise and colored noise conditions, respectively. It is shown that, because of the use of difference operation and cross-correlation matrices, even when the number of snapshots is small, the proposed method still outperforms the methods in [18,22]. Then, for the white noise, the performance of CSD is weaker than TSOD for the small number of snapshots due to the data loss caused by the difference operation between forward and backward smoothed matrices. For colored noise, CSD performs better than TSOD. In addition, by using more data information, the performance of SDMS is also less sensitive to the number of snapshots than that of CSD.
Experiment 4: Performance versus the size of subarrays. Figure 8 describes the performance versus the size of subarrays in the white noise condition. Here, we assume L = 500, S N R = 10 dB. It is easily seen that SDMS still outperforms CSD and TSOD. In addition, due to the increase of array aperture, the performance of those methods becomes better. Then, we can observe that TSOD fails with a small size of subarrays but performs better than CSD with a big size. The reason is that the data loss caused by the difference operation will become larger when the size of subarrays increases. However, the proposed method can use the cross-correlation to reduce the data loss, resulting in a better performance.

5. Conclusions

A new computationally efficient SDMS method with little information loss is proposed to suppress the effect of white noise or colored noise in LGA. The two parameters are estimated independently by using a 1D subspace-based estimation technique, the pair-matching of which is achieved by extracting the diagonal elements of URA. Simulation results show that, in LGA, the performance of the proposed method is superior to the other recently developed method in low SNR conditions and with a small number of snapshots for white noise and colored noise conditions. In the near future, based on improving the information utilization, the extended spatial differencing method with MIMO radar for joint DOA and DOD estimation will be considered carefully.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant No. 61601504) and the Natural Science Foundation Research Project of Shanxi province, China (Grant No. 2016JQ6020).

Author Contributions

Junpeng Shi, Guoping Hu and Xiaofei Zhang participated in the design of this study, and they all performed the statistical analysis and manuscript preparation. Fenggang Sun and Yu Xiao conceived and designed the experiments; Junpeng shi also wrote the paper. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Based on the fact that the rank of a matrix is unchanged by a permutation of its columns [18], the rank of SDMS in (8) can be rewritten as
r a n k D x = r a n k p = 1 P E [ z 11 z 11 H ( t ) ] J Q E [ z p 1 z p 1 H ( t ) ] * J Q .
Then, we can calculate
E [ z 11 z 11 H ( t ) ] J Q E [ z p 1 z p 1 H ( t ) ] * J Q = A x Q R s Θ Q p + 1 R s * Θ 1 + p Q A x Q H .
Since Q 2 K , the rank of D x can be repressed as
r a n k D x = r a n k p = 1 P R s Θ Q p + 1 R s * Θ 1 + p Q .
By using R s = d i a g ρ 1 , ρ 2 , , ρ K and Θ = d i a g Θ 1 , Θ 2 , , Θ K , we have
Λ p = R s Θ Q p + 1 R s * Θ 1 + p Q = b l k d i a g ρ 1 Θ 1 Q p + 1 ρ 1 * Θ 1 1 + p Q , , ρ K Θ K Q p + 1 ρ K * Θ K 1 + p Q = b l k d i a g Λ p 1 , Λ p 2 , Λ p K .
Clearly, we can get r a n k D x = r a n k p = 1 P Λ p = k = 1 K r a n k p = 1 P Λ p k . In the matrix form, we have
p = 1 P Λ p k = Γ k V Γ k H ,
where ρ k = s k s k H , Γ k = s k , Θ k Q s k * , , s k , Θ k Q P + 1 s k * and V = d i a g 1 , 1 , , 1 , 1 . From (A5), we can prove that
r a n k p = 1 P Λ p k = r a n k Γ k = r a n k s k * , Θ k s k * , Θ k P + 1 s k * = r a n k d i a g s k * Ω k = 2 ,
where Ω k = ω k , ω k * and ω k = 1 , e j π u k , , e j π ( P 1 ) u k . Therefore, in the case of Q 2 K , from (A6), we can get that the rank of D t is equal to 2K.

References

  1. Godara, L. Applications of antenna arrays to mobile communications-Part I: Performance improvement, feasibility, and system considerations. Proc. IEEE 1997, 85, 1031–1060. [Google Scholar] [CrossRef]
  2. Heidenreich, P.; Zoubir, A.; Rubsamen, M. Joint 2D DOA estimation and phase calibration for uniform rectangular arrays. IEEE Trans. Signal Process. 2012, 60, 4683–4693. [Google Scholar] [CrossRef]
  3. Lin, J.; Ma, X.; Hao, C.; Jiang, L. Direction of Arrival Estimation of Sparse Rectangular Array via Two-Dimensional Continuous Compressive Sensing. In Proceedings of the Sixth international Conference on Information Science and Technology, Dalian, China, 6–8 May 2016; Volume 5, pp. 539–543.
  4. Zhang, W.; Liu, W.; Wang, J.; Wu, S. Computationally efficient 2D DOA estimation for uniform rectangular arrays. Multidimens. Syst. Signal Process. 2014, 25, 847–857. [Google Scholar] [CrossRef]
  5. Swindlehurst, A.; Kailath, T. Azimuth/elevation direction finding using regular array geometries. IEEE Trans. Aerosp. Electron. Syst. 1993, 29, 145–155. [Google Scholar] [CrossRef]
  6. Schmidt, R.O. Multiple emitter location and signal parameter estimation. IEEE Trans. Antennas Propag. 1986, 34, 276–280. [Google Scholar] [CrossRef]
  7. 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]
  8. Yilmazer, N.; Fernandez-Recio, R.; Sarkar, T.K. Matrix pencil method for simultaneously estimating azimuth and elevation angles of arrival along with the frequency of the incoming signals. Digit. Signal Process. 2006, 16, 796–816. [Google Scholar] [CrossRef]
  9. Sen, S.; Nehorai, A. OFDM MIMO radar with mutual-Information waveform design for low-grazing angle tracking. IEEE Trans. Signal Process. 2010, 58, 3152–3162. [Google Scholar] [CrossRef]
  10. Sen, S.; Nehorai, A. OFDM MIMO radar for low-grazing angle tracking. In Proceedings of the 2009 Conference Record of the Forty-Third Asilomar Conference on Signals, Systems and Computers, Pacific Grove, CA, USA, 1–4 November 2009; Volume 11, pp. 125–129.
  11. Liu, J.; Liu, Z.; Xie, R. Low angle estimation in MIMO radar. Electron. Lett. 2010, 46, 1565–1566. [Google Scholar] [CrossRef]
  12. Shi, J.; Hu, G.; Lei, T. DOA estimation algorithms for low-angle targets with MIMO radar. Electron. Lett. 2016, 52, 652–654. [Google Scholar] [CrossRef]
  13. Shi, J.; Hu, G.; Zong, B.; Chen, M. DOA estimation using multipath echo power for MIMO radar in low-grazing angle. IEEE Sens. J. 2016, 16, 6087–6094. [Google Scholar] [CrossRef]
  14. Pillai, S.U.; Kwon, B.H. Forward/backward spatial smoothing techniques for coherent signal identification. IEEE Trans. Acoust Speech Signal Process. 1989, 37, 8–15. [Google Scholar] [CrossRef]
  15. Yeh, C.; Lee, J.; Chen, Y. Estimating two-dimensional angles of arrival in coherent source environment. IEEE Trans. Acoust Speech Signal Process. 1989, 37, 153–155. [Google Scholar] [CrossRef]
  16. Shan, T. J.; Wax, M.; Kailath, T. On spatial smoothing for direction-of-arrival estimation of coherent signals. IEEE Trans. Acoust Speech Signal Process. 1997, 33, 806–811. [Google Scholar] [CrossRef]
  17. Xu, L.; Ye, Z. Two-dimensional direction of arrival estimation by exploiting the symmetric configuration of uniform rectangular array. IET Radar Sonar Navig. 2012, 6, 307–313. [Google Scholar] [CrossRef]
  18. Liu, F.; Wang, J.; Sun, C.; Du, R. Spatial differencing method for DOA estimation under the coexistence of both uncorrelated and coherent signals. IEEE Trans. Antennas Propag. 2012, 60, 2052–2062. [Google Scholar] [CrossRef]
  19. Ma, X.; Dong, X.; Xie, Y. An improved spatial differencing method for DOA estimation with the coexistence of uncorrelated and coherent signals. IEEE Sens. J. 2016, 16, 3719–3723. [Google Scholar] [CrossRef]
  20. Hong, S.; Wan, X.; Ke, H. An improved spatial differencing method for DOA estimation with the coexistence of uncorrelated and coherent signals. Signal Process. 2015, 109, 69–73. [Google Scholar] [CrossRef]
  21. Wu, Q.; Sun, F.; Lan, P.; Ding, G. Two-dimensional direction-of-arrival estimation for co-prime planar arrays: A partial spectral search approach. IEEE Sens. J. 2016, 16, 5660–5670. [Google Scholar] [CrossRef]
  22. Wang, Y.; Lee, L.; Yang, S.; Chen, J. A tree structure one-dimensional based algorithm for estimating the two-dimensional direction of arrivals and its performance analysis. IEEE Trans. Antennas Propag. 2008, 56, 178–188. [Google Scholar] [CrossRef]
  23. Wang, H.; Liu, K. 2D spatial smoothing for multipath coherent signal separation. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 391–405. [Google Scholar] [CrossRef]
  24. Vu, D.; Renaux, A.; Boyer, R.; Marcos, S. A Cramér-Rao bounds based analysis of 3D antenna array geometries made from ULA branches. IEEE Trans. Aerosp. Electron. Syst. 2013, 24, 121–155. [Google Scholar] [CrossRef] [Green Version]
  25. Jiang, H.; Zhang, J.K.; Wong, K.M. Joint DOD and DOA Estimation for Bistatic MIMO Radar in Unknown Correlated Noise. IEEE Trans. Veh. Technol. 2015, 64, 5113–5125. [Google Scholar] [CrossRef]
Figure 1. The array geometry model for the URA.
Figure 1. The array geometry model for the URA.
Sensors 17 00470 g001
Figure 2. Rectangular subarray grouping of the URA.
Figure 2. Rectangular subarray grouping of the URA.
Sensors 17 00470 g002
Figure 3. The estimated 2D DOAs of SDMS method with 100 Monte Carlo runs.
Figure 3. The estimated 2D DOAs of SDMS method with 100 Monte Carlo runs.
Sensors 17 00470 g003
Figure 4. Performance versus SNR in the white noise condition.
Figure 4. Performance versus SNR in the white noise condition.
Sensors 17 00470 g004
Figure 5. Performance versus SNR in the colored noise condition.
Figure 5. Performance versus SNR in the colored noise condition.
Sensors 17 00470 g005
Figure 6. Performance versus the number of snapshots in the white noise condition.
Figure 6. Performance versus the number of snapshots in the white noise condition.
Sensors 17 00470 g006
Figure 7. Performance versus the number of snapshots in the colored noise condition.
Figure 7. Performance versus the number of snapshots in the colored noise condition.
Sensors 17 00470 g007
Figure 8. Performance versus the size of subarrays in the white noise condition.
Figure 8. Performance versus the size of subarrays in the white noise condition.
Sensors 17 00470 g008
Table 1. Computational complexity comparison.
Table 1. Computational complexity comparison.
FBSS-MUSICFBSS-DOAMCSDTSODProposed Method
EVDone, O ( Q 6 ) two, O ( Q 3 ) one, O ( Q 6 ) ( 2 P + 8 K + 1 ) , O ( Q 3 ) w/o
1D searchingw/otwow/otwow/o
2D searchingonew/oonew/ow/o

Share and Cite

MDPI and ACS Style

Shi, J.; Hu, G.; Zhang, X.; Sun, F.; Xiao, Y. Computationally Efficient 2D DOA Estimation with Uniform Rectangular Array in Low-Grazing Angle. Sensors 2017, 17, 470. https://doi.org/10.3390/s17030470

AMA Style

Shi J, Hu G, Zhang X, Sun F, Xiao Y. Computationally Efficient 2D DOA Estimation with Uniform Rectangular Array in Low-Grazing Angle. Sensors. 2017; 17(3):470. https://doi.org/10.3390/s17030470

Chicago/Turabian Style

Shi, Junpeng, Guoping Hu, Xiaofei Zhang, Fenggang Sun, and Yu Xiao. 2017. "Computationally Efficient 2D DOA Estimation with Uniform Rectangular Array in Low-Grazing Angle" Sensors 17, no. 3: 470. https://doi.org/10.3390/s17030470

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