The differentiation of pseudo-inverses and nonlinear least squares problems whose variables separate

Gene H Golub, Victor Pereyra
1973-04-01
Abstract:For given data , we consider the least squares fit of nonlinear models of the form \[ \eta ({\bf a},{\boldsymbol \alpha} ;t) = \sum _{j = 1}^n {a_j \varphi _j ({\boldsymbol \alpha} ;t),\qquad {\bf a} \in \mathcal{R}^n ,\qquad {\boldsymbol \alpha} \in \mathcal{R}^k .} \] For this purpose we study the minimization of the nonlinear functional It is shown that by defining the matrix , and the modified functional , it is possible to optimize first with respect to the parameters , and then to obtain, a posteriors, the optimal parameters . The matrix is the Moore–Penrose generalized inverse of . We develop formulas for the Frechet derivative of orthogonal projectors associated with and also for , under the hypothesis that is of constant (though not necessarily full) rank. Detailed algorithms are presented which make extensive use of …
What problem does this paper attempt to address?