Journal ArticleDOI
Linear least squares solutions by householder transformations
Peter A. Businger,Gene H. Golub +1 more
Reads0
Chats0
TLDR
In this paper, the euclidean norm is unitarily invariant and a vector x is determined such that x is parallel b-Ax parallel = \parallel c - QAx parallel where c denotes the first n components of c.Abstract:
Let A be a given m×n real matrix with m≧n and of rank n and b a given vector. We wish to determine a vector x such that
$$\parallel b - A\hat x\parallel = \min .$$
where ∥ … ∥ indicates the euclidean norm. Since the euclidean norm is unitarily invariant
$$\parallel b - Ax\parallel = \parallel c - QAx\parallel $$
where c=Q b and Q T Q = I. We choose Q so that
$$QA = R = {\left( {_{\dddot 0}^{\tilde R}} \right)_{\} (m - n) \times n}}$$
(1)
and R is an upper triangular matrix. Clearly,
$$\hat x = {\tilde R^{ - 1}}\tilde c$$
where c denotes the first n components of c.read more
Citations
More filters
Book
Matrix algorithms
TL;DR: This volume treats the numerical solution of dense and large-scale eigenvalue problems with an emphasis on algorithms and the theoretical background required to understand them.
Posted Content
Randomized algorithms for matrices and data
TL;DR: This monograph will provide a detailed overview of recent work on the theory of randomized matrix algorithms as well as the application of those ideas to the solution of practical problems in large-scale data analysis.
Journal ArticleDOI
Iteratively Reweighted Least Squares for Maximum Likelihood Estimation, and Some Robust and Resistant Alternatives
TL;DR: The scope of application of iteratively reweighted least squares to statistical estimation problems is considerably wider than is generally appreciated as mentioned in this paper, and it extends beyond the exponential-family-type generalized linear models to other distributions, to non-linear parameterizations, and to dependent observations.
Journal ArticleDOI
Rank revealing QR factorizations
TL;DR: An algorithm is presented for computing a column permutation Pi and a QR-factorization of an m by n (m or = n) matrix A such that a possible rank deficiency of A will be revealed in the triangular factor R having a small lower right block.
References
More filters
Journal ArticleDOI
Unitary Triangularization of a Nonsymmetric Matrix
TL;DR: This note points out that the same result can be obtained with fewer arithmetic operations, and, in particular, for inverting a square matrix of order N, at most 2(N-1) square roots are required.