English

An $\ell_p$ Variable Projection Method for Large-Scale Separable Nonlinear Inverse Problems

Numerical Analysis 2021-06-01 v1 Numerical Analysis

Abstract

The variable projection (VarPro) method is an efficient method to solve separable nonlinear least squares problems. In this paper, we propose a modified VarPro for large-scale separable nonlinear inverse problems that promotes edge-preserving and sparsity properties on the desired solution and enhances the convergence of the parameters that define the forward problem. We adopt a majorization minimization method that relies on constructing a quadratic tangent majorant to approximate a general p\ell_p regularized problem by an 2\ell_2 regularized problem that can be solved by the aid of generalized Krylov subspace methods at a relatively low cost compared to the original unprojected problem. In addition, we can use more potential general regularizers including total variation (TV), framelet, and wavelets operators. The regularization parameter can be defined automatically at each iteration by means of generalized cross validation. Numerical examples on large-scale two-dimensional imaging problems arising from blind deconvolution are used to highlight the performance of the proposed approach in both quality of the reconstructed image as well as the reconstructed forward operator.

Keywords

Cite

@article{arxiv.2105.14155,
  title  = {An $\ell_p$ Variable Projection Method for Large-Scale Separable Nonlinear Inverse Problems},
  author = {Malena Espanol and Mirjeta Pasha},
  journal= {arXiv preprint arXiv:2105.14155},
  year   = {2021}
}

Comments

23 pages, 9 figures

R2 v1 2026-06-24T02:35:30.817Z