Polarization Smoothing Generalized MUSIC Algorithm with Polarization Sensitive Array for Low Angle Estimation

Direction of Arrival (DOA) estimation of low-altitude targets is difficult due to the multipath coherent interference from the ground reflection image of the targets, especially for very high frequency (VHF) radars, which have antennae that are severely restricted in terms of aperture and height. The polarization smoothing generalized multiple signal classification (MUSIC) algorithm, which combines polarization smoothing and generalized MUSIC algorithm for polarization sensitive arrays (PSAs), was proposed to solve this problem in this paper. Firstly, the polarization smoothing pre-processing was exploited to eliminate the coherence between the direct and the specular signals. Secondly, we constructed the generalized MUSIC algorithm for low angle estimation. Finally, based on the geometry information of the symmetry multipath model, the proposed algorithm was introduced to convert the two-dimensional searching into one-dimensional searching, thus reducing the computational burden. Numerical results were provided to verify the effectiveness of the proposed method, showing that the proposed algorithm has significantly improved angle estimation performance in the low-angle area compared with the available methods, especially when the grazing angle is near zero.


Introduction
Meter wave radar is a direct and effective anti-stealth measure because the meter wavelength is located in the resonance band for most types of aircrafts, with the shape stealth technology being no longer valid in this band. Moreover, the wave absorbing efficiency of the dielectric coating is remarkably corrupted in this band [1]. However, Direction of Arrival (DOA) estimation in low-angle areas is still one of the most challenging problems, especially for the Very High Frequency (VHF) radar system, which has antennae that are severely restricted in terms of aperture and height. The poor performance of the angle estimation and the height measuring, especially in the low-angle tracking area, limits the application of meter-wave radars. The cause of this problem is the existence of the multipath effect. In low-angle estimation for the VHF radar, the direct signal and multipath signal lying in the main lobe are coherent and therefore, it is difficult to distinguish them in the spatial, temporal and Doppler domains.
Many efforts have been exerted in the recent decades to overcome the multipath effect. These efforts can be mainly divided into two broad categories: the extended monopulse methods [2][3][4][5][6][7] and the array signal processing techniques [8][9][10][11][12][13][14][15][16][17][18]. The former includes the complex indicated angle method, the double-null method and the fixed-beam method [4,5]. Although these methods have the horizontal and vertical polarized components can fix that problem by using the different phases of the fading coefficient between the horizontal and vertical polarization. In order to improve the low-angle estimation performance, especially when the grazing angle is close to zero, the Polarization Smoothing Generalized MUSIC (PS-GMUSIC) algorithm is proposed in this paper. The proposed algorithm exploits the composite searching vector instead of ordinary searching vector, making it unnecessary to solve the coherent problem. Furthermore, the proposed method based on the PSA does not have limitations in terms of the phase differences between the direct and the reflected signals. In addition, the geometry information of the symmetry multipath model is exploited to transform the two-dimensional angle search into one-dimensional angle search in order to reduce the computational burden for the proposed algorithm. The simulation results indicated that the proposed method has better low-angle estimation performance than the PS-MUSIC and PS-SS-MUSIC algorithms.
The rest of the paper is organized as follows: the multipath signal model of the polarization sensitive array and the derivation of CRB are described in Section 2. The polarization smoothing Generalized MUSIC algorithm (PS-GMUSIC) is proposed in Section 3. The simulation results and discussions are presented in Section 4. Finally, Section 5 provides our conclusions.

Multipath Model for Polarization Sensitive Array
The multipath echoes are comprised of a single specular reflection and diffuse reflections. For a rough surface, the diffuse component dominates, while the reflected signals consist mainly of the specular component for a smooth surface. The ground can be considered as a smooth surface at the VHF band, because the Rayleigh roughness criterion can be easily satisfied due to the long wavelength based on the Rayleigh roughness criterion. Thus, the diffuse reflection components are ignored and the 'double-ray' symmetry multipath model is used, which is shown in Figure 1. angle estimation performance, especially when the grazing angle is close to zero, the Polarization Smoothing Generalized MUSIC (PS-GMUSIC) algorithm is proposed in this paper. The proposed algorithm exploits the composite searching vector instead of ordinary searching vector, making it unnecessary to solve the coherent problem. Furthermore, the proposed method based on the PSA does not have limitations in terms of the phase differences between the direct and the reflected signals. In addition, the geometry information of the symmetry multipath model is exploited to transform the twodimensional angle search into one-dimensional angle search in order to reduce the computational burden for the proposed algorithm. The simulation results indicated that the proposed method has better low-angle estimation performance than the PS-MUSIC and PS-SS-MUSIC algorithms. The rest of the paper is organized as follows: the multipath signal model of the polarization sensitive array and the derivation of CRB are described in Section 2. The polarization smoothing Generalized MUSIC algorithm (PS-GMUSIC) is proposed in Section 3. The simulation results and discussions are presented in Section 4. Finally, Section V provides our conclusions.

Multipath Model for Polarization Sensitive Array
The multipath echoes are comprised of a single specular reflection and diffuse reflections. For a rough surface, the diffuse component dominates, while the reflected signals consist mainly of the specular component for a smooth surface. The ground can be considered as a smooth surface at the VHF band, because the Rayleigh roughness criterion can be easily satisfied due to the long wavelength based on the Rayleigh roughness criterion. Thus, the diffuse reflection components are ignored and the 'double-ray' symmetry multipath model is used, which is shown in Figure 1. As shown in Figure 1, the N Polarization Sensitive Array (PSA) elements are uniformly distributed on the Z-axis, while the orthogonal electric dipoles are arranged along the X-and Zdirections. The distance between elements is equal to a half-wavelength. In this figure, is the height of the PSA; is the elevation angle of the direct signal of the target; is the depression angle of the specular signal; and is the grazing angle of the target. According to the symmetry geometrical information of the symmetry multipath model (Figure 1), we can obtain that: Supposing the echo is a completely polarized wave as shown in the Figure 1, the direct electric field vector can be expressed on the polarization basis ( , ) as: where is the polarization angle; is the phase difference of the different polarized components; and are the unit vectors in the Spherical Coordinates with = [− , , 0] and = [ , , − ] ; and = and = are the electric field amplitude of the unit vectors and , respectively. Obviously, although there are three components of the completely polarized electric field, only two of them are independent. Hence, we assumed that = As shown in Figure 1, the N Polarization Sensitive Array (PSA) elements are uniformly distributed on the Z-axis, while the orthogonal electric dipoles are arranged along the Xand Z-directions. The distance between elements is equal to a half-wavelength. In this figure, H r is the height of the PSA; θ d is the elevation angle of the direct signal of the target; θ s is the depression angle of the specular signal; and ψ is the grazing angle of the target. According to the symmetry geometrical information of the symmetry multipath model (Figure 1), we can obtain that: Supposing the echo is a completely polarized wave as shown in the Figure 1, the direct electric field vector can be expressed on the polarization basis (θ, ϕ) as: where γ is the polarization angle; η is the phase difference of the different polarized components; ϕ and θ are the unit vectors in the Spherical Coordinates with e ϕ = [−sinϕ, cosϕ, 0] T and e θ = [sinθcosϕ, sinθsinϕ, −cosθ] T ; and E ϕ = cosγ and E θ = sinγe jη are the electric field amplitude of the unit vectors e ϕ and e θ , respectively. Obviously, although there are three components of the completely polarized electric field, only two of them are independent. Hence, we assumed that ϕ = π/2 in this paper as shown in Figure 1. The unit vectors of e ϕ and e θ can be simplified as e ϕ = [−1, 0, 0] T and e θ = [0, sinθ, −cosθ] T . The transformation between unit vectors in the spherical coordinates system and the Cartesian coordinate system is defined as: The polarization vector p in Cartesian coordinate system of the 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 condition can be expressed as: where s(t) is the narrow band signal reflected from the target, with the assumption that n(t) is a 2N × 1 Gaussian white noise vector with zero mean. Furthermore, ∆R = 2H r sinθ d is the path difference between the direct and the specular signals; and b(θ d , η, γ), b(θ s , η, γ) are the steering vector of the PSA for the direct and the specular signals, respectively. That can be denoted as: where denote the Kronecker product; and a(θ d ) and a(θ s ) are the steering vectors for the direct and the specular signals. These can be expressed as: The Γ(ψ) in Equation (5) is the Fresnel reflection coefficients vector that contains the horizontal and vertical polarization reflection coefficients, which can be denoted as: For the horizontal polarization, the electric field is parallel to the interface. For the vertical polarization, the electric field is perpendicular to the interface as defined in a previous study [28]. In this paper, the interface plane is XOY plane as shown in Figure 1. Hence, the electric field in θ direction is called the vertical polarization and the electric field in ϕ direction is called horizontal polarization. The Fresnel reflection coefficients for horizontal and vertical polarization are defined respectively as follows [28]: The Fresnel reflection coefficients are determined by the grazing angle ψ and ε c . where ε c = ε r − j60λσ is the complex permittivity of the smooth ground; ε r represents the relative permittivity; and σ is the conductivity of the reflective surface. When the grazing angle ψ is close to zero, we can obtain Γ h (ψ) ≈ Γ h (ψ) ≈ −1. In other words, there is very little difference between Γ h (ψ) and Γ v (ψ) when the grazing angle is close to zero. In order to simplify the received signal in Equation (4), we used c to represent the composite vectors, with c denoted as: where ρ is the fading coefficient vector of the PSA. It can be expressed as: where ρ h and ρ v are the horizontal and vertical polarization fading coefficients, respectively. These are expressed as: where ρ h and ρ v are the real parts of horizontal and vertical polarization fading coefficients, respectively; and ρ h and ρ v are the image parts of horizontal and vertical polarization fading coefficients, respectively. The fading coefficients are dependent on the grazing angle ψ, Fresnel reflection coefficient Γ h,v (ψ) and path difference ∆R. By substituting Equation (10) into Equation (5), the received signals in Equation (5) can be simplified as: In order to derive the CRB of the elevation angle θ d , we need to calculate the covariance of the received signal. The covariance of the received signal in Equation (13) is: where σ 2 s = E s(t)s(t) H and σ 2 n = E n(t)n(t) H are the covariance of the signal and noise, respectively. We defined the signal-to-noise ratio of the completely polarized wave signal as: SNR = σ 2 s /σ 2 n . The CRB is a universal lower boundary for the variance of all unbiased estimators of a set of parameters. It is defined as the inverse of the Fisher information matrix (FIM), which describes the amount of information that the data provide about unknown parameters. According to the symmetry geometrical information of the symmetry multipath model as shown in Figure 1, we can obtain θ s = −θ d . Hence, there are nine unknown parameters s , σ 2 n in the covariance matrix R. The Cramer-Rao Bound (CRB) of the elevation angle θ d for the PSA can be derived from FIM and the deduction process is included in the Appendix A. The CRB for DOA of grazing angle θ d is:

Polarization Smoothing Generalized MUSIC Algorithm
As shown in Figure 1, the received signals for the PSA under the multipath conditions for the coordinate's axis can be expressed as follows [17]: where A = [a(θ d ), a(θ s )] is the composite steering vector that contains a(θ d ) and a(θ s ). In addition, the X-direction component corresponds to the horizontal polarization and the Z-direction component corresponds to the vertical polarization. The covariance of each polarized components of the received signal in Equation (16) can be expressed as: As we can see from Equation (17), there is only one larger eigenvalue in the covariance matrix R xx and R zz due to the coherence between the direct and the specular signals. The Polarization Smoothing algorithm was proposed by Rahamim [16], which averages the data covariance matrix along the elements in the vector sensors of an array, with the averaging operation performed along the array aperture. The polarization smoothing pre-processing for VHF Polarization Sensitive Array under the multipath conditions can be summarized as: where Q is the averaged envelope covariance matrix of horizontal and vertical polarization subarrays. This can be derived as follows: As we can see from Equation (19), we can obtain that rank(Q) = 2 when ρ h = ρ v . This indicates that the polarization smoothing pre-processing could eliminate the coherence between the direct and the specular signals. Hence, there are two larger eigenvalues in the polarization smoothing average covariance R ps . The decorrelation performance of the polarization smoothing pre-processing is based on the differences of the ρ h and ρ v . The PS-MUSIC [17] algorithm is established based on the polarization smoothing technique and MUSIC algorithm for low-angle estimation problem. However, when the grazing angle is close to zero, there is very little difference between ρ h and ρ v . Therefore, the decorrelation performance of the polarization smoothing method is reduced when the grazing angle is close to zero. In order to improve the decorrelation performance, the spatial smoothing processing can be exploited after the polarization smoothing pre-processing, which can be used to achieve the PS-SS-MUSIC algorithm. The PS-SS-MUSIC algorithm suffers aperture loss when the number of overlapping subarrays for the spatial smoothing is more than one. In contrast, GMUSIC [27] can be used to solve the low-angle estimation problem, which does not need to solve the coherent problem. Hence, we developed a GMUSIC algorithm for the PSA.
As we can see from the multipath signal model in Equation (16), the covariance matrices R h and R v can be estimated with L snapshots byR h = 1 Equation (18). The steering vectors a(θ d ) and a(θ s ) that correspond to the direct and the specular incoming signals in the signal subspace are therefore orthogonal to the noise subspace U n . According to the orthogonality of the signal subspace to the noise subspace, we can derive: where P is the projection matrix that is constructed from the eigenvectors U n of the sample covariance matrixR ps corresponding to the N-2 minimum eigenvalue. By constructing the projection matrix A(θ 1 , θ 2 ) H PA(θ 1 , θ 2 ), we obtain: where θ 1 and θ 2 are the searching angles; and A = [a(θ 1 ), a(θ 2 )] is the composite steering vector. Obviously, the A(θ 1 , θ 2 ) H PA(θ 1 , θ 2 ) is singular if and only if {θ 1 = θ 2 }, {θ 1 = θ d and θ 2 = θ s }, or {θ 2 = θ d and θ 1 = θ s }. According to the geometry information of the symmetry multipath model, we assumed that θ 1 > 0 • and θ 2 < 0 • for the range of interest. Therefore, A(θ 1 , θ 2 ) H PA(θ 1 , θ 2 ) is singular only if θ 1 = θ d , θ 2 = θ s , which can obtain that: where det() denotes the determinant of the matrix. Moreover, when ρ h is equal to ρ v , A(θ d , θ s ) H PA(θ d , θ s ) also is singular, with the proof included in the Appendix B. After this, the spatial spectrum of PS-GMUSIC is formed by projecting a continuum of composite DOA vector A(θ 1 , θ 2 ) over the given range of interest onto P, with the following form: The proposed method exploits the composite steering vector A(θ 1 , θ 2 ), making it unnecessary to solve the coherent problem. However, it requires a multidimensional search, which leads to a significant increase in computational burden. Additionally, based on the geometry of symmetry multipath model (Figure 1), it can obtain the following expression for the searching angles: Substituting Equation (24) into Equation (23), the two-dimensional angle searching in Equation (23) can be transformed into one-dimensional angle searching, which is: Up to now, we have achieved the PS-GMUSIC algorithm. The procedure of the proposed method can be summarized as follows: 1.
Use the polarization smoothing pre-processing to getR ps ; 3.
Compute the noise subspace U n and Projection matrix P by means of EVD (Eigenvalue Decomposition) usingR ps ; and
The geometry information of the multipath model can be exploited to reduce the computational cost for the PS-GMUSIC algorithm.
After the transformation from the two-dimensional searching into the one-dimensional searching, the computation burden of the PS-GMUSIC is O PS−GMUSIC N 3 + 2LN 2 + n N 3 + 6N 2 + 20N + 4 , where L is the number of the snapshots and n is the number of the searching angles.

Results
In this section, some simulation results are presented to assess the performance of the proposed method in terms of CRB, RMSE and Pseudo spectra. We adopted the VHF PSA radar system with N = 10, f 0 = 300 MHZ, H r = 10 m, SNR = 10 dB, d = λ/2 and L = 256. We assumed that there is one point target over the calm seawater (ε r = 80, σ = 4 s/m) for the following simulations.

CRB Analysis
The CRB of the PSA VHF radar system for different grazing angles and polarization angles is shown in Figure 2. The polarization angle γ lies in (0 • -90 • ). The grazing angle of the target varies from 0 • to 10 • . The CRB is very high when the grazing angle is close to zero. That is still a challenging problem for low-angle estimation. There are some periodic peaks and troughs for the CRB when the polarization angle γ lies in (0 • -20 • ) and (70 • -90 • ), which is shown in Figure 2  As we can see from Figure 3, the CRBs of horizontal and vertical polarization have 'singular points' (Peak points) when the phase of fading coefficients is ±180°. For example, the grazing angle is 2.9°, 3.5°, 5.7° and 6.8°, which is shown in Figure 3a,b. It is because the SNR is reduced when the phase of fading coefficients is ±180° due to the cancellation between the direct and reflected signals. This is why the performance of GMUSIC algorithm corrupts when the phase of fading coefficients is ±180° with the traditional phase array. In contrast, the CRB of PSA could merge the 'singular points' for γ = 45° based on the different phases of fading coefficients for the horizontal and vertical polarization components in the PSA. As we can see from Figure 3, the CRBs of horizontal and vertical polarization have 'singular points' (Peak points) when the phase of fading coefficients is ±180 • . For example, the grazing angle is 2.9 • , 3.5 • , 5.7 • and 6.8 • , which is shown in Figure 3a,b. It is because the SNR is reduced when the phase of fading coefficients is ±180 • due to the cancellation between the direct and reflected signals. This is why the performance of GMUSIC algorithm corrupts when the phase of fading coefficients is ±180 • with the traditional phase array. In contrast, the CRB of PSA could merge the 'singular points' for γ = 45 • based on the different phases of fading coefficients for the horizontal and vertical polarization components in the PSA. is 2.9°, 3.5°, 5.7° and 6.8°, which is shown in Figure 3a,b. It is because the SNR is reduced when the phase of fading coefficients is ±180° due to the cancellation between the direct and reflected signals. This is why the performance of GMUSIC algorithm corrupts when the phase of fading coefficients is ±180° with the traditional phase array. In contrast, the CRB of PSA could merge the 'singular points' for γ = 45° based on the different phases of fading coefficients for the horizontal and vertical polarization components in the PSA.

RMSE of the Proposed Method
In order to assess the angle estimation performance of the proposed algorithm, the Monte Carlo simulations are presented in this part. The root-mean-square error (RMSE) is defined below: where ̂, is the i ℎ angle estimation of the DOA of grazing angle , M is the number of Monte Carlo simulations. Assume M = 1000 for the following simulations. The RMSE of the proposed method was compared with PS-MUSIC [17], PS-SS-MUSIC [16] and IBMUSIC [29] algorithms. The Spatial Smoothing utilizes the forward and backward smoothing [23] with two overlapping subarrays in the simulation. The iteration beamspace multiple signal classification (IBMUSIC) algorithm was proposed in [29] to realize unbiased estimation and reduce the calculation time for DOA of CDMA Multipath signals. This approach can be utilized for low angle estimation performance with PSA. The performance of the IBMUSIC algorithm heavily depends on initial value [29]. And assume the initial value is 1° for the following simulations. The beam number of beamspace is 5 for IBMUSIC algorithm and assume the deviation precision value as 0.001 for iteration. Figures 4a,b show the RMSE with respect to the polarization angle γ. We assumed that the direct signal arrived at an angle of 2.9° and 3.5° when the phase of fading coefficients is ±180° (Figure 3).

RMSE of the Proposed Method
In order to assess the angle estimation performance of the proposed algorithm, the Monte Carlo simulations are presented in this part. The root-mean-square error (RMSE) is defined below:  [16] and IBMUSIC [29] algorithms. The Spatial Smoothing utilizes the forward and backward smoothing [23] with two overlapping subarrays in the simulation. The iteration beamspace multiple signal classification (IBMUSIC) algorithm was proposed in [29] to realize unbiased estimation and reduce the calculation time for DOA of CDMA Multipath signals. This approach can be utilized for low angle estimation performance with PSA. The performance of the IBMUSIC algorithm heavily depends on initial value [29]. And assume the initial value is 1 • for the following simulations. The beam number of beamspace is 5 for IBMUSIC algorithm and assume the deviation precision value as 0.001 for iteration. Figure 4a,b show the RMSE with respect to the polarization angle γ. We assumed that the direct signal arrived at an angle of 2.9 • and 3.5 • when the phase of fading coefficients is ±180 • (Figure 3). Since the decorrelation performance of the polarization smoothing are determined by the polarization angle [17], the RMSE of PS-MUSIC varies according to the polarization angle as shown in  Figure 4b, the RMSE of PS-SS-MUSIC is higher than those of the PS-MUSIC algorithm due to the aperture loss. The performance of IBMUSIC algorithm is poor when γ lies in (70 • -90 • ) and (0 • -20 • ) as shown in Figure 4a. It because the decorrelation of the polarization smoothing is poor when the single polarized components dominate. Unlike the other schemes, when the single polarized components dominate the proposed method does not need to solve the coherence between the direct signal and specular signals. Hence, our proposed scheme retains more stable performance for all the polarization angles as shown in Figure 4. The RMSE of the proposed method is much higher when γ lies in (0 • -20 • ) in Figure 4a or (80 • -90 • ) in Figure 4b. This is because the SNR is reduced when the phase of fading coefficients is ±180 • for single polarized dominates. To achieve more stable angle estimation performance for all DOA of grazing angle, γ = 45 • is assumed in the following simulations.
due to the aperture loss. The performance of IBMUSIC algorithm is poor when γ lies in (70 -90 ) and (0°-20°) as shown in Figure 4a. It because the decorrelation of the polarization smoothing is poor when the single polarized components dominate. Unlike the other schemes, when the single polarized components dominate the proposed method does not need to solve the coherence between the direct signal and specular signals. Hence, our proposed scheme retains more stable performance for all the polarization angles as shown in Figure 4. The RMSE of the proposed method is much higher when γ lies in (0°-20°) in Figure 4a or (80°-90°) in Figure 4b. This is because the SNR is reduced when the phase of fading coefficients is ±180° for single polarized dominates. To achieve more stable angle estimation performance for all DOA of grazing angle, γ = 45° is assumed in the following simulations.    Figure 5). The RMSE of the proposed method still is consistent with the CRB curve, even when the grazing angle is around 0.5°. That is because the decorrelation of polarization smoothing is poor for polarization smoothing when the grazing angle is close to zero. While the proposed method doesn't need to solve the coherent problem. Hence the better performance is expected. Figure 6a,b represent the RMSE with respect to SNR with DOA of is 1° and 2° which are close to zero. As expected, the angular accuracy is improved with an increase in SNR. The proposed method has the lowest RMSE of angle estimation for different SNR in comparison with PS-MUSIC, PS-SS-MUSIC and IBMUSIC algorithms ( Figure 6). Under the low-SNR condition, the PS-MUSIC and PS-SS-MUSIC corrupts ( Figure 6). However, the proposed method still is consistent with the CRB curve, even in the low SNR conditions.   Figure 5). The RMSE of the proposed method still is consistent with the CRB curve, even when the grazing angle is around 0.5 • . That is because the decorrelation of polarization smoothing is poor for polarization smoothing when the grazing angle is close to zero. While the proposed method doesn't need to solve the coherent problem. Hence the better performance is expected. Figure 6a,b represent the RMSE with respect to SNR with DOA of θ d is 1 • and 2 • which are close to zero. As expected, the angular accuracy is improved with an increase in SNR. The proposed method has the lowest RMSE of angle estimation for different SNR in comparison with PS-MUSIC, PS-SS-MUSIC and IBMUSIC algorithms ( Figure 6). Under the low-SNR condition, the PS-MUSIC and PS-SS-MUSIC corrupts ( Figure 6). However, the proposed method still is consistent with the CRB curve, even in the low SNR conditions.

Pseudo Spectra of the Proposed Method
As we can see from Figure 7, when the DOA of grazing angle is close to zero, the pseudo spectra of PS-MUSIC and SS-PS-MUSIC algorithm do not have obvious peaks, while the PS-GMUSIC algorithm still has a peak. That indicates that the proposed method has better angle resolution and can achieve better low-angle estimation performance when the grazing angle is close to zero.

Conclusions
The PS-GMUSIC method for low-angle estimation with PSA radar system is proposed in this paper. The proposed algorithm takes advantage of the polarization smoothing method and the generalized MUSIC algorithm to achieve better low-angle estimation performance. The PS-GMUSIC algorithm utilizes the composite steering vector instead of the ordinary steering vector, which does not need to solve the coherent problem. Meanwhile, the information of the symmetry geometry of multipath model is exploited to reduce the computational burden for the proposed algorithm. The PS-GMUSIC has a more stable performance for all polarization angles. In addition, the proposed method has better angle estimation performance in low-angle area in comparison with PS-MUSIC, PS-SS-MUSIC and IBMUSIC algorithms, especially when the grazing angle is close to zero. The

Pseudo Spectra of the Proposed Method
As we can see from Figure 7, when the DOA of grazing angle θ d is close to zero, the pseudo spectra of PS-MUSIC and SS-PS-MUSIC algorithm do not have obvious peaks, while the PS-GMUSIC algorithm still has a peak. That indicates that the proposed method has better angle resolution and can achieve better low-angle estimation performance when the grazing angle is close to zero.

Pseudo Spectra of the Proposed Method
As we can see from Figure 7, when the DOA of grazing angle is close to zero, the pseudo spectra of PS-MUSIC and SS-PS-MUSIC algorithm do not have obvious peaks, while the PS-GMUSIC algorithm still has a peak. That indicates that the proposed method has better angle resolution and can achieve better low-angle estimation performance when the grazing angle is close to zero.

Conclusions
The PS-GMUSIC method for low-angle estimation with PSA radar system is proposed in this paper. The proposed algorithm takes advantage of the polarization smoothing method and the generalized MUSIC algorithm to achieve better low-angle estimation performance. The PS-GMUSIC algorithm utilizes the composite steering vector instead of the ordinary steering vector, which does not need to solve the coherent problem. Meanwhile, the information of the symmetry geometry of multipath model is exploited to reduce the computational burden for the proposed algorithm. The PS-GMUSIC has a more stable performance for all polarization angles. In addition, the proposed method has better angle estimation performance in low-angle area in comparison with PS-MUSIC, PS-SS-MUSIC and IBMUSIC algorithms, especially when the grazing angle is close to zero. The

Conclusions
The PS-GMUSIC method for low-angle estimation with PSA radar system is proposed in this paper. The proposed algorithm takes advantage of the polarization smoothing method and the generalized MUSIC algorithm to achieve better low-angle estimation performance. The PS-GMUSIC algorithm utilizes the composite steering vector instead of the ordinary steering vector, which does not need to solve the coherent problem. Meanwhile, the information of the symmetry geometry of multipath model is exploited to reduce the computational burden for the proposed algorithm. The PS-GMUSIC has a more stable performance for all polarization angles. In addition, the proposed method has better angle estimation performance in low-angle area in comparison with PS-MUSIC, PS-SS-MUSIC and IBMUSIC algorithms, especially when the grazing angle is close to zero. The strategy that proposed in this paper also can be exploited in the MIMO radar system for the low-angle estimation problem.

Appendix A
The CRB for PSA radar system is derived in this part. Fisher information Matrix (FIM) is the foundation of the Cramer-Rao Bound. The FIM of the multipath model in this paper can be denoted as: where F u i u j is the (i, j)th element of the FIM, which is defined as: There are nine unknown parameters θ d , η, γ, ρ h , ρ h , ρ v , ρ v , σ 2 s , σ 2 n in the covariance matrix R. and u i,j ∈ θ d , η, γ, ρ h , ρ h , ρ v , ρ v , σ 2 s , σ 2 n , where L is number of the snapshots of the received signal. When u i , u j ∈ σ 2 s , σ 2 n F u i u j can be simplified as: When u i , u j ∈ {θ d , η, γ, ρ h , ρ h , ρ v , ρ v }, F u i u j can be simplified as: where c is the composite steering vector of the received signal for PSA; and σ 2 s is the variance of the signal. When u i ∈ σ 2 s , σ 2 n , u j ∈ {θ d , η, γ, ρ h , ρ h , ρ v , ρ v }, F u i u j can be simplified as: In order to establish the FIM, we only need to know the partial derivatives of composite steering vector c with respect to u i ∈ {θ d , η, γ, ρ h , ρ h , ρ v , ρ v }. According to the definition of the composite steering vector c in Equation (9), the partial derivatives can be derived as: where: