scispace - formally typeset
Journal ArticleDOI

Linear least squares solutions by householder transformations

Peter A. Businger, +1 more
- 01 Jun 1965 - 
- Vol. 7, Iss: 3, pp 269-276
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.