Table of Contents
Fetching ...

Subspace Newton's Method for $\ell_0$-Regularized Optimization Problems with Box Constraint

Yuge Ye, Qingna Li

TL;DR

The paper tackles box-constrained $\ell_0$-regularized sparse optimization by introducing a $\tau$-stationary point and a subspace Newton method (BNL0R) that updates inactive components via Newton steps while updating active components with proximal gradient steps. It proves global convergence and local quadratic convergence under mild smoothness/convexity assumptions and demonstrates the method’s effectiveness through extensive compressed sensing and image-recovery experiments, where BN L0R outperforms first-order methods in iterations, time, and accuracy. The approach leverages an active/inactive index partition, an Armijo line search for descent, and a fallback projection-gradient step to ensure robust progress toward a $\tau$-stationary point or a local/global minimizer. The results establish a practical second-order framework for nonconvex, box-constrained $\ell_0$-regularized problems with strong potential for sparse recovery and signal processing applications.

Abstract

This paper investigates the box-constrained $\ell_0$-regularized sparse optimization problem. We introduce the concept of a $τ$-stationary point and establish its connection to the local and global minima of the box-constrained $\ell_0$-regularized sparse optimization problem. We utilize the $τ$-stationary points to define the support set, which we divide into active and inactive components. Subsequently, the Newton's method is employed to update the non-active variables, while the proximal gradient method is utilized to update the active variables. If the Newton's method fails, we use the proximal gradient step to update all variables. Under some mild conditions, we prove the global convergence and the local quadratic convergence rate. Finally, experimental results demonstrate the efficiency of our method.

Subspace Newton's Method for $\ell_0$-Regularized Optimization Problems with Box Constraint

TL;DR

The paper tackles box-constrained -regularized sparse optimization by introducing a -stationary point and a subspace Newton method (BNL0R) that updates inactive components via Newton steps while updating active components with proximal gradient steps. It proves global convergence and local quadratic convergence under mild smoothness/convexity assumptions and demonstrates the method’s effectiveness through extensive compressed sensing and image-recovery experiments, where BN L0R outperforms first-order methods in iterations, time, and accuracy. The approach leverages an active/inactive index partition, an Armijo line search for descent, and a fallback projection-gradient step to ensure robust progress toward a -stationary point or a local/global minimizer. The results establish a practical second-order framework for nonconvex, box-constrained -regularized problems with strong potential for sparse recovery and signal processing applications.

Abstract

This paper investigates the box-constrained -regularized sparse optimization problem. We introduce the concept of a -stationary point and establish its connection to the local and global minima of the box-constrained -regularized sparse optimization problem. We utilize the -stationary points to define the support set, which we divide into active and inactive components. Subsequently, the Newton's method is employed to update the non-active variables, while the proximal gradient method is utilized to update the active variables. If the Newton's method fails, we use the proximal gradient step to update all variables. Under some mild conditions, we prove the global convergence and the local quadratic convergence rate. Finally, experimental results demonstrate the efficiency of our method.

Paper Structure

This paper contains 12 sections, 16 theorems, 151 equations, 7 figures, 6 tables, 1 algorithm.

Key Result

Lemma 2.1

For problem eq-2-1, let $\tau > 0$ satisfy eq-tau. Then

Figures (7)

  • Figure 1: The variation of the number of non-zero components and the objective function value during the iteration process of BNL0R for E1 with $m = 0.25n$
  • Figure 2: Average recovery res and time for E1 with $m = 0.25n$
  • Figure 3: Average recovery res and time for E1 with $m = 0.15n$
  • Figure 4: Average recovery res and time for E2 with $m = 0.25n$
  • Figure 5: Average recovery res and time for E2 with $m = 0.15n$
  • ...and 2 more figures

Theorems & Definitions (35)

  • Definition 1
  • Lemma 2.1
  • proof
  • Lemma 2.2
  • proof
  • Definition 2
  • Proposition 1
  • Proposition 2
  • proof
  • Remark 1
  • ...and 25 more