Minimum-norm solutions of the non-symmetric semidefinite Procrustes problem
Abstract
Given two matrices and a set , a Procrustes problem consists in finding a matrix such that the Frobenius norm of is minimized. When is the set of the matrices whose symmetric part is positive semidefinite, we obtain the so-called non-symmetric positive semidefinite Procrustes (NSPDSP) problem. The NSPDSP problem arises in the estimation of compliance or stiffness matrix in solid and elastic structures. If has rank , Baghel et al. (Lin. Alg. Appl., 2022) proposed a three-step semi-analytical approach: (1) construct a reduced NSPDSP problem in dimension , (2) solve the reduced problem by means of a fast gradient method with a linear rate of convergence, and (3) post-process the solution of the reduced problem to construct a solution of the larger original NSPDSP problem. In this paper, we revisit this approach of Baghel et al. and identify an unnecessary assumption used by the authors leading to cases where their algorithm cannot attain a minimum and produces solutions with unbounded norm. In fact, revising the post-processing phase of their semi-analytical approach, we show that the infimum of the NSPDSP problem is always attained, and we show how to compute a minimum-norm solution. We also prove that the symmetric part of the computed solution has minimum rank bounded by , and that the skew-symmetric part has rank bounded by . Several numerical examples show the efficiency of this algorithm, both in terms of computational speed and of finding optimal minimum-norm solutions.
Cite
@article{arxiv.2406.02203,
title = {Minimum-norm solutions of the non-symmetric semidefinite Procrustes problem},
author = {Nicolas Gillis and Stefano Sicilia},
journal= {arXiv preprint arXiv:2406.02203},
year = {2025}
}
Comments
25 pages, 1 figure. The new version contains an improved algorithm, a proof of convexity of the subproblem needed to compute the minimum-norm solution, an additional experiment on a real data set