Set . Since , it is well defined, symmetric and idempotent: , . Therefore is symmetric, and
Because , this equals exactly when . Hence .
If , then and , so . For arbitrary , , giving
Thus the Householder map fixes the hyperplane orthogonal to and reverses the component in its normal direction: it is a reflection, not merely an arbitrary orthogonal map.
Choose , with , and put . The opposite sign choice guarantees , even if is already parallel to . Since ,
Consequently , yielding
The sign choice also avoids cancellation in constructing numerically.
Proceed column by column. At stage , consider the entries below the diagonal in column , indexed . If this tail is nonzero, part (c) provides a Householder matrix sending it to a multiple of its first coordinate vector. Embed this matrix as an orthogonal matrix acting identically on coordinates , and apply the similarity transformation . The tail's entries below row become zero. Symmetry is preserved, so the corresponding entries in row vanish too. If the tail is already zero, use the identity transformation.
Earlier reduced columns remain unchanged in their required zero positions: all their entries in the newly mixed tail coordinates were already zero, and the first coordinates are fixed. Thus each stage preserves the previous zeros while creating the next ones. After stages, all entries at distance greater than one from the diagonal vanish. With the product one has and
This constructive reduction uses only reflections on successive trailing subspaces.

Articles by others on the same topic (0)

There are currently no matching articles.