3.II.19D

Numerical Analysis
Part IB, 2006

(a) Define the QR factorization of a rectangular matrix and explain how it can be used to solve the least squares problem of finding an xRnx^{*} \in \mathbb{R}^{n} such that

Axb=minxRnAxb, where ARm×n,bRm,mn,\left\|A x^{*}-b\right\|=\min _{x \in \mathbb{R}^{n}}\|A x-b\|, \quad \text { where } \quad A \in \mathbb{R}^{m \times n}, \quad b \in \mathbb{R}^{m}, \quad m \geqslant n,

and the norm is the Euclidean distance y=i=1myi2\|y\|=\sqrt{\sum_{i=1}^{m}\left|y_{i}\right|^{2}}.

(b) Define a Householder transformation (reflection) HH and prove that HH is an orthogonal matrix.

(c) Using Householder reflection, solve the least squares problem for the case

A=[241121],b=[151]A=\left[\begin{array}{rr} 2 & 4 \\ 1 & -1 \\ 2 & 1 \end{array}\right], \quad b=\left[\begin{array}{l} 1 \\ 5 \\ 1 \end{array}\right]

and find the value of Axb=minxR2Axb\left\|A x^{*}-b\right\|=\min _{x \in \mathbb{R}^{2}}\|A x-b\|.