Core-Conditioned Regularized Matrix Tri-Factorization for High-Dimensional Structured Systems
For practitioners needing numerically stable low-rank approximations with diagnostic transparency, this framework offers a regularized alternative to standard decompositions, though it is incremental and limited to specific structured settings.
This paper introduces a regularized matrix tri-factorization framework that explicitly controls the conditioning of the core matrix, enabling joint optimization of reconstruction accuracy and numerical stability. The method provides competitive performance in noisy and ill-conditioned low-rank approximation tasks while offering diagnostic insights such as core condition number, though it is slower than randomized SVD for pure spectral compression and not uniformly competitive with matrix-completion baselines.
This paper studies a regularized matrix tri-factorization \(A\approx PDQ\), where \(P\) and \(Q\) are side factors and \(D\) is a central core whose conditioning can be explicitly regularized or constrained. The formulation is a structured low-rank approximation framework, not a replacement for LU, QR, Cholesky, or the singular value decomposition. In the unregularized full-data Frobenius rank-\(r\) problem, truncated SVD remains the optimal benchmark. The contribution here concerns the regularized and core-conditioned setting, where reconstruction accuracy is treated together with factor scale, numerical conditioning, perturbation behavior, and weighted approximation. The analysis establishes the algebraic scope of the \(PDQ\) representation, proves existence of minimizers under coercive regularization, identifies the non-uniqueness induced by latent-space transformations, derives well-posed block updates for the quadratic full-data objective, and gives product-level perturbation bounds. For exact alternating minimization in the full-data quadratic case, it proves descent, boundedness of iterates, and convergence to a critical point under standard Kurdyka--Łojasiewicz assumptions. A full multi-seed validation indicates competitive behavior in noisy and ill-conditioned low-rank approximation while reporting diagnostics not provided by purely spectral baselines, including the learned core condition number and block-system conditioning. The validation also clarifies the method's limits: randomized SVD remains faster for pure spectral compression, and the current weighted missing-entry variant is not uniformly competitive with matrix-completion baselines. The framework is therefore best viewed as a regularized and diagnostically transparent tri-factorization for settings where approximation quality and numerical conditioning must be controlled jointly.