Preprint
Computer Vision

Separable nonlinear least squares: the variable projection method and its applications

Gene H. Golub(Stanford University), Víctor Pereyra(Weidlinger Associates (United States))
February 14, 2003Inverse Problems768 citations

768

Citations

2

Influential Citations

Inverse Problems

Venue

2003

Year

Abstract

In this paper we review 30 years of developments and applications of the variable projection method for solving separable nonlinear least-squares problems.These are problems for which the model function is a linear combination of nonlinear functions.Taking advantage of this special structure, the method of variable projections eliminates the linear variables obtaining a somewhat more complicated function that involves only the nonlinear parameters.This procedure not only reduces the dimension of the parameter space but also results in a better-conditioned problem.The same optimization method applied to the original and reduced problems will always converge faster for the latter.We present first a historical account of the basic theoretical work and its various computer implementations, and then report on a variety of applications from electrical engineering, medical and biological imaging, chemistry, robotics, vision, and environmental sciences.An extensive bibliography is included.The method is particularly well suited for solving real and complex exponential model fitting problems, which are pervasive in their applications and are notoriously hard to solve.

Analysis

Why This Paper Matters

This paper is a seminal review of the variable projection (VP) method, a technique that has profoundly impacted how separable nonlinear least-squares problems are solved. Since its introduction in the 1970s, the VP method has been a cornerstone for problems where the model is a linear combination of nonlinear functions—a structure that appears frequently in exponential fitting, signal processing, and inverse problems. The paper consolidates 30 years of theoretical developments and practical implementations, making it an essential reference for researchers and practitioners.

The significance lies in the method's ability to reduce the dimensionality of the optimization problem by analytically eliminating linear parameters. This not only simplifies the problem but also improves its conditioning, leading to faster and more robust convergence. The review highlights that the VP method is particularly effective for real and complex exponential models, which are notoriously difficult to solve with standard techniques. By providing a unified treatment, the paper has enabled broader adoption across diverse fields, from medical imaging to robotics.

Technical Contributions

The paper's key technical contributions include:

  • Variable Projection Formulation: A rigorous derivation of the reduced objective function that depends only on nonlinear parameters, achieved by projecting out linear variables.
  • Convergence Analysis: Theoretical results showing that the reduced problem is better conditioned and converges faster than the original when using the same optimization algorithm.
  • Algorithmic Implementations: A historical account of various computer implementations, from early FORTRAN codes to more modern approaches, highlighting practical considerations.
  • Application Spectrum: A comprehensive survey of applications, including exponential fitting in chemistry, image reconstruction in medical imaging, and parameter estimation in robotics and vision.
  • Extensive Bibliography: A curated list of references that serves as a roadmap for researchers entering the field.

Results

As a review paper, it does not present new experimental results. However, it synthesizes findings from numerous studies that demonstrate the VP method's advantages. For instance, it reports that the VP method consistently outperforms full-space optimization in terms of convergence speed and numerical stability. The paper also notes that the method is particularly effective for ill-conditioned exponential fitting problems, where standard approaches often fail. The extensive bibliography provides evidence of the method's widespread success across applications.

Significance

The variable projection method has had a lasting impact on the field of optimization and inverse problems. It has become the method of choice for separable least-squares problems, especially in exponential fitting, which is pervasive in science and engineering. The paper's comprehensive review has helped standardize the approach and has inspired further research into extensions and variants. Its influence extends beyond the original domains, with applications in machine learning, where separable structures appear in neural network training and dictionary learning. The method's ability to reduce problem complexity and improve conditioning remains highly relevant in modern computational practice, making this paper a timeless reference.