1 / 26

Chem 302 - Math 252

Chem 302 - Math 252. Chapter 2 Solutions of Systems of Linear Equations / Matrix Inversion. Solutions of Systems of Linear Equations. n linear equations, n unknowns Three possibilities Unique solution No solution Infinite solutions

Download Presentation

Chem 302 - Math 252

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Chem 302 - Math 252 Chapter 2Solutions of Systems of Linear Equations / Matrix Inversion

  2. Solutions of Systems of Linear Equations • n linear equations, n unknowns • Three possibilities • Unique solution • No solution • Infinite solutions • Numerically systems that are almost singular cause problems • Range of solutions • Ill-conditioned problem Singular Systems

  3. Solutions of Systems of Linear Equations • Direct Methods • Determine solution in finite number of steps • Usually preferred • Round-off error can cause problems • Indirect Methods • Use iteration scheme • Require infinite operations to determine exact solution • Useful when Direct Methods fail

  4. Direct Methods • Cramer’s Rule • Gaussian Elimination • Gauss-Jordan Elimination • Maximum Pivot Strategy

  5. Cramer’s Rule • Write coefficient matrix (A) • Evaluate |A| • If |A|=0 then singular • Form A1 • Replace column 1 of A with answer column • Compute x1 = |A1|/|A| • Repeat 3 and 4 for other variables

  6. Cramer’s Rule Not singular: System has unique solution

  7. Cramer’s Rule

  8. Cramer’s Rule • Good for small systems • Good if only one or two variables are needed • Very slow and inefficient for large systems • n order system requires (n+1)! × & (n+1)! Additions • 2nd order 6 ×, 6 + • 10th order 3628800 ×, 3628800 + • 600th order 1.27×101408 ×, 1.27×101408 +

  9. Gaussian Elimination • Form augmented matrix • Use elementary row operations to transform the augmented matrix so that the A portion is in upper triangular form • Switch rows • Multiply row by constant • Linear combination of rows • Use back substitution to find solutions • Requires n3+n2- n ×, n3+½n2- n +

  10. Gaussian Elimination

  11. Gauss-Jordan Elimination • Form augmented matrix • Normalize 1st row • Use elementary row operations to transform the augmented matrix so that the A portion is the identity matrix • Switch rows • Multiply row by constant • Linear combination of rows • Requires ½n3+n2- 2½n+2 ×, ½n3-1½n+1 + • Can also be used to find matrix inverse

  12. Gauss-Jordan Elimination

  13. Maximum Pivot Strategy • Elimination methods can run into difficulties if one or more of diagonal elements is close to (or exactly) zero • Normalize row with largest (magnitude) element.

  14. Gauss-Jordan Elimination

  15. Comparison of Direct Methods • Small systems (n<10) not a big deal • Large systems critical Number of floating point operations

  16. Comparison of Direct Methods Time required on a 300 MFLOP computer (500 TFLOP)

  17. Indirect Methods • Jacobi Method • Gauss-Seidel Method • Use iterations • Guess solution • Iterate to self consistent • Can be combined with Direct Methods

  18. Jacobi Method • Rearrange system of equations to isolate the diagonal elements • Guess solution • Iterate until self-consistent

  19. Jacobi Method

  20. Gauss-Seidel Method • Same as Jacobi method, but use updated values as soon as they are calculated.

  21. Jacobi Method Gauss-Seidel Method

  22. Indirect Methods • Sufficient condition • Diagonally dominant • Large problems • Sparse matrix (many zeros)

More Related