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

Calculus of generalized inverses of matrices Part I. General theory

TL;DR: In this article, a generalized inverse (g-inverse) is defined and a classification of g-inverses based on their uses and their interrelationships is discussed.
Journal ArticleDOI

Overdetermined photoelastic solutions using least squares

TL;DR: In this article, the least squares method is presented for obtaining an overdetermined solution from the various relationships available in photoelastic stress analysis, and the results are illustrated in two example problems.
Proceedings ArticleDOI

Granularity issues for solving polynomial systems via globally convergent algorithms on a hypercube

TL;DR: The results of decompositions with two different granularities are presented and the experiments were conducted on an iPSC-16 hypercube using actual industrial problems.
Journal ArticleDOI

On rank-deficient pseudoinverses

TL;DR: A certain subcondition number is shown to describe the worst case rounding error growth (in some special norm) and supplies a firm theoretical basis for the application of a rank decision criterion that has already been used successfully in many real life problems.
Journal ArticleDOI

Sensor Selection With Cost Constraints for Dynamically Relevant Bases

TL;DR: In this article, the authors consider cost-constrained sparse sensor selection for full-state reconstruction, applying a well-known greedy algorithm to dynamical systems for which the usual singular value decomposition (SVD) basis may not be available or preferred.
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.