Solving Random Systems of Quadratic Equations with Tanh Wirtinger Flow

Solving quadratic systems of equations in n variables and m measurements of the form yi = |aTi x|2 , i = 1, ..., m and x ∈ Rn , which is also known as phase retrieval, is a hard nonconvex problem. In the case of standard Gaussian measurement vectors, the wirtinger flow algorithm Chen and Candes (2015) is an efficient solution. In this paper, we proposed a new form of wirtinger flow and a new spectral initialization method based on this new algorithm. We proved that the new wirtinger flow and initialization method achieve linear sample and computational complexities. We further extended the new phasing algorithm by combining it with other existing methods. Finally, we demonstrated the effectiveness of our new method in the low data to parameter ratio settings where the number of measurements which is less than information-theoretic limit, namely, m < 2n, via numerical tests. For instance, our method can solve the quadratic systems of equations with gaussian measurement vector with probability ≥ 97% when m/n = 1.7 and n = 1000, and with probability ≈ 60% when m/n = 1.5 and n = 1000.