Abstract
In some previous papers the author extended two algorithms proposed by Z. Kovarik for approximate orthogonalization of a finite set of linearly independent vectors from a Hilbert space, to the case when the vectors are rows (not necessary linearly independent) of an arbitrary rectangular matrix. In this paper we describe combinations between these two methods and the classical Kaczmarz’s iteration. We prove that, in the case of a consistent least-squares problem, the new algorithms so obtained converge to any of its solutions (depending on the initial approximation). The numerical experiments described in the last section of the paper on a problem obtained after the discretization of a first kind integral equation ilustrate the fast convergence of the new algorithms.
Similar content being viewed by others
References
A. Björck,Numerical methods for least squares problems, SIAM Philadelphia, 1996.
T. L. Boullion, P. L. Odell,Generalized inverse matrices, Wiley — Interscience, New York, 1971.
P. G. Ciarlet,Introduction à l’analyse numérique et à l’ optimisation, Masson, Paris, 1982.
H. W. Engl,Regularization methods for the stable solution of inverse problems, Surv. Math. Ind.3(1993), 71- 143.
R. A. Horn, C. R. Johnson,Theory of matrices, Cambrige Univ. Press, Cambrige, 1986.
W. J. Kammerer and M. Z. Nashed,On the convergence of the conjugate gradient method for singular linear operator equations, SIAM J. Numer. Anal.,9(1)(1972), 165 - 181.
Z. Kovarik,Some iterative methods for improving orthogonality, SIAM J. Numer. Anal.,73(1970), 386–389.
C. Popa,Characterization of the solutions set of inconsistent least-squares problems by an extended Kaczmarz algorithm, Korean Journal on Computational and Applied Mathematics,63(1999), 523 - 535.
C. Popa,Extension of an approximate orthogonalization algorithm to arbitrary rectangular matrices, Preprint 18(1998), Institute of Mathematics of The Romanian Academy, Bucharest; sent for publication in Linear Algebra and Applications, 1999.
C. Popa,On numerical solution of first kind Fredholm integral equation, to appear in Revue Roumaine Math. Pures Appl., 2000.
C. Popa,An iterative method for improving orthogonality of rows of arbitrary rectangular matrices, sent for publication in International Journal of Computer Mathematics, 1999.
K. Tanabe,Projection method for solving a singular system of linear equations and its applications, Numer. Math.,17(1971), 203 - 214.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Popa, C. A fast Kaczmarz-Kovarik algorithm for consistent least-squares problems. Korean J. Comput. & Appl. Math. 8, 9–26 (2001). https://doi.org/10.1007/BF03011619
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF03011619