Abstract
Fourier phasing is the problem of retrieving Fourier phase information from Fourier intensity data. The error-reduction (ER) algorithm consists of two projections on the subspaces generated by the Fourier magnitude constraint and the object-domain constraint. The random phase illumination (RPI) and the real-valued constraint on the object significantly reduce the complexity of the intersection of the two subspaces. In this paper, we study how to approximate the projection of the starting point onto the subspace generated by the Fourier magnitude constraint by its projection on the tangent plane and estimate the approximation error by orthogonal decompositions. Moreover, we prove that the local geometric convergence rate of the ER algorithm is less than one almost surely and can be characterized as the cosine of the angle between the two projection spaces. A theoretical estimate of the convergence rate is derived and validated by some numerical experiments.
Funding statement: This work was supported by the National Center for Theoretical Sciences (NCTS, Taiwan) and the National Science and Technology Council under grant numbers 110-2628-M-006-003-MY3 and 113-2628-M-006-002-MY3.
A Derivation of (3.22)
We first observe that
The minimum point
which implies (3.22).
B Proof of Theorem 1
Denote
and
We have the following orthogonal decompositions for
Decomposition for
where
Decomposition for
for some
The orthogonality of
where the decomposition
Lemma 1.
For any
The proof of Lemma 1 can be found in Appendix C.
Because ($C_{DEF}$) is the necessary condition for
By the method of Lagrange multipliers, the maximum in (B.2)
is attained when
where ($C_{DE}$) is a constraint obtained by substituting

Constraint ($C_{DE}$) on D and E (
The constraint ($C_{DE}$) is depicted in Figure 2, which implies that the right-hand side of (B.3) can be rewritten as
In the following,
we derive another constraint on
where
We apply (3.21) and (3.22) to write down the vectorization of
where
Because
we have
Therefore,
We start to estimate
Because
For every
and define the difference
The notation above helps us get the following estimation:
where the inequality above holds for any
By the inequality
(B.10) implies that
We claim that
whose proof can be found in Appendix D. By substituting (B.12) into (B.11),
By denoting
By combining the estimation above with (B.4) and (B.5), we obtain another constraint as follows:
The entire constraints, including ($C_{DE}$) and (${C_{E}}$), on
Constraints ($C_{DE}$) and (${C_{E}}$) on D and E (


Finally, we consider two cases as follows:
Case 1:
where
where the last inequality follows from
By (B.14), we obtain
The statement of Theorem 1 follows by substituting (B.15) into (B.13).
Case 2:
from the constraint ($C_{DE}$). Note that (B.16) implies the same conclusion as (B.15) by setting
C Proof of Lemma 1
First of all, we show that Lemma 1 follows from the estimate
for any
where the last inequality is implied by (C.1).
On the other hand, by the orthogonal decompositions of
By combining (C.2) and (C.3), we get the desired result
Now we start to prove (C.1). Denote
where Φ is a matrix consisting of the column vectors
By denoting
where
Hence, the term on the right-hand side of (C.4) can be calculated explicitly as follows:
where the minimum is taken over all diagonal matrices in
which implies that
Applying the inequality
D Proof of (B.12)
Because
Let
be an orthonormal matrix such that
Denote
where
By applying
the inequality (D.4)
with
Because
Because
which implies that
In the same way, we have
By substituting (D.7) and (D.8) into (D.6), we get
The claim (D.1) follows by (D.9) and noticing that
Acknowledgements
The author gratefully acknowledges the anonymous reviewers for their constructive comments and critiques, which significantly enhanced this manuscript. The author also wishes to thank Albert Fannjiang and Pengwen Chen for their valuable discussions and insightful feedback.
References
[1] F. Andersson and M. Carlsson, Alternating projections on nontangential manifolds, Constr. Approx. 38 (2013), no. 3, 489–525. 10.1007/s00365-013-9213-3Search in Google Scholar
[2] A. S. Bandeira, Y. Chen and D. G. Mixon, Phase retrieval from power spectra of masked signals, Inf. Inference 3 (2014), no. 2, 83–102. 10.1093/imaiai/iau002Search in Google Scholar
[3] H. H. Bauschke, J. Y. Bello Cruz, T. T. A. Nghia, H. M. Phan and X. Wang, The rate of linear convergence of the Douglas–Rachford algorithm for subspaces is the cosine of the Friedrichs angle, J. Approx. Theory 185 (2014), 63–79. 10.1016/j.jat.2014.06.002Search in Google Scholar
[4] H. H. Bauschke, P. L. Combettes and D. R. Luke, Phase retrieval, error reduction algorithm, and Fienup variants: A view from convex optimization, J. Opt. Soc. Amer. A 19 (2002), no. 7, 1334–1345. 10.1364/JOSAA.19.001334Search in Google Scholar
[5] E. J. Candès, X. Li and M. Soltanolkotabi, Phase retrieval from coded diffraction patterns, Appl. Comput. Harmon. Anal. 39 (2015), no. 2, 277–299. 10.1016/j.acha.2014.09.004Search in Google Scholar
[6] E. J. Candès, X. Li and M. Soltanolkotabi, Phase retrieval via Wirtinger flow: Theory and algorithms, IEEE Trans. Inform. Theory 61 (2015), no. 4, 1985–2007. 10.1109/TIT.2015.2399924Search in Google Scholar
[7] P. Chen and A. Fannjiang, Fourier phase retrieval with a single mask by Douglas–Rachford algorithms, Appl. Comput. Harmon. Anal. 44 (2018), no. 3, 665–699. 10.1016/j.acha.2016.07.003Search in Google Scholar PubMed PubMed Central
[8] P. Chen, A. Fannjiang and G.-R. Liu, Phase retrieval with one or two diffraction patterns by alternating projections with the null initialization, J. Fourier Anal. Appl. 24 (2018), no. 3, 719–758. 10.1007/s00041-017-9536-8Search in Google Scholar
[9] Y. Chen and E. J. Candès, Solving random quadratic systems of equations is nearly as easy as solving linear systems, Advances in Neural Information Processing Systems, Curran Associates, Red Hook (2015), 739–747. Search in Google Scholar
[10] F. H. Clarke, Y. S. Ledyaev, R. J. Stern and P. R. Wolenski, Nonsmooth Analysis and Control Theory, Grad. Texts in Math. 178, Springer, New York, 2008. Search in Google Scholar
[11] J. C. Dainty and J. R. Fienup, Phase retrieval and image reconstruction for astronomy, Image Recovery: Theory and Application, Academic Press, San Diego (1987), 231–275. Search in Google Scholar
[12] F. Deutsch, Best Approximation in Inner Product Spaces, CMS Books Math./Ouvrages Math. SMC 7, Springer, New York, 2001. 10.1007/978-1-4684-9298-9Search in Google Scholar
[13] Y. C. Eldar and S. Mendelson, Phase retrieval: Stability and recovery guarantees, Appl. Comput. Harmon. Anal. 36 (2014), no. 3, 473–494. 10.1016/j.acha.2013.08.003Search in Google Scholar
[14] A. Fannjiang, Absolute uniqueness of phase retrieval with random illumination, Inverse Problems 28 (2012), no. 7, Article ID 075008. 10.1088/0266-5611/28/7/075008Search in Google Scholar
[15] A. Fannjiang and W. Liao, Phase retrieval with random phase illumination, JOSA A 29 (2012), no. 9, 1847–1859. 10.1364/JOSAA.29.001847Search in Google Scholar PubMed
[16] A. Fannjiang and W. Liao, Fourier phasing with phase-uncertain mask, Inverse Problems 29 (2013), no. 12, Article ID 125001. 10.1088/0266-5611/29/12/125001Search in Google Scholar
[17] M. H. Hayes, The reconstruction of a multidimensional sequence from the phase or magnitude of its Fourier transform, IEEE Trans. Acoust. Speech Signal Process. 30 (1982), no. 2, 140–154. 10.1109/TASSP.1982.1163863Search in Google Scholar
[18] X. Jiang, S. Rajan and X. Liu, Wirtinger flow method with optimal stepsize for phase retrieval, IEEE Signal Process. Lett. 23 (2016), no. 11, 1627–1631. 10.1109/LSP.2016.2611940Search in Google Scholar
[19] S. Kayalar and H. L. Weinert, Error bounds for the method of alternating projections, Math. Control Signals Systems 1 (1988), no. 1, 43–59. 10.1007/BF02551235Search in Google Scholar
[20] R. Kolte and A. Ozgur, Phase retrieval via incremental truncated Wirtinger flow, preprint (2016), https://arxiv.org/abs/1606.03196. Search in Google Scholar
[21] A. S. Lewis, D. R. Luke and J. Malick, Local linear convergence for alternating and averaged nonconvex projections, Found. Comput. Math. 9 (2009), no. 4, 485–513. 10.1007/s10208-008-9036-ySearch in Google Scholar
[22] A. S. Lewis and J. Malick, Alternating projections on manifolds, Math. Oper. Res. 33 (2008), no. 1, 216–234. 10.1287/moor.1070.0291Search in Google Scholar
[23] J. Miao, P. Charalambous, J. Kirz and D. Sayre, Extending the methodology of X-ray crystallography to allow imaging of micrometre-sized non-crystalline specimens, Nature 400 (1999), no. 6742, 342–344. 10.1038/22498Search in Google Scholar
[24] R. P. Millane, Phase retrieval in crystallography and optics, JOSA A 7 (1990), no. 3, 394–411. 10.1364/JOSAA.7.000394Search in Google Scholar
[25] B. S. Mordukhovich, Variational Analysis and Generalized Differentiation. I, Grundlehren Math. Wiss. 330, Springer, Berlin, 2006. 10.1007/3-540-31246-3Search in Google Scholar
[26] P. Netrapalli, P. Jain and S. Sanghavi, Phase retrieval using alternating minimization, Advances in Neural Information Processing Systems 26, Curran Associates, Red Hook (2013), 2796–2804. Search in Google Scholar
[27] D. Noll and A. Rondepierre, On local convergence of the method of alternating projections, Found. Comput. Math. 16 (2016), no. 2, 425–455. 10.1007/s10208-015-9253-0Search in Google Scholar
[28] V. Pohl, F. Yang and H. Boche, Phase retrieval from low-rate samples, Sampl. Theory Signal Image Process. 14 (2015), no. 1, 71–99. 10.1007/BF03549588Search in Google Scholar
[29] R. Tyrrell Rockafellar and R. J.-B. Wets, Variational Analysis, Grundlehren Math. Wiss. 317, Springer, Berlin, 2009. Search in Google Scholar
[30] G. Wang, G. B. Giannakis and Y. C. Eldar, Solving systems of random quadratic equations via truncated amplitude flow, IEEE Trans. Inform. Theory 64 (2018), no. 2, 773–794. 10.1109/TIT.2017.2756858Search in Google Scholar
[31] K. Wei, Solving systems of phaseless equations via Kaczmarz methods: A proof of concept study, Inverse Problems 31 (2015), no. 12, Article ID 125008. 10.1088/0266-5611/31/12/125008Search in Google Scholar
[32] Z. Wen, C. Yang, X. Liu and S. Marchesini, Alternating direction methods for classical and ptychographic phase retrieval, Inverse Problems 28 (2012), no. 11, Article ID 115010. 10.1088/0266-5611/28/11/115010Search in Google Scholar
[33] P. Yin and J. Xin, PhaseLiftOff: An accurate and stable phase retrieval method based on difference of trace and Frobenius norms, Commun. Math. Sci. 13 (2015), no. 4, 1033–1049. 10.4310/CMS.2015.v13.n4.a10Search in Google Scholar
[34] H. Zhang, Y. Chi and Y. Liang, Provable non-convex phase retrieval with outliers: Median truncated Wirtinger flow, International Conference on Machine Learning, JMLR, New York (2016), 1022–1031. Search in Google Scholar
© 2024 Walter de Gruyter GmbH, Berlin/Boston
Articles in the same Issue
- Frontmatter
- Local convergence of the error-reduction algorithm for real-valued objects
- A rotation total variation regularization for full waveform inversion
- Stochastic data-driven Bouligand–Landweber method for solving non-smooth inverse problems
- On determining the fractional exponent of the subdiffusion equation
- Tow-parameter quasi-boundary value method for a backward abstract time-degenerate fractional parabolic problem
- Determining both leading coefficient and source in a nonlocal elliptic equation
- Discrete dynamical systems: Inverse problems and related topics
- Stability estimates for an inverse problem of determining time-dependent coefficients in a system of parabolic equations
- Numerical solutions to inverse nodal problems for the Sturm–Liouville operator and their applications
- Revisiting linear machine learning through the perspective of inverse problems
Articles in the same Issue
- Frontmatter
- Local convergence of the error-reduction algorithm for real-valued objects
- A rotation total variation regularization for full waveform inversion
- Stochastic data-driven Bouligand–Landweber method for solving non-smooth inverse problems
- On determining the fractional exponent of the subdiffusion equation
- Tow-parameter quasi-boundary value method for a backward abstract time-degenerate fractional parabolic problem
- Determining both leading coefficient and source in a nonlocal elliptic equation
- Discrete dynamical systems: Inverse problems and related topics
- Stability estimates for an inverse problem of determining time-dependent coefficients in a system of parabolic equations
- Numerical solutions to inverse nodal problems for the Sturm–Liouville operator and their applications
- Revisiting linear machine learning through the perspective of inverse problems