Abstract:We explore the connection between two problems that have arisen independently in the signal processing and related fields: the estimation of the geometric mean of a set of symmetric positive definite (SPD) matrices and their approximate joint diagonalization (AJD). Today there is a considerable interest in estimating the geometric mean of a SPD matrix set in the manifold of SPD matrices endowed with the Fisher information metric. The resulting mean has several important invariance properties and has proven very useful in diverse engineering applications such as biomedical and image data processing. While for two SPD matrices the mean has an algebraic closed form solution, for a set of more than two SPD matrices it can only be estimated by iterative algorithms. However, none of the existing iterative algorithms feature at the same time fast convergence, low computational complexity per iteration and guarantee of convergence. For this reason, recently other definitions of geometric mean based on symmetric divergence measures, such as the Bhattacharyya divergence, have been considered. The resulting means, although possibly useful in practice, do not satisfy all desirable invariance properties. In this paper we consider geometric means of covariance matrices estimated on high-dimensional time-series, assuming that the data is generated according to an instantaneous mixing model, which is very common in signal processing. We show that in these circumstances we can approximate the Fisher information geometric mean by employing an efficient AJD algorithm. Our approximation is in general much closer to the Fisher information geometric mean as compared to its competitors and verifies many invariance properties. Furthermore, convergence is guaranteed, the computational complexity is low and the convergence rate is quadratic. The accuracy of this new geometric mean approximation is demonstrated by means of simulations.

New multivariable mean from nonlinear matrix equation associated to the harmonic mean

An Adaptive Test for the Mean Vector in Large-P-small-n Problems

Bounds for the Wasserstein mean with applications to the Lie-Trotter mean

Improved Gaussian Mean Matrix Estimators In High-Dimensional Data

Inequalities for the Wasserstein mean of positive definite matrices

Geometric mean for T-positive definite tensors and associated Riemannian geometry

Weak log-majorization between the geometric and Wasserstein means

Test for mean matrix in GMANOVA model under heteroscedasticity and non-normality for high-dimensional data

Derivatives of Multilinear Functions of Matrices

Tail bounds for Multivariate Random Tensor Means

Neural Tangent Kernel of Matrix Product States: Convergence and Applications

An Efficient Iterative Approach for Hermitian Matrices Having a Fourth-Order Convergence Rate to Find the Geometric Mean

A new weighted spectral geometric mean and properties

Approximate joint diagonalization and geometric mean of symmetric positive definite matrices

Axiomatic Characterization of Nonlinear Homomorphic Means

Nonlinear asymptotic mean value characterizations of holomorphic functions

Guaranteed matrix recovery using weighted nuclear norm plus weighted total variation minimization

Geometric Mean of Partial Positive Definite Matrices with Missing Entries

High order concentrated non-negative matrix-exponential functions

Generalized Matrix Means for Semi-Supervised Learning with Multilayer Graphs