Sparse phase retrieval via Phaseliftoff

The aim of sparse phase retrieval is to recover a $k$-sparse signal $\mathbf{x}0\in \mathbb{C}^{d}$ from quadratic measurements $|\langle \mathbf{a}i,\mathbf{x}0\rangle|^2$ where $\mathbf{a}i\in \mathbb{C}^d, i=1,\ldots,m$. Noting $|\langle \mathbf{a}i,\mathbf{x}0\rangle|^{2={\text{Tr}}(AiX0)$} with $Ai=\mathbf{a}i\mathbf{a}i^*\in \mathbb{C}^{d\times d}, X0=\mathbf{x}0\mathbf{x}0^*\in \mathbb{C}^{d\times d}$, one can recast sparse phase retrieval as a problem of recovering a rank-one sparse matrix from linear measurements. Yin and Xin introduced PhaseLiftOff which presents a proxy of rank-one condition via the difference of trace and Frobenius norm. By adding sparsity penalty to PhaseLiftOff, in this paper, we present a novel model to recover sparse signals from quadratic measurements. Theoretical analysis shows that the solution to our model provides the stable recovery of $\mathbf{x}_0$ under almost optimal sampling complexity $m=O(k\log(d/k))$. The computation of our model is carried out by the difference of convex function algorithm (DCA). Numerical experiments demonstrate that our algorithm outperforms other state-of-the-art algorithms used for solving sparse phase retrieval.