Mathematical Problems in Engineering

Volume 2009, Article ID 106870, 20 pages

http://dx.doi.org/10.1155/2009/106870

## GPS/Galileo Multipath Detection and Mitigation Using Closed-Form Solutions

^{1}LMSE Laboratory, Electronics Department, The University Center of Bordj Bou Arréridj, Elanasser - Bordj Bou Arréridj 34265, Algeria^{2}LIS Laboratory, Electronics Department, University of M'sila, M'sila 28000, Algeria

Received 11 March 2009; Accepted 3 September 2009

Academic Editor: Francesco Pellicano

Copyright © 2009 Khaled Rouabah and Djamel Chikouche. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We propose an efficient method for the detection of Line of Sight (LOS) and Multipath (MP) signals in global navigation satellite systems (GNSSs) which is based on the use of virtual MP mitigation (VMM) technique. By using the proposed method, the MP signals' delay and coefficient amplitudes can be efficiently estimated. According to the computer simulation results, it is obvious that our proposed method is a solution for obtaining high performance in the estimation and mitigation of MP signals and thus it results in a high accuracy in GNSS positioning.

#### 1. Introduction

In GNSS systems such as Global Positioning System (GPS) and future Galileo the positioning accuracy is seriously degraded in the presence of MP propagation [1]. In effect, MP signals provoke tracking error in the Delay Locked Loop (DLL) [1, 2]. Therefore, it is necessary to eliminate MP errors in the DLL discriminator and track the LOS signals transmitted from satellites. Here, various techniques are proposed to mitigate the MP effect [3–10]. In [5], an analysis of feedback as well as feedforward code tracking algorithms has been done. Then, the peak tracking methods have been proposed as a combination of both feedback and feedforward structures that utilize the inherent advantages of both structures. Yet, at low signal-to-noise ratio (SNR), narrow correlator (NC) [6] is still the best choice among the considered algorithms. The MP Estimating Delay Locked Loop (MEDLL) [10] is one of the analyzed algorithms. It estimates the parameters of LOS and MP signals such as delays, amplitudes, and phases. The former MEDLL proves to have the best performance in MP environments [10, 11]. However, MEDLL requires a lot of hardware resources. Recently, VMM technique for determining the position of the LOS has been introduced [7]. This technique is shown to give superior performance on the detection of the peak location of the LOS signal in the ideal case. However, this technique proves limited due to the finite-bandwidth filter in the receiver which shows that there is an offset between the peak locations of both the virtual LOS and the actual LOS. This is due to the fact that the virtual LOS correlation function (VLOSCF) used to estimate the LOS delay is distorted. Therefore, it is necessary to estimate all the components of the received signal in order to mitigate the estimated MP signals, and locate the peak of LOS signal from the resultant signal, afterwards. Hence, in this paper, we aim at using the concept of VMM introduced in [7] differently in the sense that we use it to estimate the components of LOS and MP signals such as delays and coefficient amplitudes. In effect, we can separate all the components of the received signal. Moreover, we propose closed-form solutions to calculate the delays and the coefficient amplitudes of both LOS and MP signals. We have labeled our proposed method Reference Correlation MP Mitigation (RCMPM). The RCMPM is discussed for Both Coarse Acquisition GPS code (C/A-GPS) and Binary Offset Carrier Galileo code (BOC(1, 1)-Galileo). Since low-energy-path contribution estimators have high variance and should therefore be given weak weights or simply suppressed in combining strategies [12], it is of interest to estimate the parameters of the strongest paths only.

#### 2. Correlation Function (CF)

In the absence of MP signals, the CF of C/A-GPS and BOC(1, 1)-Galileo codes can be approached by the following expressions, respectively:

with

: LOS signal coefficient amplitude,

: code chip spacing,

: phase shift of the pseudo noise (PN) code.

Figure 1 shows the form of these CF’s in the absence of the MP signals.

#### 3. MP Signals Model

In the presence of MP signals the baseband signal model is defined as follows [7]:

with

,

: delay of LOS or MP signal (: fraction of the time delay . : portion of time delay corresponding to the integer number of chips which is determined by signal acquisition processes in the receiver),

: LOS or MP signal amplitudes,

: phase shift due to the MP signals,

: number of MP signals,

: noise.

: PN code.

Some important characteristics of the MP signals are summarised as follows:

(i) The MP signals arrive after the LOS signal because it must travel a longer propagation path.(ii) If the delay of the MP is less than one code PN chip lengths, the internally generated receiver signal partially correlates with it. (iii) The MP signals can be stronger or weaker than the LOS signal.In all the figures of the paper, the delays are normalized with respect to the LOS. In effect, “0” represents .

#### 4. Description of the Proposed Method

##### 4.1. Detection of a Single MP Signal

In presence of a single MP component, the normalized input signal with respect to is defined as follows [7]:

with:

: MP signal coefficient amplitude,

: LOS signal coefficient amplitude,

: fraction of the delay of MP signal.

###### 4.1.1. Case of C/A-GPS Signal

*(1) Presence of a Weak MP Signal*

The relationship between the amplitudes of LOS and MP signals is given as

In presence of a LOS and single weak MP signals, the receiver tries to correlate with these two components. The resulting CF is distorted as shown in Figure 2. Analytically, the LOS and MP signals may be treated separately. Thus, one may consider the CF associated with LOS (LOSCF) and the CF associated with MP signal (MCF). At any point, these two functions can be vector summed to yield the CF associated with the composite signal (CCF). To estimate each component of the received signal, we should estimate its amplitude coefficient and its delay.

The mathematical expression of the MCF can be obtained from (2.1) as follows:
The concept of the VMM, that is, the plot of CBD (VLOSCF), is illustrated in Figure 3.

The VLOSCF has peak amplitude equal to the maximum value of the CCF and has a width of two code chips on the horizontal axis.

The mathematical expression of VLOSCF can be derived from (2.1) as follows:
The mutual maximum of the VLOSCF and CCF is given as follows:
The mathematical expressions of the line segments *PA *and *PB* are given by the following expressions, respectively:
In these two equations, the slopes of and are defined as follows:
Since , then
Inequality (4.8) guarantees that the resulting CF, obtained after subtracting VLOSCF from the CCF, is aligned with MCF as illustrated in Figure 3, that is, (dashed line). The resulting CF is named “the virtual MP CF” (VMCF).

As shown in Figure 3, the VMCF is negative in the point ; this implies that the maximum of the VMCF is aligned with the MCF.

The peak location of the VMCF corresponds to the MP delay . The amplitude characterized by the point can be derived from (4.3) and (4.6) as follows:
The mathematical expressions of the line segments *P’D’ *and *P’E’* are given in the following expressions, respectively:
To estimate the delay and the coefficient amplitude of the MP signal, we use the following proposed equations:
After estimating both delay and amplitude coefficient we can subtract MCF corresponding to the MP signal from CCF. The resulting CF corresponds to only LOSCF. By using the following proposed equations, we can estimate and

*(2) Presence of a Strong MP*

In this case, the relationship between the amplitudes of the LOS and MP signals is given as

In presence of a strong MP signal, the CCF is not aligned with LOSCF but with the MCF as shown in Figure 4. This is true as long as inequality (4.8) is not satisfied. Since the maximum of CCF is aligned with MCF, we apply a virtual CF aligned with MCF “virtual MP CF (VMCF).”

As shown in Figure 5, the resulting CF “virtual LOS CF (VLOSCF),” that is, (dashed line) is aligned with the LOSCF.

The mathematical expression of VMCF can be derived from (2.1) as follows:
The amplitude characterized by the point *P’* can be derived from (4.5) and (4.13) and is given in the expression below
From Figure 5 we can verify the nonsatisfaction of inequality (4.8).

In this case, to estimate the delays and and the coefficient amplitudes and , we use the following proposed equations:

###### 4.1.2. Case of BOC(1, 1)-Galileo Signal

*(1) Presence of a Weak MP*

The relationship between the amplitudes of the LOS and MP signals is given as . Similarly, the mathematical expression of MCF and VLOSCF can be obtained from (2.2), respectively, as follows:
The CFs of the received signal and the concept of the VMM are illustrated in Figures 6 and 7, respectively.

The maximum of the VLOSCF and CCF is given as

The obtained VMCF (A1A2A3A4A5A6 in Figure 7) after subtraction of the VLOSCF is characterized by the coefficient amplitude and the delay .

To estimate the delays and and the coefficient amplitudes and , we use the following proposed equations:

*(2) Presence of a Strong MP*

With the same discussion and as the case of C/A-GPS signals, the resulting CCF and the concept of the VMM are illustrated in Figures 8 and 9, respectively.

The delays and and the coefficient amplitudes and can be estimated by the following proposed equations, respectively:

##### 4.2. Detection of Two MP Signals

In presence of two MP and a LOS signals, the received signal can be expressed by the following equation [7]:

with:

: coefficient amplitude of second MP signal (MP2),

: delay of MP2,

: phase shift due to the MP2.

The relationships of the amplitudes and the delays is:

###### 4.2.1. Case of C/A-GPS Signal

Here, the mathematical expression MCF2 can be obtained from (2.1) as follows:

Figures 10 and 11 illustrate the CF’s (LOSCF, MCF1, MCF2, and CCF) and the concept of VMM, that is, the plot of CBD, respectively.

In the presence of two MP and LOS signals, the peak position of the CCF can be located on the peak of LOSCF or MCF1 or MCF2. In this paper, we discuss only the second case. Another discussion can be done for the first and the third cases.

represents the maximum value of the CCF and the VMCF and is given as

By subtracting the VMCF from the CCF, we get the *VLOSCF *which is aligned on the maximum of the LOSCF, as illustrated in Figure 11. The values of the delays , , and can be estimated using the following proposed equations, respectively:

By solving the proposed equations’ system (4.25) and (4.26), we can estimate the amplitude of LOS and the amplitude of the MP2

Finally, by the subtraction of the two estimated CFs LOSCF and MCF2, we can get the amplitude of MP1 signal as follows:

###### 4.2.2. Case of BOC(1, 1)-Galileo Signal

With a similar discussion, Figures 12 and 13 illustrate the CFs (LOSCF, MPCF1, MPCF2, and CCF) and the concept of VMM, that is, the plot of CBD, respectively.

The values of the delays , , and can be estimated using the following proposed equations, respectively:

We can estimate the amplitudes of LOS and of MP2 signal by solving the proposed equations’ system (4.29) and (4.30).

Finally, by subtracting the two CFs (VLOSCF and MCF2), we can get the amplitude of MP1 signal as follows:

Note that the practical receiver implementation of all these closed form solutions requires that the complete correlation function of the received signal (CCF) to be measured in order to detect the shape and the distortions caused by MP. The CCF can be sampled via using a bank of correlators [10] (see the appendix).

#### 5. Simulation Results

The impact of MP on code tracking accuracy is often represented as an error envelope which represents the maximum error resulting from one single MP with a certain phase delay and amplitude. The same method of analysis applied to the MP-induced error will be used for the discriminators considered herein. It is worth noting that computing the MP-induced code tracking error envelope (CTEE) is equivalent to finding the point, where the discriminator output crosses the origin because it represents the point where the DLL will lock. It is obvious that both the finite-bandwidth filter and the correlator spacing have an influence on the envelope. So, a large correlator spacing will result in a greater susceptibility of the tracking loop with respect to MP. Usually, a narrow finite-bandwidth filter will tend to increase the MP-induced error envelope. Thus, the GNSS positioning accuracy requires a rigorous choice of these two parameters.

In order to demonstrate that our RCMPM method performs better than a single NC, two schemes have been simulated: an NC with 0.1 chip spacing and our RCMPM method. For these simulated schemes, we consider a LOS and single MP signals and a band-limited CF of 20 MHz. The errors are computed versus MP signal that has amplitude of 0.5 and a delay that varies from 0 to 1.5 chips with respect to the LOS. The MP error envelopes are calculated at the maximum points when the MP signal is at in phase or out of phase with respect to the LOS. The result is given in Figures 14 and 15, respectively.

As illustrated in Figure 14 that represents the error envelope for NC, the BOC(1, 1)-Galileo offers better resistance to the long delay MP than the C/A-GPS; however, they have exactly the same envelope for short delay MP. Since the error envelope of the NC for the Galileo receivers is inferior to that of the GPS, it suffices to compare the error envelope of our RCMPM method with that of Galileo receiver. As illustrated in Figure 15, it is clear that RCMPM method performs better than the NC in terms of the error envelope. Whatever the relative delay of the MP is, the offset error of our proposed method is always less than the NC error. In effect, our method shows the best overall MP performance, which is only sensitive for medium MP delays. The code error envelopes decrease to zero for long and short MP delays. Before ending up, we can say that there are great differences between the MP performance obtained by the NC and MP performance obtained by our RCMPM method, because the error envelope of RCMPM is greatly smaller than the NC. In other words, the bias due to the MP is reduced by 40 to 95 percent than that of the NC. Also, the band of variation of the error is completely reduced. This shows that our method has better MP rejection.

#### 6. Conclusion

An efficient method for the detection and mitigation of MP signals is proposed in this paper. This method is derived from the VMM technique [7]. The latter proves limited due to the finite-bandwidth filter in the receiver that creates an offset between the peak location of the virtual LOS and the actual LOS, leading to a code error tracking in positioning the receiver. In our proposed method we have used the concept of the VMM mitigation technique not only for the LOS delay estimation but for estimating amplitudes and delays of both LOS and MP signals as well. The estimated MP has then been subtracted from the composite signal in order to mitigate the MP effect and efficiently estimate the LOS delay. The simulation results have shown that our proposed method performs better than a single NC. Hence, the error and its band of variation are completely reduced vis-à-vis to what we observed for NC scheme.

#### Appendix

The samples of the CCF of both C/A-GPS and BOC(1, 1)-Galileo which are obtained by the use of the bank of correlators are shown in Figure 16.

In fact, we use the bank of correlators to calculate the CCF. The block diagram of the bank of correlators is shown in Figure 17. As illustrated in this figure, we use a number of correlators distributed across the received correlation function. Each correlator is positioned at unique values of *t* and it measures a single value of as shown in Figure 16.

The received signal is correlated with a number of correlators to get samples of the input correlation function .

: the samples of CCF,

: the number of correlators,

: code chip spacing,

: the estimated delay.

#### References

- E. D. Kaplan,
*Understanding GPS: Principles and Applications*, Artech House, Boston, Mass, USA, 1996. View at Zentralblatt MATH - K. Borre, D. M. Akos, N. Bertelsen, P. Rinder, and S. H. Jensen,
*A Software-Defined GPS and Galileo Receiver: A Single-Frequency Approach*, Birkhäuser, Boston, Mass, USA, 2007. View at Zentralblatt MATH - M. S. Braasch, “Performance comparison of multipath mitigating receiver architectures,” in
*Proceedings of the IEEE Aerospace Conference*, vol. 3, pp. 1309–1315, Big Sky, Mont, USA, March 2001. - E. S. Lohan, A. Lakhzouri, and M. Renfors, “Feedforward delay estimators in adverse multipath propagation for Galileo and modernized GPS signals,”
*EURASIP Journal on Applied Signal Processing*, vol. 2006, Article ID 50971, 19 pages, 2006. View at Publisher · View at Google Scholar - M. Z. H. Bhuiyan, E. S. Lohan, and M. Renfors, “Code tracking algorithms for mitigating multipath effects in fading channels for satellite-based positioning,”
*EURASIP Journal on Advances in Signal Processing*, vol. 2008, Article ID 863629, 17 pages, 2008. View at Publisher · View at Google Scholar - A. J. Van Dierendonck, P. Fenton, and T. J. Ford, “Theory and performance of narrow correlator spacing in a GPS receiver,”
*Journal of the Institute of Navigation*, vol. 39, no. 3, pp. 265–283, 1992. View at Google Scholar - Z. Zhang and C. L. Law, “Short-delay multipath mitigation technique based on virtual multipath,”
*IEEE Antennas and Wireless Propagation Letters*, vol. 4, pp. 344–348, 2005. View at Publisher · View at Google Scholar - Z. Zhang, C. L. Law, and Y. L. Guan, “BA–POC-based ranging method with multipath mitigation,”
*IEEE Antennas and Wireless Propagation Letters*, vol. 4, pp. 492–495, 2005. View at Publisher · View at Google Scholar - C.-L. Chang and J.-C. Juang, “An adaptive multipath mitigation filter for GNSS applications,”
*EURASIP Journal on Advances in Signal Processing*, vol. 2008, Article ID 214815, 10 pages, 2008. View at Publisher · View at Google Scholar - R. VanNee, “The multipath estimating delay lock loop,” in
*Proceedings of the 2nd IEEE International Symposium on Spread Spectrum Techniques and Applications (ISSTA '92)*, Yokohama, Japan, December 1992. - S. F. Matilde, A. F. Miguel, and G. A. Ana, “Performance analysis and parameter optimization of DLL and MEDLL in fading multipath environments for next generation navigation receivers,”
*IEEE Transactions on Consumer Electronics*, vol. 53, no. 4, pp. 1302–1308, 2007. View at Publisher · View at Google Scholar - O. Rabaste and T. Chonavel, “Estimation of Multipath Channels with Long Impulse Response at Low SNR via an MCMC Method,”
*IEEE Transactions on Signal Processing*, vol. 55, no. 4, pp. 1312–1325, 2007. View at Google Scholar