Primal and dual mixed-integer least-squares

distributional statistics and global algorithm

Journal Article (2024)
Author(s)

Peter J G Teunissen (Curtin University, TU Delft - Mathematical Geodesy and Positioning, University of Melbourne)

L. Massarweh (TU Delft - Mathematical Geodesy and Positioning)

Research Group
Mathematical Geodesy and Positioning
DOI related publication
https://doi.org/10.1007/s00190-024-01862-1
More Info
expand_more
Publication Year
2024
Language
English
Research Group
Mathematical Geodesy and Positioning
Issue number
7
Volume number
98
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Abstract

In this contribution we introduce the dual mixed-integer least-squares problem and study it in relation to its primal counterpart. The dual differs from the primal formulation in the order in which the integer ambiguity vector a∈Zn and baseline vector b∈Rp are estimated. As not the ambiguities, but rather the entries of b are usually the parameters of interest, the attractiveness of the dual formulation stems from its direct computation of b. It is shown that this potential advantage relies on the ease with which an implicit integer least-squares problem of the dual can be solved. For the convoluted cases, we introduce two methods of simplifying approximations. To be able to describe their quality, we provide a complete distributional analysis of their estimators, thus allowing users to judge whether or not the approximations are acceptable for their application. It is shown that this approach implicitly introduces a new class of admissible integer estimators of which we also determine the pull-in regions. As the dual function is shown to lack convexity, special care is required to be able to compute its global minimizer bˇ. Our proposed method, which has finite termination with a guaranteed ϵ-tolerance, is constructed from combining the branch-and-bound principle, with a special convex-relaxation of the dual, to which the projected-gradient-descent method is applied to obtain the required bounds. Each of the method’s three constituents are described, whereby special emphasis is given to the construction of the required continuously differentiable, convex lower bounding function of the dual.