Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2013 / Article

Research Article | Open Access

Volume 2013 |Article ID 427980 |

Fangqing Wen, Gong Zhang, "Two-Dimensional Direction-of-Arrival Estimation for Trilinear Decomposition-Based Monostatic Cross MIMO Radar", Mathematical Problems in Engineering, vol. 2013, Article ID 427980, 8 pages, 2013.

Two-Dimensional Direction-of-Arrival Estimation for Trilinear Decomposition-Based Monostatic Cross MIMO Radar

Academic Editor: Oleg V. Gendelman
Received22 Jul 2013
Revised05 Nov 2013
Accepted05 Nov 2013
Published27 Nov 2013


A low complexity monostatic cross multiple-in multiple-out (MIMO) radar scheme is proposed in this paper. The minimum-redundancy linear array (MRLA) is introduced in the cross radar to improve the efficiency of the array elements. The two-dimensional direction-of-arrival (DOA) estimation problem links to the trilinear model, which automatically pairs the estimated two-dimensional angles, requiring neither eigenvalue decomposition of received signal covariance matrix nor spectral peak searching. The proposed scheme performs better than the uniform linear arrays (ULA) configuration under the same conditions, and the proposed algorithm has less computational complexity than that of multiple signal classification (MUSIC) algorithm. Simulation results show the effectiveness of our scheme.

1. Introduction

Multiple-in multiple-out (MIMO) radar is a relative new concept in radar system that received considerable attention [1, 2]. Unlike the traditional phase-array radar, the transmit antennas of MIMO radar simultaneously transmit orthogonal waveforms and utilise multiple antennas to receive the reflected signals [3]. Angel estimation is a basic problem in MIMO radar that has been investigated in [4, 5]. Most of the previous studies are based on the uniform arrays for the simple array geometry, such as uniform linear arrays (ULA) and uniform circular-shape and L-shape arrays. The minimum-redundancy array (MRA) is an important type of the nonuniform linear geometries, with the sensor placement designed to achieve the maximum resolution for a given number of elements by reducing the redundancy of spatial correlation lags [6]. Studies show that the MRA has better angle resolution and interference suppression than the ULA for a given number of sensors [7]. Due to the capability of spanning larger aperture, a MRA configured colocated MIMO radar is proposed for one dimensional angle estimation in [8], which has achieved desired results.

Classical angle estimation algorithms such as multiple signal classification (MUSIC) and estimation of signal parameters via rotational invariance techniques (ESPRIT) require eigenvalue decomposition (EVD) or singular value decomposition (SVD) of received signal covariance matrix. In addition, MUSIC algorithm requires spectral peak searching, which would cause computation disaster in the situation of 2D angle estimation. ESPRIT algorithm is sensitive to the manifold of the arrays, making it unsuitable for MRA manifolds. Another problem that how to pair the estimated parameters arises within the ESPRIT algorithm, which requires the extra computational load, and usually fails to work in lower SNR. Furthermore, the subspace based methods suppose that the number of the sources is known, which is contrary to actual applications.

Trilinear decomposition (also called PARAFAC analysis) can be thought of as a generalization of ESPRIT and joint approximate diagonalization ideas [9], which is an iterative algorithm that does not need EVD or SVD. Trilinear decomposition algorithm has been used for detection and localisation of multiple targets in a ULA based MIMO radar system [10, 11]. Simulation in [11] shows that the trilinear decomposition algorithm can achieve better performance than that of ESPRIT with bistatic MIMO radar. However, the trilinear algorithm would be lapsed with the MRA configuration, as the phase ambiguity follows the received array signal, which would confuse the paired angles.

In this paper, the minimum-redundancy linear array (MRLA) configured monostatic cross MIMO radar is proposed for two-dimensional (2D) direction-of-arrival (DOA) estimation, which turned out to be effective than the ULA configuration for a given number of antennas. The improved trilinear model based blind angle estimation algorithm is developed for the proposed radar scheme. The DOA estimation problem links to the trilinear model. In order to eliminate the phase ambiguity, the estimated phase is adaptively compensated according to the array manifold. With this improvement, the proposed algorithm could be extended for arbitrarily array manifolds. We also analyse the complexity and error characteristics of the algorithm. Simulation results show the effectiveness of the proposed scheme.

The rest of the paper is organized as follows. In Section 2, we present the signal model for the monostatic cross MIMO radar. Section 3 specifies the DOA estimation method with trilinear decomposition. Simulation results are given in Section 4. Finally, we provide concluding remarks in Section 5.

Notation. Lowercase and capital letters in bold denote, respectively, vectors and matrices. The superscripts and represent, respectively, the transpose and Hermitian transpose of a matrix. and denote arrangement of the elements in a row and in a column, respectively.

2. Signal Model

2.1. Monostatic Cross MIMO Radar

Consider a monostatic cross MIMO radar model, as shown in Figure 1. Without loss of generality, we suppose that -element transmitter and -element receiver are symmetrically located in the -axis and -axis. The phase center of the transmits and receivers is overlapped in the origin, while the physical location of the th transmitter and the th receiver is denoted by and , respectively. It is also assumed that there are noncoherent targets appearing in the far-field of the antennas with the th target at azimuth angle and elevation angle . The transmit antennas transmit orthogonal waveforms with the same carrier frequency, and the transmit waveform of th antenna during th pulse is , where is the total pulse number, is the samples in a pulse period, and the reflection coefficient of the th target is assumed to be constant during the th pulse period. The received signal of the th receiver takes the form where is the transmit steering vector for the th target, = is the corresponding receive steering vector,  and = , with denotes the wave length of the carrier. = is the transmit wave matrix, and is the received additive white Gaussian noise of th receive antenna corresponding to the th pulse.

The received signal is correlated with a matched filter for maximizing the signal to noise ratio (SNR) in the presence of additive stochastic noise, and the output of the matched filters corresponding to the th pulse can be expressed as where accounts for the th column of .

Define as the reflection coefficients matrix, which contains the amplitude and phase information of all the targets. We suppose that the transmit waveforms fulfill Let ; then can be expressed in matrix form as where the transmit direction matrix , the receive direction matrix , and represents the interference and jamming. denotes the operation extracting the th row of its matrix argument and constructing a diagonal matrix out of it.

2.2. Minimum-Redundancy Linear Arrays

Due to the matched filters in every receiver antenna, the transmitted waveforms can be separated from each other. The virtual elements were formed through the channel of the th transmit and the th receiver, which increase the array aperture and thus improve the resolution of the MIMO radar. A uniform, linear, critically sampled array with elements is used as a reference; as depicted in Figure 2, the solid points and the hollow points represent the real elements and virtual elements, respectively. Some virtual elements overlap with real elements that are not given in the figure. Every virtual element equals a self-transmit and self-receive element, from which we can calculate that the freedom of the real aperture cross array is .

However, the referenced uniform linear array is redundant, which cannot utilize the elements of the array in the most efficient way. In the conventional methods, the minimum-redundancy linear array geometries together with appropriate matrix augmentation techniques have proven to provide better performances than the simple uniform linear array configuration in terms of ability to detect and resolve a greater number of sources [12]. Therefore, the transmit array and receive array of the cross MIMO radar are configured to the form of MRLA. Table 1 gives some MRLA configurations for a finite number of antennas.

Antenna numberThe relative locations of antennas


2.3. Trilinear Model

The data in (4) can be rearranged as where is the noise, stands for the th element of , and similarly to the others. The signal in (5) can be denoted as the trilinear slice model, which displays the reflection for three different kinds of diversity, as shown in Figure 3. The matrix can be interpreted as slicing the three dimension data in a series of slices along the spatial direction. Due to the symmetry of the trilinear model in (5), two more matrices can be interpreted as slicing the three-way data along different directions, which are denoted as

The data matrix , , and are slices along the different directions when slicing the trilinear model. Rearrange the slicing data as , , and ; the problem of cross MIMO radar 2D angle estimation is linked to trilinear decomposition, which is also known as parallel factor analysis (PARAFAC) or canonical decomposition. It will be discussed in detail in the following section.

3. The Trilinear Decomposition-Based 2D DOA Estimation Method

3.1. Trilinear Alternate Least Squares

Trilinear alternating least square (TALS) algorithm is the common data detection method for trilinear model [13]. The principle of TALS can be adopted to fit low rank trilinear models. The basic idea of TALS can be described as three major steps: (a) update one of the slicing matrices , , or each time using least squares (LS), which is conditioned on previously obtained estimates for the remaining matrices. (b) Proceed to update the other matrices using LS. (c) Repeat (a) and (b) until convergence of the LS cost function. TALS algorithm in our scheme is discussed as follows.

According to (4), the LS fitting of is where stands for the Frobenius norm. The LS update for is where stands for the pseudoinverse and and are the matrices estimated in the last iteration.

Similarly, the LS fitting of is and the LS update for is where and stand for the previously obtained estimation of and .

Finally, the LS fitting of is and the LS update for is where and stand for the previously obtained estimates of and .

Although TALS algorithm is easy to implement and guarantee to converge, it also suffers from slow convergence. In [14] the COMFAC algorithm is proposed to speed up the LS fitting of TALS. COMFAC compresses the high dimensional three-way array into the subspaces defined by the three bases of , , and . After fitting the PARAFAC model in the compressed space, the solution is recovered to the original space, which is followed by a few TALS steps in the uncompressed space. For noise-free data with no model error an exact solution can be calculated directly by posing the fitting problem as a generalized eigenvalue problem. With the presence of the noise, the COMFAC algorithm gave the LS estimation of , , and . In this paper, the COMFAC algorithm is selected for trilinear decomposition.

3.2. Uniqueness and Identifiability Analysis

Just as singular value decomposition (SVD) is a decomposition of two-way arrays into rank one two-way components, trilinear decomposition is a decomposition of three-way arrays into rank one three-way components. The fundamental difference is that the trilinear decomposition is unique under certain rank conditions, whereas bilinear decompositions into rank one components are not unique without imposing additional constraints. According to [15], the trilinear decomposition of the noiseless three-way array is unique if where , , and are -rank of the corresponding matrices and is the rank of . A matrix has -rank equal to if any subset of columns are linearly independent and there exists a subset of linearly dependent columns. When the uniqueness condition (13) is satisfied, trilinear decomposition of the received three-way array can be used to uniquely recover matrices , , and up to permutation and scaling of their columns. With the noisy received signal, any other triple , , and that construct is related to , , and via where is a permutation matrix, , , and are the corresponding estimation errors, which is closely related to the SNR. , , and are diagonal scaling matrix with .

When the matrix and are Vandermonde matrices with distinct nonzero generators, and is full -rank, (13) becomes and in general we have ; the maximal number of targets that the proposed algorithm can detect is .

3.3. 2D DOA Estimation

Once trilinear decomposition is accomplished, the estimated direction matrices and are obtained, and then the 2D DOA can be estimated with LS method. The estimated transmit vector was firstly normalized to the reference phase center. Define where with . As the phase ambiguity exits the MRLA configured receive signal, the added makes sure that is an ascending sequence. Let and define ; it is obvious that is the LS estimation of . The LS solution for is

Similarity, define where is the LS estimation of , . The LS solution for is

The azimuth angle and elevation angle can be paired automatically through the following formula:

3.4. Scheme Analysis

Traditional Capon and MUSIC algorithm requires 2D peak searching, which brings very heavy complexity. Others are effective only in condition that transmit and receive are uniform arrays. However, the proposed trilinear decomposition-based algorithm is effective for manifold of the ULA and MRLA. Compared with ESPRIT, which process the complexity of [11], the complexity of each TALS iteration is ; generally is relatively small compared to , , and , and only a few iterations are required to achieve convergence.

According to [16], the Cramér-Rao Bound (CRB) of the estimated error can be derived as where is defined as and and represent Hadamard product and Kronecker product, respectively. is the noise level after the matched filtering.   =  , = , and with and .

4. Simulation Results

To assess the angle estimation performance of the proposed algorithm, 1000 Monte Carlo simulations are presented. The root mean squared error (RMSE) of 2D-DOD is defined as follows: where and denote the estimation of azimuth and elevation angle of the th target in the th Monte Carlo trial. Note that , , , and are the number of transmit antennas, the receive antennas, the pulse, and targets, respectively. The reflection coefficients of the targets are assumed to fulfill the classic Swerling Case 2. The SNR is defined before the matched filtering. We assume that there are three uncorrelated targets located at angles , , and , respectively. and the manifold of the transmit and receive antennas are referenced in Table 1. The samples during a pulse are designated to .

Figures 4 and 5 present the estimation results of the proposed scheme with and separately and with considered in the simulation. It is observed that the proposed scheme is able to estimate azimuth angle and elevation angle for the radar system.

We compare the proposed algorithm and 2D MUSIC algorithm with , and the searching intervals of 2D MUSIC algorithm are all . Figure 6 presents the RMSE of the estimation results. We find that the RMSE performance of the proposed algorithm is very close to that of the 2D MUSIC algorithm (due to the limited searching grid, the performance of MUSIC algorithm seems better than the proposed algorithm with higher SNR). Table 2 shows the average time the algorithms run. It is obvious that the proposed algorithm processes a higher operating efficiency. The PARAFAC algorithm is iterative and good starting values are beneficial in keeping computation time low [14], but it is hard for the algorithm to achieve a perfect initial value with lower SNR. The simulation results also prove that the low SNR would affect the convergence speed of TALS.

Algorithm SNR−5 dB0 dB5 dB10 dB

2D MUSIC6.476 s6.479 s6.529 s6.482 s
Proposed algorithm0.236 s0.228 s0.227 s0.226 s

Algorithm SNR15 dB20 dB25 dB30 dB

2D MUSIC6.498 s6.559 s6.578 s6.508 s
Proposed algorithm0.224 s0.224 s0.226 s0.224 s

Figure 7 depicts the angle estimation performance comparison of the proposed MRLA based cross MIMO radar and the ULA configuration with and , respectively. Simulation results agree with what we have expected, and the performance of MRLA configuration is very close to that of the equivalent ULA configuration.

The performance of the proposed scheme under different pulse number is investigated in the simulation, as shown in Figure 8. It indicates that angle estimation performance of this algorithm is improved with snapshot number increasing, but after the snapshot reaches a level, this improvement slows down. This phenomenon can be interpreted from (22). That means that the estimated error is affected by the noise variance and more snapshots mean more accurate estimation of the noise. The estimation error would tend to stable () with sufficient snapshot number.

We also compared the RMSE performance of the cross radar with MRLA and ULA configuration, as shown in Figure 9 and is considered in the simulation. It is indicated that the angle estimation performance of our algorithm is gradually improved with the increasing number of antennas. Multiple antennas improve estimation performance because of diversity gain. Compared with the ULA configuration with the same number of antenna, the MRLA manifold brings much more virtual elements and thus achieves much better estimation performance.

5. Conclusion

In this paper, we proposed the MRLA scheme for monostatic cross MIMO radar, and the trilinear decomposition algorithm is developed for two-dimensional angle estimation. Without eigenvalue decomposition of received signal covariance matrix and peak searching, the proposed algorithm can achieve automatic pairing two-dimensional angles and works well. Furthermore, our scheme has better performance than the ULA based scheme and has a very close performance to 2D MUSIC algorithm, which has heavier computational complexity than the proposed algorithm. Simulation results proved that the proposed scheme would achieve super performance in radar parameter estimation.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This work was supported by the Chinese National Natural Science Foundation under Contract nos. 61071163, 61071164, 61271327, and 61201367 and partly funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PADA).


  1. L. Xu, J. Li, and P. Stoica, “Radar imaging via adaptive mimo techniques,” in Proceedings of the 14th European Signal Processing Conference (EUSIPCO '06), EURASIP, Florence, Italy, September 2006. View at: Google Scholar
  2. J. Li, P. Stoica, L. Xu, and W. Roberts, “On parameter identifiability of MIMO radar,” IEEE Signal Processing Letters, vol. 14, no. 12, pp. 968–971, 2007. View at: Publisher Site | Google Scholar
  3. E. Fishler, A. Haimovich, R. Blum, D. Chizhik, L. Cimini, and R. Valenzuela, “MIMO radar: an idea whose time has come,” in Proceedings of the IEEE Radar Conference, pp. 71–78, IEEE, April 2004. View at: Google Scholar
  4. 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  5. Y. Cheng, R. Yu, H. Gu, and W. Su, “Multi-svd based subspace estimation to improve angle estimation accuracy in bistatic mimo radar,” Signal Processing, 2013. View at: Google Scholar
  6. H. Duan, B. P. Ng, C. M. S. See, and J. Fang, “Broadband interference suppression performance of minimum redundancy arrays,” IEEE Transactions on Signal Processing, vol. 56, no. 6, pp. 2406–2416, 2008. View at: Publisher Site | Google Scholar | MathSciNet
  7. H. Duan, B. P. Ng, C. M. S. See, and J. Fang, “Spatial resolutions of the broadband nonredundant and minimum redundancy arrays,” IEEE Signal Processing Letters, vol. 14, no. 11, pp. 852–855, 2007. View at: Publisher Site | Google Scholar
  8. Z. Peng, D. Ben, G. Zhang, and H. Gu, “Non-uniform linear array configuration for MIMO radar,” Transactions of Nanjing University of Aeronautics and Astronautics, vol. 28, no. 2, pp. 152–156, 2011. View at: Google Scholar
  9. L. De Lathauwer, “A link between the canonical decomposition in multilinear algebra and simultaneous matrix diagonalization,” SIAM Journal on Matrix Analysis and Applications, vol. 28, no. 3, pp. 642–666, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  10. D. Nion and N. D. Sidiropoulos, “A Parafac-based technique for detection and localization of multiple targets in a MIMO radar system,” in Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '09), pp. 2077–2080, IEEE, April 2009. View at: Publisher Site | Google Scholar
  11. X. Zhang, Z. Xu, L. Xu, and D. Xu, “Trilinear decomposition-based transmit angle and receive angle estimation for multiple-input multiple-output radar,” IET Radar, Sonar and Navigation, vol. 5, no. 6, pp. 626–631, 2011. View at: Publisher Site | Google Scholar
  12. G. Gelli and L. Izzo, “Minimum-redundancy linear arrays for cyclostationarity-based source location,” IEEE Transactions on Signal Processing, vol. 45, no. 10, pp. 2605–2608, 1997. View at: Publisher Site | Google Scholar
  13. N. D. Sidiropoulos, G. B. Giannakis, and R. Bro, “Blind parafac receivers for DS-CDMA systems,” IEEE Transactions on Signal Processing, vol. 48, no. 3, pp. 810–823, 2000. View at: Publisher Site | Google Scholar
  14. R. Bro, N. D. Sidiropoulos, and G. B. Giannakis, “A fast least squares algorithm for separating trilinear mixtures,” in Workshop on Independent Component Analysis for Blind Signal Separation, pp. 289–294, 1999. View at: Google Scholar
  15. N. D. Sidiropoulos and X. Liu, “Identifiability results for blind beamforming in incoherent multipath with small delay spread,” IEEE Transactions on Signal Processing, vol. 49, no. 1, pp. 228–236, 2001. View at: Publisher Site | Google Scholar
  16. P. Stoica and A. Nehorai, “Performance study of conditional and unconditional direction-of-arrival estimation,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. 38, no. 10, pp. 1783–1795, 1990. View at: Publisher Site | Google Scholar | Zentralblatt MATH

Copyright © 2013 Fangqing Wen and Gong Zhang. 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.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.