English

Selecting Regularization Parameters for nuclear norm type minimization problems

Optimization and Control 2022-04-14 v1

Abstract

The reconstruction of low-rank matrix from its noisy observation finds its usage in many applications. It can be reformulated into a constrained nuclear norm minimization problem, where the bound η\eta of the constraint is explicitly given or can be estimated by the probability distribution of the noise. When the Lagrangian method is applied to find the minimizer, the solution can be obtained by the singular value thresholding operator where the thresholding parameter λ\lambda is related to the Lagrangian multiplier. In this paper, we first show that the Frobenius norm of the discrepancy between the minimizer and the observed matrix is a strictly increasing function of λ\lambda. From that we derive a closed-form solution for λ\lambda in terms of η\eta. The result can be used to solve the constrained nuclear-norm-type minimization problem when η\eta is given. For the unconstrained nuclear-norm-type regularized problems, our result allows us to automatically choose a suitable regularization parameter by using the discrepancy principle. The regularization parameters obtained are comparable to (and sometimes better than) those obtained by Stein's unbiased risk estimator (SURE) approach while the cost of solving the minimization problem can be reduced by 11--18 times. Numerical experiments with both synthetic data and real MRI data are performed to validate the proposed approach.

Keywords

Cite

@article{arxiv.2204.06206,
  title  = {Selecting Regularization Parameters for nuclear norm type minimization problems},
  author = {Kexin Li and Hongwei Li and Raymond H. Chan and You-wei Wen},
  journal= {arXiv preprint arXiv:2204.06206},
  year   = {2022}
}