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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

2.2. LU Factorization 101<br />

Thenextstepistotake2timesthetoprowandaddtothe bottom <strong>in</strong> the matrix on the right. To ensure that<br />

the product is unchanged, you place a −2 <strong>in</strong> the bottom left <strong>in</strong> the matrix on the left. Thus the next step<br />

yields<br />

⎡<br />

⎣<br />

1 0 0<br />

2 1 0<br />

−2 0 1<br />

⎤⎡<br />

⎦⎣<br />

1 2 3<br />

0 −1 −5<br />

0 7 4<br />

Next take 7 times the middle row on right and add to bottom row. Updat<strong>in</strong>g the matrix on the left <strong>in</strong> a<br />

similar manner to what was done earlier,<br />

⎡<br />

1 0<br />

⎤⎡<br />

0 1 2<br />

⎤<br />

3<br />

⎣ 2 1 0 ⎦⎣<br />

0 −1 −5 ⎦<br />

−2 −7 1 0 0 −31<br />

⎤<br />

⎦<br />

At this po<strong>in</strong>t, stop. You are done.<br />

♠<br />

The method just described is called the multiplier method.<br />

2.2.3 Solv<strong>in</strong>g Systems us<strong>in</strong>g LU Factorization<br />

One reason people care about the LU factorization is it allows the quick solution of systems of equations.<br />

Here is an example.<br />

Example 2.69: LU factorization to Solve Equations<br />

Suppose you want to f<strong>in</strong>d the solutions to<br />

⎡<br />

⎣<br />

1 2 3 2<br />

4 3 1 1<br />

1 2 3 0<br />

⎡<br />

⎤<br />

⎦⎢<br />

⎣<br />

x<br />

y<br />

z<br />

w<br />

⎤<br />

⎡<br />

⎥<br />

⎦ = ⎣<br />

1<br />

2<br />

3<br />

⎤<br />

⎦.<br />

Solution.<br />

Of course one way is to write the augmented matrix and gr<strong>in</strong>d away. However, this <strong>in</strong>volves more row<br />

operations than the computation of the LU factorization and it turns out that the LU factorization can give<br />

the solution quickly. Here is how. The follow<strong>in</strong>g is an LU factorization for the matrix.<br />

⎡<br />

⎣<br />

1 2 3 2<br />

4 3 1 1<br />

1 2 3 0<br />

⎤<br />

⎡<br />

⎦ = ⎣<br />

1 0 0<br />

4 1 0<br />

1 0 1<br />

⎤⎡<br />

⎦⎣<br />

1 2 3 2<br />

0 −5 −11 −7<br />

0 0 0 −2<br />

Let UX = Y and consider LY = B where <strong>in</strong> this case, B =[1,2,3] T . Thus<br />

⎡ ⎤⎡<br />

⎤ ⎡ ⎤<br />

1 0 0 y 1 1<br />

⎣ 4 1 0 ⎦⎣<br />

y 2<br />

⎦ = ⎣ 2 ⎦<br />

1 0 1 y 3 3<br />

⎤<br />

⎦.

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

Saved successfully!

Ooh no, something went wrong!