Next Article in Journal
Assessment and Prediction of Impact of Flight Configuration Factors on UAS-Based Photogrammetric Survey Accuracy
Previous Article in Journal
A Robust Sparse Imaging Algorithm Using Joint MIMO Array Manifold and Array Channel Outliers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Real Valued MUSIC Method for Height Measurement of Meter Wave Polarimetric MIMO Radar Based on Matrix Reconstruction

Air and Missile Defense College, Air Force Engineering University, Xi’an 710051, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(16), 4121; https://doi.org/10.3390/rs14164121
Submission received: 11 July 2022 / Revised: 15 August 2022 / Accepted: 15 August 2022 / Published: 22 August 2022

Abstract

:
Combining the advantages of diversity provided by polarization MIMO radar and good decoherence ability of matrix reconstruction technology, a method for height measurements based on matrix reconstruction after real valued processing is developed. To solve height measurement problem in meter wave polarization MIMO radar, we first derive the corresponding flat ground signal model; then, the received data matrix is reconstructed to eliminate the influence of multipath coherent signal on height measurements. Then, the reconstructed data matrix is transformed into a real valued matrix using a unitary matrix. In order to reduce the influence of noise on the signal subspace and reduce the data dimension, singular value decomposition technology is applied to receive the signal data. Finally, the elevation and height of the target are estimated according to the principle that the signal subspace is orthogonal to the noise subspace. The proposed method does not require prior knowledge, such as the reflection coefficient, wave path difference and polarization information. Simulation experiments show that the proposed algorithm has better estimation performance and less computational complexity than conventional algorithms.

Graphical Abstract

1. Introduction

Most existing radars are single-polarization array radars, such as MIMO radar, which can only receive one type of polarization information [1,2,3] of the electromagnetic wave signal. The vector sensor radar, namely the polarization radar, has polarization diversity characteristics and can obtain at least two types of polarization information or up to six types of polarization information [4,5,6,7,8,9,10] of the electromagnetic wave signal. Refs. [4,5,6,7] studied the angle estimation of the specially separated electromagnetic vector sensor array, and showed that the polarization information can improve the performance of angle estimation.
Refs. [8,9,10] studied the array popularity of large electric dipoles and magnetic loops and proposed corresponding angle and polarization estimation algorithms. It is well known that direction of arrival (DOA) estimation is an important part of array signal processing [11,12], in particular the height measurement of low altitude targets, which is affected by multipath and is difficult to measure accurately [13,14,15,16,17,18]. Severe multipath coherence signals exist in the low-elevation areas, which seriously affects the height measurement accuracy.
Therefore, polarimetric MIMO radar with both spatial diversity and polarization diversity has attracted the attention of many scholars [19,20,21,22]. In order to make full use of the advantages of polarimetric MIMO radar, Ref. [23] proposed a polarized smoothing generalized multiple signal classification (MUSIC) algorithm based on the generalized MUSIC algorithm, which improved the height measurement accuracy. Based on the MUSIC algorithm [24], Zheng [25] proposed the steering vector synthesis MUSIC algorithm and generalized MUSIC algorithms for the meter-wave polarimetric MIMO radar, which take full advantage of the waveform diversity and polarization diversity with good estimation accuracy and without decoherence processing.
The paper proposes the idea of matrix reconstruction as well as the operation of real-value processing to further improve the estimation accuracy and reduce the computational complexity. The work is motivated to improve the angle estimation accuracy with reduced computational complexity. The problem of the existing algorithms is that the improvement of the angle measurement accuracy and the amount of calculations are contradictory and cannot be unified. The challenges of the topic are to better apply and combine these technologies.
In order to make better use of the advantages of polarimetric MIMO radar, we first analyze the signal model of meter wave polarimetric MIMO radar in the low-elevation area and then study a method based on matrix reconstruction without any prior information. By taking full advantage of the waveform diversity and polarization diversity advantages of the polarimetric MIMO radar, the proposed method has good accuracy at the case of low snapshots or a low signal-to-noise ratio (SNR).
Finally, the proposed algorithm is further analyzed through computer simulation results, and then the correctness and advantages of the proposed algorithm are verified after comparison with the existing algorithms. The innovation of this paper is to apply the polarization MIMO radar to the height measurement of complex terrain. The matrix reconstruction in this paper does not require the prior information of the terrain, and the real value technology can improve the accuracy of height measurement and reduce the amount of calculation.
Notations: Superscript ( ) , ( ) T and ( ) H denote conjugate, transpose and conjugate transpose, respectively. denotes Kronecker product, respectively. I M is an identity matrix with dimension M . Π is the exchange matrix with ones on its anti-diagonal and zeros elsewhere. det ( ) denotes the matrix determinant.

2. Signal Model of Meter Wave Polarimetric MIMO Radar in Low Elevation Area

The schematic diagram of the polarimetric MIMO radar is shown in Figure 1. The transceiver antennas are all vertically placed with polarization sensitive antenna, and its adjust sensor spacing equals d t = d r = λ 2 , where λ is the wavelength. θ d is the direct wave angle, θ s is the reflection angle, h a is the height of the array, and h t is the height of the target.
There are M vector sensors in the polarimetric MIMO radar. The transmitted signals of polarimetric MIMO radar are assumed to be a set of orthogonal signals S = S 1 T , , S M T T C 6 M × 1 , where S m = s m , 1 , , s m , 6 T , which satisfies the following Equation (1).
s m , p s n , q = P , m = n a n d p = q 0 , o t h e r w i s e ; m , n = 1 , M ;   p , q = 1 , , 6 .
Assuming that the low elevation reflection area is smooth and flat ground, the signal of l th snapshot reaching at the target can be represented as
x ( l ) = [ b t ( θ d ) + e j δ ρ h b t ( θ s ) ] T S
where δ = 4 π h a h t R λ is the phase difference due to the range difference of wave propagation between the direct wave and reflected wave, R is the distance between the projection of the target on the ground and the radar, and ρ h is the Fresnel reflection coefficient of the horizontally polarized wave with a value of
ρ h = sin θ d ε cos 2 θ d sin θ d + ε cos 2 θ d
where ε is the surface complex dielectric constant, which can be represented by the relative dielectric constant ε r and the surface material conductivity σ e : ε = ε r j 60 λ σ e . b ( θ ) is the whole steering vector, which equals b ( θ ) = a ( θ ) g θ , where a ( θ ) is the spacial steering vector and g θ is the polarization steering vector of a single electromagnetic vector sensor, which can be expressed as:
a ( θ ) = 1 , exp ( j π sin ( θ ) ) , , exp ( j ( M 1 ) π sin ( θ ) ) T
g θ = cos ϕ cos θ sin ϕ sin ϕ cos θ cos ϕ k sin θ 0 sin ϕ cos ϕ cos θ cos ϕ sin ϕ cos θ 0 sin θ sin γ e j η cos γ
where θ represents θ d or θ s , and b r ( θ ) = b t ( θ ) . Therefore, we can obtain the values of b r ( θ d ) , b r ( θ s ) , b t ( θ d ) , b t ( θ s ) . As elevation measurement is the primary purpose of this paper, we set the azimuth as ϕ = 90 . Data for the l th snapshot received by the entire array can be expressed as
X l = [ b r ( θ d ) + e j δ ρ h b r ( θ s ) ] ξ l [ b t ( θ d ) + e j δ ρ h b t ( θ s ) ] T S + N l
where ξ l is the complex reflection coefficient, and N l is the Gaussian white noise, with a mean of zero and a variance of σ 2 . The following formula can be obtained after matching filtering using the transmitted signals to Equation (6):
Y l = [ b r ( θ d ) + e j δ ρ h b r ( θ s ) ] ξ l [ b t ( θ d ) + e j δ ρ h b t ( θ s ) ] T + N l S H l
For Equation (7), the received data under the l th snapshot can be expressed as the following formula,
Y = [ b t ( θ d ) + e j δ ρ h b t ( θ s ) ] [ b r ( θ d ) + e j δ ρ h b r ( θ s ) ] ξ l + N C 36 M N × L
According to the conclusion of Ref. [21], the noise of N is still Gaussian white noise after matching filtering and vectorization operation.

3. Real Valued MUSIC Height Measurement Method Based on Matrix Reconstruction

To solve the problem of the effect of the multi-path reflection echo signal in the low-elevation area of the polarimetric MIMO radar, this section reconstructs the received signal data Equation (8) as matrix Z l C M M × 36 , which can be expressed as the following equation:
Z l = Y 1 , 1 l , , Y 1 , M l , Y 2 , 1 l , , Y M , M l
where Y m , n l C 36 × 1 represents the data received by the n th vector sensor and the m th transmitted signal. Equation (9) can be expanded as the following formula:
Z l = A Λ ξ l G + V l = A D l + V l
where A is the steering vector after the matrix reconstruction, V l is the noise matrix of the reconstructed data, G contains all the polarization information of the received data, and Λ ξ l includes the reflection coefficient, the range difference of wave propagation and other information as follows.
A = a r θ d a t θ d a r θ d a t θ s a r θ s a t θ d a r θ s a t θ s C M M × 4
G = g θ d g θ d g θ d g θ s g θ s g θ d g θ s g θ s T C 4 × 36
Λ ξ l = d i a g ξ l , e j δ ρ h ξ l , e j δ ρ h ξ l , e j 2 δ ρ h 2 ξ l
After observation, it is not difficult to find that the column of G is linear independent and, when ξ l 0 , R D = E D l D l H , has rank min 4 , 36 = 4 . Furthermore, it is important to note that, in Equation (10), when the rank of A is 4, the rank loss generated by the multipath coherence signal has been resolved, i.e., the matrix reconstruction method eliminates the effect of the multipath coherence signal. With K noncoherent low elevation targets, the proposed method can solve the multipath coherent rank loss phenomenon for, at most, 36 4 .
However, the values in Equation (10) are complex values, and thus direct low-elevation estimation leads to greater computational complexity. To reduce the algorithm complexity, we perform the real-value processing using the unitary matrix with Equation (10) and derive the corresponding real-value MUSIC elevation estimation method.
When the received data has L snapshots, the corresponding data expression of Z L = Z 1 , , Z L C M M × 36 L can be obtained by Equation (10). To make full use of the conjugated data, Z L can be extended as
Z L , U = Z L , Π M M Z L Π 36 L C M M × 72 L
where Π N is a commutative matrix with 1 on the antidiagonal elements and the remaining position elements of 0. We can easily verify that Z L , U is the centrosymmetric Hermite matrix [26], and the real-valued matrix can be obtained by processing with Equation (14) using the unitary matrix [27].
Z R = Γ M M H Z k Γ 72 L
where Γ K is a sparse unitary matrix, whose odd and even dimensions are defined as, respectively:
Γ 2 K + 1 = 1 2 I K 0 K × 1 j I K 0 1 × K 2 0 1 × K Π K 0 K × 1 j Π K Γ 2 K = 1 2 I K j I K Π K j Π K
When the number of dimension of received data is large, the huge calculation amount will increase the practical application difficulty. In order to reduce the dimension of the receiving data, and to reduce the impact of the noise on the data, we can use Singular Value Decomposition (SVD) technology to handle Equation (15), as follows:
Z R = U s Λ s V s H + U n Λ n V n H C M M × 72 L
As there are four transmission paths in the low elevation region for one target and the reconstructed data covariance matrix rank is 4; therefore, U s C M N × 4 and U n C M N × M N 4 are composed of left singular value vectors corresponding to four large singular values and left singular values corresponding to M N 4 small singular values. Λ s and Λ n consist of four large singular values and the remaining M N 4 large singular values. Similarly, V s C 72 L × 4 and V n C 72 L × M N 4 consists of right singular value vectors corresponding to four large singular values and a vector of right singular values corresponding to M N 4 small singular values. Equation (15) right is multiplied with V s , and then we have
Z R = Z R V s C M N × 4
According to Equation (18), the dimension of matrix Z R decreases from 72 L to 4. It is not difficult to find 72 L 4 when the number of snapshots L is large, and thus it is not difficult to prove that singular value techniques can greatly reduce the matrix dimension.
In the same way, the covariance matrix of the real-valued data can be obtained by the following equation:
R Z R = 1 72 L Z R Z R H C M M × M M
Eigenvalue decomposition of matrix R Z R yields the following formula.
R Z R = E s Λ s E s H + E n Λ n E n H
Matrices E s and E n represent the signal and noise subspaces, respectively, whose eigenvalues constitute the diagonal matrices Λ s and Λ n . Since R Z R is a real-valued matrix, both E s and E n are real-valued matrices. Based on the orthogonality of the signal subspace and noise subspace, the following spectrum formula can be obtained.
P ( θ d , θ s ) = det E n H E n det A R H E n E n H A R
Equation (21) is called the generalized MUSIC algorithm. The difference between it and the traditional MUSIC algorithm is that there is no need for information of the steering vector. It is sufficient to know the steering matrix because the steering matrix is orthogonal to the noise subspace and because the angle estimation can be obtained by using this point.
As the received data is processed as real-valued data by the unitary matrix, the steering vector A R = Γ M 2 H A for the spectral peak search in Equation (21) does not contain any polarization parameter, the range difference of wave propagation and reflection coefficient information; thus, it is not difficult to prove that the proposed method can estimate the target elevation value without the polarization information.
The above equation is a two-dimensional spectral peak search, which can be reduced to a one-dimensional search according to the relationship between the angle of direct wave and the reflected wave as shown below
θ s = arctan tan θ d + 2 h a / R
Based on the distance R and elevation estimations, the height of the target can be calculated as:
H R sin θ d + h a

4. Summary and Computational Complexity of the Proposed Algorithm

The algorithm steps studied in this paper are as follows:
  • Reconstruct the received signal data using Equation (9).
  • The unitary matrix is used to process the real-value of the reconstructed received signal matrix.
  • Reduce the effect of big dimension of received data and noise using the singular value technique.
  • The signal covariance matrix was obtained by Equation (19), and then the eigenvalue was decomposed to obtain the noise subspace matrix.
  • According to the noise subspace matrix, the spectral peak search was performed as in Equation (21).
  • The target height is obtained according to the geometric relationship of Equation (23).
According to the algorithm steps, it is not difficult to find that the complexity of the proposed method is mainly divided into the following parts: (1) real-valued processing; (2) estimate the real-valued covariance matrix; (3) eigenvalue decomposition of the covariance; and (4) spectral peak search. The specific expression is as follows: (1) 24 72 M 2 L ; (2) 72 M 4 L ; (3) M 6 ; (4) n u m 8 M 4 + 36 M 2 . Therefore, the complexity of the proposed algorithm is 36 M 2 48 L + n u m + M 4 72 L + 8 n u m + M 6 .
The generalized MUSIC algorithm proposed in Ref. [28] is abbreviated as G-MUSIC; the generalized MUSIC algorithm applied in polarimetric MIMO radar is abbreviated as PG-MUSIC [23]; the height measurement method of maximum likelihood estimation proposed in Ref. [29] is abbreviated as the ML algorithm; and its extension algorithm in polarimetric MIMO radar is abbreviated as P-ML [25]. The G-MUSIC algorithm complexity is 4 M 4 L + 4 M 6 + 4 Θ 32 M 2 + 8 M 4 .
The complexity of PG-MUSIC algorithm is 4 36 M 2 3 + 4 36 M 2 2 L + 4 Θ 36 32 M 2 + 8 36 2 M 4 . The complexity of ML algorithm is 4 M 4 L + M 2 + 4 Θ 32 M 2 + 4 M 4 + M 6 . The complexity of P-ML algorithm is 4 36 2 M 4 L + 36 M 2 + 4 Θ 36 3 M 6 + 32 36 M 2 + 4 36 2 M 4 .

5. Computer Simulation Results

This section verifies the reliability and superiority of the proposed algorithm mainly through simulation experiments. In this paper, the angle search range is set as 0.5 ° 12 ° , the search angle interval is 0.01 ° , and the number of searches is 1151 . A meter-wave polarimetric MIMO radar system has the transceiver array elements M = N = 6 , and the adjusted element spacing is d t = d r = λ / 2 . In addition, the wavelength is set as λ = 1 . The height array is h a = 5 m , and number of Monte Carlo experiments is K = 500 . In simulation experiments, the number of target is set to 1, and the root mean square error is defined as
RMSE _ DOA = 1 K k = 1 K θ k θ d 2
RMSE _ H = 1 K k = 1 K H k h t 2
where θ k is the estimated elevation value obtained by the k th experiment, H k is the estimated target height obtained by the k th experiment, and K is the number of Monte Carlo experiments.
Example 1: In this experiment, the number of L = 10 , the actual elevation of the target is 3 ° . The SNR is set as 10 dB. Figure 2 is a spatial spectrum of the proposed method with 10 times. It can be seen from Figure 2 that the proposed method can correctly measure the target elevation angle, which proves the correctness of the proposed method.
Example 2: In this experiment, the number of snapshots is L = 10 , the actual elevation of the target is 3 ° , the distance between the target and the radar is R = 200 km, and the range of SNR is from −10 to 10 dB. Figure 3 shows the relationship between the RMSE of angle estimation and SNR for the five algorithms, and Figure 4 shows the RMSR of height measurement and SNR of the five algorithms.
From Figure 3 and Figure 4, the PG-MUSIC and P-ML methods show better estimation accuracy than the G-MUSIC and ML methods, which demonstrates the polarization diversity advantage of polarizing-sensitive arrays. Moreover, the proposed algorithm outperforms the PG-MUSIC and P-ML methods, which demonstrates the superiority of the proposed algorithm.
Example 3: Figure 5 shows the relationship between the computational complexity and the number of array elements of the five algorithms. It can be found from Figure 5 that the computational complexity of the PG-MUSIC and P-ML methods is significantly higher than that for the G-MUSIC and ML methods. This indicates that, although polarimetric MIMO radar has a polarization diversity advantage, it also increases the algorithmic computational complexity.
Notably, the proposed algorithm has the lowest computational complexity because of the data with real values and because the steering vector dimension for spectrum search is lower than for the PG-MUSIC and P-ML methods. This proves that the proposed method is more conducive to engineering implementations. In addition, according to the experimental results of this example and example 2, compared with the other four algorithms, the proposed algorithm reduces the computing complexity and improves the accuracy of the target elevation and height measurement.

6. Conclusions

To further improve the measurement accuracy, we studied a height measurement method based on the meter-wave polarimetric MIMO radar, which makes full use of the waveform diversity and polarization diversity. First, the receiving signal model was analyzed. According to the height measurement signal model, the received data matrix was reconstructed to eliminate the effect of the multipath reflection echo. To reduce the algorithm complexity, the reconstructed data matrix was processed with real values using a unitary matrix.
Finally, the corresponding spectrum search was presented, which requires no known polarization information. In the simulation experiment, the proposed algorithm was comprehensively compared with the state-of-the-art algorithms, which showed that the proposed algorithm had the best performance. Therefore, when selecting the estimation algorithms, we can consider using the real value processing technology in this paper to improve the diversity of samples. In addition, the decoherence method of matrix reconstruction is a good decoherence method that does not depend on the information of the terrain.

Author Contributions

Conceptualization, G.Z. and C.C.; methodology, G.Z.; software, G.Z.; writing—original draft preparation, Y.S.; writing—review and editing, C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by National Natural Science Foundation of China under Grant 61971438.

Data Availability Statement

The data can be obtained by sending an email to [email protected].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, J.; Stoica, P. MIMO Radar with Collocated Antennas. IEEE Signal Process. Mag. 2007, 24, 106–114. [Google Scholar] [CrossRef]
  2. Liu, Y.B.; Wang, C.Y.; Gong, J.; Tan, M.; Chen, G. Robust suppression of deceptive jamming with VHF-FDA-MIMO radar under multipath effects. Remote Sens. 2022, 14, 942. [Google Scholar] [CrossRef]
  3. Luo, J.; Zhang, Y.; Yang, J.; Zhang, D.; Zhang, Y.; Zhang, Y.; Huang, Y.; Jakobsson, A. Online sparse DOA estimation based on sub–aperture recursive LASSO for TDM–MIMO radar. Remote Sens. 2022, 14, 2133. [Google Scholar] [CrossRef]
  4. Zheng, G.; Wu, B.; Ma, Y.; Chen, B. DOA Estimation with a sparse uniform array of orthogonally oriented and spatially separated dipole-triads. IET Radar Sonar Navig. 2014, 8, 885–894. [Google Scholar] [CrossRef]
  5. Zheng, G. Two-dimensional DOA estimation for polarization sensitive array consisted of spatially spread crossed-dipole. IEEE Sens. J. 2018, 18, 5014–5023. [Google Scholar] [CrossRef]
  6. Wong, K.T.; Yuan, X. “Vector cross-product direction-finding” with an electromagnetic vector-sensor of six orthogonally oriented but spatially noncollocating dipoles/loops. IEEE Trans. Signal Process. 2011, 59, 160–171. [Google Scholar] [CrossRef]
  7. Song, Y.; Wong, K.T.; Yuan, X. Correction to “Vector cross-product direction-finding” with an electromagnetic vector-sensor of six orthogonally oriented but spatially noncollocating dipoles/loops. IEEE Trans. Signal Process. 2014, 62, 1028–1030. [Google Scholar] [CrossRef]
  8. Wong, K.T.; Song, Y.; Fulton, C.J.; Khan, S.; Tam, W.Y. Electrically “Long” Dipoles in a collocated/orthogonal triad—For direction finding and polarization estimation. IEEE Trans. Antennas Propag. 2017, 65, 6057–6067. [Google Scholar] [CrossRef]
  9. Khan, S.; Wong, K.T.; Song, Y.; Tam, W.Y. Electrically large circular loops in the estimation of an incident emitter’s direction-of-arrival or polarization. IEEE Trans. Antennas Propag. 2018, 66, 3046–3055. [Google Scholar] [CrossRef]
  10. Khan, S.; Wong, K.T. A six component vector comprising electrically long dipoles and large loops –to simultaneously estimate incident sources’ direction of arrival and polarizations. IEEE Trans. Antennas Propag. 2020, 68, 6355–6363. [Google Scholar] [CrossRef]
  11. Guo, Y.; Hu, X.; Feng, W.; Gong, J. Low-Complexity 2D DOA Estimation and Self-Calibration for Uniform Rectangle Array with Gain-Phase Error. Remote Sens. 2022, 14, 3064. [Google Scholar] [CrossRef]
  12. Mao, Z.; Liu, S.; Qin, S.; Haung, Y. Cramér-Rao Bound of joint DOA-range estimation for coprime frequency diverse arrays. Remote Sens. 2022, 14, 583. [Google Scholar] [CrossRef]
  13. Liu, Y.; Liu, H. Target Height Measurement under Complex Multipath Interferences without Exact Knowledge on the Propagation Environment. Remote Sens. 2022, 14, 3099. [Google Scholar] [CrossRef]
  14. Zheng, G.; Song, Y. Signal model and method for joint angle and range estimation of low-elevation target in meter-wave FDA-MIMO radar. IEEE Commun. Lett. 2022, 26, 449–453. [Google Scholar] [CrossRef]
  15. Song, Y.; Hu, G.; Zheng, G. Height measurement with meter wave MIMO radar based on precise signal model under complex terrain. IEEE Access 2021, 9, 49980–49989. [Google Scholar] [CrossRef]
  16. Chen, C.; Tao, J.; Zheng, G.; Song, Y. Meter-wave MIMO radar height measurement method based on adaptive beamforming. Digit. Signal Process. 2022, 120, 103272. [Google Scholar] [CrossRef]
  17. Chen, C.; Tao, J.; Zheng, G.; Song, Y. Beam split algorithm for height measurement with meter-wave MIMO radar. IEEE Access 2021, 9, 5000–5010. [Google Scholar] [CrossRef]
  18. Wang, H.; Zheng, G.; Liu, Y.; Song, Y. Low elevation and range joint estimation method with meter wave FDA-MIMO radar based on sparse array. IET Radar Sonar Navig. 2022, 16, 1179–1187. [Google Scholar] [CrossRef]
  19. Jiang, H.; Zhang, J.-K.; Wong, K.M. Joint DOD and DOA Estimation for Bistatic MIMO Radar in Unknown Correlated Noise. IEEE Trans. Veh. Technol. 2015, 64, 5113–5125. [Google Scholar] [CrossRef]
  20. Chintagunta, S.; Ponnusamy, P. Integrated polarisation and diversity smoothing algorithm for DOD and DOA estimation of coherent targets. IET Signal Process. 2018, 12, 447–453. [Google Scholar] [CrossRef]
  21. Wen, F.; Shi, J.; Zhang, Z. Joint 2D-DOD, 2D-DOA And Polarization Angles Estimation for Bistatic EMVS-MIMO Radar Via PARAFAC Analysis. IEEE Trans. Veh. Technol. 2020, 29, 1626–1638. [Google Scholar] [CrossRef]
  22. Wang, X.; Wan, L.; Huang, M.; Shen, C.; Zhang, K. Polarization channel estimation for circular and non-circular signals in massive MIMO systems. IEEE J. Sel. Top. Signal Process. 2019, 13, 1001–1016. [Google Scholar] [CrossRef]
  23. Tan, J.; Nie, Z. Polarization Smoothing Generalized MUSIC algorithm with PSA monostatic MIMO radar for low angle estimation. Electron. Lett. 2018, 54, 527–529. [Google Scholar] [CrossRef]
  24. Stoica, P.; Nehorai, A. MUSIC, maximum likelihood, and Cramer-Rao bound. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 720–741. [Google Scholar] [CrossRef]
  25. Zheng, G.; Song, Y.; Chen, C. Height Measurement for with wave polarimetric MIMO radar: Signal model and MUSIC algorithm. Signal Process. 2022, 190, 108344. [Google Scholar] [CrossRef]
  26. Ding, J.C.; Chen, H.W.; You, B.; Li, X.; Zhuang, Z.W. Fluctuating target detection in low-grazing angle with Multi-input Multiple-output radar. In Proceedings of the 2012 IEEE 11th International Conference on Signal Processing, Beijing, China, 21–25 October 2012; IEEE: Piscataway, NJ, USA, 2013. [Google Scholar]
  27. Li, J.; Jiang, D.; Zhang, X. DOA Estimation Based on Combined Unitary ESPRIT for Coprime MIMO Radar. IEEE Commun. Lett. 2017, 21, 96–99. [Google Scholar] [CrossRef]
  28. Haimovich, A.M.; Blum, R.S.; Cimini, L.J. MIMO Radar with Widely Separated Antennas. IEEE Signal Process. Mag. 2007, 25, 116–129. [Google Scholar] [CrossRef]
  29. Rahamim, D.; Tabrikian, J.; Shavit, R. Source localization using vector sensor array in a multipath environment. IEEE Trans. Signal Process. 2004, 52, 3096–3103. [Google Scholar] [CrossRef]
Figure 1. Signal model of meter wave polarimetric MIMO radar in a low-elevation area.
Figure 1. Signal model of meter wave polarimetric MIMO radar in a low-elevation area.
Remotesensing 14 04121 g001
Figure 2. The spatial spectrum of the proposed algorithm with 10 times.
Figure 2. The spatial spectrum of the proposed algorithm with 10 times.
Remotesensing 14 04121 g002
Figure 3. RMSE of angle estimation versus SNR.
Figure 3. RMSE of angle estimation versus SNR.
Remotesensing 14 04121 g003
Figure 4. RMSE of height measurement versus SNR.
Figure 4. RMSE of height measurement versus SNR.
Remotesensing 14 04121 g004
Figure 5. Complexity changes with the array number.
Figure 5. Complexity changes with the array number.
Remotesensing 14 04121 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zheng, G.; Chen, C.; Song, Y. Real Valued MUSIC Method for Height Measurement of Meter Wave Polarimetric MIMO Radar Based on Matrix Reconstruction. Remote Sens. 2022, 14, 4121. https://doi.org/10.3390/rs14164121

AMA Style

Zheng G, Chen C, Song Y. Real Valued MUSIC Method for Height Measurement of Meter Wave Polarimetric MIMO Radar Based on Matrix Reconstruction. Remote Sensing. 2022; 14(16):4121. https://doi.org/10.3390/rs14164121

Chicago/Turabian Style

Zheng, Guimei, Chen Chen, and Yuwei Song. 2022. "Real Valued MUSIC Method for Height Measurement of Meter Wave Polarimetric MIMO Radar Based on Matrix Reconstruction" Remote Sensing 14, no. 16: 4121. https://doi.org/10.3390/rs14164121

APA Style

Zheng, G., Chen, C., & Song, Y. (2022). Real Valued MUSIC Method for Height Measurement of Meter Wave Polarimetric MIMO Radar Based on Matrix Reconstruction. Remote Sensing, 14(16), 4121. https://doi.org/10.3390/rs14164121

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop