(i) {α - Convex Optimization
(i) {α - Convex Optimization (i) {α - Convex Optimization
Class of problems in CS (cont’d) (P σ ), (P τ ) and (P λ ) can all be cast as: (P1) : min α f 1(α) + f 2 (α) f 1 and f 2 are proper lsc convex functions. (P λ ) 1 2 ‖y − HΦα‖2 l 2 + λΨ (α) f 1 (α) = 1 2 ‖y − HΦα‖2 l 2 ,f 2 (α) =λΨ (α) (P σ ) Ψ (α) s.t. ‖y − HΦα‖ l2 ≤ σ f 1 (α) =Ψ (α) ,f 2 (α) =ı Bl2 ,σ (α) (P τ ) 1 2 ‖y − HΦα‖2 l 2 s.t. Ψ (α) ≤ τ f 1 (α) = 1 2 ‖y − HΦα‖2 l 2 ,f 2 (α) =ı Bψ,τ (α) Stanford seminar 08-8
Characterization Theorem 1 (i) Existence: (P1) possesses at least one solution if f 1 + f 2 is coercive, i.e., lim ‖α‖→+∞ f 1 (α)+f 2 (α) = +∞. (ii) Uniqueness: (P1) possesses at most one solution if f 1 + f 2 is strictly convex. This occurs in particular when either f 1 or f 2 is strictly convex. (iii) Characterization: Let α ∈ H. Then the following statements are equivalent: (a) α solves (P1). (b) α =(I + ∂(f 1 + f 2 )) −1 (α) (proximal iteration). ∂f i is the subdifferential (set-valued map), a maximal monotone operator. J ∂(f1 +f 2 ) = (I + ∂(f 1 + f 2 )) −1 is the resolvent of ∂(f 1 + f 2 ) (firmly nonexpansive operator). Stanford seminar 08-9
- Page 1 and 2: Optimization problems in compressed
- Page 3 and 4: Compressed/ive Sensing Stanford sem
- Page 5 and 6: Compressed/ive Sensing Common wisdo
- Page 7 and 8: Compressed/ive Sensing Common wisdo
- Page 9 and 10: Compressed/ive Sensing (cont’d) C
- Page 11 and 12: Compressed/ive Sensing (cont’d) C
- Page 13 and 14: Convex analysis and operator splitt
- Page 15 and 16: Class of problems in CS (cont’d)
- Page 17: Class of problems in CS (cont’d)
- Page 21 and 22: Operator splitting schemes Idea: re
- Page 23 and 24: Proximity operators Some properties
- Page 25 and 26: Example of proximity operator Stanf
- Page 27 and 28: Compressed sensing optimization pro
- Page 29 and 30: Characterizing Problem (P τ ) Stan
- Page 31 and 32: Proximity operators of Ψ Conclusio
- Page 33 and 34: DR to solve Problem (P σ ) Theorem
- Page 35 and 36: DR to solve Problem (Peq) Theorem 1
- Page 37 and 38: Pros and cons (P σ ) and (P eq ) h
- Page 39 and 40: CS reconstruction (1) H = Fourier,
- Page 41 and 42: Inpainting and CS H = Dirac, Φ = C
- Page 43 and 44: Inpainting and CS H = Dirac, Φ = C
- Page 45 and 46: Computation time CS H = Fourier, Φ
- Page 47 and 48: Ongoing and future work Beyond the
Characterization<br />
Theorem 1<br />
(i) Existence: (P1) possesses at least one solution if f 1 + f 2 is coercive, i.e.,<br />
lim ‖α‖→+∞ f 1 (α)+f 2 (α) = +∞.<br />
(ii) Uniqueness: (P1) possesses at most one solution if f 1 + f 2 is strictly convex.<br />
This occurs in particular when either f 1 or f 2 is strictly convex.<br />
(iii) Characterization: Let α ∈ H. Then the following statements are equivalent:<br />
(a) α solves (P1).<br />
(b) α =(I + ∂(f 1 + f 2 )) −1 (α) (proximal iteration).<br />
∂f i is the subdifferential (set-valued map), a maximal monotone operator.<br />
J ∂(f1 +f 2 )<br />
= (I + ∂(f 1 + f 2 )) −1 is the resolvent of ∂(f 1 + f 2 ) (firmly nonexpansive<br />
operator).<br />
Stanford seminar 08-9