06.09.2021 Views

A First Course in Linear Algebra, 2017a

A First Course in Linear Algebra, 2017a

A First Course in Linear Algebra, 2017a

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

416 Spectral Theory<br />

7.4.4 QR Factorization<br />

In this section, a reliable factorization of matrices is studied. Called the QR factorization of a matrix, it<br />

always exists. While much can be said about the QR factorization, this section will be limited to real<br />

matrices. Therefore we assume the dot product used below is the usual dot product. We beg<strong>in</strong> with a<br />

def<strong>in</strong>ition.<br />

Def<strong>in</strong>ition 7.83: QR Factorization<br />

Let A be a real m × n matrix. Then a QR factorization of A consists of two matrices, Q orthogonal<br />

and R upper triangular, such that A = QR.<br />

The follow<strong>in</strong>g theorem claims that such a factorization exists.<br />

Theorem 7.84: Existence of QR Factorization<br />

Let A be any real m × n matrix with l<strong>in</strong>early <strong>in</strong>dependent columns. Then there exists an orthogonal<br />

matrix Q and an upper triangular matrix R hav<strong>in</strong>g non-negative entries on the ma<strong>in</strong> diagonal such<br />

that<br />

A = QR<br />

The procedure for obta<strong>in</strong><strong>in</strong>g the QR factorization for any matrix A is as follows.<br />

Procedure 7.85: QR Factorization<br />

Let A be an m × n matrix given by A = [ A 1 A 2 ···<br />

]<br />

A n where the Ai are the l<strong>in</strong>early <strong>in</strong>dependent<br />

columns of A.<br />

1. Apply the Gram-Schmidt Process 4.135 to the columns of A, writ<strong>in</strong>gB i for the result<strong>in</strong>g<br />

columns.<br />

2. Normalize the B i ,tof<strong>in</strong>dC i = 1<br />

‖B i ‖ B i.<br />

3. Construct the orthogonal matrix Q as Q = [ C 1 C 2 ··· C n<br />

]<br />

.<br />

4. Construct the upper triangular matrix R as<br />

⎡<br />

⎤<br />

‖B 1 ‖ A 2 •C 1 A 3 •C 1 ··· A n •C 1<br />

0 ‖B 2 ‖ A 3 •C 2 ··· A n •C 2<br />

R =<br />

0 0 ‖B 3 ‖ ··· A n •C 3<br />

⎢<br />

⎥<br />

⎣<br />

...<br />

...<br />

...<br />

... ⎦<br />

0 0 0 ··· ‖B n ‖<br />

5. F<strong>in</strong>ally, write A = QR where Q is the orthogonal matrix and R is the upper triangular matrix<br />

obta<strong>in</strong>ed above.

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

Saved successfully!

Ooh no, something went wrong!