English

Minimum-norm solutions of the non-symmetric semidefinite Procrustes problem

Numerical Analysis 2025-09-03 v2 Numerical Analysis Optimization and Control

Abstract

Given two matrices X,BRn×mX,B\in \mathbb{R}^{n\times m} and a set ARn×n\mathcal{A}\subseteq \mathbb{R}^{n\times n}, a Procrustes problem consists in finding a matrix AAA \in \mathcal{A} such that the Frobenius norm of AXBAX-B is minimized. When A\mathcal{A} 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 XX has rank rr, Baghel et al. (Lin. Alg. Appl., 2022) proposed a three-step semi-analytical approach: (1) construct a reduced NSPDSP problem in dimension r×rr\times r, (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 rr, and that the skew-symmetric part has rank bounded by 2r2r. Several numerical examples show the efficiency of this algorithm, both in terms of computational speed and of finding optimal minimum-norm solutions.

Keywords

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

R2 v1 2026-06-28T16:52:46.543Z