Phase Retrieval from Local Measurements: Deterministic Measurement Constructions and Efficient Recovery Algorithms
|
|
- Myrtle Hampton
- 5 years ago
- Views:
Transcription
1 Phase Retrieval from Local Measurements: Deterministic Measurement Constructions and Efficient Recovery Algorithms Aditya Viswanathan Department of Mathematics and Statistics www-personal.umich.edu/ adityavv 2018 SIAM Conference on Imaging Science (IS18), Bologna, Italy June 8 th, 2018
2 Collaborators Mark Iwen Rayan Saab Brian Preskitt Yang Wang Research supported in part by NSF grant DMS
3 Motivating Application From Huang, Xiaojing, et al. Fly-scan ptychography. Scientific Reports 5 (2015). The Phase Retrieval problem arises in many molecular imaging modalities, including X-ray crystallography Ptychography Other applications can be found in optics, astronomy and speech processing. 1 / 21
4 Motivation Inverse Problem Given y, can we recover x? 2 / 21
5 Motivation Ill-Posed Inverse Problem Which is the correct reconstruction? 3 / 21
6 Mathematical Model where find 1 x C d given y i = a i, x 2 + η i i 1,..., D, y i R denotes the phaseless (or magnitude-only) measurements (D measurements acquired), a i C d are known (by design or estimation) measurement vectors, and η i R is measurement noise. 1 (upto a global phase offset) 4 / 21
7 Existing Computational Approaches Alternating projection methods [Fienup, 1978], [Marchesini et al., 2006], [Fannjiang, Liao, 2012] and many others... Methods based on semidefinite programming PhaseLift [Candes et al., 2012], PhaseCut [Waldspurger et al., 2012],... Others Frame-theoretic, graph based algorithms [Alexeev et al., 2014] (Spectral) initialization + gradient descent (Wirtinger Flow) [Candes et al., 2014] Most methods (with provable recovery guarantees) require impractical (global, random) measurement constructions. 5 / 21
8 Today... We discuss a recently introduced fast (essentially linear-time) phase retrieval algorithm based on realistic (deterministic) 2 local measurement constructions. We provide rigorous theoretical recovery guarantees and present numerical results showing the accuracy, efficiency and robustness of the method. (Time Permitting) extensions to 2D and compressive phase retrieval. 2 for a large class of real-world signals 6 / 21
9 Outline 1 Introduction 2 Solving the Phase Retrieval Problem Measurement Constructions Structured Lifting Obtaining Phase Difference Information Angular Synchronization Solving for the Individual Phases 3 Theoretical Guarantees 4 Numerical Simulations 5 Extensions
10 Local Correlation Measurements From Qian, Jianliang, et al. Efficient algorithms for ptychographic phase retrieval. Inverse Problems Appl., Contemp. Math 615 (2014). Each a i is a shift of a locally-supported vector (mask or window) m (j) C d, supp(m (j) ) = [δ] [d], j = 1,..., K 7 / 21
11 Local Correlation Measurements Each a i is a shift of a locally-supported vector (mask or window) m (j) C d, supp(m (j) ) = [δ] [d], j = 1,..., K Define the discrete circular shift operator S l : C d C d with (S l x) j = x l+j. Our measurements are then (y l ) j = x, S l m (j) 2 + η j,l, (j, l) [K] P, P {0,..., d 1} We will consider K δ and P = [d] 0 := {0,..., d 1} 7 / 21
12 Local Correlation Measurements Each a i is a shift of a locally-supported vector (mask or window) m (j) C d, supp(m (j) ) = [δ] [d], j = 1,..., K Define the discrete circular shift operator S l : C d C d with (S l x) j = x l+j. Our measurements are then (y l ) j = x, S l m (j) 2 + η j,l, (j, l) [K] P, P {0,..., d 1} We will consider K δ and P = [d] 0 := {0,..., d 1} 7 / 21
13 Local Correlation Measurements Each a i is a shift of a locally-supported vector (mask or window) m (j) C d, supp(m (j) ) = [δ] [d], j = 1,..., K Define the discrete circular shift operator S l : C d C d with (S l x) j = x l+j. Our measurements are then (y l ) j = x, S l m (j) 2 + η j,l, (j, l) [K] P, P {0,..., d 1} We will consider K δ and P = [d] 0 := {0,..., d 1} 7 / 21
14 Local Correlation Measurements Each a i is a shift of a locally-supported vector (mask or window) m (j) C d, supp(m (j) ) = [δ] [d], j = 1,..., K Define the discrete circular shift operator S l : C d C d with (S l x) j = x l+j. Our measurements are then (y l ) j = x, S l m (j) 2 + η j,l, (j, l) [K] P, P {0,..., d 1} We will consider K δ and P = [d] 0 := {0,..., d 1} 7 / 21
15 What are we measuring? Lifted System: x, S l m(j) 2 = xx, S l m(j) m (j) S l. Example: (6 6 system, δ = 2, blue denotes non-zero entries) xx, S0m (j) m (j) S 0 = xx, / 21
16 What are we measuring? Lifted System: x, S l m(j) 2 = xx, S l m(j) m (j) S l. Example: (6 6 system, δ = 2, blue denotes non-zero entries) xx, S1m (j) m (j) S 1 = xx, / 21
17 What are we measuring? Lifted System: x, S l m(j) 2 = xx, S l m(j) m (j) S l. Example: (6 6 system, δ = 2, blue denotes non-zero entries) xx, S2m (j) m (j) S 2 = xx, / 21
18 What are we measuring? Lifted System: x, S l m(j) 2 = xx, S l m(j) m (j) S l. Example: (6 6 system, δ = 2, blue denotes non-zero entries) xx, S3m (j) m (j) S 3 = xx, / 21
19 What are we measuring? Lifted System: x, S l m(j) 2 = xx, S l m(j) m (j) S l. Example: (6 6 system, δ = 2, blue denotes non-zero entries) xx, S4m (j) m (j) S 4 = xx, / 21
20 What are we measuring? Lifted System: x, S l m(j) 2 = xx, S l m(j) m (j) S l. Example: (6 6 system, δ = 2, blue denotes non-zero entries) xx, S5m (j) m (j) S 5 = xx, / 21
21 What are we measuring? Lifted System: x, S l m(j) 2 = xx, S l m(j) m (j) S l. Example: (6 6 system, δ = 2, blue denotes non-zero entries) Observation: The only entries of xx we can hope to recover (via linear inversion) are supported on a (circulant) band / 21
22 Useful Observations (I) T δ (C d d ): Let T k : C d d C d d T k (A) ij = { Aij, i j mod d < k 0, otherwise. Lifted System Revisited: x, Sl m(j) 2 = T δ (xx ), Sl m(j) m (j) S l. Bottom Line: If we can find m (j) such that Span { Sl m (j) m (j) S l }l,j = T δ(c d d ), then we can recover T δ (xx ). 9 / 21
23 Useful Observations (I) T δ (C d d ): Let T k : C d d C d d T k (A) ij = { Aij, i j mod d < k 0, otherwise. Lifted System Revisited: x, Sl m(j) 2 = T δ (xx ), Sl m(j) m (j) S l. Bottom Line: If we can find m (j) such that Span { Sl m (j) m (j) S l }l,j = T δ(c d d ), then we can recover T δ (xx ). 9 / 21
24 Useful Observations (I) T δ (C d d ): Let T k : C d d C d d T k (A) ij = { Aij, i j mod d < k 0, otherwise. Lifted System Revisited: x, Sl m(j) 2 = T δ (xx ), Sl m(j) m (j) S l. Bottom Line: If we can find m (j) such that Span { Sl m (j) m (j) S l }l,j = T δ(c d d ), then we can recover T δ (xx ). 9 / 21
25 A Highly Structured Linear System! Example linear system M q = ỹ, where q denotes the vectorized (non-zero) entries of T δ (xx ) ỹ denotes the (interleaved) measurements 10 / 21
26 Useful Observations (II) Why this is useful: (a) Diagonal entries of T δ (xx ) are x i 2. (b) Off-diagonals give the relative phases X := xx xx T δ ( X) (j,k) = e i(arg(xj) arg(x k)), j k mod d < δ Phase Synchronization: (a) The leading eigenvector (appropriately normalized) of ( ) ( ) T δ ( X) x x = diag T δ (11 ) diag x x ( ) ( ) x x = diag F ΛF diag x x Note: is the vector of phases of x. x x = [eiφ 1 e iφ 2... e iφ d] T is the (unknown) phase vector! F C d d is the discrete Fourier transform (DFT) matrix 11 / 21
27 Useful Observations (II) Why this is useful: (a) Diagonal entries of T δ (xx ) are x i 2. (b) Off-diagonals give the relative phases X := xx xx T δ ( X) (j,k) = e i(arg(xj) arg(x k)), j k mod d < δ Phase Synchronization: (a) The leading eigenvector (appropriately normalized) of ( ) ( ) T δ ( X) x x = diag T δ (11 ) diag x x ( ) ( ) x x = diag F ΛF diag x x Note: is the vector of phases of x. x x = [eiφ 1 e iφ 2... e iφ d] T is the (unknown) phase vector! F C d d is the discrete Fourier transform (DFT) matrix 11 / 21
28 Leading Eigenvector Phase Vector [ x1 2 x 1 x 2 x 2 x 1 x 2 2 x 2 x 3 x 3 x 2 x 3 2 x 3 x 4 x 4 x 3 x 4 2 x 4 x 1 x 1 x 4] T x 1 2 x 1 x 2 0 x 1 x 4 x 2 x 1 x 2 2 x 2 x x 3 x 2 x 3 2 x 3 x 4 x 4 x 1 0 x 4 x 3 x 4 2 (re-arrange) (T δ(xx ), 2δ 1 entries in band) (normalize) 1 e i(φ1 φ2) 0 e i(φ1 φ4) e i(φ2 φ1) 1 e i(φ2 φ3) 0 0 e i(φ3 φ2) 1 e i(φ3 φ4) e i(φ4 φ1) 0 e i(φ4 φ3) 1 (top eigenvector) [e iφ1 e iφ2 e iφ3 e iφ4 ] T (Signal Reconstruction) [ x1 e i φ1 x 2 e i φ2 x 3 e i φ3 x 4 e i φ4 ] T 12 / 21
29 Leading Eigenvector Phase Vector [ x1 2 x 1 x 2 x 2 x 1 x 2 2 x 2 x 3 x 3 x 2 x 3 2 x 3 x 4 x 4 x 3 x 4 2 x 4 x 1 x 1 x 4] T x 1 2 x 1 x 2 0 x 1 x 4 x 2 x 1 x 2 2 x 2 x x 3 x 2 x 3 2 x 3 x 4 x 4 x 1 0 x 4 x 3 x 4 2 (re-arrange) (T δ(xx ), 2δ 1 entries in band) (normalize) 1 e i(φ1 φ2) 0 e i(φ1 φ4) e i(φ2 φ1) 1 e i(φ2 φ3) 0 0 e i(φ3 φ2) 1 e i(φ3 φ4) e i(φ4 φ1) 0 e i(φ4 φ3) 1 (top eigenvector) [e iφ1 e iφ2 e iφ3 e iφ4 ] T (Signal Reconstruction) [ x1 e i φ1 x 2 e i φ2 x 3 e i φ3 x 4 e i φ4 ] T 12 / 21
30 Leading Eigenvector Phase Vector [ x1 2 x 1 x 2 x 2 x 1 x 2 2 x 2 x 3 x 3 x 2 x 3 2 x 3 x 4 x 4 x 3 x 4 2 x 4 x 1 x 1 x 4] T x 1 2 x 1 x 2 0 x 1 x 4 x 2 x 1 x 2 2 x 2 x x 3 x 2 x 3 2 x 3 x 4 x 4 x 1 0 x 4 x 3 x 4 2 (re-arrange) (T δ(xx ), 2δ 1 entries in band) (normalize) 1 e i(φ1 φ2) 0 e i(φ1 φ4) e i(φ2 φ1) 1 e i(φ2 φ3) 0 0 e i(φ3 φ2) 1 e i(φ3 φ4) e i(φ4 φ1) 0 e i(φ4 φ3) 1 (top eigenvector) [e iφ1 e iφ2 e iφ3 e iφ4 ] T (Signal Reconstruction) [ x1 e i φ1 x 2 e i φ2 x 3 e i φ3 x 4 e i φ4 ] T 12 / 21
31 Leading Eigenvector Phase Vector [ x1 2 x 1 x 2 x 2 x 1 x 2 2 x 2 x 3 x 3 x 2 x 3 2 x 3 x 4 x 4 x 3 x 4 2 x 4 x 1 x 1 x 4] T x 1 2 x 1 x 2 0 x 1 x 4 x 2 x 1 x 2 2 x 2 x x 3 x 2 x 3 2 x 3 x 4 x 4 x 1 0 x 4 x 3 x 4 2 (re-arrange) (T δ(xx ), 2δ 1 entries in band) (normalize) 1 e i(φ1 φ2) 0 e i(φ1 φ4) e iφ1 e iφ1 e i(φ2 φ1) 1 e i(φ2 φ3) 0 e iφ2 0 e i(φ3 φ2) 1 e i(φ3 φ4) e iφ3 = (2δ 1) e iφ2 e iφ3 e i(φ4 φ1) 0 e i(φ4 φ3) 1 e iφ4 e iφ4 (leading eigenvector) [e iφ1 e iφ2 e iφ3 e iφ4 ] T (Signal Reconstruction) [ x1 e i φ1 x 2 e i φ2 x 3 e i φ3 x 4 e i φ4 ] T 12 / 21
32 Leading Eigenvector Phase Vector [ x1 2 x 1 x 2 x 2 x 1 x 2 2 x 2 x 3 x 3 x 2 x 3 2 x 3 x 4 x 4 x 3 x 4 2 x 4 x 1 x 1 x 4] T x 1 2 x 1 x 2 0 x 1 x 4 x 2 x 1 x 2 2 x 2 x x 3 x 2 x 3 2 x 3 x 4 x 4 x 1 0 x 4 x 3 x 4 2 (re-arrange) (T δ(xx ), 2δ 1 entries in band) (normalize) 1 e i(φ1 φ2) 0 e i(φ1 φ4) e iφ1 e iφ1 e i(φ2 φ1) 1 e i(φ2 φ3) 0 e iφ2 0 e i(φ3 φ2) 1 e i(φ3 φ4) e iφ3 = (2δ 1) e iφ2 e iφ3 e i(φ4 φ1) 0 e i(φ4 φ3) 1 e iφ4 e iφ4 (leading eigenvector) [e iφ1 e iφ2 e iφ3 e iφ4 ] T (Signal Reconstruction) [ x1 e i φ1 x 2 e i φ2 x 3 e i φ3 x 4 e i φ4 ] T 12 / 21
33 Recovery Algorithm Define the map A : C d d C D A(Z) (l,j) = Z, Sl m (j) m (j) S l (l,j). and its restriction A Tδ (C d d ) to our subspace. 13 / 21
34 Recovery Algorithm Define the map A : C d d C D A(Z) (l,j) = Z, Sl m (j) m (j) S l (l,j). and its restriction A Tδ (C d d ) to our subspace. In the noisy setting: Step 1: Step 2: Estimate T δ (xx ) by the banded matrix Z = T δ (Z) := ( A 1 T δ (C d d ) y ) ( + A 1 2 T δ (C d d ) y ). 2 Estimate ( the ) phase by computing the leading eigenvector Z of T δ. Z Step 3: Combine phase with of diagonal entries of T δ (Z) to estimate x. 13 / 21
35 Recovery Algorithm Define the map A : C d d C D A(Z) (l,j) = Z, Sl m (j) m (j) S l (l,j). and its restriction A Tδ (C d d ) to our subspace. In the noisy setting: Step 1: Step 2: Estimate T δ (xx ) by the banded matrix Z = T δ (Z) := ( A 1 T δ (C d d ) y ) ( + A 1 2 T δ (C d d ) y ). 2 Estimate ( the ) phase by computing the leading eigenvector Z of T δ. Z Step 3: Combine phase with of diagonal entries of T δ (Z) to estimate x. 13 / 21
36 Recovery Algorithm Define the map A : C d d C D A(Z) (l,j) = Z, Sl m (j) m (j) S l (l,j). and its restriction A Tδ (C d d ) to our subspace. In the noisy setting: Step 1: Step 2: Estimate T δ (xx ) by the banded matrix Z = T δ (Z) := ( A 1 T δ (C d d ) y ) ( + A 1 2 T δ (C d d ) y ). 2 Estimate ( the ) phase by computing the leading eigenvector Z of T δ. Z Step 3: Combine phase with of diagonal entries of T δ (Z) to estimate x. 13 / 21
37 Recovery Algorithm Define the map A : C d d C D A(Z) (l,j) = Z, Sl m (j) m (j) S l (l,j). and its restriction A Tδ (C d d ) to our subspace. In the noisy setting: Step 1: Estimate T δ (xx ) by Cost: O(d δ 3 + δ d log d) flops Step 2: Z = T δ (Z) := ( A 1 T δ (C d d ) y ) ( + A 1 2 T δ (C d d ) y ). 2 Estimate ( the ) phase by computing the leading eigenvector Z of T δ. Cost: O(δ 2 d log d) flops Z Step 3: Combine phase with of diagonal entries of T δ (Z) to estimate x. Total Cost: O(δ 2 d log d + d δ 3 ) flops 13 / 21
38 Outline 1 Introduction 2 Solving the Phase Retrieval Problem Measurement Constructions Structured Lifting Obtaining Phase Difference Information Angular Synchronization Solving for the Individual Phases 3 Theoretical Guarantees 4 Numerical Simulations 5 Extensions
39 Well-Conditioned Linear Systems Theorem (Iwen, V., Wang 2015) Choose entries of the measurement mask (m (i) ) as follows: { { } (m (i) e 4 l/a ) l = 2δ 1 e 2πi i l 2δ 1, l δ 0, l > δ, a := max 4, δ 1 2, i = 1, 2,..., N. Then, the resulting system matrix for the phase differences (step 1), A Tδ, has condition number } κ(a Tδ ) < max {144e 2, 9e2 (δ 1)2. 4 Deterministic (windowed DFT-type) measurement masks! δ is typically chosen to be c log 2 d with c small (2 3). Extensions: oversampling, random masks / 21
40 Recovery Guarantee Theorem (Iwen, Preskitt, Saab, V. 2016) Let x min := min j x j be the smallest magnitude of any entry in x. Then, the estimate z produced by the proposed algorithm satisfies min θ [0,2π] x e iθ z ( ) ( ) 2 x d 2 C x 2 κ η 2 + Cd 1 4 κ η 2, min δ where C R + is an absolute universal constant. This result yields a deterministic recovery result for any signal x which contains no zero entries. A randomized result can be derived for arbitrary x by right multiplying the signal x with a random flattening matrix. (this is also useful for performing sparse phase retrieval!) 15 / 21
41 Outline 1 Introduction 2 Solving the Phase Retrieval Problem Measurement Constructions Structured Lifting Obtaining Phase Difference Information Angular Synchronization Solving for the Individual Phases 3 Theoretical Guarantees 4 Numerical Simulations 5 Extensions
42 Executime Time (secs.) Efficiency FFT time phase retrieval 10 4 Execution Time (D = d2d log 2 (d)e measurements) 10 2 O(d 3 ) 10 0 O(d log 2 d) 10-2 PhaseLift HIO+ER (300 iters.) BlockPR (w/ 60 iters. HIO+ER) BlockPR (Alg. 1) Problem Size, d 16 / 21
43 Robustness to Measurement Errors / 21
44 Robustness to Measurement Errors / 21
45 Reconstruction Error vs. No. of Measurements / 21
46 Summary and Current/Future Research Directions Today Phase retrieval is an immensely challenging problem seen in important applications such as x-ray crystallography. Proposed mathematical framework: Essentially linear-time robust phase retrieval from deterministic local correlation measurement constructions with rigorous recovery guarantee. Current and Future Directions More robust measurement constructions Compressive phase retrieval Extensions to 2D and Ptychographic datasets Continuous problem formulation 19 / 21
47 Outline 1 Introduction 2 Solving the Phase Retrieval Problem Measurement Constructions Structured Lifting Obtaining Phase Difference Information Angular Synchronization Solving for the Individual Phases 3 Theoretical Guarantees 4 Numerical Simulations 5 Extensions
48 Extension 2D Phase Retrieval Preliminary results for 2D masks with tensor product structure Results from 1D extend to 2D; 2D linear system is a tensor product of the 1D linear system (up to row permutations) Eigenvector-based phase synchronization also works calculation of spectral gap and error analysis pending Test Image ( pixels) Recon. (Rel. error ) 20 / 21
49 Extension Compressive Phase Retrieval Model find x C d given Mx 2 + n = y R D where x is k-sparse, with k d, is entry-wise absolute value, and M is a measurement matrix. Measurement Design Assume M = PC where P C D d is an admissible phase retrieval matrix with an associated recovery algorithm Φ P : R D C d, and C C d d is an admissible compressive sensing matrix with an associated recovery algorithm C : C d C d. Recovery Algorithm (Two-stage) C Φ P : R D C d Performance Metrics No. of measurements required is O(k ln(d/k)) Computational cost (sub-linear) is O(k ln c k ln d) 21 / 21
50 Pubs./Preprints/Code (see www-personal.umich.edu/~adityavv) M. Iwen, B. Preskitt, R. Saab and A. Viswanathan. Phase Retrieval from Local Measurements: Improved Robustness via Eigenvector-Based Angular Synchronization. arxiv: , M. Iwen, A. Viswanathan, and Y. Wang. Fast Phase Retrieval from Local Correlation Measurements. SIAM J. Imag. Sci., Vol. 9(4), pp , Oct Compressive Phase Retrieval M. Iwen, A. Viswanathan, and Y. Wang. Robust Sparse Phase Retrieval Made Easy. Appl. Comput. Harmon. Anal., Vol. 42(1), pp , Jan D Phase Retrieval Mark Iwen, Brian Preskitt, Rayan Saab and A. Viswanathan. Phase Retrieval from Local Measurements in Two Dimensions., Proc. SPIE 10394, Wavelets and Sparsity XVII, X, Aug Code
51 Questions?
Fast Angular Synchronization for Phase Retrieval via Incomplete Information
Fast Angular Synchronization for Phase Retrieval via Incomplete Information Aditya Viswanathan a and Mark Iwen b a Department of Mathematics, Michigan State University; b Department of Mathematics & Department
More informationPhase Retrieval from Local Measurements in Two Dimensions
Phase Retrieval from Local Measurements in Two Dimensions Mark Iwen a, Brian Preskitt b, Rayan Saab b, and Aditya Viswanathan c a Department of Mathematics, and Department of Computational Mathematics,
More informationFast and Robust Phase Retrieval
Fast and Robust Phase Retrieval Aditya Viswanathan aditya@math.msu.edu CCAM Lunch Seminar Purdue University April 18 2014 0 / 27 Joint work with Yang Wang Mark Iwen Research supported in part by National
More informationFast Compressive Phase Retrieval
Fast Compressive Phase Retrieval Aditya Viswanathan Department of Mathematics Michigan State University East Lansing, Michigan 4884 Mar Iwen Department of Mathematics and Department of ECE Michigan State
More informationRecovery of Compactly Supported Functions from Spectrogram Measurements via Lifting
Recovery of Compactly Supported Functions from Spectrogram Measurements via Lifting Mark Iwen markiwen@math.msu.edu 2017 Friday, July 7 th, 2017 Joint work with... Sami Merhi (Michigan State University)
More informationPHASE RETRIEVAL FROM LOCAL MEASUREMENTS: IMPROVED ROBUSTNESS VIA EIGENVECTOR-BASED ANGULAR SYNCHRONIZATION
PHASE RETRIEVAL FROM LOCAL MEASUREMENTS: IMPROVED ROBUSTNESS VIA EIGENVECTOR-BASED ANGULAR SYNCHRONIZATION MARK A. IWEN, BRIAN PRESKITT, RAYAN SAAB, ADITYA VISWANATHAN Abstract. We improve a phase retrieval
More informationNonlinear Analysis with Frames. Part III: Algorithms
Nonlinear Analysis with Frames. Part III: Algorithms Radu Balan Department of Mathematics, AMSC, CSCAMM and NWC University of Maryland, College Park, MD July 28-30, 2015 Modern Harmonic Analysis and Applications
More informationPhase recovery with PhaseCut and the wavelet transform case
Phase recovery with PhaseCut and the wavelet transform case Irène Waldspurger Joint work with Alexandre d Aspremont and Stéphane Mallat Introduction 2 / 35 Goal : Solve the non-linear inverse problem Reconstruct
More informationarxiv: v1 [math.na] 16 Apr 2019
A direct solver for the phase retrieval problem in ptychographic imaging arxiv:1904.07940v1 [math.na] 16 Apr 2019 Nada Sissouno a,b,, Florian Boßmann d, Frank Filbir b, Mark Iwen e, Maik Kahnt c, Rayan
More informationPHASE RETRIEVAL FROM STFT MEASUREMENTS VIA NON-CONVEX OPTIMIZATION. Tamir Bendory and Yonina C. Eldar, Fellow IEEE
PHASE RETRIEVAL FROM STFT MEASUREMETS VIA O-COVEX OPTIMIZATIO Tamir Bendory and Yonina C Eldar, Fellow IEEE Department of Electrical Engineering, Technion - Israel Institute of Technology, Haifa, Israel
More informationProvable Alternating Minimization Methods for Non-convex Optimization
Provable Alternating Minimization Methods for Non-convex Optimization Prateek Jain Microsoft Research, India Joint work with Praneeth Netrapalli, Sujay Sanghavi, Alekh Agarwal, Animashree Anandkumar, Rashish
More informationarxiv: v1 [cs.it] 26 Oct 2015
Phase Retrieval: An Overview of Recent Developments Kishore Jaganathan Yonina C. Eldar Babak Hassibi Department of Electrical Engineering, Caltech Department of Electrical Engineering, Technion, Israel
More informationStatistical Issues in Searches: Photon Science Response. Rebecca Willett, Duke University
Statistical Issues in Searches: Photon Science Response Rebecca Willett, Duke University 1 / 34 Photon science seen this week Applications tomography ptychography nanochrystalography coherent diffraction
More informationThe Iterative and Regularized Least Squares Algorithm for Phase Retrieval
The Iterative and Regularized Least Squares Algorithm for Phase Retrieval Radu Balan Naveed Haghani Department of Mathematics, AMSC, CSCAMM and NWC University of Maryland, College Park, MD April 17, 2016
More informationSolving Corrupted Quadratic Equations, Provably
Solving Corrupted Quadratic Equations, Provably Yuejie Chi London Workshop on Sparse Signal Processing September 206 Acknowledgement Joint work with Yuanxin Li (OSU), Huishuai Zhuang (Syracuse) and Yingbin
More informationSelf-Calibration and Biconvex Compressive Sensing
Self-Calibration and Biconvex Compressive Sensing Shuyang Ling Department of Mathematics, UC Davis July 12, 2017 Shuyang Ling (UC Davis) SIAM Annual Meeting, 2017, Pittsburgh July 12, 2017 1 / 22 Acknowledgements
More informationRapid, Robust, and Reliable Blind Deconvolution via Nonconvex Optimization
Rapid, Robust, and Reliable Blind Deconvolution via Nonconvex Optimization Shuyang Ling Department of Mathematics, UC Davis Oct.18th, 2016 Shuyang Ling (UC Davis) 16w5136, Oaxaca, Mexico Oct.18th, 2016
More informationGlobally Convergent Levenberg-Marquardt Method For Phase Retrieval
Globally Convergent Levenberg-Marquardt Method For Phase Retrieval Zaiwen Wen Beijing International Center For Mathematical Research Peking University Thanks: Chao Ma, Xin Liu 1/38 Outline 1 Introduction
More informationA new method on deterministic construction of the measurement matrix in compressed sensing
A new method on deterministic construction of the measurement matrix in compressed sensing Qun Mo 1 arxiv:1503.01250v1 [cs.it] 4 Mar 2015 Abstract Construction on the measurement matrix A is a central
More informationStrengthened Sobolev inequalities for a random subspace of functions
Strengthened Sobolev inequalities for a random subspace of functions Rachel Ward University of Texas at Austin April 2013 2 Discrete Sobolev inequalities Proposition (Sobolev inequality for discrete images)
More informationBM3D-prGAMP: Compressive Phase Retrieval Based on BM3D Denoising
BM3D-prGAMP: Compressive Phase Retrieval Based on BM3D Denoising Chris Metzler, Richard Baraniuk Rice University Arian Maleki Columbia University Phase Retrieval Applications: Crystallography Microscopy
More informationPhase Retrieval using the Iterative Regularized Least-Squares Algorithm
Phase Retrieval using the Iterative Regularized Least-Squares Algorithm Radu Balan Department of Mathematics, AMSC, CSCAMM and NWC University of Maryland, College Park, MD August 16, 2017 Workshop on Phaseless
More informationPHASE RETRIEVAL OF SPARSE SIGNALS FROM MAGNITUDE INFORMATION. A Thesis MELTEM APAYDIN
PHASE RETRIEVAL OF SPARSE SIGNALS FROM MAGNITUDE INFORMATION A Thesis by MELTEM APAYDIN Submitted to the Office of Graduate and Professional Studies of Texas A&M University in partial fulfillment of the
More informationNorms and embeddings of classes of positive semidefinite matrices
Norms and embeddings of classes of positive semidefinite matrices Radu Balan Department of Mathematics, Center for Scientific Computation and Mathematical Modeling and the Norbert Wiener Center for Applied
More informationSparse Phase Retrieval Using Partial Nested Fourier Samplers
Sparse Phase Retrieval Using Partial Nested Fourier Samplers Heng Qiao Piya Pal University of Maryland qiaoheng@umd.edu November 14, 2015 Heng Qiao, Piya Pal (UMD) Phase Retrieval with PNFS November 14,
More informationCapacity-Approaching PhaseCode for Low-Complexity Compressive Phase Retrieval
Capacity-Approaching PhaseCode for Low-Complexity Compressive Phase Retrieval Ramtin Pedarsani, Kangwook Lee, and Kannan Ramchandran Dept of Electrical Engineering and Computer Sciences University of California,
More informationNon-convex Robust PCA: Provable Bounds
Non-convex Robust PCA: Provable Bounds Anima Anandkumar U.C. Irvine Joint work with Praneeth Netrapalli, U.N. Niranjan, Prateek Jain and Sujay Sanghavi. Learning with Big Data High Dimensional Regime Missing
More informationof Orthogonal Matching Pursuit
A Sharp Restricted Isometry Constant Bound of Orthogonal Matching Pursuit Qun Mo arxiv:50.0708v [cs.it] 8 Jan 205 Abstract We shall show that if the restricted isometry constant (RIC) δ s+ (A) of the measurement
More informationGauge optimization and duality
1 / 54 Gauge optimization and duality Junfeng Yang Department of Mathematics Nanjing University Joint with Shiqian Ma, CUHK September, 2015 2 / 54 Outline Introduction Duality Lagrange duality Fenchel
More informationFourier phase retrieval with random phase illumination
Fourier phase retrieval with random phase illumination Wenjing Liao School of Mathematics Georgia Institute of Technology Albert Fannjiang Department of Mathematics, UC Davis. IMA August 5, 27 Classical
More informationCombining Sparsity with Physically-Meaningful Constraints in Sparse Parameter Estimation
UIUC CSL Mar. 24 Combining Sparsity with Physically-Meaningful Constraints in Sparse Parameter Estimation Yuejie Chi Department of ECE and BMI Ohio State University Joint work with Yuxin Chen (Stanford).
More informationAn algebraic perspective on integer sparse recovery
An algebraic perspective on integer sparse recovery Lenny Fukshansky Claremont McKenna College (joint work with Deanna Needell and Benny Sudakov) Combinatorics Seminar USC October 31, 2018 From Wikipedia:
More informationLarge-Scale L1-Related Minimization in Compressive Sensing and Beyond
Large-Scale L1-Related Minimization in Compressive Sensing and Beyond Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, U.S.A. Arizona State University March
More informationCoherent imaging without phases
Coherent imaging without phases Miguel Moscoso Joint work with Alexei Novikov Chrysoula Tsogka and George Papanicolaou Waves and Imaging in Random Media, September 2017 Outline 1 The phase retrieval problem
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationA Geometric Analysis of Phase Retrieval
A Geometric Analysis of Phase Retrieval Ju Sun, Qing Qu, John Wright js4038, qq2105, jw2966}@columbiaedu Dept of Electrical Engineering, Columbia University, New York NY 10027, USA Abstract Given measurements
More informationIntroduction to Compressed Sensing
Introduction to Compressed Sensing Alejandro Parada, Gonzalo Arce University of Delaware August 25, 2016 Motivation: Classical Sampling 1 Motivation: Classical Sampling Issues Some applications Radar Spectral
More informationSemidefinite Programming Based Preconditioning for More Robust Near-Separable Nonnegative Matrix Factorization
Semidefinite Programming Based Preconditioning for More Robust Near-Separable Nonnegative Matrix Factorization Nicolas Gillis nicolas.gillis@umons.ac.be https://sites.google.com/site/nicolasgillis/ Department
More informationReshaped Wirtinger Flow for Solving Quadratic System of Equations
Reshaped Wirtinger Flow for Solving Quadratic System of Equations Huishuai Zhang Department of EECS Syracuse University Syracuse, NY 344 hzhan3@syr.edu Yingbin Liang Department of EECS Syracuse University
More informationSample Optimal Fourier Sampling in Any Constant Dimension
1 / 26 Sample Optimal Fourier Sampling in Any Constant Dimension Piotr Indyk Michael Kapralov MIT IBM Watson October 21, 2014 Fourier Transform and Sparsity Discrete Fourier Transform Given x C n, compute
More informationSuper-resolution by means of Beurling minimal extrapolation
Super-resolution by means of Beurling minimal extrapolation Norbert Wiener Center Department of Mathematics University of Maryland, College Park http://www.norbertwiener.umd.edu Acknowledgements ARO W911NF-16-1-0008
More informationRobust Phase Retrieval Algorithm for Time-Frequency Structured Measurements
Robust Phase Retrieval Algorithm for Time-Frequency Structured easurements Götz E. Pfander Palina Salanevich arxiv:1611.0540v1 [math.fa] Oct 016 July 17, 018 Abstract We address the problem of signal reconstruction
More informationInformation and Resolution
Information and Resolution (Invited Paper) Albert Fannjiang Department of Mathematics UC Davis, CA 95616-8633. fannjiang@math.ucdavis.edu. Abstract The issue of resolution with complex-field measurement
More informationSuper-resolution via Convex Programming
Super-resolution via Convex Programming Carlos Fernandez-Granda (Joint work with Emmanuel Candès) Structure and Randomness in System Identication and Learning, IPAM 1/17/2013 1/17/2013 1 / 44 Index 1 Motivation
More informationMismatch and Resolution in CS
Mismatch and Resolution in CS Albert Fannjiang, UC Davis Coworker: Wenjing Liao, UC Davis SPIE Optics + Photonics, San Diego, August 2-25, 20 Mismatch: gridding error Band exclusion Local optimization
More informationDISTRIBUTION A: Distribution approved for public release.
AFRL-AFOSR-VA-TR-2016-0105 BRI) Fault Tolerant Paradigms BENJAMIN ONG MICHIGAN STATE UNIV EAST LANSING 02/26/2016 Final Report Air Force Research Laboratory AF Office Of Scientific Research (AFOSR)/ RTA2
More informationFrequency Subspace Amplitude Flow for Phase Retrieval
Research Article Journal of the Optical Society of America A 1 Frequency Subspace Amplitude Flow for Phase Retrieval ZHUN WEI 1, WEN CHEN 2,3, AND XUDONG CHEN 1,* 1 Department of Electrical and Computer
More informationPhase Retrieval using Alternating Minimization
Phase Retrieval using Alternating Minimization Praneeth Netrapalli Department of ECE The University of Texas at Austin Austin, TX 787 praneethn@utexas.edu Prateek Jain Microsoft Research India Bangalore,
More informationLow-rank Matrix Completion with Noisy Observations: a Quantitative Comparison
Low-rank Matrix Completion with Noisy Observations: a Quantitative Comparison Raghunandan H. Keshavan, Andrea Montanari and Sewoong Oh Electrical Engineering and Statistics Department Stanford University,
More informationRecovering overcomplete sparse representations from structured sensing
Recovering overcomplete sparse representations from structured sensing Deanna Needell Claremont McKenna College Feb. 2015 Support: Alfred P. Sloan Foundation and NSF CAREER #1348721. Joint work with Felix
More informationSensing systems limited by constraints: physical size, time, cost, energy
Rebecca Willett Sensing systems limited by constraints: physical size, time, cost, energy Reduce the number of measurements needed for reconstruction Higher accuracy data subject to constraints Original
More informationRecovery of Sparse Signals from Noisy Measurements Using an l p -Regularized Least-Squares Algorithm
Recovery of Sparse Signals from Noisy Measurements Using an l p -Regularized Least-Squares Algorithm J. K. Pant, W.-S. Lu, and A. Antoniou University of Victoria August 25, 2011 Compressive Sensing 1 University
More informationDictionary Learning Using Tensor Methods
Dictionary Learning Using Tensor Methods Anima Anandkumar U.C. Irvine Joint work with Rong Ge, Majid Janzamin and Furong Huang. Feature learning as cornerstone of ML ML Practice Feature learning as cornerstone
More informationStatistical Geometry Processing Winter Semester 2011/2012
Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian
More informationSparse Subspace Clustering
Sparse Subspace Clustering Based on Sparse Subspace Clustering: Algorithm, Theory, and Applications by Elhamifar and Vidal (2013) Alex Gutierrez CSCI 8314 March 2, 2017 Outline 1 Motivation and Background
More informationRapidly Computing Sparse Chebyshev and Legendre Coefficient Expansions via SFTs
Rapidly Computing Sparse Chebyshev and Legendre Coefficient Expansions via SFTs Mark Iwen Michigan State University October 18, 2014 Work with Janice (Xianfeng) Hu Graduating in May! M.A. Iwen (MSU) Fast
More informationDeterministic Compressed Sensing for Efficient Image Reconstruction
Chirp for Efficient Image University of Idaho Pacific Northwest Numerical Analysis Seminar 27 October 2012 Chirp Collaborators: Connie Ni (HRL Lab, Malibu, CA) Prasun Mahanti (Arizona State University)
More informationNumerical Approximation Methods for Non-Uniform Fourier Data
Numerical Approximation Methods for Non-Uniform Fourier Data Aditya Viswanathan aditya@math.msu.edu 2014 Joint Mathematics Meetings January 18 2014 0 / 16 Joint work with Anne Gelb (Arizona State) Guohui
More informationAPPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES
APPROXIMATING THE INVERSE FRAME OPERATOR FROM LOCALIZED FRAMES GUOHUI SONG AND ANNE GELB Abstract. This investigation seeks to establish the practicality of numerical frame approximations. Specifically,
More informationsparse and low-rank tensor recovery Cubic-Sketching
Sparse and Low-Ran Tensor Recovery via Cubic-Setching Guang Cheng Department of Statistics Purdue University www.science.purdue.edu/bigdata CCAM@Purdue Math Oct. 27, 2017 Joint wor with Botao Hao and Anru
More informationCoSaMP. Iterative signal recovery from incomplete and inaccurate samples. Joel A. Tropp
CoSaMP Iterative signal recovery from incomplete and inaccurate samples Joel A. Tropp Applied & Computational Mathematics California Institute of Technology jtropp@acm.caltech.edu Joint with D. Needell
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationNew Coherence and RIP Analysis for Weak. Orthogonal Matching Pursuit
New Coherence and RIP Analysis for Wea 1 Orthogonal Matching Pursuit Mingrui Yang, Member, IEEE, and Fran de Hoog arxiv:1405.3354v1 [cs.it] 14 May 2014 Abstract In this paper we define a new coherence
More informationFourier Phase Retrieval: Uniqueness and Algorithms
Fourier Phase Retrieval: Uniqueness and Algorithms Tamir Bendory 1, Robert Beinert 2 and Yonina C. Eldar 3 arxiv:1705.09590v3 [cs.it] 6 Nov 2017 1 The Program in Applied and Computational Mathematics,
More informationSparse Legendre expansions via l 1 minimization
Sparse Legendre expansions via l 1 minimization Rachel Ward, Courant Institute, NYU Joint work with Holger Rauhut, Hausdorff Center for Mathematics, Bonn, Germany. June 8, 2010 Outline Sparse recovery
More informationRecent Developments in Compressed Sensing
Recent Developments in Compressed Sensing M. Vidyasagar Distinguished Professor, IIT Hyderabad m.vidyasagar@iith.ac.in, www.iith.ac.in/ m vidyasagar/ ISL Seminar, Stanford University, 19 April 2018 Outline
More informationConstructing Approximation Kernels for Non-Harmonic Fourier Data
Constructing Approximation Kernels for Non-Harmonic Fourier Data Aditya Viswanathan aditya.v@caltech.edu California Institute of Technology SIAM Annual Meeting 2013 July 10 2013 0 / 19 Joint work with
More informationNonconvex Methods for Phase Retrieval
Nonconvex Methods for Phase Retrieval Vince Monardo Carnegie Mellon University March 28, 2018 Vince Monardo (CMU) 18.898 Midterm March 28, 2018 1 / 43 Paper Choice Main reference paper: Solving Random
More informationRapid, Robust, and Reliable Blind Deconvolution via Nonconvex Optimization
Rapid, Robust, and Reliable Blind Deconvolution via Nonconvex Optimization Xiaodong Li, Shuyang Ling, Thomas Strohmer, and Ke Wei June 15, 016 Abstract We study the question of reconstructing two signals
More informationStatistics of the Stability Bounds in the Phase Retrieval Problem
Statistics of the Stability Bounds in the Phase Retrieval Problem Radu Balan Department of Mathematics, AMSC, CSCAMM and NWC University of Maryland, College Park, MD January 7, 2017 AMS Special Session
More informationCovariance Sketching via Quadratic Sampling
Covariance Sketching via Quadratic Sampling Yuejie Chi Department of ECE and BMI The Ohio State University Tsinghua University June 2015 Page 1 Acknowledgement Thanks to my academic collaborators on some
More informationGREEDY SIGNAL RECOVERY REVIEW
GREEDY SIGNAL RECOVERY REVIEW DEANNA NEEDELL, JOEL A. TROPP, ROMAN VERSHYNIN Abstract. The two major approaches to sparse recovery are L 1-minimization and greedy methods. Recently, Needell and Vershynin
More informationGossip algorithms for solving Laplacian systems
Gossip algorithms for solving Laplacian systems Anastasios Zouzias University of Toronto joint work with Nikolaos Freris (EPFL) Based on : 1.Fast Distributed Smoothing for Clock Synchronization (CDC 1).Randomized
More informationPhase Retrieval with Random Illumination
Phase Retrieval with Random Illumination Albert Fannjiang, UC Davis Collaborator: Wenjing Liao NCTU, July 212 1 Phase retrieval Problem Reconstruct the object f from the Fourier magnitude Φf. Why do we
More informationCompressibility of Infinite Sequences and its Interplay with Compressed Sensing Recovery
Compressibility of Infinite Sequences and its Interplay with Compressed Sensing Recovery Jorge F. Silva and Eduardo Pavez Department of Electrical Engineering Information and Decision Systems Group Universidad
More informationAn Adaptive Sublinear Time Block Sparse Fourier Transform
An Adaptive Sublinear Time Block Sparse Fourier Transform Volkan Cevher Michael Kapralov Jonathan Scarlett Amir Zandieh EPFL February 8th 217 Given x C N, compute the Discrete Fourier Transform (DFT) of
More informationSpectral Processing. Misha Kazhdan
Spectral Processing Misha Kazhdan [Taubin, 1995] A Signal Processing Approach to Fair Surface Design [Desbrun, et al., 1999] Implicit Fairing of Arbitrary Meshes [Vallet and Levy, 2008] Spectral Geometry
More informationMultiresolution Analysis
Multiresolution Analysis DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Frames Short-time Fourier transform
More informationTensor Low-Rank Completion and Invariance of the Tucker Core
Tensor Low-Rank Completion and Invariance of the Tucker Core Shuzhong Zhang Department of Industrial & Systems Engineering University of Minnesota zhangs@umn.edu Joint work with Bo JIANG, Shiqian MA, and
More informationSparse and Low Rank Recovery via Null Space Properties
Sparse and Low Rank Recovery via Null Space Properties Holger Rauhut Lehrstuhl C für Mathematik (Analysis), RWTH Aachen Convexity, probability and discrete structures, a geometric viewpoint Marne-la-Vallée,
More informationA New Estimate of Restricted Isometry Constants for Sparse Solutions
A New Estimate of Restricted Isometry Constants for Sparse Solutions Ming-Jun Lai and Louis Y. Liu January 12, 211 Abstract We show that as long as the restricted isometry constant δ 2k < 1/2, there exist
More informationSparsity Regularization
Sparsity Regularization Bangti Jin Course Inverse Problems & Imaging 1 / 41 Outline 1 Motivation: sparsity? 2 Mathematical preliminaries 3 l 1 solvers 2 / 41 problem setup finite-dimensional formulation
More informationCS 229r: Algorithms for Big Data Fall Lecture 19 Nov 5
CS 229r: Algorithms for Big Data Fall 215 Prof. Jelani Nelson Lecture 19 Nov 5 Scribe: Abdul Wasay 1 Overview In the last lecture, we started discussing the problem of compressed sensing where we are given
More informationCompressed Sensing and Robust Recovery of Low Rank Matrices
Compressed Sensing and Robust Recovery of Low Rank Matrices M. Fazel, E. Candès, B. Recht, P. Parrilo Electrical Engineering, University of Washington Applied and Computational Mathematics Dept., Caltech
More informationTowards a Mathematical Theory of Super-resolution
Towards a Mathematical Theory of Super-resolution Carlos Fernandez-Granda www.stanford.edu/~cfgranda/ Information Theory Forum, Information Systems Laboratory, Stanford 10/18/2013 Acknowledgements This
More informationOn the recovery of measures without separation conditions
Norbert Wiener Center Department of Mathematics University of Maryland, College Park http://www.norbertwiener.umd.edu Applied and Computational Mathematics Seminar Georgia Institute of Technology October
More informationFast, Sample-Efficient Algorithms for Structured Phase Retrieval
Fast, Sample-Efficient Algorithms for Structured Phase Retrieval Gauri jagatap Electrical and Computer Engineering Iowa State University Chinmay Hegde Electrical and Computer Engineering Iowa State University
More informationRobust Sparse Recovery via Non-Convex Optimization
Robust Sparse Recovery via Non-Convex Optimization Laming Chen and Yuantao Gu Department of Electronic Engineering, Tsinghua University Homepage: http://gu.ee.tsinghua.edu.cn/ Email: gyt@tsinghua.edu.cn
More informationMitigating the Effect of Noise in Iterative Projection Phase Retrieval
Mitigating the Effect of Noise in Iterative Projection Phase Retrieval & David Hyland Imaging and Applied Optics: Optics and Photonics Congress Introduction Filtering Results Conclusion Outline Introduction
More information8.1 Concentration inequality for Gaussian random matrix (cont d)
MGMT 69: Topics in High-dimensional Data Analysis Falll 26 Lecture 8: Spectral clustering and Laplacian matrices Lecturer: Jiaming Xu Scribe: Hyun-Ju Oh and Taotao He, October 4, 26 Outline Concentration
More informationInterpolation via weighted l 1 minimization
Interpolation via weighted l 1 minimization Rachel Ward University of Texas at Austin December 12, 2014 Joint work with Holger Rauhut (Aachen University) Function interpolation Given a function f : D C
More informationRich Tomography. Bill Lionheart, School of Mathematics, University of Manchester and DTU Compute. July 2014
Rich Tomography Bill Lionheart, School of Mathematics, University of Manchester and DTU Compute July 2014 What do we mean by Rich Tomography? Conventional tomography reconstructs one scalar image from
More informationSparse analysis Lecture VII: Combining geometry and combinatorics, sparse matrices for sparse signal recovery
Sparse analysis Lecture VII: Combining geometry and combinatorics, sparse matrices for sparse signal recovery Anna C. Gilbert Department of Mathematics University of Michigan Sparse signal recovery measurements:
More informationA Quantum Interior Point Method for LPs and SDPs
A Quantum Interior Point Method for LPs and SDPs Iordanis Kerenidis 1 Anupam Prakash 1 1 CNRS, IRIF, Université Paris Diderot, Paris, France. September 26, 2018 Semi Definite Programs A Semidefinite Program
More informationNoisy Signal Recovery via Iterative Reweighted L1-Minimization
Noisy Signal Recovery via Iterative Reweighted L1-Minimization Deanna Needell UC Davis / Stanford University Asilomar SSC, November 2009 Problem Background Setup 1 Suppose x is an unknown signal in R d.
More informationSolution Recovery via L1 minimization: What are possible and Why?
Solution Recovery via L1 minimization: What are possible and Why? Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, U.S.A. Eighth US-Mexico Workshop on Optimization
More informationOrthogonal Nonnegative Matrix Factorization: Multiplicative Updates on Stiefel Manifolds
Orthogonal Nonnegative Matrix Factorization: Multiplicative Updates on Stiefel Manifolds Jiho Yoo and Seungjin Choi Department of Computer Science Pohang University of Science and Technology San 31 Hyoja-dong,
More informationarxiv: v2 [stat.ml] 12 Jun 2015
Phase Retrieval using Alternating Minimization Praneeth Netrapalli Prateek Jain Sujay Sanghavi June 15, 015 arxiv:1306.0160v [stat.ml] 1 Jun 015 Abstract Phase retrieval problems involve solving linear
More informationReconstruction of Block-Sparse Signals by Using an l 2/p -Regularized Least-Squares Algorithm
Reconstruction of Block-Sparse Signals by Using an l 2/p -Regularized Least-Squares Algorithm Jeevan K. Pant, Wu-Sheng Lu, and Andreas Antoniou University of Victoria May 21, 2012 Compressive Sensing 1/23
More informationECE 8201: Low-dimensional Signal Models for High-dimensional Data Analysis
ECE 8201: Low-dimensional Signal Models for High-dimensional Data Analysis Lecture 7: Matrix completion Yuejie Chi The Ohio State University Page 1 Reference Guaranteed Minimum-Rank Solutions of Linear
More information