Biomechanics_3

# Biomechanics_3 - Lecture 4 BMEn 3001 Biomechanics 10...

This preview shows pages 1–10. Sign up to view the full content.

Lecture 4 BMEn 3001 Biomechanics 10 September 2008

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
Last Lecture • Linear Algebraic Systems • LU Decomposition • Pivoting
Topics for Today • LU (Fnish) • Pivoting • Matrix Inversion • The Eigenproblem • Condition Number

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
An LU Decomposition • Note that we can get back the original matrix by multiplying the lower and upper parts. • Note also that A = 1 0 0 2 1 0 4 0 1 1 1 1 0 2 4 0 0 2 = 1 1 1 2 0 2 4 4 6 A L + U
Use LU to solve. .. 1 1 1 2 0 2 4 4 6 x 1 x 2 x 3 = 12 2 62

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
Storage The L and U matrices only overlap on the diagonal, and L has one's on the diagonal, so we really don't need to store them. In fact, if we don't need to get the original matrix back, we can store L and U over it. U 11 U 12 U 13 L U 1 N L 21 U 22 U 23 L U 2 N L 31 L 32 U 33 L U 3 N M M M O M L N 1 L N 2 L N 3 L U NN
Pivoting • Let's solve the matrix problem at right by Gauss elimination. • First, we'll subtract the frst row ±rom the second row, right? 1 2 0 1 2 3 0 1 1 x y z = 7 10 3

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
Pivoting • All looks good, but our next step should be to divide the second row by its diagonal element, which is zero! 1 2 0 0 0 3 0 1 1 x y z = 7 3 3 • The situation looks dire, but the determinant of the original matrix is not zero, so there is hope. (1•2•1+2•3•0+0•1•1-0•2•0-2•1•1-1•3•1 = -3) 1 2 0 1 2 3 0 1 1 x y z = 7 10 3
Pivoting • Fortunately, the

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
This is the end of the preview. Sign up to access the rest of the document.

## This note was uploaded on 09/16/2008 for the course BMEN 3001 taught by Professor Victorbarocas during the Spring '08 term at Minnesota.

### Page1 / 35

Biomechanics_3 - Lecture 4 BMEn 3001 Biomechanics 10...

This preview shows document pages 1 - 10. Sign up to view the full document.

View Full Document
Ask a homework question - tutors are online