06.09.2021 Views

Linear Algebra, Theory And Applications, 2012a

Linear Algebra, Theory And Applications, 2012a

Linear Algebra, Theory And Applications, 2012a

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

15.2. THE QR ALGORITHM 389<br />

Proof: Let<br />

Q =(q 1 , ··· , q n ) ,Q (k) = ( q k 1, ··· , q k )<br />

n<br />

where the q are the columns. Also denote by rij k the ijth entry of R k . Thus<br />

⎛<br />

Q (k) R k = ( q k 1, ··· , q k ) ⎜<br />

n ⎝<br />

r k 11<br />

. ..<br />

∗<br />

0 r k nn<br />

⎞<br />

⎟<br />

⎠<br />

It follows<br />

and so<br />

Therefore,<br />

Next consider the second column.<br />

r k 11q k 1 → q 1<br />

r11 k = ∣ r<br />

k<br />

11 q k ∣<br />

1 → 1<br />

q k 1 → q 1 .<br />

r k 12q k 1 + r k 22q k 2 → q 2<br />

Taking the inner product of both sides with q k 1 it follows<br />

(<br />

lim<br />

k→∞ rk 12 = lim q2 · q k )<br />

1 =(q2 · q 1 )=0.<br />

k→∞<br />

Therefore,<br />

lim<br />

k→∞ rk 22q k 2 = q 2<br />

and since r k 22 > 0, it follows as in the first part that r k 22 → 1. Hence<br />

Continuing this way, it follows<br />

for all i ≠ j and<br />

lim<br />

k→∞ rk jj =1,<br />

lim<br />

k→∞ qk 2 = q 2 .<br />

lim<br />

k→∞ rk ij =0<br />

lim<br />

k→∞ qk j = q j .<br />

Thus R k → I and Q (k) → Q. This proves the first part of the lemma.<br />

The second part follows immediately. If QR = Q ′ R ′ = A where A −1 exists, then<br />

Q ∗ Q ′ = R (R ′ ) −1<br />

and I need to show both sides of the above are equal to I. The left side of the above is<br />

unitary and the right side is upper triangular having positive entries on the diagonal. This<br />

is because the inverse of such an upper triangular matrix having positive entries on the<br />

main diagonal is still upper triangular having positive entries on the main diagonal and<br />

the product of two such upper triangular matrices gives another of the same form having<br />

positive entries on the main diagonal. Suppose then that Q = R where Q is unitary and R<br />

is upper triangular having positive entries on the main diagonal. Let Q k = Q and R k = R.<br />

It follows<br />

IR k → R = Q<br />

and so from the first part, R k → I but R k = R and so R = I. Thus applying this to<br />

Q ∗ Q ′ = R (R ′ ) −1 yields both sides equal I.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!