Efficient iterative solutions to complex-valued nonlinear least-squares problems with mixed linear and antilinear operators
Year: 2021
Type: article
Abstract: We consider a setting in which it is desired to find an optimal complex vector $${\mathbf {x}}\in {\mathbb {C}}^N$$ x ∈ C N that satisfies $${\mathcal {A}}({\mathbf {x}}) \approx {\mathbf {b}}$$ A ( x ) ≈ b in a least-squares sense, where $${\mathbf {b}} \in {\mathbb {C}}^M$$ b ∈ C M is a data vector (possibly noise-corrupted), and $${\mathcal {A}}(\cdot ): {\mathbb {C}}^N \rightarrow {\mathbb {C}}^M$$ A ( · ) : C N → C M is a measurement operator. If $${\mathcal {A}}(\cdot )$$ A ( · ) were linear, this reduces to the classical linear least-squares problem, which has a well-known analytic solution as well as powerful iterative solution algorithms. However, instead of linear least-squares, this work considers the more complicated scenario where $${\mathcal {A}}(\cdot )$$ A ( · ) is nonlinear, but can be represented as the summation and/or composition of some operators that are linear and some operators that are antilinear. Some common nonlinear operations that have this structure include complex conjugation or taking the real-part or imaginary-part of a complex vector. Previous literature has shown that this kind of mixed linear/antilinear least-squares problem can be mapped into a linear least-squares problem by considering $${\mathbf {x}}$$ x as a vector in $${\mathbb {R}}^{2N}$$ R 2 N instead of $${\mathbb {C}}^N$$ C N . While this approach is valid, the replacement of the original complex-valued optimization problem with a real-valued optimization problem can be complicated to implement, and can also be associated with increased computational complexity. In this work, we describe theory and computational methods that enable mixed linear/antilinear least-squares problems to be solved iteratively using standard linear least-squares tools, while retaining all of the complex-valued structure of the original inverse problem. An illustration is provided to demonstrate that this approach can simplify the implementation and reduce the computational complexity of iterative solution algorithms.
Source: Optimization and Engineering
Authors: Tae Hyung Kim, Justin P. Haldar
Institution: University of Southern California
Language: English
Subfield: Computational Mechanics
Field: Engineering
Domain: Physical Sciences
Open Access status: hybrid
Estimated APC paid: $2,780
Funders: National Science Foundation, Kwanjeong Educational Foundation, National Institutes of Health, Annenberg School for Communication and Journalism, University of Southern California, National Institute of Mental Health, National Cancer Institute, National Institute of Neurological Disorders and Stroke, National Institute of Biomedical Imaging and Bioengineering
Awards: G1066422085, CAREER: Low-Rank Matrix Modeling for Constrained Reconstruction from Noisy and Sparsely-Sampled Data, Area B: Precise DCE-MRI Assessment of Brain Tumors, G3048631115, Faster MRI with Sparse Sampling and Low-Rank Modeling, G9000024958, G9000050521, G9000156017, Next generation in-vivo diffusion imaging at submillimeter resolution
OpenAlex ID: W3042574396
