Meter Wave Polarization-Sensitive Array Radar for Height Measurement Based on MUSIC Algorithm

Obtaining good measurement performance with meter wave radar has always been a difficult problem. Especially in low-elevation areas, the multipath effect seriously affects the measurement accuracy of meter wave radar. The generalized multiple signal classification (MUSIC) algorithm is a well-known measurement method that dose not require decorrelation processing. The polarization-sensitive array (PSA) has the advantage of polarization diversity, and the polarization smoothing generalized MUSIC algorithm demonstrates good angle estimation performance in low-elevation areas when based on a PSA. Nevertheless, its computational complexity is still high, and the estimation accuracy and discrimination success probability need to be further improved. In addition, it cannot estimate the polarization parameters. To solve these problems, a polarization synthesis steering vector MUSIC algorithm is proposed in this paper. First, the MUSIC algorithm is used to obtain the spatial spectrum of the meter wave PSA. Second, the received data are properly deformed and classified. The Rayleigh–Ritz method is used to decompose the angle to realize the decoupling of polarization and the direction of the arrival angle. Third, the geometric relationship and prior information of the direct wave and the reflected wave are used to continue dimension reduction processing to reduce the computational complexity of the algorithm. Finally, the geometric relationship is used to obtain the target height measurement results. Extensive simulation results illustrate the accuracy and superiority of the proposed algorithm.


Introduction
In recent years, with the emergence of stealth aircraft, anti-radiation missiles, and other weapons, meter wave radar, which can detect stealth targets and anti-radiation missiles, has received extensive attention [1,2]. However, due to meter wave radar's wide beam, low band, and long wavelength, the ground-reflected echo cannot be ignored in low-elevation areas, which leads to the echo signal having a low signal-to-noise ratio and serious multipath coherence, thereby reducing the target detection ability of meter wave radar [3,4]. It is known that the essence of using radar to measure the target height is to estimate the target elevation angle first and to then calculate the target height according to the geometric relationship. Therefore, meter wave radar's low-elevation estimation problem represents one of the key difficulties in the field of radar measurement.
During direction of arrival (DOA) estimation with low-angle targets, there is a serious multipath coherent signal in low-elevation areas. Additionally, there is a correlation and even coherence between the signal and the direct wave signal, which leads to the rank deficiency of the covariance matrix of the received data, destroys the orthogonality between the signal subspace and the noise subspace, and greatly reduces the accuracy of DOA estimation. Therefore, the conventional eigen-subspace super-resolution algorithm cannot accurately estimate the DOA of the source when there is a coherent source in space, and it requires decoherence preprocessing. The spatial smoothing algorithm is a commonly tions of complex terrain, there is a broadband radar method based on a super-resolution algorithm [19]. Compared to conventional narrowband radar, this method has a better effect. The above decoherence algorithm can be applied to the meter wave radar height measurement of the polarization-sensitive array with appropriate deformation. For a direct height measurement model of a meter wave PSA, there are three main approaches: First, there is the polarization smoothing MUSIC algorithm [3] with the classical multipath signal model, which analyzes measurement performance. Second, there is a method that applies both the polarization smoothing algorithm and the spatial smoothing algorithm [20] for decoherence processing, solving the coherence between the direct wave and the reflected wave in low-elevation regions to a certain extent. Third, Tan proposed a method for elevation estimation using the generalized MUSIC algorithm after polarization smoothing [21]. It achieves better low-elevation angle estimation, but its computational complexity is still high, and the estimation accuracy and discrimination success probability need to be further improved. In addition, it cannot estimate the polarization parameters.
In order to further improve the height measurement accuracy and discrimination success probability of meter wave PSA radar in low-elevation areas, to reduce the computational complexity and to estimate the polarization parameters, the polarization synthesis steering vector MUSIC (P-SSV-MUSIC) algorithm is proposed in this paper. The proposed Sensors 2022, 22, 7298 3 of 20 algorithm synthesizes the steering vector formed by the reflected wave to the direct waveguide vector and does not need to solve the coherence problem. Then, the MUSIC algorithm is used to obtain the elevation results of the target in a low-elevation area. The simulation results indicate that the proposed algorithm can achieve higher height measurement accuracy and discrimination success probability, can estimate the polarization parameters, and has lower computational complexity and a lower resolution threshold.
The rest of the paper is organized as follows: The multipath signal model of the meter wave PSA is given in Section 2. The polarization smoothing generalized MUSIC (PS-GMUSIC) algorithm is introduced in Section 3. The modified polarization smoothing generalized MUSIC (MPS-GMUSIC) algorithm is described in Section 4. The P-SSV-MUSIC algorithm is discussed in Section 5. The calculations of three the algorithms are given in Section 6. The correctness and effectiveness of the algorithm according to computer simulation results are proved in Section 7. Finally, Section 8 provides our conclusions.

Multipath Signal Model of Meter Wave PSA
Assume that a meter wave PSA radar system uses the classical multipath receiving signal model. The low-elevation reflection area is a smooth and flat surface, f 0 denotes the signal frequency, λ denotes the wavelength, and the PSA is composed of N biorthogonal dipoles that are arranged along the Z-axis and that are parallel to the X-axis and the Z-axis. As shown in Figure 1, h a and h t are the height of the PSA antenna and the height of the target, respectively; R is the horizontal distance between the point perpendicular to the ground and the radar antenna; θ d is the incident angle of the target direct wave signal; θ s is the incident angle of the target reflection multipath signal; and d is the element spacing.
In order to further improve the height measurement accuracy and discrimination success probability of meter wave PSA radar in low-elevation areas, to reduce the computational complexity and to estimate the polarization parameters, the polarization synthesis steering vector MUSIC (P-SSV-MUSIC) algorithm is proposed in this paper. The proposed algorithm synthesizes the steering vector formed by the reflected wave to the direct waveguide vector and does not need to solve the coherence problem. Then, the MUSIC algorithm is used to obtain the elevation results of the target in a low-elevation area. The simulation results indicate that the proposed algorithm can achieve higher height measurement accuracy and discrimination success probability, can estimate the polarization parameters, and has lower computational complexity and a lower resolution threshold.
The rest of the paper is organized as follows: The multipath signal model of the meter wave PSA is given in Section 2. The polarization smoothing generalized MUSIC (PS-GMUSIC) algorithm is introduced in Section 3. The modified polarization smoothing generalized MUSIC (MPS-GMUSIC) algorithm is described in Section 4. The P-SSV-MUSIC algorithm is discussed in Section 5. The calculations of three the algorithms are given in Section 6. The correctness and effectiveness of the algorithm according to computer simulation results are proved in Section 7. Finally, Section 8 provides our conclusions.

Multipath Signal Model of Meter Wave PSA
Assume that a meter wave PSA radar system uses the classical multipath receiving signal model. The low-elevation reflection area is a smooth and flat surface, 0 f denotes the signal frequency, λ denotes the wavelength, and the PSA is composed of N biorthogonal dipoles that are arranged along the Z-axis and that are parallel to the X-axis and the Z-axis. As shown in Figure 1 Suppose that the target echo is a fully polarized electromagnetic wave and that ϕ is the target azimuth angle. We assumed that 2 ϕ π = , which means that the target is incident in the YOZ plane; [ ) , η π π ∈ − is the polarization phase difference; and is the polarization auxiliary angle. The polarization vector p in the Cartesian coordinate system of the meter wave PSA can be denoted as Suppose that the target echo is a fully polarized electromagnetic wave and that ϕ is the target azimuth angle. We assumed that ϕ = π/2, which means that the target is incident in the YOZ plane; η ∈ [−π, π) is the polarization phase difference; and γ ∈ [0, π/2) is the polarization auxiliary angle. The polarization vector p in the Cartesian coordinate system of the meter wave PSA can be denoted as where θ is the grazing angle of the direct signals and θ = θ d (Figure 1). The received signal for the PSA under multipath conditions can be expressed as where s is the echo signal vector after target scattering; n is the noise signal vector; α = 4πh a h t /Rλ is the phase difference generated by the delay difference between the reflected wave and the direct wave [22]; and a(θ d , η, γ) and a(θ s , η, γ) are the polarization spatial joint guidance vectors of the PSA corresponding to the direct wave and the ground-reflected wave, respectively, and are denoted as where ⊗ denotes the Kronecker product, and b(θ d ) and b(θ s ) are the steering vectors for the direct wave and the ground-reflected wave. These can be expressed as ρ is the Fresnel reflection coefficient's vector matrix. which contains the horizontal and vertical polarization reflection coefficients and can be denoted as where ρ h and ρ v are the Fresnel reflection coefficients for horizontal and vertical polarization, respectively. These can be defined as follows [23]: where ε is the complex permittivity of the smooth ground. It can be expressed by ε r and σ e as follows [23]: where ε r represents the relative permittivity, and σ e is the conductivity of the reflective surface. As we can see from Equation (6), the Fresnel reflection coefficients are determined by θ d and ε. Hence, we can obtain ρ h ≈ ρ v ≈ −1 when the grazing angle θ is very low. In other words, there is very little difference between ρ h and ρ v when the grazing angle is very low. We use vector C to define the synthetic steering vector of the PSA for simplicity, with C being denoted as where ψ is the fading coefficient vector of the received signal. It can be expressed as By substituting Equation (8) into Equation (2), the received signal in Equation (2) can be simplified as Additionally, the covariance of the received signal in Equation (10) is: where the superscript (·) H denotes conjugate transpose; I N is an N × N identity unit matrix; and σ 2 s = E ss H and σ 2 n = E nn H are the covariance of the signal and noise, respectively. We defined the signal-to-noise ratio of the received signal as SNR = σ 2 s /σ 2 n .

PS-GMUSIC Algorithm
The PS-GMUSIC algorithm combines polarization smoothing and the generalized MUSIC algorithm for PSA via the following procedure [21]: Firstly, according to the horizontal polarization component signal and the vertical polarization component signal received by the PSA [3], the received signals for the PSA on the coordinate's axis can be expressed as follows: where n h and n v are noise signals; According to Equation (11) and Equation (12), the covariance matrices of the multipath signals received via vertical polarization and horizontal polarization can be expressed as Second, the covariance matrix requires polarization smoothing. Polarization smoothing is an effective decoherence algorithm that can remove the correlation between the direct wave and the ground-reflected wave. The process of polarization smoothing [4] can be summarized as where Q ps is the signal envelope matrix of the covariance matrix R ps after polarization smoothing. This can be derived as follows: In Equation (15), R h and R v can be estimated as follows where T is the number of received signal snapshots. Third, we need singular value decomposition for covariance matrix B after polarization smoothing, which has the following form: where λ 1 > λ 2 ≥ λ 3 ∼ = λ 4 · · · λ N is the eigenvalue; U ps s is the eigenvector corresponding to two larger eigenvalues, representing the signal subspace; and U ps n is the eigenvector group corresponding to (N − 2) smaller eigenvalues, representing the noise subspace. Therefore, we can obtain the noise subspace projection matrix P, which is Next, we can obtain the spatial spectrum of the PS-GMUSIC algorithm as follows: where det[·] denotes the determinant of the matrix. As we can see from Equation (20), it requires a two-dimensional angle search, which increases the computational complexity of the algorithm. In order to reduce the computational complexity, we can use the prior information from the multipath model ( Figure 1): the geometric relationship [24] between the direct wave θ d and the reflected wave θ s : Substituting Equation (21) into Equation (20), the two-dimensional angle search in Equation (20) can be simplified to a one-dimensional angle search: Finally, we need to convert the estimation of the low-elevation angle to the height of the target after conducting the spectral peak search. The conversion formula is as follows: where R r denotes the distance between the target and the radar.

Proposed MPS-GMUSIC Algorithm
As we can see from Equation (16), the decoherence ability of polarization smoothing in low-elevation regions is related to the difference in the ground reflection coefficient between vertical polarization and parallel polarization. The rank of the signal envelope matrix Q ps after polarization smoothing is 2 when ψ h = ψ v , which indicates that the polarization smoothing algorithm can remove the coherence between the direct wave and the ground reflection wave. However, the rank of the signal envelope matrix Q ps after polarization smoothing is 1 when ψ h = ψ v . Additionally, the difference in the ground reflection coefficients between vertical polarization and parallel polarization is very small in low-elevation areas, and their values are all approximately −1. Therefore, the polarization smoothing decoherence ability is poor at this time, resulting in the meter wave PSA having poor angle measurement accuracy in low-elevation areas. In order to improve the decoherence performance of the polarization smoothing algorithm, we propose the MPS-GMUSIC algorithm, which is modified after polarization smoothing. The essence of the modified processing is forward and backward spatial smoothing processing with the number of subarrays being 1 [24]. The process of the MPS-GMUSIC algorithm is as follows: where the superscript (·) * denotes complex conjugation, and I v denotes the inverse diagonal unit matrix, which is: According to the definition of the composite steering vector A and I v , the following relation equations can be obtained: By substituting Equation (26) into Equation (24), R psm in Equation (24) can be simplified as: where Q psm is the signal envelope matrix of the modified polarization smoothing covariance matrix. It is defined as follows: where ψ h is the real part of the horizontal polarization multipath attenuation coefficient ψ h , and ψ v is the real part of the vertical polarization multipath attenuation coefficient ψ v . It can be seen from Equation (28) that when ψ h = ψ v , the rank of the signal envelope matrix Q psm after modified polarization smoothing is 2, which means that the modified polarization smoothing algorithm can also effectively remove the coherence between the direct wave and the ground reflection wave. Additionally, the modified polarization smoothing covariance matrix can be estimated as follows: whereR psm is the polarization smoothing covariance matrix. Then, we need to conduct singular value decomposition forR psm via the following form: is the eigenvector corresponding to two larger eigenvalues, representing the signal subspace; and U psm n is the eigenvector group corresponding to (N − 2) smaller eigenvalues, representing the noise subspace. Therefore, we can obtain the noise subspace projection matrix P psm , which is Next, we can obtain the spatial spectrum of the MPS-GMUSIC algorithm as follows: Substituting Equation (21) into Equation (32), we can obtain: As we can see from Equation (33), the two-dimensional spectral function becomes a one-dimensional spectral function, which greatly reduces the amount of calculation.
In the end, we need to use Equation (23) to convert the estimated low-elevation angle to the height of the target after the peak search. From the above two methods, it can be seen that the advantages of the PS-GMUSIC algorithm and the MPS-GMUSIC algorithm are that they are independent of the reflection coefficient and have strong robustness to the position. However, the disadvantages are that the accuracy is not high enough, and the polarization parameters cannot be estimated. Therefore, we propose the P-SSV-MUSIC algorithm, which has higher estimation accuracy and can estimate both the incident angle of the target direct wave signal θ d and the polarization parameters γ, η. It also has low algorithm complexity.

Proposed P-SSV-MUSIC Algorithm
As we can see from Equation (8), the signal synthesis steering vector C is onedimensional, and its rank is 1, expressing that the reflected wave polarization airspace joint steering vector is synthesized into the direct wave polarization airspace joint steering vector. Additionally, a rank of 1 means that the array-receiving signal model only has one incident signal source. Therefore, many conventional super-resolution DOA estimation algorithms can be directly applied to this signal model. It is necessary to demonstrate that the classical MUSIC algorithm has good directional resolution characteristics for incoherent signals as well as good measurement accuracy. In this paper, the conventional MUSIC algorithm is applied to the meter wave polarization-sensitive array elevation model, and dimension reduction is also carried out. The whole process can be derived as follows: First, we need to decompose the eigenvalues of Equation (10) and to divide its eigenvector. Then, the only eigenvector corresponding to a large eigenvalue is used to constitute the signal subspace E s , and its 2M × (2M − 1)-dimension eigenvector is used to constitute the noise subspace E n . After this, according to the conventional MUSIC algorithm, the spatial spectrum of the meter wave polarization-sensitive array can be obtained as follows: As we can see from Equation (34), it contains six unknowns and requires six-dimensional search processing, which is not suitable for practical engineering applications. Therefore, we need to reduce the number of dimensions. We divided the dimensionality reduction into two stages. In the first stage, we decoupled the polarization information and the DOA information to reduce the dimensionality. In the second stage, we used the relationship between the direct wave and reflected wave and the relationship between the reflection coefficient and the direct wave to reduce the number of dimensions.
In the first stage of dimension reduction, the polarization information and DOA information need to be decoupled to reduce the dimensionality, so the signal model in Section 2 needs to be deformed and classified. First, we rewrite Equation (1) in the following form: where E(θ) and g(γ, η) are defined as follows: Therefore, Equation (8) can be transformed into Next, we classify the reflection coefficient and the wave path difference into one class, the steering vector into another class, and the polarization information into the same category. After this, Equation (38) can be further deformed to obtain Since the reflection coefficient ρ h , ρ v and range difference α are functions of θ d , θ s , we can obtain the following definition: According to Equations (40)-(42), Equation (39) can be simplified as follows: Then, according to Equation (43), we can define the MUSIC cost function as: It is not difficult to find g(γ, η) H g(γ, η) = 1. Therefore, Equation (44) can be further transformed into Through the observation of Equation (45), it can be found that Equation (45) conforms to the criteria for the maximum and minimum Rayleigh quotients [25]. Therefore, we can obtain λ min (θ d , θ s ) = min where λ min (θ d , θ s ) denotes the minimum eigenvalue obtained by the eigendecomposition . Therefore, we can obtain the fourdimensional search equation, which is The eigenvector corresponding to the minimum eigenvalue is the estimated value of the polarization parameterĝ(γ, η), so the polarization parameter can also be obtained by the following formula: As we can see from Equation (47), the polarization information and DOA information have been decoupled, which reduces the two search dimensions. However, the number of calculations required for a four-dimensional search element is still large. The second stage of dimension reduction is carried out below: First, by substituting Equation (21) into Equation (47), Equation (47) can be reduced to a three-dimensional angle search.
Secondly, as we can see from Equation (6), the reflection coefficient ρ h , ρ v is determined by the incident angle θ d , the relative dielectric constant ε r , and the surface material conductivity σ e . Additionally, the relative dielectric constant ε r and the surface material conductivity σ e are different and are in different positions, such as in water, land, and vegetation. Fortunately, previous scholars have measured and summarized specific values under different scenarios [23], as shown in Table 1. According to Table 1, the relative dielectric constant ε r and the surface material conductivity σ e are known. Therefore, there is only one unknown in Equation (47), reducing it to a one-dimensional search, and the MUSIC dimension reduction process is completed. Finally, we use Equation (23) to convert the estimated low-elevation angle of the target after the spectral peak search into the height H of the target.
It should be noted that relative dielectric constant ε r and surface material conductivity σ e cannot be accurately measured in a complex position scene, and there are certain errors. In this situation, a three-dimensional search process is needed, and an alternative search method [26] can be used to alleviate the above problems, but this is not the core content of this paper, so it is not expanded upon here.

Computational Complexity
In this section, we analyze the computational complexity of the three algorithms and compare them with each other. The computational complexity of the different algorithms is shown in Table 2, where N represents the number of antenna elements, T represents the number of snapshots of the received signals, and n represents the number of DOA angle searches. As we can see from Table 2, the computational complexity of the MPS-GMUSIC algorithm is only 2N 2 more than that of the PS-GMUSIC algorithm, and the computational complexity of the P-SSV-MUSIC algorithm is at least nN 2 less than the other two algorithms. Table 2. Computational complexity of three algorithms.

Algorithm Name
Computational Complexity

Simulation
As described in this section, some simulations we executed in terms of spectrum estimation, computational complexity, RSME, discrimination success probability and the tracking measurement results of the simulated track to evaluate the proposed algorithm. We adopted the meter wave PSA radar system with N = 13, d = λ/2, f 0 = 150 MHz, λ = 2 m, γ = 85 • , η = 170 • , h a = 10 m, and R r = 100 km. Additionally, we assumed that the reflection dielectric is seawater. As shown in Table 1, the surface material conductivity and relative permittivity can be set as σ e = 0.5 and ε r = 75, respectively. Figure 2 demonstrates the spectrum estimation of the three methods. As we can see in Figure 2, the proposed MPS-GMUSIC and P-SSV-MUSIC algorithms can estimate the low-elevation angle of the target, and the overall effect is better than that of the PS-GMUSIC algorithm. However, the elevation estimation results of the MPS-GMUSIC algorithm are only closer to the actual low-elevation angle value; the estimated angle of the P-SSV-MUSIC algorithm is the same as the actual low-elevation angle value, and its spectral peak is sharper, indicating that the method has a better angular resolution, making it more accurate in practical applications.

Computational Complexity of the Proposed Method
In order to compare the computational complexity of the three algorithms mo tuitively, we set the number of antenna elements between 2 and 30. In addition, sinc search interval set by Simulation 7.1 is 0.01  , the number of DOA angle search n 1000 = . Figure 3 shows the computational complexity of the three algorithms wi spect to number of antenna elements. It can be seen that the computational complex the P-SV-MUSIC algorithm is significantly lower than that of the other two algorithm

Computational Complexity of the Proposed Method
In order to compare the computational complexity of the three algorithms more intuitively, we set the number of antenna elements between 2 and 30. In addition, since the search interval set by Simulation 7.1 is 0.01 • , the number of DOA angle searches is n = 1000. Figure 3 shows the computational complexity of the three algorithms with respect to number of antenna elements. It can be seen that the computational complexity of the P-SV-MUSIC algorithm is significantly lower than that of the other two algorithms.
In order to compare the computational complexity of the three algorithms more intuitively, we set the number of antenna elements between 2 and 30. In addition, since the search interval set by Simulation 7.1 is 0.01  , the number of DOA angle searches is n 1000 = . Figure 3 shows the computational complexity of the three algorithms with respect to number of antenna elements. It can be seen that the computational complexity of the P-SV-MUSIC algorithm is significantly lower than that of the other two algorithms.

RMSE of the Proposed Method
In order to assess the angle estimation performance of the proposed algorithm, this section presents the RMSE simulations.

RMSE with Respect to the SNR
We assumed that the SNR between 10 dB − and10 dB . Figure 4 shows the RMSE of the elevation with respect to the SNR when d θ is 1  and 2  . As we can see in Figure 4, the estimation accuracy of the MPS-GMUSIC algorithm and P-SSV-MUSIC algorithm is improved when there is an increase in the SNR, which is consistent with the expected results. The accuracy of the MPS-GMUSIC algorithm cannot exceed that of the PS-GMUSIC algorithm when d θ is 1  . However, the P-SSV-MUSIC algorithm is superior to the other two algorithms in terms of the angle measurement accuracy and height measurement accuracy within the given SNR range, regardless of whether d θ is 1  or 2  ,

RMSE of the Proposed Method
In order to assess the angle estimation performance of the proposed algorithm, this section presents the RMSE simulations.

RMSE with Respect to the SNR
We assumed that the SNR between −10 dB and 10 dB. Figure 4 shows the RMSE of the elevation with respect to the SNR when θ d is 1 • and 2 • . As we can see in Figure 4, the estimation accuracy of the MPS-GMUSIC algorithm and P-SSV-MUSIC algorithm is improved when there is an increase in the SNR, which is consistent with the expected results. The accuracy of the MPS-GMUSIC algorithm cannot exceed that of the PS-GMUSIC algorithm when θ d is 1 • . However, the P-SSV-MUSIC algorithm is superior to the other two algorithms in terms of the angle measurement accuracy and height measurement accuracy within the given SNR range, regardless of whether θ d is 1 • or 2 • , and it is an order of magnitude higher. This proves that the P-SSV-MUSIC algorithm has excellent angle measurement accuracy.   Figure 5 shows the RMSE of the elevation with respect to the SNR of the two-polarization parameters and its CRB using the P-SSV-MUSIC algorithm, which is indicates that the algorithm accurately estimates the two-dimensional polarization parameters and can be close to the optimal estimation performance. The CRB corresponding to the direct wave angle and the two-polarization parameters and the deduction process are included in Appendix A.   Figure 5 shows the RMSE of the elevation with respect to the SNR of the twopolarization parameters and its CRB using the P-SSV-MUSIC algorithm, which is indicates that the algorithm accurately estimates the two-dimensional polarization parameters and can be close to the optimal estimation performance. The CRB corresponding to the direct wave angle and the two-polarization parameters and the deduction process are included in Appendix A. Figure 5 shows the RMSE of the elevation with respect to the SNR of the two-polarization parameters and its CRB using the P-SSV-MUSIC algorithm, which is indicates that the algorithm accurately estimates the two-dimensional polarization parameters and can be close to the optimal estimation performance. The CRB corresponding to the direct wave angle and the two-polarization parameters and the deduction process are included in Appendix A.

RMSE with Respect to the Snapshot Number
We assumed that the snapshot number is between 10 and 250 . Figure 6 shows the RMSE of the elevation with respect to the snapshot number when d θ is 1  and 2  . As we can see in Figure 6, the estimation accuracy of the MPS-GMUSIC algorithm and P-SSV-MUSIC algorithm improves when there is an increase in the snapshot number, which is consistent with the expected results. The accuracy of the MPS-GMUSIC algorithm cannot exceed that of the PS-GMUSIC algorithm when d θ is 1  . However, the P-SSV-MUSIC algorithm is superior to the other two algorithms in terms of angle measurement accuracy and height measurement accuracy within the given snapshot number range, regardless of whether d θ is 1  or 2  , and it is an order of magnitude higher.

RMSE with Respect to the Snapshot Number
We assumed that the snapshot number is between 10 and 250 . Figure 6 shows the RMSE of the elevation with respect to the snapshot number when θ d is 1 • and 2 • . As we can see in Figure 6, the estimation accuracy of the MPS-GMUSIC algorithm and P-SSV-MUSIC algorithm improves when there is an increase in the snapshot number, which is consistent with the expected results. The accuracy of the MPS-GMUSIC algorithm cannot exceed that of the PS-GMUSIC algorithm when θ d is 1 • . However, the P-SSV-MUSIC algorithm is superior to the other two algorithms in terms of angle measurement accuracy and height measurement accuracy within the given snapshot number range, regardless of whether θ d is 1 • or 2 • , and it is an order of magnitude higher.

RMSE with Respect to the Phase Aberration
We assumed that the phase aberration is between 0  and 45  . Figure 7 shows the RMSE of elevation with respect to the phase aberration when d θ is 1  and 2  . As we can see in Figure 7, as the phase aberration increases, the angle RSME of the three algorithms increases, and the performance of the algorithm decreases. However, the RSME of the P-SV-MUSIC algorithm is lower than that of the PS-GMUSIC algorithm and the MPS-GMUSIC algorithm, regardless of whether d θ is 1  or 2  , which proves that it has higher measurement accuracy.

RMSE with Respect to the Phase Aberration
We assumed that the phase aberration is between 0 • and 45 • . Figure 7 shows the RMSE of elevation with respect to the phase aberration when θ d is 1 • and 2 • . As we can see in Figure 7, as the phase aberration increases, the angle RSME of the three algorithms increases, and the performance of the algorithm decreases. However, the RSME of the P-SV-MUSIC algorithm is lower than that of the PS-GMUSIC algorithm and the MPS-GMUSIC algorithm, regardless of whether θ d is 1 • or 2 • , which proves that it has higher measurement accuracy.
We assumed that the phase aberration is between 0  and 45  . Figure 7 shows the RMSE of elevation with respect to the phase aberration when d θ is 1  and 2  . As we can see in Figure 7, as the phase aberration increases, the angle RSME of the three algorithms increases, and the performance of the algorithm decreases. However, the RSME of the P-SV-MUSIC algorithm is lower than that of the PS-GMUSIC algorithm and the MPS-GMUSIC algorithm, regardless of whether d θ is 1  or 2  , which proves that it has higher measurement accuracy.

Discrimination Success Probability of the Proposed Method
In this simulation, we set the search angle to between 0.1  and 4  . In addition, we consider that the discrimination is successful when the spectrum estimation produces obvious peaks and when the estimated angle is less than or equal to 0.05  . Figure 8 shows the discrimination success probability of the three algorithms with respect to the search angle when the SNR is 15 dB and 20 dB . The simulation results show that the P-SSV-MUSIC algorithm has a higher discrimination success probability and a lower discrimination threshold.

Discrimination Success Probability of the Proposed Method
In this simulation, we set the search angle to between 0.1 • and 4 • . In addition, we consider that the discrimination is successful when the spectrum estimation produces obvious peaks and when the estimated angle is less than or equal to 0.05 • . Figure 8 shows the discrimination success probability of the three algorithms with respect to the search angle when the SNR is 15 dB and 20 dB. The simulation results show that the P-SSV-MUSIC algorithm has a higher discrimination success probability and a lower discrimination threshold.

Tracking Measurement Results of Simulated Track
In this simulation, we set the target's range to be between 50 km and 100 km. The elevation measurement results of the target for the three algorithms are shown in Figure  9. The height measurement results of target for the three algorithms are shown in Figure  10. Compared to the real elevation angle and the real height of the target, we found that the angle and height estimation results of the P-SSV-MUSIC algorithm are almost consistent with the actual value. The angle and height estimation results of the PS-GMUSIC algorithm and MPS-GMUSIC algorithm begin to deviate from the actual values after 85 km, but the MPS-GMUSIC algorithm demonstrates less deviation than the PS-GMUSIC algorithm. Figure 11 shows the error results of the target elevation measurements for the three algorithms. We can see intuitively in Figure 11 that the error results of the target elevation measurements of the three algorithms are within a small range when the horizontal distance lies before 85 but that large errors begin to occur after 85. It increases as the distance of the target increases, which is consistent with the expected results. Additionally, the error results of the target elevation measurements of the P-SSV-MUSIC algorithm are

Tracking Measurement Results of Simulated Track
In this simulation, we set the target's range to be between 50 km and 100 km. The elevation measurement results of the target for the three algorithms are shown in Figure 9. The height measurement results of target for the three algorithms are shown in Figure 10. Compared to the real elevation angle and the real height of the target, we found that the angle and height estimation results of the P-SSV-MUSIC algorithm are almost consistent with the actual value. The angle and height estimation results of the PS-GMUSIC algorithm and MPS-GMUSIC algorithm begin to deviate from the actual values after 85 km, but the MPS-GMUSIC algorithm demonstrates less deviation than the PS-GMUSIC algorithm. Figure 11 shows the error results of the target elevation measurements for the three algorithms. We can see intuitively in Figure 11 that the error results of the target elevation measurements of the three algorithms are within a small range when the horizontal distance lies before 85 but that large errors begin to occur after 85. It increases as the distance of the target increases, which is consistent with the expected results. Additionally, the error results of the target elevation measurements of the P-SSV-MUSIC algorithm are within the range of ±0.01 • when the horizontal distance is 85-100. The error results of the target elevation measurements of the MPS-GMUSIC algorithm are within the range of ±0.02 • when the horizontal distance is 85-100. The error results of the target elevation measurements of the PS-GMUSIC algorithm are within the range of ±0.12 • when the horizontal distance is 85-100. Figure 12 shows the error results of the target height measurements for the three algorithms. It is clear that the error results of the target height measurements for the P-SSV-MUSIC algorithm are within 25 m and that those of the MPS-GMUSIC algorithm are within 50 m. However, the maximum error of the PS-GMUSIC algorithm is close to 200 m. These indicate that the P-SSV-MUSIC algorithm has better angle and height estimation performance and that it can realize precise height measurements of meter wave radar.