Wireless Communications and Mobile Computing
Volume 2021 (2021), Article ID 9973253, 10 pages
https://doi.org/10.1155/2021/9973253

Research Article

Source Localization Using Distributed Electromagnetic Vector Sensors

Tingping Zhang,¹ Di Wan,¹ Xinhai Wang,² and Fangqing Wen³

¹School of Information Science & Engineering, Chongqing Jiaotong University, Chongqing 400074, China
²Nanjing Marine Radar Institute, Nanjing 211153, China
³Electronic and Information School of Yangtze University, Jingzhou 434023, China

Correspondence should be addressed to Tingping Zhang; ztp@cqjtu.edu.cn

Received 21 March 2021; Accepted 27 April 2021; Published 11 May 2021

Academic Editor: Jian Feng Li

Copyright © 2021 Tingping Zhang 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.

Electromagnetic vector sensor (EVS) array has drawn extensive attention in the past decades, since it offers two-dimensional direction-of-arrival (2D-DOA) estimation and additional polarization information of the incoming source. Most of the existing works concerning EVS array are focused on parameter estimation with special array architecture, e.g., uniform manifold and sparse arrays. In this paper, we consider a more general scenario that EVS array is distributed in an arbitrary geometry, and a novel estimator is proposed. Firstly, the covariance tensor model is established, which can make full use of the multidimensional structure of the array measurement. Then, the higher-order singular value decomposition (HOSVD) is adopted to obtain a more accurate signal subspace. Thereafter, a novel rotation invariant relation is exploited to construct a normalized Poynting vector, and the vector cross-product technique is utilized to estimate the 2D-DOA. Based on the previous obtained 2D-DOA, the polarization parameter can be easily achieved via the least squares method. The proposed method is suitable for EVS array with arbitrary geometry, and it is insensitive to the spatially colored noise. Therefore, it is more flexible than the state-of-the-art algorithms. Finally, numerical simulations are carried out to verify the effectiveness of the proposed estimator.

1. Introduction

Direction-of-arrival (DOA) estimation using sensor array is one of the most popular methods for source localization, and it has struck a series of technical elevations in wireless communications [1–3], radars, sonars, etc. As a nonlinear issue, many efforts have been devoted and many superresolution methodologies have been developed, e.g., the estimation method of signal parameters via rotational invariance technique (ESPRIT) [4, 5], maximum-likelihood [6], multiple signal classification (MUSIC) [7–9], sparsity-aware strategy [10, 11], and tensor approach [12–14]. Nevertheless, most of the current works are concentrating on one-dimensional (1D) DOA estimation. A small number of works have paid attention to the two-dimensional (2D) DOA estimation problem, by cooperating with nonlinear scaler arrays [15, 16], e.g., L-shaped array, circular array, and rectangular array. A common characteristic of these arrays is that sensors must be strictly nonlinear. Besides, the sensors’ locations must be known priors, otherwise algorithms will fail to work. As an alternative of nonlinear scaler array, electromagnetic vector sensor (EVS) array has drawn more and more interest in the past decades. Unlike scaler arrays, a signal EVS consists of six collocated components, i.e., three orthogonally oriented dipoles and three orthogonally oriented loops, which measure the electric field and magnetic field, respectively. EVS array can provide not only 2D-DOA estimation but also the polarization status of the incoming sources, which maybe particularly important in detecting weak signal [17].

Many attempts have been done with respect to 2D-DOA estimation using EVS array, and various estimators have been proposed. The most popular method is the vector cross product [18], which is based on the Poynting vector of the polarization steering vector. The vector cross-product technique is suitable for a single EVS. Another important branch is combining the traditional subspace method with vector cross-product. In [19], MUSIC-like algorithm was presented. However, it is computationally inefficient owing to high-dimensional spectrum grid search. To avoid such drawback, the ESPRIT-like approach was investigated in [20]. It firstly obtains the signal subspace via eigendecomposition, and then, a normalized Poynting vector is constructed; finally, the vector cross product is performed to achieve 2D-DOA estimation. ESPRIT-like approach is suitable for arbitrary array geometry, and it offers closed-form solution for parameter estimation. Similar works concerning subspace methods and vector cross-product technique have been done in [21–24]. To further lower the computational burden, the propagator method was carried out in [25], which replaces the eigendecomposition with least squares. To exploit the tensor nature, the tensor-aware estimators were derived in [26, 27]. Since the tensor algebra can suppress the noise much better than the matrix-based tools, tensor methods provide more accurate estimation performance than the matrix algorithms. Some efforts have been devoted to further improve the estimation accuracy or reduce the mutual coupling effort. In [28], the coprime array (see [29] and the references therein) geometry was introduced into EVS, which occupies much larger aperture than the traditional uniform manifold. In [30, 31], the EVS architecture has been redesigned. The component number in the improved EVS is much smaller than six as to release the mutual coupling. In addition, the EVS array has been extended to MIMO radar in [32–36], which bring new insights into parameter estimation to MIMO systems. It should be emphasized that most of the existing algorithms are only effective to specific array geometry, e.g., uniform linear array and coprime array. In the presence of limited space, these array geometries may be not practical. Besides, none of the existing works has taken the spatially colored noise into consideration.

In this paper, we generalize the issue of source localization using the EVS array. We consider that the EVSs are displaced in arbitrary geometry. Moreover, spatially colored noise appears in the array. An improved higher-order singular value decomposition (HOSVD) estimator is proposed. Firstly, the temporal cross-correlation scheme is utilized to eliminate the spatially colored noise. Thereafter, the crosscovariance measurement is arranged into a fourth-order tensor. The HOSVD is adopted to obtain an enhanced signal subspace. Then, a normalized Poynting vector is constructed, and 2D-DOAs are obtained via the shift invariance property of the previously constructed Poynting vector. Finally, the polarization parameters are achieved via the least squares technique. Our algorithm is insensitive to the sensors’ position and spatially colored noise. Numerical examples are given to show the superiority of the proposed algorithm.

1.1. Notations

Throughout the paper, lowercase letters represent vectors and uppercase letters denote matrices, respectively; is a identity matrix; , , and denote the Kronecker product, the Khatri-Rao product, and the vector cross product, respectively; the superscript , , , and indicate the operations of transpose, Hermitian transpose, inverse, and pseudoinverse, respectively; denotes a diagonal matrix with the diagonal entities which are the th row of ; accounts for the th entity of the vector ; accounts for the Frobenius norm; is to get the phase; and returns the mathematical expectation.

2. Preliminaries and Problem Formulation

2.1. Tensor Preliminaries

Let us first introduce some preliminaries concerning tensors. Interested readers are recommend to refer to [37] for more details.

Definition 1 (tensor unfolding). Let be an -order tensor. The -mode tensor unfolding of is denoted by , which is a matrix. The th entity of is the th element of , , and.

Definition 2 (mode- tensor-matrix product). The mode- product of by a matrix , is denoted by , where and .

Definition 3 (important principles of the mode- tensor-matrix product). The following principles are important for mode- tensor-matrix product:

Definition 4 (HOSVD). The HOSVD of an -order tensor is denoted by

where is the so-called core tensor and is the unitary matrix. In mode- tensor unfolding format, (2) is rewritten as

2.2. Preliminaries of EVS

A single EVS consists of six-component electromagnetic antennas. For an incoming source that impinges on an isolated EVS, the electromagnetic response of the EVS is given by [20]

where , , , and are the elevation angle, azimuth angle, auxiliary polarization angle, and polarization phase difference of the source, respectively; and are measurements that sense the electric field and magnetic field, respectively. is the direction parameter-depend only matrix, and is the polarization parameter-depend only vector.

An important characteristic is that and are mutually orthogonal to each other. Moreover, the normalized Poynting-vector fulfills [20]

2.3. Problem Formulation

In this paper, we consider an -element EVS array with arbitrary geometry, as illustrated in Figure 1. Suppose that the th EVS is located at coordinate . Assume that there are far-field source signals that impinge on the EVS array. Let be an angle pair that parameterized the th () source signal. The received array signal from the EVS array can be written as [20]

where is the spatial steering vector of the th source signal, is the carrier wavelength, and . is the associated polarization steering vector, is the th source signal, and is the source signal vector. is the array noise, which is assumed to be spatially correlated, temporally independent Gaussian process, i.e.,

where is the covariance matrix of the noise, which is a Hermitian matrix. In the presence of Gaussian white noise, .

Figure 1: Illustration of the EVS array with arbitrary geometry.

3. The Proposed Algorithm

3.1. Spatially Colored Noise Suppression

Suppose the source signals are uncorrelated and is uncorrelated with ; then, the covariance of is given by

where is the covariance matrix of . In the presence of spatially colored noise, the signal subspace will be corrupted by the noise subspace; thus, the performance of the traditional subspace methods would be decreased or even fail to work. Therefore, it is necessary for us to eliminate the effect of the spatially colored noise. As illustrated in (7), the noise between different snapshots is uncorrelated; thus, we can construct the following cross-correlation matrix:

where is the cross-correlation matrix of the source signal. In the presence of temporally correlated source signal, e.g., sinusoidal signal, is a diagonal matrix. When snapshots are available, can be estimated via

Obviously, the noise is eliminated in . In practice, can be approximated by its truncated eigenvalue decomposition as

where is the dominate eigenvalue matrix, consist of the K dominate eigenvectors.

3.2. HOSVD

Actually, can be arranged as a fourth-order tensor as [35]

The HOSVD of is given by

where is the core tensor, , , , and are the unitary matrices, which consist of the left singular vectors of the mode-1, mode-2, mode-3, and mode-4 unfolding of , respectively. Since the rank of is , we can approximate by the truncated HOSVD as [35]

where , , , and consist of the eigenvectors associated to the dominate eigenvalues. is called the signal component of , which is given by

Inserting into (14) yields

By Hermitian unfolding , we can construct an enhanced covariance matrix as

Combined with the result of (11), we can obtain

Since is Hermitian, we have and . Performing eigendecomposition on , we can obtain a signal subspace as

Since the virtual response matrix spans the same subspace with , we have

where is a nonsingular matrix.

3.3. 2D-DOA Estimation

Let . We can get the following relationship:

with

where . Define the following selection matrices:

where is the th row of . Therefore, the relationships in (21) becomes

Inserting (24) into (20) establishes

Equivalently,

According to the first line of (26), we can perform eigendecomposition of , the eigenvalues of which reveal the estimation of , and the associated eigenvectors form the estimation of . Compute the left parts of the reminder equations, and then, we can get the estimation of , , , and , respectively.

In fact, can be written in the following normalized format:

Then, we have

Since is a constant, it can be easily removed from the normalized operation. Let and . According to the property in (5), we have

Accordingly, we can get

Denote the estimates of as ; then, we can get the estimation of , , and , which are denoted by , , and , respectively. Then, 2D-DOA estimation can be achieved via

Obviously, the estimated 2D-DOA are automatically paired.

3.4. Polarization Parameter Estimation

It is easy to find

So the th column of is

According to (4), we should have

where and are the direction-only and polarization-only variables. Define

Then, can also be expressed as

Therefore, can be estimated via

the th column of which is denoted by . After the 2D-DOA have been estimated, we can construct the matrix according to (4); then, we can compute the following least squares:

Obviously, is the estimation of . Finally, the polarization parameters can be achieved via

Also, the polarization parameters are paired automatically.

4. Algorithm Analysis

Now, we have accomplished the estimation algorithm for distributed EVS array. To help the reader better understand the proposed algorithm, we summarize its algorithmic steps in Table 1.

Table 1: Algorithmic steps of the proposed algorithm.

Step No.
Operation

Step 1
Estimate the noiseless covariance matrix through (10)

Step 2
Rearrange the covariance measurement into a fourth-order tensor via (12), and perform HOSVD on

Step 3
Construct from (17), and then, get via the Hermitian unfolding of

Step 4
Perform eigendecomposition on to get the signal subspace

Step 5
Calculate , and perform eigendecomposition on it to get and

Step 6
To get , , , and via (26)

Step 7
Estimate via (30), and then, get the 2D-DOA through (31)

Step 8
Estimate via (37), and compute through (38). Finally, get the polarization parameters via (39)

4.1. Complexity

Next, we summarize the main complexity (number of complex multiplication that is required) of typical algorithms. The main complexity of ESPRIT in [20] is eigendecomposition, which is on the order . The main computational burden of the PARAFAC estimator in [28] is the iteration, which requires complex multiplications. In the proposed algorithm, the main complexity is the HOSVD, which is the same as that of the eigendecomposition. We listed the main complexity of the above algorithms in Table 2. From this perspective, the proposed approach is less efficient than PARAFAC. However, it should be noticed, which will be shown in the simulation section, that the proposed approach is suitable for colored noise, while the others are sensitive to the colored noise.

Table 2: Comparison of the complexity.

Method
Complexity

ESPRIT

PARAFAC

Proposed

4.2. Identifiability

The identifiability of the proposed approach is equal to the maximum , which is also the same as that of ESPRIT. According to (26), the proposed approach is effective if . Therefore, the proposed algorithm can identify sources at most. In comparison, PARAFAC can identify sources. Usually, is much larger than ; thus, PARAFAC provides much better identifiability than ESPRIT and the proposed approach.

5. Simulation Results

Numerical simulations are carried out to show the superiority of the proposed algorithm. In the following simulations, we assumed an EVS array with elements, the coordinates of the EVSs fulfill uniform distribution in . Assume that there are uncorrelated far-field sinusoidal signals with frequency which is 0.2 MHz, 0.4 MHz, and 0.85 MHz, respectively. The associated parameters are , , , and , respectively. samples are collected. The signal-to-noise ratio (SNR) in the simulation is defined as the ratio of signal power to noise power in (6). The spatially correlated noise is modeled as an autoregressive (AR) process with the coefficients . All the results rely on 200 independent trials.

Example 1. Scatter results of the proposed algorithm. In such simulation, we consider , , and L =1000 and SNR is set to 10 dB. The scatter results of 2D-DOA estimation and polarization parameter estimation are depicted in Figure 2. It is obvious to us that all the parameters can be accurately estimated and correctly paired. It is evident that the proposed algorithm is suitable for distributed EVSs.

Example 2. The root mean square error (RMSE) performance versus SNR. Herein, RMSE for a parameter is defined as

where is the estimation of in the th trial. In this simulation, we consider , , and L =1000. For comparison purposes, we add the performance of the traditional ESPRIT algorithm [20] and the PARAFAC algorithm [28]. Figure 3 presents the average RMSE on 2D-DOA estimation (marked with suffix “-d”) and average RMSE on polarization parameter estimation (marked with suffix “-p”) in the presence of white Gaussian noise. It is seen that the proposed algorithm provides better RMSE than ESPRIT, but a slight worse RMSE performance than PARAFAC. Figure 4 displays the RMSE results with colored noise. Obviously, all the algorithms achieve better RMSE with higher SNR. As expected, both ESPRIT and PARAFAC perform worse at low SNR regions since they are affected by the spatially colored noise. On the other hand, the proposed method offers much better performance than ESPRIT and PARAFAC at low SNR regions. However, all the algorithms provide very similar RMSE performance at high SNR regions (SNR is larger than 10 dB).

Example 3. The RMSE performance versus . In this simulation, we consider and , and SNR is fixed at 0 dB. The RMSE curves are illustrated in Figure 5. It seems that RMSE can be slightly improved with increasing. The ESPRIT method provides the worst RMSE performance. Since the tensor nature has been exploited, the PARAFAC estimator offers better RMSE than ESPRIT. However, the proposed method outperforms all the compared methods since it is free from the spatially colored noise.

Example 4. The RMSE performance versus colored parameter . In this simulation, we consider and and SNR is set to 0 dB. Figure 6 gives the RMSE curves. It should be pointed that when , the noise can be regarded as spatially uncorrelated (or with very weak relevance). Under such condition, all the algorithms offer very close RMSE performance. However, RMSE associated with ESPRIT and PARAFAC become worse when , since they are sensitive to the spatially colored noise. By comparison, the proposed algorithm offers much better RMSE than the compared methods. It is evident that the proposed method is insensitive to the spatially colored noise.

(a)

(b)

Figure 2: Scatter results of the proposed algorithm: (a) 2D-DOA estimation, (b) polarization parameter estimation.

(a)

(b)

Figure 3: Average RMSE performance comparison versus SNR with Gaussian white noise: (a) average RMSE on 2D-DOA estimation; (b) average RMSE on polarization parameter estimation.

(a)

(b)

Figure 4: Average RMSE performance comparison versus SNR with colored noise: (a) average RMSE on 2D-DOA estimation; (b) average RMSE on polarization parameter estimation.

(a)

(b)

Figure 5: Average RMSE performance comparison versus with colored noise: (a) average RMSE on 2D-DOA estimation; (b) average RMSE on polarization parameter estimation.

(a)

(b)

Figure 6: Average RMSE performance comparison versus : (a) average RMSE on 2D-DOA estimation; (b) average RMSE on polarization parameter estimation.

6. Conclusion

In this paper, a novel HOSVD algorithm is proposed for distributed EVS array, the core of which is to estimate 2D-DOA by the vector cross product of the normalized Poynting vector, which is achieved from cross-correlation tensor signal subspace. The proposed algorithm is better than the state-of-the-art algorithms since it is insensitive to the spatially colored noise and sensor position. Numerical simulation results verify the effectiveness of the proposed algorithm. It should have a bright future in radar, sonar, and wireless communication as well as Internet of Things.

Data Availability

There is no available data for this paper.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

This work was supported by the Humanities and Social Sciences projects of the Ministry of Education under grant 20YJAZH132 and Natural Science Foundation project of CQCST under grant cstc2018jcyjAX0398.

References

G. Gui, M. Liu, F. Tang, N. Kato, and F. Adachi, “6G: opening new horizons for integration of comfort, security and intelligence,” IEEE Wireless Communications Magazine, vol. 27, no. 5, pp. 126–132, 2020.

L. Sun, L. Wan, and X. Wang, “Learning-based resource allocation strategy for industrial IoT in UAV-enabled MEC systems,” IEEE Transactions on Industrial Informatics, vol. 17, no. 7, pp. 5031–5040, 2021.

L. Wan, L. Sun, K. Liu, X. Wang, Q. Lin, and T. Zhu, “Autonomous vehicle source enumeration exploiting non-cooperative UAV in software defined internet of vehicles,” IEEE Transactions on Intelligent Transportation Systems, 2020.

J. Chen, H. Gu, and W. Su, “Angle estimation using ESPRIT without pairing in MIMO radar,” Electronics Letters, vol. 44, no. 24, pp. 1422-1423, 2008.

J. Li, D. Jiang, and X. Zhang, “DOA estimation based on combined unitary ESPRIT for coprime MIMO radar,” IEEE Communications Letters, vol. 21, no. 1, pp. 96–99, 2017.

B. Tang, J. Tang, Y. Zhang, and Z. Zheng, “Maximum likelihood estimation of DOD and DOA for bistatic MIMO radar,” Signal Processing, vol. 93, no. 5, pp. 1349–1357, 2013.

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.

J. Li, D. Li, D. Jiang, and X. Zhang, “Extended-aperture unitary root MUSIC-based DOA estimation for coprime array,” IEEE Communications Letters, vol. 22, no. 4, pp. 752–755, 2018.

F. Wen, J. Wang, J. Shi, and G. Gui, “Auxiliary vehicle positioning based on robust DOA estimation with unknown mutual coupling,” IEEE Internet of Things Journal, vol. 7, no. 6, pp. 5521–5532, 2020.

J. Cong, X. Wang, M. Huang, and L. Wan, “Robust DOA estimation method for MIMO radar via deep neural networks,” IEEE Sensors Journal, vol. 21, no. 6, pp. 7498–7507, 2021.

X. Wang, L. T. Yang, D. Meng, M. Dong, K. Ota, and H. Wang, “Multi-UAV cooperative localization for marine targets based on weighted subspace fitting in SAGIN environment,” IEEE Internet of Things Journal, 2021.

F. Wen and C. Liang, “Improved tensor-mode based direction-of-arrival estimation for massive MIMO systems,” IEEE Communications Letters, vol. 19, no. 12, pp. 2182–2185, 2015.

F. E. D. Raimondi and P. Comon, “Tensor DOA estimation with directional elements,” IEEE Signal Processing Letters, vol. 24, no. 5, pp. 648–652, 2017.

A. Weiss, “Blind direction-of-arrival estimation in acoustic vector-sensor arrays via tensor decomposition and Kullback-Leibler divergence covariance fitting,” IEEE Transactions on Signal Processing, vol. 69, pp. 531–545, 2021.

J. Li, Y. Li, and X. Zhang, “Two-dimensional off-grid DOA estimation using unfolded parallel coprime array,” IEEE Communications Letters, vol. 22, no. 12, pp. 2495–2498, 2018.

P. Ma, J. Li, G. Zhao, and X. Zhang, “CAP-3 coprime array for DOA estimation with enhanced uniform degrees of freedom and reduced mutual coupling,” IEEE Communications Letters, p. 1, 2021.

L. Wan, K. Liu, Y. C. Liang, and T. Zhu, “DOA and polarization estimation for non-circular signals in 3-D millimeter wave polarized massive MIMO systems,” IEEE Transactions on Wireless Communications, p. 1, 2020.

A. Nehorai and E. Paldi, “Vector-sensor array processing for electromagnetic source localization,” IEEE Transactions on Signal Processing, vol. 42, no. 2, pp. 376–398, 1994.

J. Li, “Direction and polarization estimation using arrays with small loops and short dipoles,” IEEE Transactions on Antennas and Propagation, vol. 41, no. 3, pp. 379–387, 1993.

K. T. Wong and M. D. Zoltowski, “Closed-form direction finding and polarization estimation with arbitrarily spaced electromagnetic vector-sensors at unknown locations,” IEEE Transactions on Antennas and Propagation, vol. 48, no. 5, pp. 671–681, 2000.

K. T. Wong and M. D. Zoltowski, “Uni-vector-sensor ESPRIT for multisource azimuth, elevation, and polarization estimation,” IEEE Transactions on Antennas and Propagation, vol. 45, no. 10, pp. 1467–1474, 1997.

K. T. Wong and M. D. Zoltowski, “Self-initiating music-based direction finding and polarization estimation in spatio-polarizational beamspace,” IEEE Transactions on Antennas and Propagation, vol. 48, no. 8, pp. 1235–1245, 2000.

M. D. Zoltowski and K. T. Wong, “ESPRIT-based 2-D direction finding with a sparse uniform array of electromagnetic vector sensors,” IEEE Transactions on Signal Processing, vol. 48, no. 8, pp. 2195–2204, 2000.

S. Miron, N. Le Bihan, and J. Mars, “Vector-sensor MUSIC for polarized seismic sources localization,” EURASIP Journal on Advances in Signal Processing, vol. 2005, no. 1, 2005.

J. He and Z. Liu, “Computationally efficient 2D direction finding and polarization estimation with arbitrarily spaced electromagnetic vector sensors at unknown locations using the propagator method,” Digital Signal Processing, vol. 19, no. 3, pp. 491–503, 2009.

X. Zhang and D. Xu, “Blind PARAFAC signal detection for polarization sensitive array,” EURASIP Journal on Advances in Signal, vol. 2007, no. 1, 2007.

M. Y. Cao, X. Mao, X. Long, and L. Huang, “Tensor approach to DOA estimation of coherent signals with electromagnetic vector-sensor array,” Sensors, vol. 18, no. 12, p. 4320, 2018.

A. Tanveer, X. Zhang, and W. Zheng, “DOA estimation for coprime EMVS arrays via minimum distance criterion based on PARAFAC analysis,” IET Radar, Sonar & Navigation, vol. 13, no. 1, pp. 65–73, 2019.

J. Shi, F. Wen, and T. Liu, “Nested MIMO radar: coarrays, tensor modeling and angle estimation,” IEEE Transactions on Aerospace and Electronic Systems, vol. 57, no. 1, pp. 573–585, 2021.

J. He, L. Li, and T. Shu, “Sparse nested arrays with spatially spread square acoustic vector sensors for high accuracy underdetermined direction finding,” IEEE Transactions on Aerospace and Electronic Systems, p. 1, 2021.

J. He, L. Li, and T. Shu, “Sparse nested arrays with spatially spread orthogonal dipoles: high accuracy passive direction finding with less mutual coupling,” IEEE Transactions on Aerospace and Electronic Systems, p. 1, 2021.

S. Chintagunta and P. Palanisamy, “2D-DOD and 2D-DOA estimation using the electromagnetic vector sensors,” Signal Processing, vol. 147, pp. 163–172, 2018.

F. Wen, J. Shi, and Z. Zhang, “Joint 2D-DOD, 2D-DOA, and polarization angles estimation for bistatic EMVS-MIMO radar via PARAFAC analysis,” IEEE Transactions on Vehicular Technology, vol. 69, no. 2, pp. 1626–1638, 2020.

F. Wen, J. Shi, and Z. Zhang, “Fast direction finding for bistatic EMVS-MIMO radar without pairing,” Signal Processing, vol. 173, article 107512, 2020.

X. Wang, M. Huang, and L. Wan, “Joint 2D-DOD and 2D-DOA estimation for coprime EMVS-MIMO radar,” Circuits, Systems, and Signal Processing, 2021.

F. Wen, J. Shi, and Z. Zhang, “Closed-form estimation algorithm for EMVS-MIMO radar with arbitrary sensor geometry,” vol. 186, Article ID 108117, 2021.

M. Haardt, F. Roemer, and G. Del Galdo, “Higher-order SVD-based subspace estimation to improve the parameter estimation accuracy in multidimensional harmonic retrieval problems,” IEEE Transactions on Signal Processing, vol. 56, no. 7, pp. 3198–3213, 2008.