International Journal of Antennas and Propagation
Volume 2018 (2018), Article ID 4370637, 7 pages
https://doi.org/10.1155/2018/4370637

Research Article

Joint DOA and DOD Estimation Based on Tensor Subspace with Partially Calibrated Bistatic MIMO Radar

Junxiang Wang,1,2 Ping Huang,¹ and Dingjie Xu²

¹College of Automation, Harbin Engineering University, Harbin, Heilongjiang 150001, China
²College of Electrical Engineering and Automation, Harbin Institute of Technology University, Harbin 150001, China

Correspondence should be addressed to Junxiang Wang; wangjunxiang@hrbeu.edu.cn

Received 2 August 2017; Accepted 18 December 2017; Published 22 April 2018

Academic Editor: Piotr Samczynski

Copyright © 2018 Junxiang Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A joint direction-of-departure (DOD) and direction-of-arrival (DOA) estimation algorithm based on tensor subspace approach for partially calibrated bistatic multiple-input multiple-output (MIMO) radar is proposed. By exploiting the multidimensional structure of the received data, a third-order measurement tensor is constructed. Consequently, the tensor-based signal subspace is achieved using the higher-order singular value decomposition (HOSVD). To achieve accurate DOA estimation with partially calibrated array, a closed-form solution is provided to estimate the gain-phase uncertainties of the transmit and receive arrays by modeling the imperfections of the arrays. Simulation results demonstrate the effectiveness of the proposed calibration algorithm.

1. Introduction

Array signal processing, which includes parameter estimation [1, 2] and beamforming [3, 4] has been widely investigated in past decades [5] which has found various applications in radar, sonar, and satellite navigation system. In recent years, multiple-input multiple-output (MIMO) radar [6, 7] has drawn increasing attention due to its significant advantages, such as higher angular resolution and more degrees-of-freedom (DOFs). MIMO radar using multiple antennas to simultaneously emit orthogonal waveforms and multiple antennas to receive the reflected signals has recently drawn increasing attention owing to its superior performance over to the conventional phased array. Direction-of-departure (DOD) and direction-of-arrival (DOA) estimation is an important objective of bistatic MIMO radar system. Various algorithms for DOD and DOA estimation have been proposed in literatures [8–12]. In [8], the rotational invariance property of the MIMO radar is exploited which enables the DOD and DOA estimation with bistatic MIMO radar. However, this method determines the DODs and DOAs using two independent estimation of signal parameters via rotational invariance technique (ESPRIT) and requires an extra procedure to achieve the paired DOD and DOA estimation. To avoid the pair matching procedure, the relationship between the two independent ESPRIT is exploited in [9] to estimate the DODs and DOAs. In [10], a reduced-dimension MUltiple SIgnal Classification (MUSIC) algorithm is proposed which only requires one-dimensional peak search. This method achieves a similar performance compared with the two-dimensional MUSIC algorithm and has a lower computational complexity. Based on the double polynomial root-finding procedure, a polynomial root-finding technique for joint DOD and DOA estimation is proposed in [11]. This algorithm allows an efficient estimation with automatic pairing. In [12], the property of the Kronecker product is utilized to formulate a cost function in the quadratic form. Then the DODs can be estimated from the constraint function. The DOAs can be achieved from the eigenvector which is related to the corresponding DODs.

However, the algorithms mentioned above are based on the assumption that arrays are well calibrated which is not true in practice. Therefore, the above algorithms will suffer from performance degradation or even fail to achieve accurate DOD and DOA estimation in many scenarios. A number of algorithms have been proposed to deal with the array calibration problem [13–17]. An iterative algorithm based on the MUSIC technique is proposed in [13], which can simultaneously achieve the angle and gain-phase uncertainties estimation. Due to the iterative procedure, this algorithm has a high computational complexity. In addition, this method is based on the assumption that the gain-phase uncertainties are small, otherwise it will suffer from suboptimal convergence. The algorithms proposed in [15, 16] do not require iteration and exploit the ESPRIT-like technique to achieve accurate angle estimations. To achieve higher DOFs, the method proposed in [17] utilizes the property of the quasistationary signals to solve the underdetermined DOA estimation problem with partly calibrated arrays. Moreover, a kind of calibration methods based on the sparse reconstruction is discussed in [18, 19].

However, the multiple-dimensional property of the received signal is not exploited in the abovementioned algorithms. The inherently multiple-dimensional structure of the received signal can be utilized by expressing the received signal as a third-order tensor [20]. Thus, the higher-order singular value decomposition (HOSVD) can be used to achieve a tensor-based signal subspace which improves the estimation performance especially with low signal to noise ratio (SNR).

In this paper, we exploit a calibration algorithm based on tensor approach for partially calibrated bistatic MIMO radar. We first exploit the multidimensional structure of the received data, and then, the tensor-based signal subspace is calculated using HOSVD. By formulating the gain-phase imperfections in the generalized manifold matrix, a closed-form solution is achieved to estimate the gain-phase uncertainties of the transmit and receive arrays. With the estimated gain-phase uncertainties, the arrays can be calibrated and the DODs and DOAs of the sources can be estimated accurately. Simulation results demonstrate the effectiveness of the proposed joint DOD and DOA estimation algorithm for partially calibrated bistatic MIMO radar.

The rest of this paper is organized as follows. In Section 2, we briefly introduce the tensor algebra basics and the tensor-based signal model of the bistatic MIMO radar. In Section 3, we first describe the partially calibrated MIMO radar model. Then, the proposed method is introduced. Numerical simulation results are presented in Section 4. The conclusions are made in Section 5.

Notations used in this paper are as follows: lowercase (uppercase) bold characters are used to denote vectors (matrices). stands for the entry of . denotes the identity matrix. and denote the transpose and Hermitian operations of a matrix or vector, respectively. stands for the vectorization operator. denotes a diagonal matrix with the elements of constituting the diagonal entries. denotes the Kronecker product, denotes the Hadamard product, and denotes the Khatri-Rao product. In addition, denotes the -norm, whereas and express an complex matrix or vector set.

2. Tensor Algebra Basic and Signal Model

In this section, we first introduce tensor operations used in this paper and provide the bistatic MIMO radar signal model.

2.1. Tensor Algebra Basic

In this subsection, the tensor algebra basic is introduced for readers’ convenience [20, 21].

Definition 1 (-mode unfolding). Denote a -order tensor as , the corresponding -mode unfolding is to convert into a matrix which is expressed as . The entry of the tensor maps to the entry of , where

Definition 2 (-mode tensor product). The -mode matrix product of the tensor with a matrix is denoted by with size and

Moreover, the properties of the mode tensor product are provided below:

Definition 3 (HOSVD). The HOSVD of the tensor can be written as

where is the core tensor of which can be calculated by

In addition, , , is a unitary , which is the left singular matrix of .

2.2. Bistatic MIMO Radar Signal Model

Consider a bistatic MIMO radar system where the transmit and receive arrays are both uniform linear arrays (ULAs) with and colocated antennas, respectively. In addition, both the transmit and receive antennas are arranged with half-wavelength interelement space. The transmit antennas simultaneously emit orthogonal waveforms where is the number of samples during pulse period. Assume that there are uncorrelated narrowband far-field source signals in the range bin of interest. The baseband received signal at the pulse can be expressed as:

where and are the receive and transmit manifold matrix, respectively. and are the steering vectors of the source, respectively. is the pair of the DOA and DOD of the source. with accounts for the Doppler effect and the radar reflection coefficient fading which follows the Swerling I model [22]. Finally, is zero mean Gaussian noise term. Under the fact that the transmitted waveforms are orthogonal to each other, that is,

Then, the received signal can be processed by matched filtering using . Thus, the output of the matched filter can be written as

where is additive zero mean Gaussian noise term with unknown variance after matched filtering.

Usually, the received signal matrix can be achieved by vectorizing and stacking the pulses received signal, that is,

where and . However, this expression ignores the multidimensional structure inherent of the received signals. In order to exploit the multidimensional structure of the received signal, we collect pulses and gather as a third-order tensor given by

where is the identity tensor. is the noise tensor. According to Definition 1, the relationship between (8) and (10) can be expressed as .

3. Proposed Method

In this section, the partially calibrated MIMO radar signal model is first provided. Then, a robust DOA estimation algorithm based on tensor algebra is presented.

3.1. Partially Calibrated MIMO Radar Model

Consider the case where only a part of the transmit and receive antennas is well-calibrated. Without loss of generality, we assume that the first and antennas of the transmit and receive arrays are calibrated, respectively, whereas the rest of the sensors are uncalibrated with unknown gain-phase uncertainties. Moreover, we assume that the uncertainties are direction-independent. Thus, the transmit manifold matrix and the receive manifold matrix can be, respectively, expressed as

where and with gain-phase uncertainties and , respectively. Thus, the received signal matrix can be expressed as

For the case of partially calibrated array, the conventional DOD and DOA estimation algorithm will suffer from severe performance degradation. In order to achieve satisfactory angle estimation result, the partially calibrated problem should be concerned.

3.2. ESPRIT-Like Calibration Method

The HOSVD of the received signal (12) can be written as

where , , , and . Denote the signal subspace of as

where , contains the first dominant singular vectors of with . Thus, according to Definition 3, we have

Then, the signal subspace in matrix form can be calculated from the -mode unfolding of , that is, . Note that the signal subspace expands the same space with the manifold matrix, that is,

where is a full rank matrix.

The generalized manifold matrix can be rewritten as

where . Thus, by selecting the first rows of , we have

where and are the first and last rows of , respectively. and denote the first and last entries of , respectively. Obviously, and has the rotational invariance property, that is, with . Thus, we have the following expression

where and

with and denoting the first and last entries of , respectively. Once the matrix is calculated, the DODs of the sources can be estimated. Thus, we construct the following optimization problem

where is used to ensure that the first sensors of the transmit array are well calibrated.

From (19), we can achieve the least square solution of as

By substituting (22) into (21), we have

where .

Utilizing the fact that , , and for matrices , , and vectors with proper dimensions, the objective function in (23) is further simplified as

Thus, the optimization in (23) is expressed as

This problem can be solved through the Lagrange multiplier method and the optimal solution is expressed as:

where . Thus, can be estimated through (22). Note that , we can use eigenvalue decomposition to achieve the estimation of and the DOD estimation of the source can be calculated as

where is the eigenvalue of .

In addition, the element of also can be estimated as

where , is the elements of .

We select a transformation matrix which makes , then the DOAs of the sources can be estimated by using a similar manner with the estimation of DODs. Thus, the detailed procedure is omitted in this paper. Furthermore, the DODs and DOAs can be paired using the method proposed in [9].

4. Simulation Results

In this section, several experimental simulations are carried out to demonstrate the performance of the proposed method.

A bistatic MIMO radar with transmit antennas and receive antennas is used in the following simulations. In addition, both the transmit and receiver arrays are assumed to be partially calibrated. The corresponding gain-phase uncertainties are and . Thus, we have and . The number of samples during one pulse period is , and the number of pulses is set as unless otherwise specified.

In the first simulation, 100 realizations of the proposed method are provided in Figure 1. In this simulation, the SNR is set as 5 dB. From Figure 1, it is clear that the proposed method achieves accurate DOD and DOA estimation which are paired correctly. The estimation of gain-phase uncertainties are provided in Table 1. The mean and standard deviation (STD) of the estimations are calculated from independent Monte Carlo trials.

Figure 1: 100 realizations of the proposed method. SNR = 5 dB.

Table 1: Performance of gain and phase estimation.

(a) Gain estimation results of

Index
True value
Mean
STD

0.8000
0.7999
0.0020

1.2500
1.2502
0.0045

0.7500
0.7503
0.0037

(b) Phase estimation results of (radian)

Index
True value
Mean
STD

0.6283
0.6285
0.0030

−1.0472
−1.0468
0.0042

0.7854
−0.7859
0.0057

(c) Gain estimation results of

Index
True value
Mean
STD

1.3712
1.3722
0.0053

1.1199
1.1212
0.0068

1.0000
1.0016
0.0081

0.9730
0.9744
0.0086

1.0000
1.0018
0.0119

0.9600
0.9622
0.0132

(d) Phase estimation results of (radian)

Index
True value
Mean
STD

−0.2675
−0.2677
0.0037

0.0676
0.0682
0.0064

0
0.0000
0.0084

−0.2509
−0.2500
0.0094

0
−0.0011
0.0123

0.0381
0.0392
0.0147

To evaluate the performance of the proposed method, we utilize the root mean square error (RMSE) result in the following simulations. In this paper, the RMSE is calculated from the following formula

where is the number of Monte Carlo trials and and are the estimations of and of the trial, respectively. In the following experiments, we use independent trails to calculate the RMSE of the angle estimation results.

In Figure 2, the RMSEs of the conventional ESPRIT algorithm, tensor-based ESPRIT algorithm, method in [16] (denoted as Guo’s method), and the proposed algorithm are compared with different SNRs. Due to the partially calibrated arrays, the conventional ESPRIT algorithm and the tensor-based ESPRIT algorithm fail to achieve accurate angle estimations. The Guo’s method [16] is robust under the partially calibrated array case, however, it has a worse RMSE performance. On the other hand, the proposed method achieves the best estimation among the existing algorithms especially for the case of low SNR.

Figure 2: RMSE versus SNR.

The probability of successful detection is provided in Figure 3. The successful detection means that the absolute error of all the DOA and DOD estimations is within . As mentioned above, both the conventional ESPRIT and the tensor-based ESPRIT algorithms exhibit a failed detection due to the partially calibrated arrays used in this simulation. However, the proposed method achieves a higher probability of successful detection especially for SNR greater than  dB. This simulation result verifies the improved performance of the proposed algorithm over existing algorithms.

Figure 3: Probability of successful detection.

In Figure 4, we evaluate the performance of the algorithms under consideration for the case of varying number of the pulses. It is observed that the existing algorithms experience a heavily degraded estimation performance when the number of pulses is smaller than . In addition, both the proposed algorithm and Guo’s method show improvement with the increased number of pulses. However, the proposed algorithm converges fast under low SNR and it always outperforms the other algorithms under consideration.

Figure 4: RMSE versus pulses.

5. Conclusion

In this paper, we proposed a novel joint DOD-DOA estimation algorithm based on tensor subspace approach for the case of partially calibrated bistatic MIMO radar. We first develop a third-order measurement tensor of the received signal. The HOSVD is then exploited to evaluate the signal subspace of the tensor-based received signal matrix. A close-form expression for estimating the gain-phase uncertainties of the transmit and receive arrays is derived, which enables us to estimate the joint DOA-DOA effectively. Compared with the existing algorithms, the proposed algorithm achieves a significant improvement in estimation performance. Finally, the simulation results demonstrate the effectiveness of the proposed algorithm.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work of Ping Huang is supported by the National Natural Science Foundation of China (61403091), the Heilongjiang Province Natural Science Foundation (QC2015049), and the launch of the postdoctoral research in Heilongjiang Province (LBH-Q16065).

References

R. Schmidt, “Multiple emitter location and signal parameter estimation,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276–280, 1986.

Z. Shi, C. Zhou, Y. Gu, N. A. Goodman, and F. Qu, “Source estimation using coprime array: a sparse reconstruction perspective,” IEEE Sensors Journal, vol. 17, no. 3, pp. 755–765, 2017.

Y. Gu, N. A. Goodman, S. Hong, and Y. Li, “Robust adaptive beamforming based on interference covariance matrix sparse reconstruction,” Signal Process, vol. 96, pp. 375–381, 2014.

Y. Gu and A. Leshem, “Robust adaptive beamforming based on interference covariance matrix reconstruction and steering vector estimation,” IEEE Transactions on Signal Processing, vol. 60, no. 7, pp. 3881–3885, 2012.

H. Krim and M. Viberg, “Two decades of array signal processing research: the parametric approach,” IEEE Signal Processing Magazine, vol. 13, no. 4, pp. 67–94, 1996.

E. Fishler, A. Haimovich, R. S. Blum, L. J. Cimini, D. Chizhik, and R. A. Valenzuela, “MIMO radar: an idea whose time has come,” in Proceedings of the 2004 IEEE Radar Conference (IEEE Cat. No.04CH37509), pp. 71–78, Philadelphia, PA, USA, 2004.

J. Li and P. Stoica, “MIMO radar with colocated antennas,” IEEE Signal Processing Magazine, vol. 24, no. 5, pp. 106–114, 2007.

C. Duofang, C. Baixiao, and Q. Guodong, “Angle estimation using ESPRIT in MIMO radar,” Electronics Letters, vol. 44, no. 12, p. 770, 2008.

C. Jinli, G. Hong, and S. Weimin, “Angle estimation using ESPRIT without pairing in MIMO radar,” Electronics Letters, vol. 44, no. 24, p. 1422, 2008.

X. Zhang, L. Xu, L. Xu, and D. Xu, “Direction of departure (DOD) and direction of arrival (DOA) estimation in MIMO radar with reduced-dimension MUSIC,” IEEE Communications Letters, vol. 14, no. 12, pp. 1161–1163, 2010.

M. L. Bencheikh, Y. Wang, and H. He, “Polynomial root finding technique for joint DOA DOD estimation in bistatic MIMO radar,” Signal Processing, vol. 90, no. 9, pp. 2723–2730, 2010.

R. Xie, Z. Liu, and J. Wu, “Direction finding with automatic pairing for bistatic MIMO radar,” Signal Processing, vol. 92, no. 1, pp. 198–203, 2012.

B. Friedlander and A. Weiss, “Eigenstructure methods for direction finding with sensor gain and phase uncertainties,” in ICASSP-88., International Conference on Acoustics, Speech, and Signal Processing, pp. 2681–2684, New York, NY, USA, April 1988.

A. Liu, G. Liao, C. Zeng, Z. Yang, and Q. Xu, “An eigenstructure method for estimating DOA and sensor gain-phase errors,” IEEE Transactions on Signal Processing, vol. 59, no. 12, pp. 5944–5956, 2011.

B. Liao and S. C. Chan, “Direction finding with partly calibrated uniform linear arrays,” IEEE Transactions on Antennas and Propagation, vol. 60, no. 2, pp. 922–929, 2012.

Y. D. Guo, Y. S. Zhang, and N. N. Tong, “ESPRIT-like angle estimation for bistatic MIMO radar with gain and phase uncertainties,” Electronics Letters, vol. 47, no. 17, p. 996, 2011.

B. Wang, W. Wang, Y. Gu, and S. Lei, “Underdetermined DOA estimation of quasi-stationary signals using a partly-calibrated array,” Sensors, vol. 17, no. 4, p. 702, 2017.

B. Wang, Y. D. Zhang, and W. Wang, “Robust DOA estimation in the presence of miscalibrated sensors,” IEEE Signal Processing Letters, vol. 24, no. 7, pp. 1073–1077, 2017.

B. Wang, S. Qin, Y. D. Zhang, and M. G. Amin, “Robust DOA estimation in the presence of mis-calibrated sensors,” in 2017 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), pp. 3076–3080, New Orleans, LA, USA, March 2017.

D. Nion and N. D. Sidiropoulos, “Tensor algebra and multidimensional harmonic retrieval in signal processing for MIMO radar,” IEEE Transactions on Signal Processing, vol. 58, no. 11, pp. 5693–5705, 2010.

T. G. Kolda and B. W. Bader, “Tensor decompositions and applications,” SIAM Review, vol. 51, no. 3, pp. 455–500, 2009.

M. Skolnik, Introduction to Radar Systems, McGraw-Hill, New York, NY, USA, 3rd edition, 2002.