We consider the problem of approximate joint triangularization of a set of noisy jointly diagonalizable real matrices. Approximate joint triangularizers are commonly used in the estimation of the joint eigenstructure of a set of matrices, with applications in signal processing, linear algebra, and tensor decomposition. By assuming the input matrices to be perturbations of noise-free, simultaneo...

The computation of spherical harmonic coefficients of the Earth’s gravity field from satellite-to-satellite tracking (SST) data and satellite gravity gradiometry (SGG) data is considered. As long as the functional model related to SST data contains nuisance parameters (e.g. unknown initial state vectors), assembling of the corresponding normal matrix must be supplied with the back-substitution ...