Global Navigation Satellite System (GNSS) observations are subject to various errors during their propagation process. A reasonable correction of these errors can improve the positioning, navigation, and timing (PNT) service capability. The impact of multipaths on pseudorange observations can reach a decimeters or even meters level. However, their mechanism is complex and there is currently no universally accepted high-precision correction model. The correlation between the pseudorange multipaths (MP) of BDS-2 satellites and satellite elevation has been confirmed, while there have been fewer analyses of the MP characteristics for different frequencies of BDS-3 satellites. The broadcasting of multi-frequency observations in BDS-3 should theoretically make the extracted MP more accurate compared to traditional methods. Based on this, in this contribution, a multi-frequency MP extraction algorithm based on the least squares principle is proposed, which can simultaneously eliminate the influence of higher-order ionospheric delay. The analytical expression for only eliminating first-order ionospheric delay is successfully derived. Subsequently, the characteristics of the MPs extracted from different frequency combinations and the impact of combination noise on the extraction accuracy are discussed. The influence of second-order ionospheric delay on the MPs for each frequency under different combination noises, as well as the periodic behavior exhibited in long-term observations of the BDS-3 medium earth orbit (MEO) and inclined geosynchronous orbit (IGSO) satellites, are also analyzed. Finally, the correlations between the MPs of each frequency of BDS satellite and elevation are quantitatively analyzed based on observations from 35 stations. Overall, this work has positive implications for the study of the MP characteristics of BDS-3 and subsequent modeling efforts.
Keywords: BDS-3; analytical expression; combination noise; correlation analysis; ionospheric delay; multi-frequency observations; pseudorange multipath.