International Journal of Antennas and Propagation

Volume 2018, Article ID 4370637, 7 pages

https://doi.org/10.1155/2018/4370637

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

^{1}College of Automation, Harbin Engineering University, Harbin, Heilongjiang 150001, China^{2}College of Electrical Engineering and Automation, Harbin Institute of Technology University, Harbin 150001, China

Correspondence should be addressed to Junxiang Wang; nc.ude.uebrh@gnaixnujgnaw

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.