Next Article in Journal
Automatic Non-Destructive Growth Measurement of Leafy Vegetables Based on Kinect
Next Article in Special Issue
Fusion of an Ensemble of Augmented Image Detectors for Robust Object Detection
Previous Article in Journal
Electrode Coverage Optimization for Piezoelectric Energy Harvesting from Tip Excitation
Previous Article in Special Issue
Lightweight Active Object Retrieval with Weak Classifiers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Target Angle Tracking Algorithm for Bistatic Multiple-Input Multiple-Output (MIMO) Radar Based on the Elements of the Covariance Matrix

National University of Defense Technology, Hefei 230037, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(3), 805; https://doi.org/10.3390/s18030805
Submission received: 22 January 2018 / Revised: 28 February 2018 / Accepted: 1 March 2018 / Published: 7 March 2018
(This article belongs to the Special Issue Sensor Signal and Information Processing)

Abstract

:
In this paper, we consider the problem of tracking the direction of arrivals (DOA) and the direction of departure (DOD) of multiple targets for bistatic multiple-input multiple-output (MIMO) radar. A high-precision tracking algorithm for target angle is proposed. First, the linear relationship between the covariance matrix difference and the angle difference of the adjacent moment was obtained through three approximate relations. Then, the proposed algorithm obtained the relationship between the elements in the covariance matrix difference. On this basis, the performance of the algorithm was improved by averaging the covariance matrix element. Finally, the least square method was used to estimate the DOD and DOA. The algorithm realized the automatic correlation of the angle and provided better performance when compared with the adaptive asymmetric joint diagonalization (AAJD) algorithm. The simulation results demonstrated the effectiveness of the proposed algorithm. The algorithm provides the technical support for the practical application of MIMO radar.

1. Introduction

In recent years, the multiple-input multiple-output (MIMO) radar has been proposed as a new system radar [1]. Multiple array elements of the MIMO radar can transmit mutually orthogonal waveforms, which have a high degree of freedom [2,3,4]. Compared with the phased array radar, the MIMO radar has more accurate performance in target detection, identification, parameter estimation, and tracking. According to the array element configuration, the MIMO radar is divided into the statistical MIMO radar and coherent MIMO radar [5]. The array elements of the statistical MIMO radar are far away from each other, therefore, the statistical MIMO radar obtains spatial diversity gain, which can effectively improve the estimate performance of the scintillation target. The transmit and receive elements of the coherent MIMO radar are closely spaced, which can effectively improve the estimation accuracy of the target parameters, increase the number of the maximum identification targets, and so on. In the coherent MIMO radar, the bistatic MIMO radar is an important structure. The bistatic MIMO radar, which combines the advantages of MIMO radar and bistatic radar, effectively reduces the requirement of the three synchronizations (time, space, frequency). Therefore, the bistatic MIMO radar was used as the research object in this paper.
The existing parameter estimation algorithms in the bistatic MIMO radar are mostly aimed at stationary targets [6,7,8,9,10,11,12,13,14,15,16], which contain the estimation of signal parameters via rotational invariance technique (ESPRIT) algorithms [6,7,8,9], the Capon algorithm [10], the multiple signal classification (MUSIC) algorithms [11,12,13,14,15,16], and so on.
However, when the target is moving, the performance of the algorithm in [6,7,8,9,10,11,12,13,14,15,16] will decrease or even fail, which cannot be applied for tracking moving targets. The algorithms described above are all based on feature subspaces and require eigendecomposition. The computational complexity of the eigendecomposition process for the N dimension square matrix is generally o ( N 3 ) . The covariance matrix obtained in bistatic MIMO radar has a large dimension, which is the product of the number of transmitters and receivers. Therefore, the computational complexity required for eigendecomposition is large and the project realization is more difficult. In addition, these eigendecomposition and eigensubspace methods are a class of batch processing methods. Obviously, these algorithms cannot be applied to time-varying signals, but in the actual battlefield environment, the goal often moves. The target angle tracking is the key problem that restricts the practical application of bistatic MIMO radar. Therefore, this paper studies the angle tracking problem in bistatic MIMO radar.
There are some published studies on MIMO radar tracking. The monostatic MIMO radar tracking algorithm is given in [17], which has low complexity, but the cost is the reduction of tracking performance. In [18], Kalman was introduced into the projective approximation subspace tracking with deflation (PASTd) algorithm. The Kalman filter was used to realize data association, and the algorithm converges quickly. In [19], a low-complexity angle tracking algorithm in monostatic MIMO radar was proposed. The studies are all about target tracking for the monostatic MIMO radar.
Bistatic MIMO radar is different from monostatic MIMO radar and consists of a transmitting and receiving base. The corresponding direction of arrivals (DOA) and direction of departure (DOD) are not equal. Therefore, the joint steering vector is more complicated. The above algorithms in [17,18,19] cannot solve the problem of target tracking in bistatic MIMO radar.
In [20], the PASTd algorithm in the array-signal-processing was introduced to the bistatic MIMO radar, and the target tracking problem of the MIMO radar was successfully solved. However, it requires an additional data correlation operation and cannot track the targets of the same DOD or DOA. In order to solve the deficiency of [20], reference [21] proposed a low complexity tracking algorithm in bistatic MIMO radar. The algorithm deduces the formula of the covariance matrix difference of an adjacent moment, and then achieves the target angle tracking. However, the algorithm from [21] uses only partial covariance matrix information, and the tracking performance is low. A target tracking algorithm based on Adaptive Asymmetric Joint Diagonalization (AAJD) was proposed in [22]. The algorithm does not need an additional correlation operation and can track the targets whose DOD or DOA is the same. Nevertheless, the performance of the algorithm is reduced by reusing the estimation angle of the last time.
In this paper, we consulted the monostatic MIMO radar angle tracking idea in [19] to propose a DOD and DOA tracking algorithm suitable for bistatic MIMO radar. The covariance matrix, constructed by the monostatic MIMO radar signal, satisfied the Toeplitz form in [19]. However, the DOD and DOA were different and the covariance matrix did not satisfy the Toeplitz structure in the bistatic MIMO radar. Therefore, we proposed an improved DOA and DOD tracking method based on three approximations. Then, we found that the covariance matrix satisfied the approximate Toeplitz property, whose partial elements in a straight line paralleled to the principal diagonal were equal. On this basis, our algorithm used the approximate Toeplitz property to take the average operation, which was equivalent to improve the signal-to-noise ratio (SNR). The proposed algorithm could realize automatic matching and association of DOD and DOA. Error analysis was also derived in this paper. Finally, the simulation results were presented to verify the effectiveness of the proposed algorithm.
There were some differences between the algorithm in [19] and the proposed algorithm. (1) Reference [19] proposed a low-complexity angle tracking algorithm in monostatic MIMO radar. However, we solved the DOD and DOA tracking problems of the bistatic MIMO radar. (2) The monostatic MIMO radar only needed to estimate DOA, but the algorithm in this paper needed to solve the DOD and DOA estimation, therefore the problem is more complicated. In this paper, the algorithm extended the approximate idea in [19] and solved the target angle by three approximation operations. (3) The monostatic MIMO radar angle tracking algorithm requires that the covariance matrix satisfy the Toeplitz form. Since the DOD and DOA are different, the covariance matrix is more complicated and does not satisfy the Toeplitz form in bistatic MIMO radar. Therefore, the algorithm in [19] cannot be used to solve the problem directly in this paper. So we used the approximate Toeplitz properties to improve the tracking performance.
Thus, the algorithm cannot only be seen as an extension of the work in [19] but is also an improved algorithm. The simulation results showed that the proposed algorithm had a better tracking performance than the angle tracking algorithm in [19,20,21,22].
The rest of this paper is organized as follows. In Section 2, the signal model of bistatic MIMO radar is presented. Section 3 establishes our angle tracking algorithm based on the elements of the covariance matrix of the receive signal. Section 4 compares the performance of the algorithm in [19,20,21,22] and our algorithm. The simulation results verify the effectiveness of the proposed algorithm. Finally, Section 5 concludes the paper.
Notations: ( ) T , ( ) H , ( ) + , and ( ) 1 denote the transpose, Hermitian transpose, pseudoinverse and inverse operations, respectively. I K is an K × K identity matrix; vec ( ) is the vectorization of a matrix; diag ( υ ) stands for diagonal matrix whose diagonal is a vector v ; and are the Kronecker product and Hadamard product, respectively.

2. Signal Model

In this paper, a bistatic MIMO radar was used to observe the moving targets in the air. The distances between the targets and bases are far, so the target satisfies the point target model. The bistatic MIMO radar is composed of M transmit antennas and N receive antennas. The space between the transceiver antennas is the same and half of the wavelength. The configuration of the bistatic MIMO radar is shown in Figure 1.
It is assumed that there is a P far-field moving point target in the air, and the DOD and DOA at time t is [ ( φ t , 1 , θ t , 1 ) , ( φ t , 2 , θ t , 2 ) , ( φ t , P , θ t , P ) ] , respectively. The velocity of the target i is v i , and the angles between the moving direction and the DOD and DOA directions are φ i and θ i , respectively.
The transmit signal radiates to P targets, and the signal that arrives at the receive elements after scattering is
x ¯ ( t ) = i = 1 P a r ( θ t , i ) ε i a t T ( φ t , i ) s ( t τ i ) exp ( j ω i t ) + n ¯ ( t )
where ω i = 2 π v i ( cos θ i + cos φ i ) λ ; f i = v i ( cos θ i + cos φ i ) λ is the doppler shift. The receive steering vector is a r ( θ t , i ) = [ 1 , e j π sin θ t , i , , e j π ( N 1 ) sin θ t , i ] T ; and the transmit steering vector is a t ( φ t , i ) = [ 1 , e j π sin φ t , i , e j π ( M 1 ) sin φ t , i ] T . ε 1 , , ε P is the scattering coefficient of the observation target and satisfies the Swerling II model, which is invariable within the pulse time. s ( t τ i ) = [ s 1 ( t τ i ) exp [ j 2 π f c ( t τ i ) ] , , s M ( t τ i ) exp [ j 2 π f c ( t τ i ) ] ] T , n ¯ ( t ) is Gaussian additive white noise.
The signal x ¯ ( t ) is filtered through a set of matched filters, let each filter match only one transmit signal. Let the delay of matched filter of i th target signal is τ i , and τ i = τ i .
The output of the signal in Equation (2) after matched filtering is
x ( t ) = i = 1 P a r ( θ t , i ) ε i a t T ( φ t , i ) s ( t τ i ) exp ( j ω i t ) s ( t τ i ) + n ¯ ( t ) s ( t τ i ) = i = 1 P a r ( θ t , i ) ε i exp ( j ω i t ) a t T ( φ t , i ) + n ( t )
ε i exp ( j ω i t ) and ε i has the same statistical properties, ε = [ ε 1 exp ( j ω 1 t ) , , ε P exp ( j ω P t ) ] T still satisfy the Swerling II model. The mean and variance of n ( t ) is 0 and σ 2 .
Further simplify Equation (3),
x ( t ) = A r ( θ ) d i a g ( ε ) A t T ( φ ) + n ( t )
Consider the vectorization of x ( t ) in (4),
y ( t ) = A t ( φ ) A r ( θ ) vec ( d i a g ( ε ) ) + n ( t ) = W t ( φ , θ ) ε + n ( t )
where W t ( φ , θ ) = [ a r ( θ t , 1 ) a t ( φ t , 1 ) , a r ( θ t , 2 ) a t ( φ t , 2 ) , , a r ( θ t , P ) a t ( φ t , P ) ] denotes M N × P joint steering vector.
Consider
R t = E [ y ( t ) y H ( t ) ] = W t R ε W t H + R n ( t )
where R ε = E [ ε ( t ) ε H ( t ) ] = d i a g ( | ε 1 | 2 | ε 2 | 2 | ε P | 2 ) , R n ( t ) = E [ n ( t ) n ( t ) H ] .
Assuming that DOD and DOA change slowly, and consider that the DOD and DOA of the target are the same in the time interval [ ( k 1 ) T s , k T s ] , ε i and ε j are uncorrelated for the different targets and all targets are in the same range bin. During the interval [ ( k 1 ) T s , k T s ] , φ t , P , θ t , P remains constant and L snapshots of sensor data are available for the signal processing.

3. Angle Tracking Algorithm Description

The angle tracking algorithm in [19] requires that the steering vector satisfies the Vandermonde form. The joint steering vector in Equation (4) does not satisfy the Vandermonde form, so the angle tracking algorithm in [19] cannot be applied directly to bistatic MIMO radar. We improved the algorithm in [19] and proposed an angle tracking algorithm suitable for bistatic MIMO radar.

3.1. Estimation of the Covariance Matrix Difference and the Angle Difference

At the time t , the DOD and DOA of the P targets are recorded as γ t = [ φ t , 1 , φ t , 2 , , φ t , P , θ t , 1 , θ t , 2 , , θ t , P ] . Similarly, the DOD and DOA at time t + 1 is recorded as γ t + 1 = [ φ t + 1 , 1 , φ t + 1 , 2 , , φ t + 1 , P , θ t + 1 , 1 , θ t + 1 , 2 , , θ t + 1 , P ] .
We define
Δ γ t = γ t + 1 γ t
where Δ γ t = [ Δ φ t , 1 , Δ φ t , 2 , Δ φ t , P , Δ θ t , 1 , Δ θ t , 2 , Δ θ t , P ] is the angle difference between t and t + 1 , Δ φ t , i = φ t + 1 , i φ t , i and Δ θ t , i = θ t + 1 , i θ t , i .
We define R t + 1 as the covariance matrices of the signal at time t + 1 . The covariance matrices are
R t + 1 = E [ y ( t + 1 ) y H ( t + 1 ) ] = W t + 1 R ε W t + 1 H + R n ( t + 1 )
then we can obtain
Δ R t = R t + 1 R t = ( W t + 1 R ε W t + 1 H W t R ε W t H ) + ( R n ( t + 1 ) R n ( t ) ) .
Supposing that the noise covariance matrix at time t + 1 is approximately equal to that at time t , then we have
Δ R t W t + 1 R ε W t + 1 H W t R ε W t H .
It can be seen that the covariance matrix difference of adjacent moment is caused by the angle difference of adjacent moment, so there is a relationship between the two. Therefore, by deriving the relationship between the two, the angle difference of adjacent moment can be obtained.

3.2. Estimation of DOD and DOA

In Section 3.1, we obtained the covariance matrix difference Δ R t and the angle difference Δ γ t . This section will deduce the linear relationship between the covariance matrix difference and the angle difference.
We first analyzed the properties of the elements of the covariance matrix. The literature [19] proved that Δ R t in Equation (9) can be expressed as
Δ R t = [ 0 b 1 , 0 b M 1 , 0 b 0 , 1 b M 1 , 1 b M 1 , N 1 b 1 , 0 * 0 0 b M 1 , 0 * 0 b 0 , 1 * 0 0 b M 1 , 1 * 0 0 b M 1 , N 1 * 0 ]
where b m , n = R ( 1 , m + n * M + 1 ) , b m , n = R ( m + n * N + 1 , 1 ) .
b m , n = i = 1 P ρ i ( e j π n sin ( θ t , i + Δ θ t , i ) e j π m sin ( φ t , i + Δ φ t , i ) e j π n sin θ t , i e j π m sin φ t , i )
where ρ i is the ( i , i ) element of the matrix R ε and ρ i = | ε i | 2 , m = 0 , 1 , M 1 ; n = 0 , 1 , N 1 . b m , n in Equation (11) can be expanded as
b m , n = i = 1 P ρ i ( e j π n [ sin θ t , i cos Δ θ t , i + cos θ t , i sin Δ θ t , i ] e j π m [ sin φ t , i cos Δ φ t , i + cos φ t , i sin Δ φ t , i ] e j π n sin θ t , i e j π m sin φ t , i ) .
From Equation(12), it can be see that b m , n is related to the angle difference and the angle of the previous moment. Δ θ t , i and Δ φ t , i are the parameters to be estimated. Considering that Δ θ t , i and Δ φ t , i is very small,
sin ( θ t , i + Δ θ t , i ) = sin θ t , i cos Δ θ t , i + cos θ t , i sin Δ θ t , i sin θ t , i + Δ θ t , i cos θ t , i
sin ( φ t , i + Δ φ t , i ) = sin φ t , i cos Δ φ t , i + cos φ t , i sin Δ φ t , i sin φ t , i + Δ φ t , i cos φ t , i
substituting Equations (13) and (14) into Equation(12), then b m , n in Equation (12) can be denoted as
b m , n i = 1 P ρ i ( e j π n [ sin θ t , i + Δ θ cos θ t , i ] e j π m [ sin φ t , i + Δ φ t , i cos φ t , i ] e j π n sin θ t , i e j π m sin φ t , i ) = i = 1 P ρ i ( e j π n sin θ t , i e j π m sin φ t , i [ e j π n Δ θ t , i cos θ t , i e j π m Δ φ t , i cos φ t , i 1 ] ) .
Considering that x is very small, e x 1 x . Then, b m , n can be rewritten as
b m , n i = 1 P ρ i e j π n sin θ t , i e j π m sin φ t , i [ j π n Δ θ t , i cos θ t , i j π m Δ φ t , i cos φ t , i ] .
By Equation(16), we can construct the following equation
V t Δ γ t = b
where
b = [ b 1 , 0 b 2 , 0 b M 1 , 0 b 0 , 1 b M 1 , 1 b M 1 , N 1 b 1 , 0 * b 2 , 0 * b M 1 , 0 b 0 , 1 b M 1 b M 1 , N 1 * ] T
Δ γ t = [ Δ φ t , 1 Δ φ t , 2 Δ φ t , P Δ θ t , 1 Δ θ t , 2 Δ θ t , P ] T
To give the V t , we define ξ φ i = e j π sin φ i , ξ θ i = e j π sin θ i , β φ i = j π cos φ i , β θ i = j π cos θ i .
V t = [ ρ 1 ξ φ 1 β φ 1 ρ P ξ φ P β φ P 0 0 ρ 1 ξ φ 1 2 2 β φ 1 ρ P ξ φ P 2 2 β φ P 0 0 ρ 1 ξ φ 1 M 1 ( M 1 ) β φ 1 ρ P ξ φ P M 1 ( M 1 ) β φ P 0 0 0 0 ρ 1 ξ θ 1 β θ 1 ρ P ξ θ P β θ P ρ 1 ξ φ 1 M 1 ξ θ 1 ( M 1 ) β φ 1 ρ P ξ φ P M 1 ξ θ P ( M 1 ) β φ P ρ 1 ξ φ 1 M 1 ξ θ 1 β θ 1 ρ P ξ φ P M 1 ξ θ P β θ P 0 0 ρ 1 ξ θ 1 N 1 ( N 1 ) β θ 1 ρ P ξ θ P N 1 ( N 1 ) β θ P ρ 1 ξ φ 1 M 1 ξ θ 1 N 1 ( M 1 ) β φ 1 ρ P ξ φ P M 1 ξ θ P N 1 ( M 1 ) β φ P ρ 1 ξ φ 1 M 1 ξ θ 1 N 1 ( N 1 ) β θ 1 ρ P ξ φ P M 1 ξ θ P N 1 ( N 1 ) β θ P ρ 1 ξ φ 1 1 ( β φ 1 ) ρ P ξ φ P 1 ( β φ P ) 0 0 ρ 1 ξ φ 1 2 ( 2 ) β φ 1 ρ P ξ φ P 2 ( 2 ) β φ P 0 0 ρ 1 ξ φ 1 ( M 1 ) ( ( M 1 ) ) β φ 1 ρ P ξ φ P ( M 1 ) ( ( M 1 ) ) β φ P 0 0 0 0 ρ 1 ξ θ 1 β θ 1 ρ P ξ θ P β θ P ρ 1 ξ φ 1 ( M 1 ) ξ θ 1 1 ( ( M 1 ) ) β φ 1 ρ P ξ φ P ( M 1 ) ξ θ P 1 ( ( M 1 ) ) β φ P ρ 1 ξ φ 1 ( M 1 ) ξ θ 1 1 β θ 1 ρ P ξ φ P ( M 1 ) ξ θ P 1 β θ P 0 0 ρ 1 ξ θ 1 ( N 1 ) ( ( N 1 ) ) β θ 1 ρ P ξ θ P ( N 1 ) ( ( N 1 ) ) β θ P ρ 1 ξ φ 1 ( M 1 ) ξ θ 1 ( N 1 ) ( ( M 1 ) ) β φ 1 ρ P ξ φ P ( M 1 ) ξ θ P ( N 1 ) ( ( M 1 ) ) β φ P ρ 1 ξ φ 1 ( M 1 ) ξ θ 1 ( N 1 ) ( ( N 1 ) ) β θ 1 ρ P ξ φ P ( M 1 ) ξ θ P ( N 1 ) ( ( N 1 ) ) β θ P ]
Using the least square method to estimate Δ γ t in Equation (17), we get
Δ γ t = ( V t H V t ) 1 V t H b .
The final estimate of the angle is
γ t + 1 = γ t + Δ γ t

3.3. Covariance Element Average Operation

The algorithm in [19] makes full use of the Toeplitz matrix property to improve performance. Since the joint steering vector of bistatic MIMO radar does not satisfy the Vandermonde form, the covariance matrix difference does not satisfy the Toeplitz property. Therefore, we first analyzed the structure of the steering vector.
W t ( φ , θ ) = [ 1 , e j π sin φ t , i , e j π ( M 1 ) sin φ t , i , e j π sin θ t , i , e j π sin θ t , i e j π sin φ t , i , , e j π sin θ t , i e j π ( M 1 ) sin φ t , i , ,          e j π ( N 1 ) sin θ t , i , e j π ( N 1 ) sin θ t , i e j π sin φ t , i , , e j π ( N 1 ) sin θ t , i e j π ( M 1 ) sin φ t , i ] T .
It can be found that some elements in the steering vector satisfied the Vandermonde structure in Equation (23). This structure is called the approximate Vandermonde structure in this paper.
We further analyzed the structure of the covariance matrix, taking R t as an example.
R t = [ r 1 , 1 r 1 , 2 r 1 , M r 1 , M + 1 r 1 , M + 2 r 1 , 2 M r 1 , ( N 1 ) M + 1 r 1 , ( N 1 ) M + 2 r 1 , N M r 1 , 2 r 1 , 1 r 1 , M 1 r 1 , M r 1 , M 1 r 1 , 1 r 1 , M + 1 r 1 , 1 r 1 , 2 r 1 , M r 1 , M + 2 r 1 , 2 r 1 , 1 r 1 , M 1 r 1 , 2 M r 1 , M r 1 , M 1 r 1 , 1 r 1 , ( N 1 ) M + 1 r 1 , 1 r 1 , 2 r 1 , M r 1 , ( N 1 ) M + 2 r 1 , 2 r 1 , 1 r 1 , M 1 r 1 , 1 r 1 , N M r 1 , M r 1 , M 1 r 1 , 1 ]
Since the steering vector satisfied the approximate Vandermonde structure, the sub-matrices in a straight line parallel to the principal diagonal of R t are the same. R t + 1 has a similar structure.
From the above analysis, we can see that Δ R t is an approximate Toeplitz matrix whose sub-matrices in a straight line parallel to the principal diagonal are the same. We used the average operation to estimate b m , n and b m , n to eliminate the noise. The following steps can be used to update the b m , n and b m , n .
b ^ m , n = 1 ( M m ) ( N n ) k = 1 M m k k = n + 1 N Δ R t ( k + ( k k n 1 ) M , k + m + ( k k 1 ) M ) b ^ m , n = 1 ( M m ) ( N n ) k = 1 M m k k = n + 1 N Δ R t ( k + m + ( k k 1 ) M , k + ( k k n 1 ) M )
where Δ R t ( i , j ) is the ( i , j ) element of the matrix Δ R t . Substituting Equation (25) into Equation (18), then b can be rewritten as
b = [ b ^ 1 , 0 b ^ 2 , 0 b ^ M 1 , 0 b ^ 0 , 1 b ^ M 1 , 1 b ^ M 1 , N 1 b ^ 1 , 0 * b ^ 2 , 0 * b ^ M 1 , 0 b ^ 0 , 1 b ^ M 1 b ^ M 1 , N 1 * ] T
The proposed algorithm fully uses the approximate Toeplitz matrix property and more receiving information to eliminate the noise, and improve angle tracking performance. Now, we discuss the performance comparison between our algorithm and the algorithms in the literature [19,20,21,22].
The AAJD and PASTd algorithm sestimate the target angle through optimizing the function. Because it is difficult to find the optimal solution of the optimization function, the performance of the algorithm is low. The algorithm in [19,21] and our algorithm obtain angle via the difference between the previous and current covariance matrix of the receiving signal. So the performance of the algorithm in [19,21] and our algorithm is better than the AAJD and PASTd algorithms.
The algorithm in [19] can be improved to solve the angle tracking problem of bistatic MIMO radar, but it only uses the ( M 2 + M ) ( N 1 ) + M 2 M elements of the covariance matrix ( M is the number of transmit antennas, and N is the number of receive antennas). The algorithm in [21] and our algorithm uses the 2 ( M N 1 ) and ( M 2 + M ) ( N 2 N ) 2 + ( M 2 M ) N elements of the covariance matrix to track the target, respectively. Owning to ( M 2 + M ) ( N 2 N ) 2 + ( M 2 M ) N > ( M 2 + M ) ( N 1 ) + M 2 M > 2 ( M N 1 ) , the proposed algorithm used more covariance matrix information than the algorithm in [19,21]. So, the performance of our algorithm was better than that of [19,21].
In summary, the performance of our algorithm was the best.
Until now, we show the major steps of the angle tracking algorithm in bistatic MIMO radar as follows.
Step 1. Calculate the covariance matrix R t and R t + 1 via Equations (5) and (7).
Step 2. Calculate the covariance matrix difference Δ R t via Equation (10).
Step 3. The vector b is obtained via Equation (26), and the vector V t is obtained via Equation (20).
Step 4. We estimate Δ γ t via Equation (21), and the angle at time t + 1 is γ t + 1 = γ t + Δ γ t = γ 1 + i = 1 t Δ γ i .
Step 5. Repeat steps 1 to 4 to estimate the angle of the next moment.
Note 1.
This paper assumes that the number of targets in the bistatic MIMO radar is known. If we do not know the number in advance, we can use the existing target-number estimation algorithm in [23] to estimate the number of targets.
Note 2.
The algorithm in this paper only obtains the angle difference of adjacent time, and therefore we need to get the initial DOD and DOA of the target. The initial DOD and DOA can be obtained using the MUSIC algorithm or another angle estimation algorithm.
Note 3.
The algorithm in this paper was valid effectively when the target velocity was low and the DOD and DOA changed slowly. When the target moves faster, the performance of the algorithm in this paper will be reduced or even invalidated. When the target was far from the transceiver base, the angle difference was generally small, so this algorithm is suitable for tracking the long-distance target.
Note 4.
The noise covariance matrices of adjacent time can be approximately assumed to be equal. No matter what kind of noise, the noise component in Equation (8) can be eliminated. The algorithm is still effective under the colored noise conditions.

3.4. Computational Complexity Analysis and Advantages of the Proposed Algorithm

For the proposed algorithm, the calculation of the covariance matrix needs O ( ( M N ) 2 L ) , and the computation of ( V t H V t ) 1 V t H b requires O ( 2 P 2 ( M N 1 ) + P 3 + 2 P ( M N 1 ) + P 2 ) . The main computational complexity of the proposed algorithm is O ( ( M N ) 2 L + 2 P 2 ( M N 1 ) + P 3 + 2 P ( M N 1 ) + P 2 ) .
The advantages of this algorithm are listed as follows:
(1)
The proposed algorithm does not need eigenvalue decomposition of the covariance matrix, so the complexity is lower.
(2)
The proposed algorithm not only introduces the tracking algorithm in [19], but also expands it.
(3)
The proposed algorithm makes full use of the elements in the covariance matrix to improve the tracking performance. The performance of this algorithm is better than the AAJD algorithm.
(4)
The algorithm in this paper can automatically match and associate the angles of adjacent moment and reduce the computational complexity.

4. Error Analysis

In this section, we deduce the variance of DOD and DOA tracking. We assume that the observed noise variances are nearly the same at the adjacent time. When estimating DOD and DOA, we used approximate calculations such as e x 1 x and sin x x when x was smaller. This leads to a slight difference from the real value. Consider
sin x = x x 3 3 ! + Λ
cos x = 1 x 2 2 ! + Λ
e x 1 = x + Λ
where Λ , Λ , and Λ are the high-order expansion terms.
According to Equations (13), (14), (27) and (28), we have
b m , n = i = 1 P ρ i ( e j π n [ sin θ t , i cos Δ θ t , i + cos θ t , i sin Δ θ t , i ] e j π m [ sin φ t , i cos Δ φ t , i + cos φ t , i sin Δ φ t , i ] e j π n sin θ t , i e j π m sin φ t , i ) = i = 1 P ρ i ( e j π n [ sin θ t , i ( 1 ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( Δ θ t , i ( Δ θ t , i 3 / 3 ! ) + Λ ) ] e j π m [ sin φ t , i ( 1 ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( Δ φ t , i ( Δ φ t , i 3 / 3 ! ) + Λ ) ] e j π n sin θ t , i e j π m sin φ t , i ) = i = 1 P ρ i ( e j π n [ sin θ t , i + cos θ t , i Δ θ t , i ] e j π n [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ] e j π m [ sin φ t , i + cos φ t , i Δ φ t , i ] e j π m [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] e j π n sin θ t , i e j π m sin φ t , i ) = i = 1 P ρ i e j π ( n sin θ t , i + m sin φ t , i ) ( e j π ( n cos θ t , i Δ θ t , i + m cos φ t , i Δ φ t , i ) e j π n [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ] e j π m [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] 1 ) = i = 1 P ρ i e j π ( n sin θ t , i + m sin φ t , i ) ( e j π ( n cos θ t , i Δ θ t , i + m cos φ t , i Δ φ t , i ) j π n [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ] j π m [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] 1 )
By Equations (15), (16), and (29), then
b m , n = i = 1 P ρ i e j π ( n sin θ t , i + m sin φ t , i ) ( j π ( n cos θ t , i Δ θ t , i + m cos φ t , i Δ φ t , i ) j π n [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ] j π m [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] + j Λ ) = i = 1 P ρ i e j π ( n sin θ t , i + m sin φ t , i ) j π ( n cos θ t , i Δ θ t , i + m cos φ t , i Δ φ t , i ) + i = 1 P ρ i e j π ( n sin θ t , i + m sin φ t , i ) ( j π n [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ] j π m [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] + j Λ )
b m , n is the estimation error of b m , n , and b m , n can be shown as follows:
b m , n = i = 1 P ρ i e j π ( n sin θ t , i + m sin φ t , i ) ( j π n [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ] j π m [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] + j Λ )
According to Equations (21) and (32), the variance of Δ b m , n is denoted by
E [ | b m , n | 2 ] = 1 ( M m ) ( N n ) i = 1 P ρ i 2 ( π n [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ]                 + π m [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] + Λ ) 2
Then, we have
E [ b m , n b r , s ] = 0 , m r , n s
E [ b m , n 2 ] = 0
We define that b is the estimation error of b .
b = [ b ^ 1 , 0 b ^ 2 , 0 b ^ M 1 , 0 b ^ 0 , 1 b ^ M 1 , 1 b ^ M 1 , N 1 b ^ 1 , 0 * b ^ 2 , 0 * b ^ M 1 , N 1 * ] T
Then the variance of b is
E [ | b | 2 ] = [ E [ | b 1 , 0 | 2 ] E [ | b M 1 , 0 | 2 ] E [ | b 0 , 1 | 2 ] E [ | b M 1 , 1 | 2 ] E [ | b M 1 , N 1 | 2 ] E [ | b 1 , 0 | 2 ] E [ | b M 1 , N 1 | 2 ] ] = [ 1 ( M 1 ) N i = 1 P ρ i 2 ( π [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] + Λ ) 2 1 N i = 1 P ρ i 2 ( π ( M 1 ) [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] + Λ ) 2 1 M ( N 1 ) i = 1 P ρ i 2 ( π [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ] + Λ ) 2 1 ( N 1 ) i = 1 P ρ i 2 ( π [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ] + π ( M 1 ) [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] + Λ ) 2 i = 1 P ρ i 2 ( π ( N 1 ) [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ] + π ( M 1 ) [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] + Λ ) 2 1 ( M 1 ) N i = 1 P ρ i 2 ( π [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] + Λ ) 2 i = 1 P ρ i 2 ( π ( N 1 ) [ sin θ t , i ( ( Δ θ t , i 2 / 2 ! ) + Λ ) + cos θ t , i ( ( Δ θ t , i 3 / 3 ! ) + Λ ) ] + π ( M 1 ) [ sin φ t , i ( ( Δ φ t , i 2 / 2 ! ) + Λ ) + cos φ t , i ( ( Δ φ t , i 3 / 3 ! ) + Λ ) ] + Λ ) 2 ]
According to [24], we get the variance of Δ γ t , i
var [ Δ γ t , i ] = E [ V t , i + b ] + Re ( E [ V t , i + b ] 2 ) 2
where V t , i + denotes the i th row of V t + and V t + is the pseudoinverse of V t . According to Equations (33)–(35) and (37), we get
var [ Δ γ t , i ] = V t , i + d i a g ( E [ | b | 2 ] V t , i + H + Re ( V t , i + E [ b b T ] V t , i + T ) ) 2 = V t , i + d i a g ( E [ | b | 2 ] ) V t , i + H 2
where E [ | b | 2 ] is shown in Equation (37).
From Equations (37) and (39), we can obtain an effective conclusion where the theoretical variance of the proposed algorithm is gradually decreased with the number of transmit/receive antennas increases. Multiple transmit/receive antennas improve the angle tracking performance.

5. Simulation Results

Assuming that both the transmit and receive arrays of the bistatic MIMO radar are linearly configured, the spacing of the array elements is half wavelength. The carrier frequency of the array element is 1   GHz , the pulse width is 10   μ s , and the pulse repetition rate is 10   kHz . The emission waveform uses the Hadamard Code Pulse (HCP) signal, and the number of the transmit and receive array elements is M = N = 5 (except experiments 6,7). We defined the root-mean square error (RMSE) as R M S E ( θ ) = 1 F m = 1 F 1 P k = 1 P 1 T t = 1 T [ ( θ ^ k , m , t θ k , m , t ) 2 + ( φ ^ k , m , t φ k , m , t ) 2 ] , where θ ^ k , m , t and φ ^ k , m , t is the estimate of angle θ k , m , t and φ k , m , t ; F are the times of the Monte Carlo trial. The targets are tracked over an interval of 6 s, during each 0.1 s interval, L snapshots of sensor data are generated and used to estimate angles.
Figure 2 shows the result of the tracking angle of the targets of uniform speed for P = 2 , L = 100 , and S N R = 15   dB . The simulation results showed that the estimated trajectory coincided with the real trajectory, which proved the effectiveness of the algorithm.
Figure 3 depicts the tracking result of the proposed algorithm for non-uniform moving targets and showed that our algorithm could successfully track the moving target at a non-uniform speed. The DOD and DOA could be automatically associated. The estimated angle trajectory coincided with the true trajectory, indicating the effectiveness and robustness of the proposed algorithm.
Figure 4 depicts the tracking result comparison between the proposed algorithm and the angle tracking algorithm in [19] with S N R = 10   dB . From Figure 4, it can be seen that the tracking result of our algorithm had a high degree of coincidence with the real trajectory of the target, and the performance was better than the tracking algorithm in [19].
To better verify the performance of the proposed algorithm, Figure 5 shows the tracking performance comparison with P = 2 , L = 100 , F = 200 , and S N R = 5 10   dB , where we compared the proposed algorithm against the angle tracking algorithm in [19,21], the PASTd algorithm in [20], and AAJD algorithm in [22]. At the same time, we gave the theoretical error caused by the approximate operation. It can be seen from Figure 5 that the RMSE of our algorithm was lower than that of the angle tracking algorithm in [19,21], PASTd algorithm, and AAJD algorithm, which showed that the tracking performance of our algorithm was the best and the correctness of the theoretical analysis was verified. This was because the proposed algorithm made full use of the elements of covariance matrix to eliminate the noise and improved the estimation performance. The theoretical variance was lower than the actual variance, because the noise was ignored when the theoretical variance is derived.
Figure 6 and Figure 7 display the performance of angle tracking via our algorithm in the condition of P = 2 , L = 100 , F = 200 , S N R = 10   dB , and variable numbers of M / N . It was clearly shown that the angle tracking performance of the proposed algorithm gradually improved with the increased number of transmit/receive antennas. Multiple transmit/receive antennas improved the angle tracking performance because of diversity gain. The correctness of the conclusions drawn from the theoretical error was validated.

6. Conclusions

In this paper, we proposed a moving multi-target angle tracking algorithm for bistatic MIMO radar. The proposed algorithm obtained the linear relationship between the covariance matrix difference and the angle difference through the three approximate processes. The proposed algorithm reduced the computational complexity and realized the automatic association of DOA and DOD. The proposed algorithm made full use of the elements of the covariance matrix by taking the average method, eliminating the noise, and improving the tracking performance. The research in this paper provides technical support for the practical application of the MIMO radar. In future work, we will analyze wideband signal processing to improve performance [25,26] and study signal processing in complex backgrounds to increase the robustness of the algorithm.

Author Contributions

Zhengyan Zhang and Jianyun Zhang participated in the design of this study. Zhengyan Zhang and Qingsong Zhou conceived and designed the experiments and performed the experiments; Zhengyan Zhang and Xiaobo Li analyzed the result of the experiments; Zhengyan Zhang wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fishler, E.; Haimovich, A.; Blum, R.; Chizhik, D.; Cimini, L.; Valenzuela, R. MIMO radar: An idea whose timehas come. In Proceedings of the IEEE Radar Conference, Philadelphia, PA, USA, 26–29 April 2004; pp. 71–78. [Google Scholar]
  2. Li, J.; Stoica, P. MIMO Radar with Colocated Antennas. IEEE Signal Process. Mag. 2007, 24, 106–114. [Google Scholar] [CrossRef]
  3. Haimovich, A.M.; Blum, R.S.; Cimini, L.J. MIMO Radar with Widely Separated Antennas. IEEE Signal Process. Mag. 2008, 25, 116–129. [Google Scholar] [CrossRef]
  4. Wang, X.P.; Huang, M.X.; Wu, X.Q.; Bi, G.A. Direction of arrival estimation for MIMO radar via unitary nuclear norm minimization. Sensor 2017, 17, 939. [Google Scholar] [CrossRef] [PubMed]
  5. Tang, B.; Tang, J. Joint Design of Transmit Waveforms and Receive Filters for MIMO Radar Space-Time Adaptive Processing. IEEE Trans. Signal Process. 2016, 64, 4707–4722. [Google Scholar] [CrossRef]
  6. 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]
  7. Wang, W.; Wang, X.; Song, H.; Ma, Y. Conjugate ESPRIT for DOA estimation in monostatic MIMO radar. Signal Process. 2013, 93, 2070–2075. [Google Scholar] [CrossRef]
  8. Yang, X.; Zheng, G.M.; Tang, J. ESPRIT algorithm for coexistence of circular and noncircular signals in bistatic MIMO radar. In Proceedings of the 2016 IEEE Radar Conference (RadarConf), Philadelphia, PA, USA, 2–6 May 2016; pp. 1–4. [Google Scholar]
  9. Li, J.F.; Jiang, D.F.; Zhang, X.F. DOA estimation based on combined unitary ESPRIT for coprime MIMO radar. IEEE Commun. Lett. 2017, 21, 96–99. [Google Scholar] [CrossRef]
  10. Zhang, X.; Huang, Y.; Chen, C.; Li, J.; Xu, D. Reduced complexity Capon for direction of arrival estimation in a monostatic multiple-input multiple-output radar. IET Radar Sonar Navig. 2012, 6, 796–801. [Google Scholar] [CrossRef]
  11. Belfiori, F.; Rossum, W.V.; Hoogeboom, P. Coherent MUSIC technique for range/angle information retrieval: Application to a frequency-modulated continuous wave MIMO radar. IET Radar Sonar Navig. 2014, 8, 75–83. [Google Scholar] [CrossRef]
  12. Tan, J.; Nie, Z.P.; Wen, D.B. Low complexity MUSIC-based direction-of-arrival for monostatic MIMO radar. Electron. Lett. 2017, 53, 275–277. [Google Scholar] [CrossRef]
  13. Zahernia, A.; Dehahani, M.J.; Javidan, R. MUSIC algorithm for DOA estimation using MIMO arrays. In Proceedings of the 2011 Sixth International Conference on Telecommunication Systems, Services, and Application, Bali, Indonesia, 20–21 October 2011; pp. 149–153. [Google Scholar]
  14. Devaney, A.J. Time reversal imaging of obscured targets from multistatic data. IEEE Trans. Antennas Propag. 2005, 53, 1600–1610. [Google Scholar] [CrossRef]
  15. Ciuonzo, D.; Romano, G.; Solimene, R. Performance analysis of time-reversal MUSIC. IEEE Trans. Signal Process. 2015, 63, 2650–2662. [Google Scholar] [CrossRef]
  16. Ciuonzo, D.; Salvo, R.P. Noncolocated time-reversal MUSIC: High-SNR distribution of null spectrum. IEEE Signal Process. Lett. 2017, 24, 397–401. [Google Scholar] [CrossRef]
  17. Wu, H.L.; Zhang, X.F. DOA tracking in monostatic MIMO radar using PARAFAC-RLST algorithm. In Proceedings of the 3rd International Conference on Information Science and Engineering (ICISE’11), Yangzhou, China, 29 September–1 October 2011; pp. 958–961. [Google Scholar]
  18. Zhang, X.F.; Li, J.F.; Feng, G.P.; Wu, H.L. Kalman-PASTd based DOA tracking algorithm for monostatic MIMO radar. In Proceedings of the International Conference on Information, Services and Management Engineering, Beijing, China, 26–28 December 2011; pp. 220–224. [Google Scholar]
  19. Yu, H.X.; Zhang, X.F.; Chen, X.Q.; Wu, H.L. Computationally efficient DOA tracking algorithm in monostatic MIMO radar with automatic association. Int. J. Antennas Propag. 2014, 2014. [Google Scholar] [CrossRef]
  20. Wu, H.; Zhang, X. DOD and DOA tracking algorithm for bistatic MIMO radar using PASTD without additional angles pairing. In Proceedings of the IEEE 5th International Conference on Advanced Computational Intelligence (ICACI’12), Nanjing, China, 18–20 October 2012; pp. 1132–1136. [Google Scholar]
  21. Zhang, Z.Y.; Li, X.B.; Xu, X.Y.; Huang, C. Target angle rapid tracking algorithm for bistatic MIMO radar. J. Signal Process. 2016, 32, 701–706. (In Chinese) [Google Scholar]
  22. Zhang, W.T.; Lou, S.T.; Li, X.J. Tracking multiple target in MIMO radar Via adaptive asymmetric joint diagonalization. IEEE Trans. Signal Process. 2016, 64, 2880–2893. [Google Scholar] [CrossRef]
  23. Xin, J.M.; Zheng, N.N.; Sano, A. Simple and efficient nonparametric method for estimating the number of signals without eigendecomposition. IEEE Trans. Signal Process. 2007, 55, 1405–1420. [Google Scholar] [CrossRef]
  24. 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]
  25. Shi, G.; Nehorai, A. Maximum Likelihood Estimation of Point Scatterers for Computational Time-reversal Imaging. Commun. Inf. Syst. 2005, 5, 227–256. [Google Scholar]
  26. Ciuonzo, D. On Time-Reversal Imaging by Statistical Testing. IEEE Signal Process. Lett. 2017, 24, 1024–1028. [Google Scholar] [CrossRef]
Figure 1. Bistatic multiple-input multiple-output (MIMO) radar transceiver element configuration.
Figure 1. Bistatic multiple-input multiple-output (MIMO) radar transceiver element configuration.
Sensors 18 00805 g001
Figure 2. Angle tracking results of uniform moving target at S N R = 15   dB : (a) The direction of departure (DOD); (b) The direction of arrival (DOA); (c) DOD and DOA trajectory.
Figure 2. Angle tracking results of uniform moving target at S N R = 15   dB : (a) The direction of departure (DOD); (b) The direction of arrival (DOA); (c) DOD and DOA trajectory.
Sensors 18 00805 g002
Figure 3. Angle tracking results of non-uniform moving targets with S N R = 15   dB : (a) The direction of departure; (b) The direction of arrival; (c) DOD and DOA trajectory.
Figure 3. Angle tracking results of non-uniform moving targets with S N R = 15   dB : (a) The direction of departure; (b) The direction of arrival; (c) DOD and DOA trajectory.
Sensors 18 00805 g003aSensors 18 00805 g003b
Figure 4. Angle tracking result comparison with S N R = 10   dB : (a) The direction of departure; (b) The direction of arrival; (c) DOD and DOA trajectory comparison.
Figure 4. Angle tracking result comparison with S N R = 10   dB : (a) The direction of departure; (b) The direction of arrival; (c) DOD and DOA trajectory comparison.
Sensors 18 00805 g004
Figure 5. RMSE changes with signal-to-noise ratio (SNR).
Figure 5. RMSE changes with signal-to-noise ratio (SNR).
Sensors 18 00805 g005
Figure 6. Angle tracking with different values of M t .
Figure 6. Angle tracking with different values of M t .
Sensors 18 00805 g006
Figure 7. Angle tracking with different values of M r .
Figure 7. Angle tracking with different values of M r .
Sensors 18 00805 g007

Share and Cite

MDPI and ACS Style

Zhang, Z.; Zhang, J.; Zhou, Q.; Li, X. Multi-Target Angle Tracking Algorithm for Bistatic Multiple-Input Multiple-Output (MIMO) Radar Based on the Elements of the Covariance Matrix. Sensors 2018, 18, 805. https://doi.org/10.3390/s18030805

AMA Style

Zhang Z, Zhang J, Zhou Q, Li X. Multi-Target Angle Tracking Algorithm for Bistatic Multiple-Input Multiple-Output (MIMO) Radar Based on the Elements of the Covariance Matrix. Sensors. 2018; 18(3):805. https://doi.org/10.3390/s18030805

Chicago/Turabian Style

Zhang, Zhengyan, Jianyun Zhang, Qingsong Zhou, and Xiaobo Li. 2018. "Multi-Target Angle Tracking Algorithm for Bistatic Multiple-Input Multiple-Output (MIMO) Radar Based on the Elements of the Covariance Matrix" Sensors 18, no. 3: 805. https://doi.org/10.3390/s18030805

APA Style

Zhang, Z., Zhang, J., Zhou, Q., & Li, X. (2018). Multi-Target Angle Tracking Algorithm for Bistatic Multiple-Input Multiple-Output (MIMO) Radar Based on the Elements of the Covariance Matrix. Sensors, 18(3), 805. https://doi.org/10.3390/s18030805

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