In this paper, we solve the ℓ2-ℓ1 sparse recovery problem by transforming the objective function of this problem into an unconstrained differentiable function and applying a limited-memory trust-region method. Unlike gradient projection-type methods, which uses only the current gradient, our approach uses gradients from previous iterations to obtain a more accurate Hessian approximation. Numerical experiments show that our proposed approach eliminates spurious solutions more effectively while improving computational time.
Reconstructing high-dimensional sparse signals from low-dimensional low-count photon observations is a challenging nonlinear optimization problem. In this paper, we build upon previous work on minimizing the Poisson log-likelihood and incorporate recent work on the generalized nonconvex Shannon entropy function for promoting sparsity in solutions. We explore the effectiveness of the proposed approach using numerical experiments.