A Multipath Processing Technology Based on Multiparameter-Combined Observation in GNSS

In the Global Navigation Satellite System (GNSS), themultipath error affected bymany aspects is themain error source that affects satellite navigation, and it is difficult to establish a much more accurate model to analyze it. Based on the ground multipath reflection model, it firstly deeply studies the influence of GNSS satellite orbit parameters on multipath fading frequency and establishes a multipath signal model related to satellite orbit parameters. Secondly, the influence of carrier phase cycle slip, receiver clock adjustment, and GNSS satellite orbit on multiparameter(MP-) combined observations is analyzed in detail based on the measured data. Finally, aiming at the common phenomenon of code-carrier divergence in the Beidou system, the elevation-based pseudorange correction model and a sidereal filtering are built to correct the MP errors; experiments with measured data show that there is a fluctuation range reduction of 35.7% after sidereal filtering when the receiver reaches a steady state.


Introduction
e satellite navigation system has become one of the most important technical means in modern navigation and positioning measurement [1,2]. To pursue more accurate positioning results, it is necessary to analyze the reasons that affect the accuracy of a satellite navigation system and take targeted measures to reduce and correct the relevant errors and obtain a better positioning performance [3,4]. ere are many factors which will affect the positioning error of satellite navigation, among which ionospheric error and tropospheric error can be eliminated by differential technique or some mathematical methods [5,6]. However, the multipath error is affected by many aspects, and it is difficult to establish a more accurate model for it and eliminate it by the general method, so multipath error has become the main error source that affects satellite navigation [7,8]. As early as the GPS demonstration test, the influence of the multipath effect on satellite ranging performance and positioning results was considered by the designers. e mathematical relationship between the multipath effect and pseudorandom code tracking derived by Hagerman in 1973 is the basis of multipath effect analysis. In the late 1970s, the researchers tested the multipath effect and its effect on the receiver [9,10]. e influence of the multipath effect on short baseline carrier phase differential positioning was studied by MIT and Draper laboratories in the early 1980s. For GPS systems, the effect of the multipath effect on carrier phase differential positioning can be reduced to several centimeters by averaging the observation results greater than 1 hour [11,12]. Since then, with the deepening of the research, the research on multipath effect in the field of satellite positioning and navigation mainly focuses on the environment and channel modelling, multipath measurement technology, antimultipath antenna, anti-multipath receiver digital signal processing technology, anti-multipath data postprocessing method, and multipath research under the new signal system and evaluation of anti-multipath performance. Based on the traditional antenna design, an antenna technology with a standing wave ratio below 1.5 m is studied by using L feed and double-layer structure in [4]. A technique of adding parasitic circular metal patch as the top-loaded antenna is studied based on a prefeed network in [11].

Literature Review
Multipath is the main systematic error of the GNSS short baseline positioning. Multipath error cannot be eliminated by the double-differenced technique and is difficult to parameterize, which severely restrict the high-precision GNSS positioning application. A correlation reconstruction algorithm for energy compensation under the BOC system is studied in [13]. A time-division multiplexing dual reference waveform multipath suppression algorithm for TMBOC signals is studied in [14]. Based on Bayesian estimation, a multipath estimation technique of the particle filter algorithm is studied in [15]. Based on maximum likelihood estimation, an off-line mesh search and online iteration method is studied in [16] to suppress short multipath delay. In [17], the authors propose a new multipath mitigation strategy in the coordinate domain that shakes off the formation mechanism of multipath; a convolutional neural network (CNN) long short-term memory (LSTM) method is used to mine the deep multipath features in GNSS coordinate series. An advanced receiver autonomous integrity monitoring (ARAIM) method is proposed to evaluate GNSS integrity in [18], which mainly focuses on validating whether the multipath error model in the ARAIM algorithm is conservative enough for GEO satellites. In [19], the authors address the problem of GNSS multipath mitigation using antenna arrays, in which a new data-dependent beamforming technique is proposed, which is based on the wellknown Capon beamformer. is technique aims to avoid the typical cancellation phenomenon between signal and correlated multipaths, by exploiting the known power of the direct signal at the receiver. To this effect, a measure of the correlation between the signal and multipaths is obtained in the matrix form, and it is then subtracted from the spatial correlation matrix of the received signal. is results in a new spatial correlation matrix that is used for the final Power-Based Capon beamformer. e behaviour of this technique is justified mathematically, and it is supported by several numerical results.
Here, this paper establishes a much more accurate multipath signal model, and the basic ideas can be concluded as follows: firstly, the influences of GNSS satellite orbit parameters on multipath fading frequency are studied and analyzed, and a multipath fading frequency model related to satellite orbit parameters is established. Secondly, the MPcombined observations are extracted by the code subtraction carrier phase, and the influences of carrier phase cycle slip, receiver clock correction, and GNSS satellite orbit to MP observations are analyzed by measured data. Finally, the elevation correlation pseudorange correction model and sidereal filtering are used to correct the common codecarrier divergence phenomenon in the Beidou system. e rest of the paper is organized as follows. In Section 3, the mathematical model of multipath signals is established. In Section 4, the MP-combined observations are analyzed. In Section 5, the suppression of MP signals is discussed. Finally, in Section 6, we conclude our work.

Fading Frequency of the Multipath Signal.
Considering generality, assume that the satellite signal reaches the receiver at the same time as the reflected signal, and the received synthetic signal s(t) can be described as where p(t) is the signal amplitude (or spread spectrum code or message), |Γ| is the attenuation coefficient, Δt is the multipath delay (s), ω 0 is the angular frequency of the direct signal and contains the influence of the doppler shift Δω 0 � 2πΔf 0 , Δϕ m (t) � Δφ m + (Δw m − Δw 0 )t is the phase difference between the reflected signal and the direct signal, where Δφ m is the initial phase of the reflected signal, (Δw m − Δw 0 ) is the Doppler frequency difference between the reflected signal and the direct signal, and |Γ|, Δt, and Δϕ m (t) describe the basic characteristics of the synthesized signal.
Due to the relative motion between the antenna of the ground receiver and the satellite, the multipath delay Δt and phase difference Δϕ m (t) will change randomly with time, which will cause the change of the carrier frequency of the reflected wave and produce an increment Δf m (t) on the carrier nominal frequency, which is called multipath fading frequency. According to the relationship between the frequency and phase of the sinusoidal signal, Δf m (t) can be expressed in [10]: It can be seen from equation (2) that Δf m (t) depends on the Doppler frequency difference (Δw m − Δw 0 ); when it changes rapidly with time, the relative speed between the ground receiving antenna and the satellite will be relatively fast, and Δf m (t) will be larger.

Multipath Reflection Model.
e multipath reflection model (see Figure 1), in which θ s (t) is the elevation angle of the satellite at t time in the Earth-centered Earth-fixed (ECEF) coordinate system. A is the phase center of the user receiver, whose height from the ground is h, R is the incident point of the reflected signal, and L LOS (t) and L NLOS (t) represent the propagation path length of the direct signal and reflected signal, respectively.
Since the satellite is far from the satellite, it is generally considered that the direct signal is parallel to the reflected signal, and the vertical line of the direct signal and the reflected signal is made at R; with the vertical foot is D, the multipath geometric distance delay of the reflected signal relative to the direct signal L(t) can be expressed as follows: According to equation (3), L(t) is a function of antenna height h and satellite elevation angle θ s (t). According to the theory of electromagnetic wave propagation, the relative phase of multipath Δϕ m (t) is explained in [11]: where λ is the signal wavelength. Combine equations (4) with (2), the multipath fading frequency Δf m (t) can be rewritten as

Satellite Orbital Parameters.
In ECEF coordinates, the geometric relationship between the satellite orbits and Earth (see Figure 2, in which, R e is the radius of the Earth whose center is O, and the GNSS receiver antenna is placed at A, a point on the ground (which can fall at any point in the sphere; for easy drawing, it is placed at the north pole, the vector OA ⇀ point to zenith). P is the level plane where the user A is, S represents the spatial position of the satellite relative to the Earth at t time, whose projection point on the horizontal P is D, r(t) � ‖AS ⇀ ‖ is the line of sight distance between the user and the satellite, R s � ‖OS ⇀ ‖ represents the geometric distance between the geocentric center O and the satellite, and the satellite elevation angle θ s (t) is the angle between AS ⇀ and A D ⇀ at time t. According to the geometric relation shown in Figure 2 and the cosine theorem, it can be obtained that Combine (5) with (6), then Δf m (t) can be represented as where h, λ, R e , θ s (t), and r(t) can be regarded as known, that is, Δf m (t) mainly depends on the unknown r ′ (t), which can be understood as the relative velocity in the direction of the signal propagation. Let K � − (2h/λ)((1/R e ) + (sin θ s (t)/r(t))), equation (7) can be simplified as For a ground static receiver, let v , that is [10], where β is the angle between v ⇀ s and I ⇀ s and satisfies where f s is the true near point angle of the satellite and c is the angle between the satellite observation vector and geocentric distance vector. Based on the transformation theory of the coordinate system, v s can be expressed as where where α is the inclination angle of the satellite's orbit, w e is the rate of the Earth's self-rotation, w is the distance of the orbit near the Earth's angle, and w s is the angular rate of the satellite.

MP-Combined Observations
where M p and M ϕ are the specular reflection of pseudorange and carrier phase, respectively, v P and v ϕ are the diffuse reflection of pseudorange and carrier phase, respectively, HW is the time delay of the processor, and URE is the deviation of the satellite orbit. Ignore the carrier phase errors ε ϕ , v ϕ , and M ϕ with very small values; the difference between pseudorange and carrier phase observation at the same frequency, also named as multiparameter(MP)-combined observation, can be obtained as follows: where the ionospheric delay 2cI can be modified with a double-frequency ionospheric correction method described in [13]. e carrier multipath error can be ignored, for its maximum multipath error is much smaller than that of code multipath error [14]. Hence, the MP-combined observation is usually used to describe the multipath signal. According to the method mentioned above, the modified MP f i of f i frequency with dual-frequency ionospheric method can be described as Since the transformation factor between time and distance is a constant c, there is no distinction between the two different physical quantities representing time and distance without confusion in the GNSS positioning field. erefore, equation (18) can be further simplified as follows:

Effect of Carrier Phase Cycle Slip on MP Observations.
If there are cycle slips in the carrier phase of the observations, the calculated carrier phase ionospheric delay is inaccurate, which will, in turn, affect the MP value; in particular, a continuous small cycle slip will cause fluctuations in the MP values. Hence, it is necessary to analyze and detect the influence of the carrier phase slip on MP values. e dual-frequency carrier phase method and ionospheric residual method are generally used to detect the cycle slip [15]. Here, the double-frequency pseudorange observations are used for cycle slip detection.
For the convenience of analysis, the smaller error term in the phase measurement equation is ignored, and the carrier phase measurement equation described in equation (16) can be simplified as follows: e phase observations of two adjacent epochs at times t and t + 1 are calculated; most of the ionospheric and tropospheric delays can be eliminated without the carrier phase slip [16], and the carrier phase difference Δϕ can be expressed as e two-frequency phase observations for GPS L1 and L2 satisfies Let Technically, the cycle slip search of MP observations can be carried out in two steps: firstly, the severe cycle slips and signal tracking interruption are detected with a threshold of 5σ or above [20], where σ is the mean square error (MSE) of carrier phase difference, marked in the output results and repaired in sections. According to the all-day data collected from Beidou on January 3rd in 2018, the cycle slips and their influences on MP of SV5 B1 are shown in Figure 3.
In Figure 3, the X axis represents epoch time (s) or samples and (a) and (b) show the first difference of ionospheric residuals before and after cycle slip repair, respectively, whose Y axis represents the ionospheric residual (m). e red horizontal line in (a) is the severe cycle slip search threshold, while the blue pulse is the severe cycle slip and tracking interruption that needs to be repaired. e green horizontal line in (b) is the position where the small cycle slips appear. e threshold is chosen according to [21]. (c) is the ionospheric residual curve, in which the Y axis represents the ionospheric residual value, the red asterisk is the location of the severe cycle slip and tracking interrupt, and the green circle is the location where the small cycle slips occur. (d) is the final MP observations, where the blue and red curves correspond to B1 and B2 frequencies, respectively.  Mobile Information Systems

Effect of Receiver Clock Adjustment on MP Observations.
ere are usually two methods for the receiver to deal with the local clock [21]: one is to use an interval clock adjustment when the cumulative error reaches a certain threshold (such as 0.5 ms), the local clock is set back or forward twice the threshold (1 ms), and the other is to use a continuous clock  In Figure 5, the X-axis represents epoch time (s) or samples, (a) and (b), Y-axis represents the ionospheric residual (m), which is the first difference of ionospheric residuals before and after the cycle slip repairment, respectively, where the red and green horizontal lines are the cycle slip detection thresholds. (c) is the ionospheric residual curve, where the red mark indicates the location of the cycle slip, and (d) is the final MP observations. It can be seen from Figures 4 and 5 that there is no serious cycle slip with the ionospheric residuals, and the high order difference results of the carrier phase are stable within ±20 cycles.

MP Observations of IGSO/MEO Satellites.
ere exists code and carrier divergence in the MP observations of Beidou IGSO/MEO satellites which is related to the elevation angle. In the process of satellite transits, the MP decreases gradually with the increase of the elevation angle, and conversely, the MP value increases with the falling of the elevation angle [21]. According to the MP measurement theory, the reason for the decrease of the MP is that the pseudorange measurements are small or the carrier phase    All the data in Figures 6 and 7 are collected with the same receiver.
e analysis of the IGS MEGX data shows that there are no obvious changes related to the elevation angle in the MP observations of QZSS, which is independent of the multipath environment and its surrounding area and more inclined to come from the satellite itself.

MP Observations of GEO Satellites.
MP observations of GEO satellites are more complex than that of IGSO; the GEO satellites of the Satellite-Based Augmentation System (SBAS) also emit L1 and L5 dual-frequency signals. e MP change law of the SBAS system is different from that of Beidou GEO satellites. MP observations of SBAS S127 collected at CUT0, JFNG, and MAYG stations with an elevation of 16°, 17°, and 71°, respectively, are shown in Figure 8.
e MP values after the 2500 epoch in Figure 8 change strongly for both L1 and L5, and the fluctuation period of L5 is 2 hours at the same time in the receiver of the three different stations, which is caused mainly by the satellite signals. e fluctuation of L1 decreases with the increase of the elevation angle.

Pseudorange Correction Model Based on Elevation.
MP observations of Beidou satellites show a change related to elevation, which is called the code-carrier divergence phenomenon [21]. In Section 4, a mathematical model of pseudorange multipath error is established, which can be applied to MEO, IGSO, and GEO satellites. But it does not consider the code-carrier divergence phenomenon of the Beidou satellite, so a statistical model based on elevation can be used to describe the code and carrier divergence of Beidou IGSO and MEO satellites. Instead of considering the cause of code and carrier divergence, the model only repairs it statistically. e GPS SVN49 satellite launched in March 2009 also showed a code and carrier divergence; different from that of Beidou, the code and carrier divergence of the L1 frequency of SVN49 is not obvious when the elevation is less than 40°; otherwise, it increases with the elevation; however, the code and carrier divergence of L2 frequency will decrease first and then increase when the elevation exceeds 40°. Hauschild et al. [22] collected the orthogonal signal of GPS SVN49 with a high gain antenna, separated the delay, relative amplitude, and phase difference of the multipath signal on L1 and L2 by impulse response and chip shape, explained the origin of the multipath by the GPS satellite antenna parameters and radiofrequency circuit structure, and, finally, reproduced the MP fluctuation and code-carrier divergence with these parameters.
Using the above model, the simulation results of MP fluctuation of the GEO satellite are illustrated (see Figure 9).
Based on the statistical analysis of the measured MP, an elevation-based pseudorange correction model is established with B1, B2, and B3, three frequencies, using the continuous observations for several cycles, which is shown in detail (see Table 1).
Because of the ambiguity of carrier phase observation, the MP value is a relative quantity, and its average value is artificially set to zero. erefore, the absolute value of the pseudorange correction cannot be obtained when the elevation pseudorange correction model is established based on the MP observations. at is, if a constant is added to the pseudorange correctionΔp in the model, the modified MP value remains unchanged, but the location results may change.

Sidereal Filtering Based on Satellite Elevation.
Usually, the multipath error period of the GEO satellite is between 86184 and 86201 seconds, which is very close to a stellar day period (86164 seconds). erefore, to weaken the influence of multipath error on positioning, a sidereal filtering is often used to correct the pseudorange observation.  Mobile Information Systems e main principle of the sidereal filtering is to use the strong correlation of the two consecutive days' data to correct the multipath error. Usually, the data features of the first day are extracted to correct the data for the second day. Due to the difference of 236 seconds between a stellar day and a solar day, it is necessary to translate the previous day's data when using a sidereal filtering [23]. e principle of the sidereal filtering is shown in Figure 10.
e pseudorange MP fluctuation of Beidou GEO and IGSO satellites can be alleviated with a sidereal filtering. In  the elevation-based pseudorange correction model, the mean value of the all-day pseudorange measurement must remain unchanged; the pseudorange measurement noises of the two consecutive days are independent, even if the lowpass filter has filtered out the high-frequency component of the noise; the cost of the sidereal filtering is the increase of the measurement noise. Data of Beidou C06 collected with the Trimble BD930 receiver on October 7th and 8th in 2018 is used for analysis; the MP observations before and after the correction of sidereal filtering are shown in Figure 11, where the MP observations range within (− 0.4, 1) m before sidereal filtering, but after sidereal filtering, the MP observations range within (− 0.3, 0.2) m with a fluctuation range reduction of 35.7% when the receiver reaches a steady state.

Conclusions
In GNSS systems, satellite multipath interference is a nonnegligible factor restricting the overall performance of the satellite navigation systems. Based on the ground multipath reflection model, it deeply studies the influence of GNSS satellite orbit parameters on multipath fading frequency and establishes a multipath signal model related to satellite orbit parameters. rough collecting the measured data of different stations and satellites, the influence of carrier phase cycle slip, receiver clock adjustment, and GNSS satellite orbit on MP observation is analyzed in detail. e statistical analysis of the measured data shows that the multipath fading characteristics of the IGSO satellite and the MEO satellite are related to the satellite elevation, and there is a code-carrier divergence phenomenon, which is not obvious in the GEO satellite. Finally, aiming at the common phenomenon of code-carrier divergence in the Beidou system, the elevation-based pseudorange correction model and a sidereal filtering are built to correct the MP errors, and finally, the measured data are used to verify the model. e MP model built in this paper explains the source of the periodic fluctuation of the MP value but does not explain the phenomenon of Beidou satellite codecarrier divergence. Also, the elevation model and sidereal filtering can be used to repair code-carrier divergence in precise single point positioning, but its practical application is very limited. So it is necessary to study the cause of code-carrier divergence and establish a widely used MP model in the future.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.