Abstract The Partial Total Least Squares (PTLS) subroutine solves the Total Least Squares (TLS) problem AX ≈ B by using a Partial Singular Value Decomposition (PSVD), hereby improving the computational efficiency of the classical TLS algorithm (a factor 2 approximately). The computational improvement of PTLS is based on the observation that the TLS solution can be computed from any orthogonal basis of the right singular subspace corresponding to the smallest singular values of the augmented matrix [A;B]. The dimension of this subspace may be determined by the rank of the TLS approximation of [A;B] or depend on a given upper bound (e.g. a perturbation level) for those smallest singular values. The routine PTLS allows more than one vector in B. It also solves determined and underdetermined sets of equations by computing the minimum norm solution.
[1]
J. Vandewalle,et al.
Analysis and solution of the nongeneric total least squares problem
,
1988
.
[2]
Sabine Van Huffel,et al.
An efficient and reliable algorithm for computing the singular subspace of a matrix, associated with its smallest singular values
,
1987
.
[3]
Sabine Van Huffel.
Analysis of the Total Least Squares Problem and its Use in Parameter Estimation
,
1987
.
[4]
V. Klema.
LINPACK user's guide
,
1980
.
[5]
Tony F. Chan,et al.
An Improved Algorithm for Computing the Singular Value Decomposition
,
1982,
TOMS.
[6]
Charles L. Lawson,et al.
Basic Linear Algebra Subprograms for Fortran Usage
,
1979,
TOMS.