Smooth Bias Estimation for Multipath Mitigation Using Sparse Estimation

Abstract

Multipath remains the main source of error when using global navigation satellite systems (GNSS) in constrained environment, leading to biased measurements and thus to inaccurate estimated positions. This paper formulates the GNSS navigation problem as the resolution of an overdetermined system, which depends nonlinearly on the receiver position and linearly on the clock bias and drift, and possible biases affecting GNSS measurements. The extended Kalman filter is used to linearize the navigation problem whereas sparse estimation is considered to estimate multipath biases. We assume that only a part of the satellites are affected by multipath, i.e., that the unknown bias vector is sparse in the sense that several of its components are equal to zero. The natural way of enforcing sparsity is to introduce an $\ell_1$ regularization associated with the bias vector. This leads to a least absolute shrinkage and selection operator (LASSO) problem that is solved using a reweighted-$\ell_1$ algorithm. The weighting matrix of this algorithm is designed carefully as functions of the satellite carrier to noise density ratio and the satellite elevations. The smooth variations of multipath biases versus time are enforced using a regularization based on total variation. An experiment conducted on real data allows the performance of the proposed method to be appreciated.

Publication
Proceedings of IEEE International Conference on Information Fusion (FUSION)

Matlab folder contains:

  • the same positions creation script, weights compuation, and $\ell_1$ algorithms as in IEEE Trans. Intell. Transp. Syst code
  • create_measurements_fusion.m, which creates the measurements from the measurement equation (19)
  • main_fusion.m, the main file to launch
  • plot_results_fusion.m, which plots the estimated biases and various RMSE over time
  • position_filtering_fusion.m, which filters the state vector over time thanks to the measurements and the sparse bias estimation according to Section II.B) steps 1) to 5) and Section IV.A) step 3)
  • solveLassoSmooth.m, which solves problem (35)
  • sparse_smooth_gnss.m, which estimates the sparse smooth multipath bias
Julien LESOUPLE
Julien LESOUPLE
Lecturer/Researcher

My research interests include statistical signal processing applied to navigation.