Robust Sparse Recovery via Non-Convex Optimization

Transcription

1 Robust Sparse Recovery via Non-Convex Optimization Laming Chen and Yuantao Gu Department of Electronic Engineering, Tsinghua University Homepage: August 23, 2014

2 Contents 1 Introduction 2 Preliminary 3 Main Contribution 4 Simulation 5 Summary

3 Contents 1 Introduction 2 Preliminary 3 Main Contribution 4 Simulation 5 Summary

4 Introduction Sparse Recovery Problem Sparse recovery problem is expressed as argmin x 0, subject to y = Ax x where A R M N is a fat sensing matrix. Computationally intractable Sparse Recovery Algorithms Greedy pursuits (OMP, CoSaMP, SP, etc) Convex optimization based algorithms Non-convex optimization based algorithms Better performance tends to be derived Still lacking complete convergence analysis

5 Introduction Current Convergence Analysis For p-irls (0 < p < 1), the convergence is guaranteed in a sufficiently small neighborhood of the sparse signal. For MM subspace algorithm, it is shown that the generated sequence will converge to a critical point. For SL0, the complete convergence analysis is done due to the local convexity of the penalties, and the algorithm needs to solve a sequence of optimization problems rather than a single one to guarantee convergence. Objective To provide theoretical convergence guarantees of a non-convex approach for sparse recovery from the initial solution to the global optimal.

6 Contents 1 Introduction 2 Preliminary 3 Main Contribution 4 Simulation 5 Summary

7 Preliminary Problem Setup Non-convex optimization problem is formulated as argminj(x) = x N F (x i ), subject to y = Ax i=1 F ( ) belongs to a class of sparseness measures (proposed by Rémi Gribonval et al.) No. 1 No. 2 No. 3 No. 4 No. 5 No

8 Preliminary Null Space Property Define γ(j, A, K) as the smallest quantity such that J(z S ) γ(j, A, K)J(z S c) holds for any set S {1, 2,..., N} with #S K and for any z N (A). Performance analysis with γ(j, A, K) (provided by Rémi Gribonval et al.) If γ(j, A, K) < 1, for any K-sparse x and y = Ax, the optimization problem returns x ; If γ(j, A, K) > 1, there exists a K-sparse x and y = Ax such that the problem returns a signal differs from x ; γ(l 0, A, K) γ(j, A, K) γ(l 1, A, K).

9 Preliminary Weak Convexity F ( ) is ρ-convex if and only if ρ is the largest quantity such that there exists convex H( ) and F (t) = H(t) + ρt 2. F ( ) is strongly convex if ρ > 0; convex if ρ = 0; weakly convex (also known as semi-convex) if ρ < 0. Generalized gradient (a generalization of subgradient for convex functions) for weakly convex functions can be seen as any element in the set F (t) := H(t) + 2ρt.

10 Contents 1 Introduction 2 Preliminary 3 Main Contribution 4 Simulation 5 Summary

11 Main Contribution Sparsity-inducing Penalty F ( ) is a weakly convex sparseness measure; Vast majority of non-convex sparsity-inducing penalties are formed by weakly convex sparseness measures; There exists α such that F (t) α t. Define ρ/α as the non-convexity of the penalty (a measure of how quickly the generalized gradient of F ( ) decreases) ρ/α = 1/3 ρ/α = 1/2 ρ/α =

12 Main Contribution Projected Generalized Gradient (PGG) Method Initialization: Calculate A, x(0) = A y, n = 0; Repeat: x(n + 1) = x(n) κ J(x(n)); x(n + 1) = x(n + 1) + A (y A x(n + 1)); n = n + 1; Until: Stopping criterion satisfied; Gradient Update x* x(0) y=ax Projection

13 Main Contribution Performance of PGG Theorem 1. Assume y = Ax + e and x 0 K. For any tuple (J, A, K) with J( ) formed by weakly convex sparseness measure F ( ) and γ(j, A, K) < 1, and for any positive constant M 0, if the non-convexity of J( ) satisfies ρ α 1 M 0 1 γ(j, A, K) 5 + 3γ(J, A, K), the recovered solution ˆx by PGG satisfies ˆx x 2 4α2 N C 1 κ + 8C 2 e 2 provided that x(0) x 2 M 0.

14 Main Contribution Approximate PGG (APGG) Method To reduce the computational burden of calculating A = A T (AA T ) 1, assume we adopt A T B to approximate A, and let ζ = I AA T B 2. Theorem 2. Under the same assumptions as Theorem 1, if the non-convexity of J( ) satisfies ρ α 1 M 0 1 γ(j, A, K) 5 + 3γ(J, A, K), and ζ < 1, the recovered solution ˆx by APGG satisfies ˆx x 2 2C 3 (ζ)κ + 2C 4 (ζ) e 2 provided that x(0) x 2 M 0.

15 Contents 1 Introduction 2 Preliminary 3 Main Contribution 4 Simulation 5 Summary

16 Simulation Parameter Setting Matrix A: ; independent and identically distributed Gaussian entries; Vector x : nonzero entries independently satisfy Gaussian distribution; normalized to have unit l 2 norm; The approximate A is calculated using an iterative method (introduced by Ben-Israel et al.);

17 Simulation First Experiment Recovery performance of the PGG method is tested in the noiseless scenario with different sparsity-inducing penalties and different choices of non-convexity. Kmax No.1 No.2 No.3 No.4 No.5 No Non convexity

18 Simulation Second Experiment The recovery performance of the (A)PGG method is compared in the noiseless scenario with some typical sparse recovery algorithms Successful Recovery Probability OMP 0.3 L1 RL1 0.2 ISL0 IRLS 0.1 PGG APGG K

19 Simulation Third Experiment The recovery precisions of the (A)PGG method are simulated under different settings of step size and measurement noise. Average RSNR PGG APGG MSNR = 20dB MSNR = 30dB MSNR = 40dB MSNR = 50dB MSNR = infinity Step Size κ 10 6

20 Contents 1 Introduction 2 Preliminary 3 Main Contribution 4 Simulation 5 Summary

21 Summary Conclusion A class of weakly convex sparseness measures is adopted to constitute the sparsity-inducing penalties; The convergence analysis of the (A)PGG method reveals that when the non-convexity is below a threshold, the recovery error is linear in both the step size and the noise term; As for the APGG method, the influence of the approximate projection is reflected in the coefficients instead of an additional error term.

22 Related Works [1] Laming Chen and Yuantao Gu, The Convergence Guarantees of a Non-convex Approach for Sparse Recovery, IEEE Transactions on Signal Processing, 62(15): , [2] Laming Chen and Yuantao Gu, The Convergence Guarantees of a Non-convex Approach for Sparse Recovery Using Regularized Least Squares, IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP), , May 4-9, 2014, Florence, Italy. [3] Laming Chen and Yuantao Gu, Robust Recovery of Low-Rank Matrices via Non-Convex Optimization, International Conference on Digital Signal Processing (DSP), , Aug , 2014, Hong Kong, China.

23 Thanks!

24 Supplementary Weakly Convex Sparseness Measures with Parameter ρ No. F (t) ρ 1. t 0 2. (p 1)σ p 2 t ( t +σ) 1 p 3. 1 e σ t σ 2 /2 4. ln(1 + σ t ) σ 2 /2 5. atan(σ t ) 3 3σ 2 /16 6. (2σ t σ 2 t 2 )X t 1 σ + X t > 1 σ σ 2