Next Article in Journal
Combining SDAE Network with Improved DTW Algorithm for Similarity Measure of Ultra-Weak FBG Vibration Responses in Underground Structures
Next Article in Special Issue
Low Complexity Beamspace Super Resolution for DOA Estimation of Linear Array
Previous Article in Journal
Coherent Integration Method Based on Radon-NUFFT for Moving Target Detection Using Frequency Agile Radar
Previous Article in Special Issue
ECCM Schemes against Deception Jamming Using OFDM Radar with Low Global PAPR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

2D-DOD and 2D-DOA Estimation for a Mixture of Circular and Strictly Noncircular Sources Based on L-Shaped MIMO Radar

1
School of Faculty of Information Science and Engineering, Ningbo University, Ningbo 315211, China
2
Key Laboratory of Intelligent Perception and Advanced Control of State Ethnic Affairs Commission, Dalian 116600, China
3
School of Software, Qufu Normal University, Qufu 273165, China
*
Author to whom correspondence should be addressed.
Submission received: 9 March 2020 / Revised: 10 April 2020 / Accepted: 11 April 2020 / Published: 12 April 2020

Abstract

:
In this paper, a joint diagonalization based two dimensional (2D) direction of departure (DOD) and 2D direction of arrival (DOA) estimation method for a mixture of circular and strictly noncircular (NC) sources is proposed based on an L-shaped bistatic multiple input multiple output (MIMO) radar. By making full use of the L-shaped MIMO array structure to obtain an extended virtual array at the receive array, we first combine the received data vector and its conjugated counterpart to construct a new data vector, and then an estimating signal parameter via rotational invariance techniques (ESPRIT)-like method is adopted to estimate the DODs and DOAs by joint diagonalization of the NC-based direction matrices, which can automatically pair the four dimensional (4D) angle parameters and solve the angle ambiguity problem with common one-dimensional (1D) DODs and DOAs. In addition, the asymptotic performance of the proposed algorithm is analyzed and the closed-form stochastic Cramer–Rao bound (CRB) expression is derived. As demonstrated by simulation results, the proposed algorithm has outperformed the existing one, with a result close to the theoretical benchmark.

1. Introduction

A multiple input multiple output (MIMO) radar can provide increased degrees of freedom by exploiting waveform diversity, with an enhanced performance for spatial resolution, parameter estimation, and target detection [1,2,3,4,5,6,7]. In MIMO radar, by focusing on both directions of departure (DODs) and directions of arrival (DOAs), target localization [8,9,10] is an important issue that has drawn significant attention in recent years.
For MIMO radar systems based on one dimensional (1D) uniform linear arrays (ULAs), in reference [11], by employing the property of Kronecker product, a reduced-dimension multiple signal classification (MUSIC) method was developed, and only 1D search was required to locate the DOD and DOA of the target. A double polynomial root MUSIC method was proposed to jointly estimate the DOA and DOD in [12]. To avoid an exhaustive search over the whole angle space, joint DOA and DOD estimation methods were proposed based on the computationally efficient estimating signal parameter via rotational invariance techniques (ESPRIT) method with pairing required in reference [13] and pairing-free in reference [14], respectively. Furthermore, joint diagonalization-based ESPRIT method was presented in reference [15], where a closed-form expression for both DOA and DOD was obtained and automatically paired. For MIMO radar systems based on planar arrays, a method for joint estimation of 2D-DOD and 2D-DOA was presented in references [16,17,18] by transforming the four dimensional (4D) angle estimates into four 1D estimates without any pairing procedures. In addition, the method in reference [17] can also work well when common 1D angles are presented by joint diagonalization, but the method in reference [16] cannot. Based on the electromagnetic vector sensors (EVSs), 2D-DODs and 2D-DOAs were estimated using an ESPRIT-based method for the bistatic MIMO radar, with an additional optimization function for pair matching [18].
However, the above-mentioned methods for MIMO radar failed to exploit some potential information of radar signals, such as noncircularity [19,20,21,22,23], which could lead to significant improvement in its performance. With the aid of noncircular (NC) property of the signals, a series of angle estimation methods [24,25,26] for bistatic MIMO radar were proposed, which can improve the accuracy of angle estimation and detect more signals. In reference [24], with NC incoming signals, a combined ESPRIT and MUSIC approach was applied to MIMO radar for joint estimation of DOA and DOD by decoupling the 2D direction findings into two 1D ones. In reference [25], a kind of ESPRIT algorithm of low complexity was proposed with real-valued computation by Euler’s formula, for MIMO radar with NC signals employed, but it requires additional cost function to avoid ambiguity. For a more general situation with the coexistence of noncircular and circular signals, ESPRIT and unitary ESPRIT were employed in reference [26] for MIMO radar; however, no theoretical error performance analysis was provided for the proposed method, and the derived stochastic Cramer–Rao bound (CRB) does not have a closed-form expression.
In this paper, a NC ESPRIT-like method based on joint diagonalization is proposed to estimate 2D-DOD and 2D-DOA for an L-shaped bistatic MIMO radar. The main contributions of the work are given as follows.
(1)
A general model including a mixture of circular and strictly noncircular sources is built for the L-shaped bistatic MIMO radar by stacking received data vector and its conjugated counterpart. Four NC-based direction matrices are then constructed and by joint diagonalization an ESPRIT-like algorithm is developed employing four block selection matrices.
(2)
The proposed algorithm can work in the case of common 1D DODs and DOAs, and automatically pair the 4D angle parameters.
(3)
The asymptotic performance of the proposed algorithm is analyzed, and the stochastic Cramer–Rao bound (CRB) for the problem is derived with a closed-form expression to serve as the performance benchmark.
The rest of this paper is organized as follows. Section 2 introduces the general mixed signal model for MIMO radar. The proposed algorithm is described in detail in Section 3. The asymptotic performance of the proposed algorithm and the closed-form stochastic CRB are analyzed in Section 4. Simulation results are presented in Section 5, and conclusions are drawn in Section 6.
Notations: ( ) , ( ) T , ( ) 1 , and ( ) H denote conjugate, transpose, inverse, and conjugate transpose, respectively. E ( ) and var ( ) are the expectation and variance operations, respectively; Re ( ) and Im ( ) denote the real and imaginary parts; d i a g ( ) denotes the diagonal matrix; b l k d i a g ( ) represents the generation of a block diagonal matrix; and ⊙ are the Kronecker and Hadamard products, respectively; I k denotes the k -dimensional identity matrix; γ k represents the k -dimensional exchange matrix; 0 k × l and 1 k × l denote the k × l zero matrix and all-one matrix, respectively; arg ( ) is the phase operation; and t r ( ) represents the trace of a matrix.

2. General Signal Model

Consider a bistatic MIMO radar system with an L-shaped antenna array for signal transmission and a second L-shaped antenna array for signal reception, as shown in Figure 1. It is assumed that the target fluctuates according to the Swerling II model [1,2,3,4], i.e., the reflection coefficient changes from pulse to pulse. The transmit array has a total number of M = M 1 + M 2 1 antennas, with M 1 and M 2 antennas located on the X and Y axes, respectively, and the receive array has N = N 1 + N 2 1 antennas, of which N 1 and N 2 elements are located on the X and Y axes, respectively. The four subarrays are all uniform linear arrays (ULAs) with omnidirectional antennas and a half-wavelength inter-element spacing d . The M transmitted waveforms are supposed to be circular (QPSK)- or strictly noncircular (BPSK)-modulated. The targets that slowly move are far field with their directions parameterized as ( θ k 1 , θ k 2 , θ k 3 , θ k 4 ) , where ( θ k 1 , θ k 2 ) is the 2D-DOD of the k th target and ( θ k 3 , θ k 4 ) is its 2D-DOA. The received signals reflected by K targets at the receive array can be written as
r ( l , t ) = k = 1 K α k ( t ) a ( θ k 1 , θ k 2 ) b T ( θ k 3 , θ k 4 ) u ( l , t ) + w ( l , t )
where a k = a ( θ k 1 , θ k 2 ) and b k = b ( θ k 3 , θ k 4 ) are the M × 1 transmit array and N × 1 receive array manifold vectors, with a k = [ e j η d M 2 cos θ k 2 , , 1 , , e j η d M 1 cos θ k 1 ] T and b k = [ e j η d N 2 cos θ k 4 , , 1 , , e j η d N 1 cos θ k 3 ] T , η = 2 π λ 1 , and u ( τ , t ) = [ u 1 T ( τ , t ) , , u M T ( τ , t ) ] T ; u m ( τ , t ) is the mth transmitter antenna signal that is supposed to be circular (QPSK)- or strictly noncircular (BPSK)-modulated, α k ( t ) is the reflection coefficient of the kth target depending on the target radar cross section (RCS), and w ( τ , t ) is the additive white Gaussian noise vector with zero mean and variance σ n 2 . τ and t indicate the time within pulse (fast time) and the index of radar pulse (slow time), respectively. Thus, the output of the matched filters at the receive array can be expressed as
y m ( t ) = k = 1 K a k b k T α k ( t ) r k ( t ) [ 0 0 1 0 0 ] + w m ( t ) = k = 1 K a k b k , m α k ( t ) r k ( t ) + w m ( t ) = k = 1 K a k b k , m s k ( t ) + w m ( t )
where b k , m denotes the m th element of the transmitter steering vector, s k ( t ) = α k ( t ) r k ( t ) is circular or strictly noncircular baseband signal, and w m ( t ) is the noise vector after matched filter. As for strictly noncircular baseband signal, the s k ( t ) can also be written as s k ( t ) = s n , k ( t ) e j φ k / 2 [19,20,21,22,23], where s n , k ( t ) is real-value and φ k / 2 is arbitrary phase shifts that can be different for each signal but are constant with time.
Let x ( t ) = [ y 1 T ( t ) , y M T ( t ) ] T be the output of all the received signal, which is shown as
x ( t ) = C ( θ k 1 , θ k 2 , θ k 3 , θ k 4 ) s ( t ) + n ( t )
where x ( t ) is the M N × 1 data vector, C = [ c 1 , c 2 , , c K ] T is the M N × K extended virtual array manifold matrix, and c k = b k a k is the M N × 1 extended virtual array manifold vector. n ( t ) = [ n 1 ( t ) , , n M N ( t ) ] T is the M N × 1 additive white Gaussian noise vector with zero mean and variance σ n 2 . s ( t ) = [ s 1 ( t ) , , s K ( t ) ] T is the K × 1 mixed signal vector, which contains K n strictly noncircular signals s n , k ( t ) , k = 1 , 2 , , K n and K c circular signals s c , k ( t ) , k = 1 , 2 , , K c , satisfying K = K n + K c . As shown in [27,28], each of the circular signals can be separated into two uncorrelated strictly noncircular signals. Thus, s ( t ) can be rewritten as
s ( t ) = [ Φ 1 0 0 [ I K c j I K c ] ] [ s n ( t ) s c r ( t ) s c q ( t ) ] = Φ s ˜ ( t )
where Φ 1 = d i a g ( e j φ 1 / 2 , , e j φ K n / 2 ) is the K n × K n arbitrary phase matrix corresponding to the strictly noncircular signals s n ( t ) ; furthermore, Φ is of size K × K with K = K n + 2 K c , and the K × 1 real-valued vector s ˜ ( t ) contains the symbols of the K n strictly noncircular signals s n ( t ) cum the K c real parts s c r ( t ) and K c imaginary parts s c q ( t ) of the circular signals s c ( t ) . Therefore, the extended virtual array manifold matrix C can be rewritten as
C = [ C n C c ]
where C n and C c represent the M N × K n and M N × K c array manifold matrix related to strictly noncircular and circular signals, respectively.
According to Equations (4) and (5), the data vector of Equation (3) can be expressed as
x ( t ) = C Φ s ˜ ( t ) + n ( t ) = C n Φ 1 s n ( t ) + C c s c ( t ) + n ( t )
For notional convenience, the angle pair ( θ k 1 , θ k 2 , θ k 3 , θ k 4 ) and time t will be omitted in the following sections.

3. The Proposed Algorithm

In order to utilize the noncircularity characteristic of the strictly noncircular signals and the virtual noncircularity characteristic of the circular signals, a new data matrix is constructed by stacking the original data matrix X = [ x ( 1 ) , , x M N ( T ) ] ( T is the number of snapshots) and its corresponding conjugated counterpart as
Y = [ X γ M N X ] = [ C Φ S ˜ γ M N C Φ S ˜ ] + [ N γ M N N ] = C ˜ S ˜ + N ˜
where
C ˜ = [ C Φ γ M N C Φ ] = [ C n Φ 1 C c [ I K c j I K c ] γ M N C n Φ 1 γ M N C c [ I K c j I K c ] ]
is the 2 M N × K extended array manifold matrix, N ˜ = [ N γ M N N ] is the 2 M N × T noise matrix with N = [ n ( 1 ) , , n M N ( T ) ] , and S ˜ = S ˜ with S ˜ = [ s ˜ ( 1 ) , , s ˜ M N ( T ) ] .
Performing singular value decomposition (SVD) on Y , one can get
Y = ( U s U n ) ( Σ s 0 0 Σ n ) ( V s H V n H ) = U s Σ s V s H + U n Σ n V n H
where the 2 M N × K matrix U s and the T × K matrix V s are the left and right singular signal subspace associated with corresponding left and right singular values matrices Σ s = d i a g ( λ 1 , λ 2 , , λ K ) and Σ n = d i a g ( λ K + 1 , λ K + 2 , , λ 2 M N ) , respectively, while the 2 M N × ( 2 M N K ) matrix U n and the T × ( 2 M N K ) matrix V n are the left and right singular noise subspace, respectively.
By defining a new matrix E s as E s = U s Σ s , and the following selection matrices
J 1 a = [ 0 ( a 1 ) × ( M a ) 0 ( a 1 ) × 1 I ( a 1 ) ] , a = M 1 , N 1
J 2 a = [ 0 ( a 1 ) × ( M a ) I ( a 1 ) 0 ( a 1 ) × 1 ] , a = M 1 , N 1
J 1 b = [ I ( b 1 ) 0 ( b 1 ) × 1 0 ( b 1 ) × ( M b ) ] , b = M 2 , N 2
J 2 b = [ 0 ( b 1 ) × 1 I ( b 1 ) 0 ( b 1 ) × ( M b ) ] , b = M 2 , N 2
the selection matrices displayed in Figure 2 for θ k l ( l = 1 , 2 ) of the mixed strictly noncircular and circular signals can be expressed as
K l 1 = b l k d i a g ( J l 1 , γ ( M l 1 ) M J l 2 γ M N ) , l = 1 , 2
K l 2 = b l k d i a g ( J l 2 , γ ( M l 1 ) M J l 1 γ M N ) , l = 1 , 2
where J 11 = I N J 1 M 1 , J 12 = I N J 2 M 1 , J 21 = I N J 1 M 2 , and J 22 = I N J 2 M 2 . Similarly, as shown in Figure 3, the selection matrices for θ k l ( l = 3 , 4 ) of the mixed signals can be expressed as
K l 1 = b l k d i a g ( J l 1 , γ ( N l 1 ) N J l 2 γ N M ) , l = 3 , 4
K l 2 = b l k d i a g ( J l 2 , γ ( N l 1 ) N J l 1 γ N M ) , l = 3 , 4
where J 31 = J 1 N 1 I M , J 32 = J 2 N 1 I M , J 41 = J 1 N 2 I M , and J 42 = J 2 N 2 I M .
Instead of complex peak-seeking methods [29,30,31,32], following the principle of the ESPRIT algorithm [33,34,35], we define the direction matrices G l related to θ k l ( l = 1 , 2 , 3 , 4 ) as follows
G l = ( K l 2 E s ) + K l 1 E s = E Θ l E H , l = 1 , 2 , 3 , 4
where
Θ l = d i a g ( e j 2 π λ 1 d cos θ 1 l , , e j 2 π λ 1 d cos θ K l ) = d i a g ( η 1 l , , η K l )
is a diagonal matrix and E is the K × K unitary matrix.
It can be seen that G l in Equation (18) satisfies the joint diagonalization condition. Then, we define a set G = { G 1 , G 2 , G 3 , G 4 } and use the joint diagonalization method in [28,36,37] to obtain the unitary matrix E = [ e 1 , e 2 , , e K ] , where e k is the eigenvector of G . It should be mentioned that the proposed method does not require the 4D angle pairing process, as the eigenvalues of G maintain a one-to-one correspondence in the joint diagonalization process. Then, the eigenvalues of G can be computed as
η k l = e k H G l e k , k = 1 , 2 , , K , l = 1 , 2 , 3 , 4
Thus, it can be easily obtained that
θ k l = arccos ( λ arg ( η k l ) 2 π d ) , k = 1 , 2 , , K , l = 1 , 2 , 3 , 4
It should be noted that each circular signal is treated as two strictly noncircular signals, and K angle estimates are obtained for the mixed targets. However, only K actual 2D-DODs and 2D-DOAs are present, so the circular and strictly noncircular signals can be discriminated according to the number of repetitions of the angle estimates. Then, the two estimated angles of the circular signal are reliable, and θ k l , c ( l = 1 , 2 , 3 , 4 ) can be obtained by calculating the average of two identical estimates
θ k l , c = θ k l , c 1 + θ k l , c 2 2 , k = 1 , , K c
Till now, the proposed method has provided closed-form of 2D-DOA and 2D-DOD angle estimates that are automatically paired and summarized in Table 1.
Remark 1. The major computational effort the proposed algorithm contains SVD of Y ^ , pseudo inverse operation for G l in Equation (18), and joint diagonalization of the set G ^ . Performing SVD of Y ^ requires the amount of complex multiplications of O ( ( 2 M N ) 2 L ) , the pseudo inverse in Equation (18) costs O ( 4 ( 2 M N ) 3 ) , and jointly diagonalizing the set G ^ is of O ( 4 ( K ) 3 ) . The total computational complexity of the proposed algorithm is about O ( ( 2 M N ) 2 L + 4 ( 2 M N ) 3 + 4 ( K ) 3 ) .
Remark 2. Compared to the Xia’s algorithm [17], the proposed algorithm exploits the redundancy existing in the noncircular signals, which improves the array virtual aperture. Additionally, the maximum numbers of detectable signals by the proposed algorithm is based on the new data vetor in Equation (7) as well as the matrices K l 1 and K l 2 , l = 1 , 2 in Equations (14) and (15) for 2D DODs, K l 1 and K l 2 , l = 3 , 4 in Equations (16) and Equation (17) for 2D DOAs, which are shown in Table 2 compared to Xia’s algorithm. Obviously, the proposed algorithm can distinguish more mixed signals than Xia’s algorithm.

4. Performance Analysis

4.1. Asymptotic Performance Analysis

In this section, the asymptotic performance of the proposed algorithm is derived, which is consistent with the first-order analysis done by Rao [38] and the backward error analysis of Li [39]. For the ESPRIT-like subspace algorithm, we need to analyze subspace perturbation as a criterion for evaluation. Therefore, we can perform SVD on the noiseless extended observation model Y ˜ as follows:
Y ˜ = ( U s U n ) ( Σ s 0 0 Σ n ) ( V s H V n H )
In line with the first-order approximation principle [38,39] of eigenvalues in Equation (18), we get
δ η k l e k H δ G l e k , k = 1 , 2 , , K , l = 1 , 2 , 3 , 4 = e k H ( K l 2 U y ) + ( K l 1 η k l K l 2 ) δ U y e k
where
U y = U s Σ s
δ U y = δ U s Σ s
δ U s = U n U n H N ˜ V s Σ s 1
According to Equations (26) and (27), Equation (24) can be rewritten as
δ η k l e k H ( K l 2 U y ) + ( K l 1 U n η k l K l 2 U n ) U n H N ˜ V s e k
By performing the first-order Taylor series expansion on Equation (21), the perturbation of θ k l can be expressed as
δ θ k l = ξ k l Im ( δ η k l η k l )
where ξ k l = λ ( 2 π d sin θ k l ) 1 . The error-variances of the estimated 2D-DODs and 2D-DOAs of the mixed sources are
var ( δ θ k l ) = ξ k l 2 var { Im ( δ η k l η k l ) }
It is worth noting here that Equation (30) can only calculate the mean-squared error for the strictly noncircular signals, while the variances of the four estimated angles for the k th k = 1 , 2 , , K c circular signal are calculated as
var ( δ θ k l , c ) = var ( δ θ k l , c 1 + δ θ k l , c 2 2 ) = var ( δ θ k l , c 1 ) + var ( δ θ k l , c 2 ) + 2 E ( δ θ k l , c 1 δ θ k l , c 2 ) 4

4.2. Stochastic Cramer–Rao Bound

The CRB, which has a lower bound on the variance of any unbiased estimator, is often adopted for the performance benchmark. In reference [20], the CRB is analyzed by assuming that the set of incident sources are all strictly noncircular, while the CRB analyzed in reference [40], provided that the incident sources are all circular. However, when considering a scenario that both the strictly noncircular and circular sources coexist, the two signal models mentioned above are not applicable. In this section, the stochastic closed-form CRB is first derived for the estimates of 2D-DOA and 2D-DOA of the mixed strictly noncircular and circular sources based on the L-shaped bistatic MIMO radar.
Here, we focus on a real-valued vector of the interest parameters ω = [ θ n 1 T , θ n 2 T , θ n 3 T , θ n 4 T , ψ n T , θ c 1 T , θ c 2 T , θ c 3 T , θ c 4 T ] T with θ n 1 = [ θ n 1 , 1 , , θ n 1 , K n ] T , θ n 2 = [ θ n 2 , 1 , , θ n 2 , K n ] T , θ n 3 = [ θ n 3 , 1 , , θ n 3 , K n ] T , θ n 4 = [ θ n 4 , 1 , , θ n 4 , K n ] T , ψ n = [ φ n , 1 , , φ n , K n ] T , θ c 1 = [ θ c 1 , K n + 1 , , θ c 1 , K ] T , θ c 2 = [ θ c 2 , K n + 1 , , θ c 2 , K ] T , θ c 3 = [ θ c 3 , K n + 1 , , θ c 3 , K ] T , and θ c 4 = [ θ c 4 , K n + 1 , , θ c 4 , K ] T . Then, followed by references [40,41,42,43], the (p, q)th entry of the ( 5 K n + 4 K c ) × ( 5 K n + 4 K c ) CRB matrix for the parameter ω estimates is given by
[ C R B 1 ( ω ) ] p , q = 2 L σ n 2 Re { t r [ A ω p P A A ω q Q ] }
where ω k is the kth element of ω ,
A = [ C n Φ 1 C c ]
P A = I A ( A H A ) 1 A H
Q = P s A H R 1 A P s = [ Q 1 Q 2 Q 3 Q 4 ]
and Q 1 , Q 2 , Q 3 , and Q 4 are sizes of K n × K n , K n × K c , K c × K n , and K c × K c matrices, respectively,
P s = E ( [ s n s c ] [ s n s c ] H )
R = E [ x x H ]
Define
D n = [ D θ n 1 , D θ n 2 , D θ n 3 , D θ n 4 , j 2 C n Φ 1 , 0 M N × 4 K c ]
D c = [ 0 M N × 5 K n , D θ c 1 , D θ c 2 , D θ c 3 , D θ c 4 ]
with D θ n 1 = [ C n θ n 1 , 1 Φ 1 , , C n θ n 1 , K n Φ 1 ] , D θ n 2 = [ C n θ n 2 , 1 Φ 1 , , C n θ n 2 , K n Φ 1 ] , D θ 3 = [ C n θ n 3 , 1 Φ 1 , , C n θ n 3 , K n Φ 1 ] , D θ n 4 = [ C n θ n 4 , 1 Φ 1 , , C n θ n 4 , K n Φ 1 ] , D θ c 1 = [ C c θ c 1 , K n + 1 , , C c θ c 1 , K ] , D θ c 2 = [ C c θ c 2 , K n + 1 , , C c θ c 2 , K ] , D θ c 3 = [ C c θ c 3 , K n + 1 , , C c θ c 3 , K ] , and D θ c 4 = [ C c θ c 4 , K n + 1 , , C c θ c 4 , K ] .
After some simplifications, we obtain the closed-form expression for the CRB ω as
C R B ( ω ) = σ n 2 2 L [ C R B 1 ( ω ) C R B 2 ( ω ) C R B 3 ( ω ) C R B 4 ( ω ) ] 1
where
C R B 1 ( ω ) = Re { J n ( D n H P A D n ) J n T ( 1 5 1 5 T Q 1 T ) } 1
C R B 2 ( ω ) = Re { J n ( D n H P A D c ) J c T ( 1 4 1 5 T Q 3 T ) } 1
C R B 3 ( ω ) = Re { J c ( D c H P A D n ) J n T ( 1 5 1 4 T Q 2 T ) } 1
C R B 4 ( ω ) = Re { J c ( D c H P A D c ) J c T ( 1 4 1 4 T Q 4 T ) } 1
with J n = [ I 5 K n 0 5 K n × 4 K c ] , J c = [ 0 4 K c × 5 K n I 4 K c ] , 1 5 = [ 1 , 1 , 1 , 1 , 1 ] and 1 4 = [ 1 , 1 , 1 , 1 ] .
It should be noted that the C R B ( ω ) will degenerate into C R B 1 ( ω ) and C R B 4 ( ω ) , which corresponds to the incident sources are all strictly noncircular and circualr sources, respectively, namely, C R B ( ω ) = C R B 1 ( ω ) and C R B ( ω ) = C R B 4 ( ω ) .

5. Simulation Results

In this part, we evaluate the effectiveness of the proposed method in terms of several simulations. The proposed algorithm is compared with Xia’s algorithm [17], asymptotic performance analysis (Proposed asy.) in Equations (30) and (31), and the derived stochastic CRB in Equation (40). We use the root mean square error (RMSE) given by R M S E _ D O D = 1 K M c k = 1 K m = 1 M c [ ( θ ^ k 1 θ k 1 ) 2 + ( θ ^ k 2 θ k 2 ) 2 ] and R M S E _ D O A = 1 K M c k = 1 K m = 1 M c [ ( θ ^ k 3 θ k 3 ) 2 + ( θ ^ k 4 θ k 4 ) 2 ] as the performance criterion, where M c is the number of Monte-Carlo trials. The first experiment is based on an L-shaped MIMO array with N 1 = 3 , M 1 = M 2 = N 2 = 2 , and for the next two experiments, M 1 = M 2 = N 1 = N 2 = 3 , d is half wavelength.
Experiment 1. In the first experiment, we verify that the mixture of strictly noncircular and circular signals can be estimated successfully by the proposed method with the increasing of the maximum number of detectable signal in comparison with Xia’s method. Here, we consider six uncorrelated targets in the experiment, which have five BPSK signals with direction pairs ( 70 ° , 105 ° , 70 ° , 105 ° ) , ( 60 ° , 80 ° , 60 ° , 80 ° ) , ( 85 ° , 90 ° , 85 ° , 90 ° ) , ( 100 ° , 70 ° , 100 ° , 70 ° ) , and ( 110 ° , 100 ° , 110 ° , 100 ° ) , and one QPSK signal with direction pairs ( 120 ° , 85 ° , 120 ° , 85 ° ) . Then, we have K n + 2 K c = 7 , min { 2 N ( M 1 1 ) , 2 N ( M 2 1 ) } D O D = 8 , and min { 2 M ( N 1 1 ) , 2 M ( N 2 1 ) } D O A = 6 . The signal-to-noise ratio (SNR) is set at 25 dB, the number of snapshots is 500, and M c = 100 . Figure 4 and Figure 5 show the 2D-DODs and 2D-DOAs scattergram of six mixed signals, respectively. It can be seen that the proposed algorithm can estimate the 2D-DODs and 2D-DOAs of six targets correctly, while the algorithm in reference [17] fails to work, because the former can detect more signals with available noncircular information. It should also be noted that the 2D-DOA estimation for QPSK signal is slightly inaccurate, because the number of mixed targets has exceeded the maximum number that the 2D-DOA can detect, but it still can roughly estimate the 2D-DOA of QPSK signal.
Experiment 2. In the second experiment, the performance of the proposed algorithm is investigated with SNR varying from −5 dB to 15 dB. We consider four uncorrelated mixed signals with direction pairs ( 60 ° , 50 ° , 60 ° , 50 ° ) , ( 70 ° , 50 ° , 80 ° , 70 ° ) , ( 70 ° , 60 ° , 80 ° , 90 ° ) , and ( 80 ° , 70 ° , 100 ° , 90 ° ) . We consider the cases of one, two, three, and four BPSK signals, and the remaining signals are QPSK, respectively. The number of snapshots is 300 and M c = 2000. In Figure 6 and Figure 7, the estimation performance of the proposed algorithm is shown to be superior to Xia’s algorithm for both 2D-DOD and 2D-DOA estimation in all four cases. From case 1 to case 4, the performance of the proposed algorithm is improved in turn, as more noncircular information is available. In addition, the RMSEs of the proposed algorithm vary almost in accordance with their asymptotic error-variances, and both of them are close to the CRBs, especially for case 3 and case 4.
Experiment 3. In the third experiment, we investigate the performance with respect to a varying number of snapshots ranging from 50 to 950. The SNR is set at 5Db, and the other parameters are the same as Experiment 2. As shown in Figure 8 and Figure 9, we can draw similar conclusions as Experiment 2, that the proposed algorithm has better performance with the number of snapshots increases in all four cases, again outperforms the Xia’s method, and is close to the theoretical benchmark.

6. Conclusions

Based on the joint diagonalization technique, a 2D-DOD and 2D-DOA estimation algorithm for mixed strictly noncircular and circular signals in L-shaped bistatic MIMO radar is proposed in this paper. It utilizes the noncircularity characteristic to construct a virtual array, and then derives the joint diagonalization-based NC-ESPRIT method to achieve automatic pairing and the identification of the estimated 4D angles of mixed signals. The asymptotic performance of the proposed method as well as the stochastic CRB for the mixed signals scenario is also derived. Simulation results show that the proposed algorithm has a better angle estimation performance than the algorithm without noncircularity characteristics.

Author Contributions

J.F., Y.L. (Yonghong Liu), and H.C. did the overall analysis and writing; Y.J. and Y.L. (Yang Lu) helped carry out the data statistics; and Z.Z. and L.W. did the review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by Key Laboratory of Intelligent Perception and Advanced Control of State Ethnic Affairs Commission under Grant MD-IPAC-2019102, and by Zhejiang Provincial Natural Science Foundation of China under Grant LQ19F010002, and by Natural Science Foundation of Ningbo Municipality under Grant 2018A610094, and by K. C. Wong Magna Fund in Ningbo University, and by K. C. Wong Magna Fund in Ningbo University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, W.; Liu, W.; Wang, J.; Wu, S. Joint Transmission and Reception Diversity Smoothing for Direction Finding of Coherent Targets in MIMO Radar. IEEE J. Sel. Top. Signal Process. 2014, 8, 115–124. [Google Scholar] [CrossRef] [Green Version]
  2. Shi, J.; Hu, G.; Zhang, X.; Sun, F. Sparsity-Based DOA Estimation of Coherent and Uncorrelated Targets With Flexible MIMO Radar. IEEE Trans. Veh. Technol. 2019, 68, 5835–5848. [Google Scholar] [CrossRef]
  3. Wen, F.; Shi, J.; Zhang, Z. Joint 2D-DOD, 2D-DOA, and Polarization Angles Estimation for Bistatic EMVS-MIMO Radar via PARAFAC Analysis. IEEE Trans. Veh. Technol. 2020, 69, 1626–1638. [Google Scholar] [CrossRef]
  4. Liao, B. Fast Angle Estimation for MIMO Radar with Nonorthogonal Waveforms. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 2091–2096. [Google Scholar] [CrossRef]
  5. Zhou, C.; Gu, Y.; Fan, X.; Shi, Z.; Mao, G.; Zhang, Y.D. Direction-of-Arrival Estimation for Coprime Array via Virtual Array Interpolation. IEEE Trans. Signal Process. 2018, 66, 5956–5971. [Google Scholar] [CrossRef]
  6. Wen, F.; Zhang, Z.; Zhang, G. Joint DOD and DOA Estimation for Bistatic MIMO Radar: A Covariance Trilinear Decomposition Perspective. IEEE Access 2019, 7, 53273–53283. [Google Scholar] [CrossRef]
  7. Shi, Z.; Zhou, C.; Gu, Y.; Goodman, N.A.; Qu, F. Source Estimation Using Coprime Array: A Sparse Reconstruction Perspective. IEEE Sens. J. 2017, 17, 755–765. [Google Scholar] [CrossRef]
  8. Li, J.; Jiang, D.; Zhang, X. DOA Estimation Based on Combined Unitary ESPRIT for Coprime MIMO Radar. IEEE Commun. Lett. 2017, 21, 96–99. [Google Scholar] [CrossRef]
  9. Zhou, C.; Gu, Y.; He, S.; Shi, Z. A Robust and Efficient Algorithm for Coprime Array Adaptive Beamforming. IEEE Trans. Veh. Technol. 2018, 67, 1099–1112. [Google Scholar] [CrossRef]
  10. Wen, F. Computationally Efficient DOA Estimation Algorithm for MIMO Radar with Imperfect Waveforms. IEEE Commun. Lett. 2019, 23, 1037–1040. [Google Scholar] [CrossRef]
  11. Zhang, X.; Xu, L.; Xu, L.; Xu, D. Direction of Departure (DOD) and Direction of Arrival (DOA) Estimation in MIMO Radar with Reduced-Dimension MUSIC. IEEE Commun. Lett. 2010, 14, 1161–1163. [Google Scholar] [CrossRef]
  12. Bencheikh, M.L.; Wang, Y.; He, H. Polynomial root finding technique for joint DOA DOD estimation in bistatic MIMO radar. Signal Process. 2010, 90, 2723–2730. [Google Scholar] [CrossRef]
  13. Duofang, C.; Baixiao, C.; Guodong, Q. Angle estimation using ESPRIT in MIMO radar. Electron. Lett. 2008, 44, 770–771. [Google Scholar] [CrossRef]
  14. Zheng, G.; Chen, B.; Yang, M. Unitary ESPRIT algorithm for bistatic MIMO radar. Electron. Lett. 2012, 48, 179–181. [Google Scholar] [CrossRef]
  15. Xia, T.Q. Joint diagonalization based DOD and DOA estimation for bistatic MIMO radar. Signal Process. 2015, 108, 159–166. [Google Scholar] [CrossRef]
  16. Li, J.; Zhang, X. Closed-form blind 2D-DOD and 2D-DOA estimation for MIMO radar with arbitrary arrays. Wirel. Pers. Commun. 2013, 69, 175–186. [Google Scholar] [CrossRef]
  17. Xia, T.Q. Joint diagonalization based 2D-DOD and 2D-DOA estimation for bistatic MIMO radar. Signal Process. 2015, 116, 7–12. [Google Scholar] [CrossRef]
  18. Chintagunta, S.; Ponnusamy, P. 2D-DOD and 2D-DOA estimation using the electromagnetic vector sensors. Signal Process. 2018, 147, 163–172. [Google Scholar] [CrossRef]
  19. Chen, H.; Liu, W.; Zhu, W.P.; Swamy, M.N.S.; Wang, Q. Mixed rectilinear sources localization under unknown mutual coupling. J. Frankl. Inst. (Elsevier) 2019, 356, 2372–2394. [Google Scholar] [CrossRef] [Green Version]
  20. Abeida, H.; Delmas, J.P. Direct derivation of the stochastic CRB of DOA estimation for rectilinear sources. IEEE Signal Process. Lett. 2017, 24, 1522–1526. [Google Scholar] [CrossRef] [Green Version]
  21. Wang, X.; Wan, L.; Huang, M.; Shen, C.; Zhang, K. Polarization Channel Estimation for Circular and Non-Circular Signals in Massive MIMO Systems. IEEE J. Sel. Top. Signal Process. 2019, 13, 1001–1016. [Google Scholar] [CrossRef]
  22. Chen, H.; Wang, W.; Liu, W. Joint DOA, Range, and Polarization Estimation for Rectilinear Sources with a COLD Array. IEEE Wirel. Commun. Lett. 2019, 8, 1398–1401. [Google Scholar] [CrossRef]
  23. Wang, X.; Wan, L.; Huang, M.; Shen, C.; Han, Z.; Zhu, T. Low-Complexity Channel Estimation for Circular and Noncircular Signals in Virtual MIMO Vehicle Communication Systems. IEEE Trans. Veh. Technol. 2020. [Google Scholar] [CrossRef]
  24. Bencheikh, M.L.; Wang, Y. Non circular ESPRIT-Root MUSIC joint DOA-DOD estimation in bistatic MIMO radar. In Proceedings of the International Workshop on Systems, Signal Processing and their Applications, Tipaza, Algeria, 9–11 May 2011; pp. 9–11. [Google Scholar]
  25. Guo, Y.D.; Zhang, Y.S.; Gong, J.; Zheng, G.M. Direction Finding with Real-Valued ESPRIT for Noncircular Signal in Bistatic MIMO Radar. Wirel. Pers. Commun. 2017, 95, 3321–3332. [Google Scholar] [CrossRef]
  26. Zheng, G.M.; Tang, J.; Yang, X. ESPRIT and Unitary ESPRIT Algorithms for Coexistence of Circular and Noncircular Signals in Bistatic MIMO Radar. IEEE Access. 2016, 4, 7232–7240. [Google Scholar] [CrossRef]
  27. Steinwandt, J.; Roemer, F.; Haardt, M. ESPRIT-type algorithms for a received mixture of circular and strictly non-circular signals. In Proceedings of the IEEE Int. Conference on Acoustics, Speech, and Signal Processing (ICASSP), Brisbane, Australia, 19–24 April 2015; pp. 2809–2813. [Google Scholar]
  28. Chen, H.; Hou, C.P.; Zhu, W.P.; Liu, W.; Dong, Y.Y.; Peng, Z.J.; Wang, Q. ESPRIT-like two-dimensional direction finding for mixed circular and strictly noncircular sources based on joint diagonalization. Signal Process. (Elsevier) 2017, 141, 48–56. [Google Scholar] [CrossRef] [Green Version]
  29. Zhou, C.; Gu, Y.; Zhang, Y.D.; Shi, Z.; Jin, T.; Wu, X. Compressive Sensing based Coprime Array Direction-of-Arrival Estimation. IET Commun. 2017, 11, 1719–1724. [Google Scholar] [CrossRef]
  30. Chen, H.; Zhu, W.P.; Liu, W.; Swamy MN, S.; Li, Y.; Wang, Q.; Peng, Z. RARE-based localization for mixed near-field and far-field rectilinear sources. Digit. Signal Process. 2019, 85, 54–61. [Google Scholar] [CrossRef]
  31. Zhou, C.; Gu, Y.; Shi, Z.; Zhang, Y.D. Off-Grid Direction-of-Arrival Estimation Using Coprime Array Interpolation. IEEE Signal Process. Lett. 2018, 25, 1710–1714. [Google Scholar] [CrossRef]
  32. Zhou, C.; Shi, Z.; Gu, Y.; Shen, X. DECOM: DOA Estimation with Combined MUSIC for Coprime Array. In Proceedings of the International Conference on Wireless Communications and Signal Processing (WCSP), Hangzhou, China, 24–26 October 2013; pp. 1–5. [Google Scholar]
  33. Zheng, Z.; Wang, W.; Meng, H.; So, H.C.; Zhang, H. Efficient beamspace-based algorithm for two-dimensional DOA estimation of incoherently distributed sources in massive MIMO systems. IEEE Trans. Veh. Technol. 2018, 67, 11776–11789. [Google Scholar] [CrossRef]
  34. Zhou, C.; Zhou, J. Direction-of-arrival estimation with coarray ESPRIT for coprime array. Sensors 2017, 17, 1779. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Chen, H.; Wang, W.; Liu, W. Augmented Quaternion ESPRIT-Type DOA Estimation with a Crossed-Dipole Array. IEEE Commun. Lett. 2020, 24, 548–552. [Google Scholar] [CrossRef]
  36. Belouchrani, A.; Abed-Meraim, K.; Cardoso, J.F.; Moulines, E. A blind source separation technique using second-order statistics. IEEE Trans. Signal Process. 1997, 45, 434–444. [Google Scholar] [CrossRef] [Green Version]
  37. Dong, Y.; Dong, C.; Zhu, Y.; Zhao, G.; Liu, S. Two-dimensional DOA estimation for L-shaped array with nested subarrays without pair matching. IET Signal Process. 2016, 10, 1112–1117. [Google Scholar] [CrossRef]
  38. Rao, B.D.; Hari, K.V.S. Performance analysis of ESPRIT and TAM in determining the direction of arrival of plane waves in noise. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 1990–1995. [Google Scholar] [CrossRef]
  39. Li, F.; Vaccaro, R.J.; Tufts, D.W. Performance analysis of the state-space realization (TAM) and ESPRIT algorithms for DOA estimation. IEEE Trans. Antennas Propag. 1991, 39, 418–423. [Google Scholar] [CrossRef]
  40. Stoica, P.; Larsson, A.G.; Gershman, A.B. The stochastic CRB for array processing: A textbook derivation. IEEE Signal Process. Lett. 2001, 8, 148–150. [Google Scholar] [CrossRef]
  41. Wen, F.; Wang, J.; Shi, J.; Gui, G. Auxiliary Vehicle Positioning Based on Robust DOA Estimation with Unknown Mutual Coupling. IEEE Internet Things J. 2020. [Google Scholar] [CrossRef]
  42. He, J.; Swamy, M.N.S.; Ahmad, M.O. Efficient Application of MUSIC Algorithm Under the Coexistence of Far-Field and Near-Field Sources. IEEE Trans. Signal Process. 2012, 60, 2066–2070. [Google Scholar] [CrossRef]
  43. Wen, F.; Zhang, X.; Zhang, Z. CRBs for direction-of-departure and direction-of-arrival estimation in collocated MIMO radar in the presence of unknown spatially coloured noise. IET Radar Sonar Navig. 2019, 13, 530–537. [Google Scholar] [CrossRef]
Figure 1. A diagram for the L-shaped multiple input multiple output (MIMO) array structure.
Figure 1. A diagram for the L-shaped multiple input multiple output (MIMO) array structure.
Sensors 20 02177 g001
Figure 2. Block selection matrices for estimating ( θ k 1 , θ k 2 ) .
Figure 2. Block selection matrices for estimating ( θ k 1 , θ k 2 ) .
Sensors 20 02177 g002
Figure 3. Block selection matrices for estimating ( θ k 3 , θ k 4 ) .
Figure 3. Block selection matrices for estimating ( θ k 3 , θ k 4 ) .
Sensors 20 02177 g003
Figure 4. Two dimensional (2D)-direction of departure (DOD) scattergram of six targets for mixed signals.
Figure 4. Two dimensional (2D)-direction of departure (DOD) scattergram of six targets for mixed signals.
Sensors 20 02177 g004
Figure 5. 2D-DOA scattergram of six mixed signals.
Figure 5. 2D-DOA scattergram of six mixed signals.
Sensors 20 02177 g005
Figure 6. Root mean square error (RMSE) of 2D-DOD for mixed signals versus signal-to-noise ratio (SNR) (a–d).
Figure 6. Root mean square error (RMSE) of 2D-DOD for mixed signals versus signal-to-noise ratio (SNR) (a–d).
Sensors 20 02177 g006
Figure 7. RMSE of 2D-DOA for mixed signals versus SNR (a–d).
Figure 7. RMSE of 2D-DOA for mixed signals versus SNR (a–d).
Sensors 20 02177 g007
Figure 8. RMSE of 2D-DOD for mixed signals versus snapshots (a–d).
Figure 8. RMSE of 2D-DOD for mixed signals versus snapshots (a–d).
Sensors 20 02177 g008
Figure 9. RMSE of 2D-DOA for mixed signals versus snapshots (a–d).
Figure 9. RMSE of 2D-DOA for mixed signals versus snapshots (a–d).
Sensors 20 02177 g009
Table 1. Summary of the proposed method.
Table 1. Summary of the proposed method.
Input:   { x ^ ( t ) } t = 1 , T : T snapshots of the new constructed array vector.
Output:   { θ ^ k l } k = 1 , K , l = 1 , 2 , 3 , 4 : pair-free 2D-DODs and 2D-DOAs of K mixed signals
Step 1: Perform SVD on Y ^ to get U ^ s , and then compute E ^ s = U ^ s Σ ^ s ;
Step 2: Define a set G ^ = { G ^ 1 , G ^ 2 , G ^ 3 , G ^ 4 } according to Equation (18)
Step 3: Implement the joint diagonalization to the set G ^ to obtain the unitary matrix E ^ by a series of Givens rotations;
Step 4: Compute the eigenvalues η ^ k l according to Equation (20), and then compute θ ^ k l according to Equation (21);
Step 5: Compute the 2-D DODs and 2-D DOAs of circular signals according to Equation (22).
Table 2. Maximum number of detection signals.
Table 2. Maximum number of detection signals.
AlgorithmAngleMaximum Number
Proposed algorithmDOD K n + 2 K c = min { 2 N ( M 1 1 ) , 2 N ( M 2 1 ) }
DOA K n + 2 K c = min { 2 M ( N 1 1 ) , 2 M ( N 2 1 ) }
Xia’s algorithmDOD K n + K c = min { N ( M 1 1 ) , N ( M 2 1 ) }
DOA K n + K c = min { M ( N 1 1 ) , M ( N 2 1 ) }

Share and Cite

MDPI and ACS Style

Fang, J.; Liu, Y.; Jiang, Y.; Lu, Y.; Zhang, Z.; Chen, H.; Wang, L. 2D-DOD and 2D-DOA Estimation for a Mixture of Circular and Strictly Noncircular Sources Based on L-Shaped MIMO Radar. Sensors 2020, 20, 2177. https://0-doi-org.brum.beds.ac.uk/10.3390/s20082177

AMA Style

Fang J, Liu Y, Jiang Y, Lu Y, Zhang Z, Chen H, Wang L. 2D-DOD and 2D-DOA Estimation for a Mixture of Circular and Strictly Noncircular Sources Based on L-Shaped MIMO Radar. Sensors. 2020; 20(8):2177. https://0-doi-org.brum.beds.ac.uk/10.3390/s20082177

Chicago/Turabian Style

Fang, Jiaxiong, Yonghong Liu, Yifang Jiang, Yang Lu, Zehao Zhang, Hua Chen, and Laihua Wang. 2020. "2D-DOD and 2D-DOA Estimation for a Mixture of Circular and Strictly Noncircular Sources Based on L-Shaped MIMO Radar" Sensors 20, no. 8: 2177. https://0-doi-org.brum.beds.ac.uk/10.3390/s20082177

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