next up previous


Householder QR Factorization


Apply the Householder reflection to the matrix tex2html_wrap_inline7395 to obtain the QR factorization.

Example 2.3.1. Suppose tex2html_wrap_inline7399 and assume that the Householder matrices H1 and H2 have been computed so that
displaymath7323
Concentrating on the highlighted vector tex2html_wrap_inline7405 , we determine a Householder matrix tex2html_wrap_inline7407 such that
displaymath7324
Choosing tex2html_wrap_inline7409we get
displaymath7325
Next consider the highlighted vector tex2html_wrap_inline7411 and determine tex2html_wrap_inline7413 such that
displaymath7326
Choosing tex2html_wrap_inline7415 we have


displaymath7327
By setting Q=H1H2H3H4, we obtain QR=H1H2H3H4H4H3H2H1A=A.

Proposition 2.3.1. If tex2html_wrap_inline7421then there exist the Householder matrices tex2html_wrap_inline7423 such that
displaymath7328

displaymath7329
and
displaymath7330
where tex2html_wrap_inline7425 is orthogonal and tex2html_wrap_inline7427 is upper triangular.

Example 2.3.2. Find the Householder QR factorization for
displaymath7331
In example 2.1.1 there has been found the Householder matrix for the transformation of the first column vector tex2html_wrap_inline7431 of A:
displaymath7332
Find that
displaymath7333
To find tex2html_wrap_inline7435 we compute the according Householder vector
displaymath7334
Hence
displaymath7335
and
displaymath7336
and also
displaymath7337

displaymath7338
Find also the orthogonal matrix
displaymath7339

displaymath7340
and check the result

displaymath7341

Example 2.3.3*. Find the Householder QR factorization of
displaymath7342

The vector that has to be transformed is tex2html_wrap_inline7441where tex2html_wrap_inline7443 Construct the vector
displaymath7343
Choose the sign minus for the coefficient of tex2html_wrap_inline7205 and take into account that H depends only on the direction of tex2html_wrap_inline7153:
displaymath7344
Find the Householder matrix
displaymath7345

displaymath7346
Verify that H1 annihilates all the elements of the first column of A but the first one. Indeed,
displaymath7347
Further we transform the vector tex2html_wrap_inline7455where tex2html_wrap_inline7457 Find the Householder vector according to tex2html_wrap_inline7177
displaymath7348
Choose the sign minus for coefficient of tex2html_wrap_inline7205:
displaymath7349
We obtain the Householder matrix according to this vector
displaymath7350

displaymath7351

displaymath7352
and find that
displaymath7353
Thus,
displaymath7354
and
displaymath7355
Let us check the result:

displaymath7356

Exercise 2.3.1.* Find the QR factorization of A if
displaymath7357

displaymath7358


next up previous