Mathematical Problems in Engineering
Volume 2020 (2020), Article ID 5720189, 8 pages
https://doi.org/10.1155/2020/5720189

Research Article

Tensor-Based Angle and Range Estimation Method in Monostatic FDA-MIMO Radar

Tengxian Xu,¹ Yongqin Yang,¹ Mengxing Huang,¹ Han Wang,² Di Wu,¹ and Qu Yi³

¹State Key Laboratory of Marine Resource Utilization in South China Sea, School of Information and Communication Engineering, Haikou 570228, China
²College of Physical Science and Technology, Yichun University, Yichun 336000, China
³College of Physics and Electronic Engineering, Hainan Normal University, Haikou 571158, China

Correspondence should be addressed to Yongqin Yang; yangyq@hainanu.edu.cn and Mengxing Huang; huangmx09@163.com

Received 3 July 2020; Accepted 15 July 2020; Published 12 August 2020

Guest Editor: Wang Zheng

Copyright © 2020 Tengxian Xu 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.

In the paper, joint angle and range estimation issue for monostatic frequency diverse array multiple-input multiple-output (FDA-MIMO) is proposed, and a tensor-based framework is addressed to solve it. The proposed method exploits the multidimensional structure of matched filters in FDA-MIMO radar. Firstly, stack the received data to form a third-order tensor so that the multidimensional structure information of the received data can be acquired. Then, the steering matrices contain the angle and rang information are estimated by using the parallel factor (PARAFAC) decomposition. Finally, the angle and range are achieved by utilizing the phase characteristic of the steering matrices. Due to exploiting the multidimensional structure of the received data to further suppress the effect of noise, the proposed method performs better in angle and range estimation than the existing algorithms based on ESPRIT, simulation results can prove the proposed method’s effectiveness.

1. Introduction

Multiple-input multiple-output (MIMO) radar was first proposed in [1–3], which is a key research point in today’s radar field. In MIMO radar, all the antennas are omnidirectional and the transmitted signals are orthogonal to each other, which can achieve the waveform diversity gain in the receiver side. Currently, there are two main types of MIMO radars, namely, collocated MIMO radar [4–6] and statistical MIMO radar [7, 8]. The statistical MIMO radar is composed with separated transmit and receive antennas for obtaining both the waveform and spatial diversity gain. In contrast, in order to improve the estimation performance, the collocated MIMO radar closely places the transmitting and receiving arrays to form a virtual array that has large aperture. In the past few years, parameter estimation, including direction of arrival (DOA), direction of departure (DOD), and Doppler frequency, became a hot topic and investigated by a lot of researchers. In the common subspace-based methods, such as multiple signal classification algorithm (MUSIC) and estimation of signal parameters via rotational invariance techniques (ESPRIT), MIMO radar has used them to achieve DOD and DOA estimation [9–14], which can obtain the desired performance with reasonable SNR and snapshots. In addition, a joint angle and Doppler frequency estimation method is investigated, which can be used to track the targets. However, it is noted that these methods for MIMO radar with narrow-band signals cannot achieve the rang information which is very important for target localization in practice.

FDA radar can get the distance information of the target [15]. There is one difference between the FDA radar and the traditional phased array radar, and for the FDA radar, the transmitting frequency of each transmitting antenna has a frequency increment so that it can draw the relevant beam pattern, improve the degree of freedom in space, and then achieve the joint angle and distance estimation about target in the FDA radar. However, the angle and range of the target in the beam domain of the FDA radar are mixed, so it cannot directly obtain the angle and distance information of the target through spectral peak search. At present, there are some solutions to this problem. In [16], a nonlinear FDA radar is proposed, but this method needs to change the array element position in actual operation, which is difficult to achieve. In [17, 18], a method based on subarrays to achieve angle and distance decoupling is proposed, but this method has the problem of distance ambiguity. In addition, in [19], the target positioning method of the coprime array FDA radar is studied, and in [20, 21], the nested FDA parameter estimation method is studied.

In [18], FDA-MIMO radar with both frequency and waveform diversity is proposed for the first time, which makes use of the unique advantages of FDA radar and MIMO radar, which not only has excellent target detection performance in target detection but also has a high freedom of the spatial degree; it can realize the joint estimation of the target angle and range [17, 22–24]. At present, there are many kinds of FDA-MIMO radar target positioning algorithms, such as MUSIC methods and ESPRIT methods, but the MUSIC method involves a spectral peak search and results in a high algorithm complexity. Although the ESPRIT method can be applied to reduce the computational complexity and obtain the angle and range information of the target through rotation invariance, this method is not suitable for the phase ambiguity situation.

The advantages of FDA are developed based on the characteristics of the multidimensional structure of the signal. The abovementioned traditional matrix processing method cannot make good use of the multidimensional characteristics, which limits their performance [25]. To take advantage of the multidimensional structural characteristics of the signal, the PARAFAC decomposition is used for angle estimation in unknown target localization by modelling tensor signals [26–28]. At present, there are few studies on tensor-based FDA MIMO target estimation methods. The paper proposes the algorithm based on parallel factorization in monostatic FDA-MIMO radar. This method first establishes a third-order tensor signal model and uses parallel factorization to obtain the transmission direction matrix and the reception direction matrix [29, 30]. Then, obtain the angle and distance information of the target by taking the phase of the direction matrix and finally realize the angle and distance estimation.

The proposed method compares the achievement with traditional estimation of signal parameters via rotational invariance techniques (ESPRIT) method [31], Unitary ESPRIT method [32], and Cramér–Rao Bound (CRB). The paper structure is divided into multiple parts. Section 2 introduces the tensor-based data model. Section 3 deduces the method of angle estimation and distance estimation. The CRB of the target estimation method in FDA-MIMO radar is shown in Section 4. The Section 5 shows the simulation results.

Notation: , , and represent the Kronecker product, Khatri–Rao product, and Hadamard product, respectively. , , , , and stand for transpose, conjugate-transpose, inverse, conjugate, and pseudoinverse, respectively. represents a diagonal matrix, which composed by the th row of . represents the Frobenius norm. indicates a identity matrix. angle and represent the phase angle of each element and the minimum element in the array. denotes the real part operator. is denoted as a third-order tensor.

2. Tensor-Based Data Model of Monostatic FDA-MIMO Radar

The paper is based on a narrowband monostatic FDA-MIMO radar. It is composed by -element transmitting antennas and -element receiving antennas. Figure 1 shows the structure of monostatic FDA-MIMO radar. The transmitting array and receiving array are placed together so the direction of arrival (DOA) is the same as the direction of departure (DOD). In general, the transmitting array and receiving array are uniform linear arrays (ULA) and are separated by half a wavelength. In the paper, there are far-field targets by default which are independent of each other. Based on the first element of the transmitting antenna, the frequency of the th element transmitting antenna is derived as

where the reference frequency is denoted by and denotes frequency increase of the transmitting antenna. Moreover, .

Figure 1: Monostatic FDA-MIMO radar.

After the matched filter, the receiving data can be rearranged as

where denotes an noise vector that is the additional white Gaussian noise vector, and stand for the transmit steering vector and the receive steering vector, respectively, and stands for the received data of samples, which can be expressed as

where , the transmitting direction matrix is defined as , and the receiving direction matrix is defined as . stands for Khatri–Rao product. So, is expressed as

Under the interference of noise, the signal model is updated as

where stands for the received noise of the th slice. By the characteristics of the tensor model, equation (4) is expanded by

where is defined as the element of the matrix and is defined as the element of the matrix . Another term for the PARAFAC decomposition of a third-order tensor is trilinear decomposition. Figure 2 shows the specific decomposition process.

Figure 2: Schematic diagram of trilinear decomposition.

Under the definition of PARAFAC decomposition, equation (6) can be regarded as a trilinear model with symmetrical properties. So, we can get and . By the symmetry of the trilinear decomposition, the matrices and are constructed as

We apply the tensor signal model and parallel factor decomposition to the monostatic FDA-MIMO radar and derive the PARAFAC-based angle and range estimation method. In the following section, we will derive the method in more detail.

3. The Proposed Method

3.1. Parallel Factor Decomposition

In this section, according to the trilinear alternating least square (TALS) method, we can estimate the transmitting direction matrix and the receiving direction matrix. The TALS method can be explored for data analysis in trilinear models.

According to [33], the estimation of the direction matrices can be obtained by the parallel factor decomposition. The least squares (LS) update of equation (3) is written by

where the noise signal is denoted by . The LS update of matrix is written by

where and stand for the estimates of and which have previously obtained, respectively. The LS fitting of equation (7) can be written as

where the noise signal is denoted by . The LS update of matrix is written by

where and on behalf of the estimates of and which have previously obtained, respectively. Similarly, the LS fitting of equation (8) can be written as

where the noise signal is denoted by . The LS update of matrix is

where and on behalf of the estimates of and which have previously obtained, respectively.

According to equations (10), (12), and (14), , , and are given by the LS. The LS update does not stop until convergence, and the constraint condition can be expressed as .

For the received noise signals, according to the trilinear decomposition, we can get the estimated parameter matrices: , , and , where is a permutation matrix, , , and stand for the matrices in correspondence to estimation errors, and , , and denote the diagonal scaling matrices and are subject to . The inherent scale ambiguity and permutation of trilinear decomposition and normalization can be used to eliminate scale ambiguity.

3.2. Angle Estimation

From the trilinear decomposition in the previous section, we can get estimates of the direction matrices and by the TALS method.

It can be seen from the signal model that the steer vector of is only related to the direction of arrival (DOA), and we can get the estimated value of the angle through .

In the paper, the receive array can be preset as half-wavelengths spaced uniform linear arrays (ULA), so the receive steering vector of is

Let denote the phase of , which can be expressed as

After the receiving steering vector is well calibrated, it is represented as and utilizes normalization to dispel scale ambiguity. So, can be obtained from equation (16). The estimate of can be calculated by LS principle. We can construct the LS fitting as

where the and stand for the estimated vector and the estimated phase of steering vector, is defined as

is the LS solution of by equation (17), which stands for the estimated vector and can be expressed as

The receive angle is derived from

where represents the element in the second row of the vector.

3.3. Range Estimation

For the tensor-based data model, the range and transmit angle of the uncorrelated target are both included in the transmit steering vector. So, the transmit steering vector can be expressed as

where and denote the transmitting angle steering vector and range steering vector, respectively.

The transmit array is preset as half-wavelengths spaced ULA, so the transmitting steer vector of can be reduced as

In the previous section, the direction of arrival (DOA) has been obtained. Referring to the derivation process in the previous section and letting represent the phase of , which can be expressed as

In the previous section, we have obtained the estimated receive vector . Similarly, we can get by equation (17), which stands for the estimated transmit receive vector and can be expressed as

The range of target is derived from

where represents the element in the second row of the vector .

So far, we already get the angle and range of the uncorrelated target in the monostatic FDA-MIMO radar. We can recap the key processes of the proposed method as follows:

(i)Step 1: according to equation (6), the received data is transformed into a third-order tensor model.(ii)Step 2: apply the trilinear decomposition principle; we can obtain the transmitting direction matrix and the receiving direction matrix .(iii)Step 3: take the phase and of and by equations (16) and (23), respectively.(iv)Step 4: the phase of the receive direction matrix was plugged into equation (17) and obtain by LS fitting; then, the receive angle can be calculated by equation (20).(v)Step 5: similarly, we can obtain; then, the range of target can be calculated by substituting the angle into equation (25).

4. CRB of FDA-MIMO Radar

The input signal spectrum is defined as

where is the unknown signal spectrum matrix. CRB formula can be expressed as [17, 32]

where , and .

When is large enough, and are not close to singularity and the signals are unknown; the above formula can be approximated as an approximate CRB:

5. Simulation Results

In this section, the effectiveness and advantages of our proposed method can be proved by some numerical simulations. The ESPRIT method and Unitary-ESPRIT method can be utilized to contrast with our proposed method. We default the monostatic FDA-MIMO radar with transmitting antennas and receiving antennas in this paper.

Unless otherwise specified, it is supposed that there are uncorrelated targets which are objected to far-field. The three targets are placed at , , and . The performance of our proposed method can be appraised by the root mean square error (RMSE), which can be expressed as

where and are the estimation of DOA and range of the th target for the th Monte Carlo trials, respectively, denotes the total amount of Monte Carlo trials, and is preset in this simulation.

In addition, the probability of successful detection is another metric used to appraise the achievement of our method, which is defined by

where represents the number of successful estimates, and the criterion for the success of the angle and range experiments is that the absolute value of all experimental numerical errors are less than the minimum and , respectively.

We preset in the first simulation; Figure 3 shows the estimation results of our proposed method. Figure 3 shows the estimation of the range and angle of the uncorrelated targets are correct, and the landing points are highly concentrated. It directly proves the stability and accuracy of the proposed method.

Figure 3: The range and angle estimation performance of the proposed method with .

And then, denotes the number of snapshots, and it is preset to in this simulation. We first investigate the liaison between RMSE and SNR of range estimation and angle estimation in the second simulation. We use two comparison methods, which are the ESPRIT method [31] and the Unitary ESPRIT method [32], respectively. After introducing CRB, they are compared with the proposed method. Figures 4 and 5 correspond to the simulation. From the two figures, it can be demonstrated that the estimated performance of the proposed method is better than the ESPRIT method and the Unitary ESPRIT method. In addition, the RMSE of the proposed method is closer to CRB.

Figure 4: RMSE of angle estimation versus SNR.

Figure 5: RMSE of range estimation versus SNR.

From the following, we explore the relationship between RMSE and snapshots of range estimation and angle estimation in the third simulation, and the results are shown in Figures 6 and 7, respectively. The is preset in this simulation. Similar to the first simulation, after introducing CRB, we use the contrast methods to compare with the proposed algorithm. As the number of snapshots increase, we can know from the figures that the RMSE of three methods and CRB all decrease. So, we conclude that the accuracy of the methods is improving, and the proposed method has the best performance.

Figure 6: RMSE of angle estimation versus the total amount of snapshots.

Figure 7: RMSE of range estimation versus the total amount of snapshots.

In the fourth simulation, the relationship between probability and SNR of angle and range probability successful detection is obtained. The number of snapshots is preset to . Figures 8 and 9 correspond to the fourth simulation. From two figures, we can know the success rate of three methods can reach 100% with the increase of SNR. Moreover, with the increase of SNR, the probability of successful detection also increases. It can reach a detection success rate of , when the SNR reaches a sufficiently high level, this SNR is generally called the SNR threshold. Obviously, the threshold of the proposed method is the lowest among the three methods, which also illustrates the superiority of the proposed method.

Figure 8: Probability of angle successful detection versus SNR.

Figure 9: Probability of range successful detection versus SNR.

6. Conclusions

In the paper, we proposed a tensor-based range and angle estimation method in monostatic FDA-MIMO radar. The proposed method uses the trilinear model to obtain the direction matrices through PARAFAC decomposition and extracts the phase from the direction matrix to estimate the distance and angle. This method uses the multidimensional information of the received data. Compared with the subspace methods such as ESPRIT and Unitary-ESPRIT methods, the proposed method has the best performance. The superiority of the proposed method can be verified by simulation.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Hainan Provincial Natural Science Foundation of China (no. 2018CXTD336) and National Natural Science Foundation of China (no. 61864002).

References

E. Fishler, A. Haimovich, R. Blum et al., “MIMO radar: an idea whose time has come,” in Proceedings of the IEEE Radar conference, vol. 7, pp. 71–78, Philadelphia, PA, USA, April 2004.

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. Blum et al., “Performance of MIMO radar systems: advantages of angular diversity,” in Proceeding of the Thirty-Eighth Asilomar Conference on Signals, Systems and Computer, vol. 1, pp. 305–309, Pacific Grove, CA, USA, November 2004.

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

A. A. Gorji, R. Tharmarasa, T. Kirubarajan et al., “Optimal antenna allocation in MIMO radars with collocated antennas,” IEEE Transactions on Aerospace and Electronic Systems, vol. 50, no. 1, pp. 542–558, 2014.

A. Hassanien and S. A. Vorobyov, “Phased-MIMO radar: a tradeoff between phased-array and MIMO radars,” IEEE Transactions Signal Processing, vol. 58, no. 6, pp. 3137–3151, 2009.

E. Fishler, A. Haimovich, R. S. Blum, L. J. Cimini, D. Chizhik, and R. A. Valenzuela, “Spatial diversity in radars-models and detection performance,” IEEE Transactions on Signal Processing, vol. 54, no. 3, pp. 823–838, 2006.

A. Haimovich, R. Blum, L. Cimini et al., “MIMO radar with widely separated antennas,” IEEE Signal Processing Magazine, vol. 25, no. 1, pp. 116–129, 2008.

A. J. Fenn, D. H. Temme, W. P. Delaney et al., “The development of phased-array radar technology,” Lincoln Laboratory Journal, vol. 12, no. 2, pp. 321–340, 2000.

Y. Hua, T. K. Sarkar, D. D. Weiner et al., “An L-shaped array for estimating 2-D directions of wave arrival,” IEEE Transactions on Antennas and Propagation, vol. 39, no. 2, pp. 143–146, 1991.

X. Wang, L. Wan, M. Huang, C. Shen, and K. Zhang, “Polarization channel estimation for circular and non-circular signals in massive MIMO systems,” IEEE Journal of Selected Topics in Signal Processing, vol. 13, no. 5, pp. 1001–1016, 2019.

H. Wang, L. Wan, M. Dong, K. Ota, and X. Wang, “Assistant vehicle localization based on three collaborative base stations via SBL-based robust DOA estimation,” IEEE Internet of Things Journal, vol. 6, no. 3, pp. 5766–5777, 2019.

X. Wang, D. Meng, M. Huang, and L. Wan, “Reweighted regularized sparse recovery for DOA estimation with unknown mutual coupling,” IEEE Communications Letters, vol. 23, no. 2, pp. 290–293, 2019.

D. Meng, X. Wang, M. Huang, L. Wan, and B. Zhang, “Robust weighted subspace fitting for DOA estimation via block sparse recovery,” IEEE Communications Letters, vol. 24, no. 3, pp. 563–567, 2020.

P. Antonik, M. C. Wicks, H. D. Griffiths et al., “Frequency diverse array radars,” in Proceedings of the 2006 IEEE Radar Conference, pp. 215–217, Verona, NY, USA, April 2006.

P. F. Sammartino, C. J. Baker, H. D. Griffiths et al., “Frequency diverse MIMO techniques for radar,” IEEE Transactions on Aerospace and Electronic Systems, vol. 49, no. 1, pp. 201–222, 2013.

W.-Q. Wang and H. C. So, “Transmit subaperturing for range and angle estimation in frequency diverse array radar,” IEEE Transactions on Signal Processing, vol. 62, no. 8, pp. 2000–2011, 2014.

W.-Q. Wang, “Subarray-based frequency diverse array radar for target range-angle estimation,” IEEE Transactions on Aerospace and Electronic Systems, vol. 50, no. 4, pp. 3057–3067, 2014.

S. Qin, Y. D. Zhang, M. G. Amin, and F. Gini, “Frequency diverse coprime arrays with coprime frequency offsets for multitarget localization,” IEEE Journal of Selected Topics in Signal Processing, vol. 11, no. 2, pp. 321–335, 2017.

W.-Q. Wang and C. Zhu, “Nested array receiver with time-delayers for joint target range and angle estimation,” IET Radar, Sonar & Navigation, vol. 10, no. 8, pp. 1384–1393, 2016.

S. Qin and Y. M. D. Zhang, “Frequency diverse array radar for target range-angle estimation,” International Journal of Computations and Mathematics in Electrical, vol. 35, no. 3, pp. 1257–1270, 2016.

P. F. Sammartino, C. J. Baker, H. D. Griffiths et al., “Range-angle dependent waveform,” in Proceedindgs of the IEEE Radar Conference, pp. 511–515, Washington, DC, USA, May 2010.

Y.-Q. Yang, H. Wang, H.-Q. Wang, S.-Q. Gu, D.-L. Xu, and S.-L. Quan, “Optimization of sparse frequency diverse array with time-invariant spatial-focusing beampattern,” IEEE Antennas and Wireless Propagation Letters, vol. 17, no. 2, pp. 351–354, 2018.

A. M. Yao, W. Wu, D. G. Fang et al., “Frequency diverse array radar with time-dependent frequency offset,” IEEE Antennas and Wireless Propagation Letters, vol. 13, pp. 758–761, 2014.

D. Nion and N. D. Sidiropoulos, “Adaptive algorithms to track the PARAFAC decomposition of a third-order tensor,” IEEE Transactions on Signal Processing, vol. 57, no. 6, pp. 2299–2310, 2009.

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.

B. Xu, Y. Zhao, Z. Cheng, H. Li et al., “A novel unitary PARAFAC method for DOD and DOA estimation in bistatic MIMO radar,” Signal Processing, vol. 138, no. 11, pp. 273–279, 2017.

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 Speech and Signal & Processing, pp. 2077–2080, Taipei, Taiwan, April 2009.

F. Wen, X. Xiong, and Z. Zhang, “Angle and mutual coupling estimation in bistatic MIMO radar based on PARAFAC decomposition,” Digital Signal Processing, vol. 65, pp. 1–10, 2017.

J. F. Li and X. F. Zhang, “A method for joint angle and array gain-phase error estimation in bistatic multiple-input multiple-output non-linear arrays,” IET Signal Processing, vol. 8, no. 2, pp. 131–137, 2014.

B. Li, W. Bai, G. Zheng et al., “Successive ESPRIT algorithm for joint DOA-range-polarization estimation with polarization sensitive FDA-MIMO radar,” IEEE Access, vol. 6, pp. 36376–36382, 2018.

F. Liu, X. Wang, M. Huang et al., “A novel unitary ESPRIT algorithm for monostatic FDA-MIMO radar,” Sensors, vol. 20, no. 3, 827 pages, 2020.

X. Zhang, Z. Xu, and D. Xu, “Trilinear decomposition-based transmit angle and receive angle estimation for multiple-input multiple-output radar,” IET Radar, Sonar & Navigation, vol. 5, no. 6, pp. 626–631, 2011.