Scribd
Upload
151 views10 pages

Introduction To Numerical Analysis For Engineers: - Systems of Linear Equations Mathews

This document summarizes numerical methods for solving systems of linear equations, including Gaussian elimination, LU factorization, and iterative methods. It discusses solving tridiagonal systems that arise from finite difference approximations. Special matrices are covered, such as banded, symmetric, and positive definite matrices, and how their structure can be exploited using techniques like skyline storage and Cholesky factorization.

Uploaded by

Rija Hossain
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
151 views10 pages

Introduction To Numerical Analysis For Engineers: - Systems of Linear Equations Mathews

This document summarizes numerical methods for solving systems of linear equations, including Gaussian elimination, LU factorization, and iterative methods. It discusses solving tridiagonal systems that arise from finite difference approximations. Special matrices are covered, such as banded, symmetric, and positive definite matrices, and how their structure can be exploited using techniques like skyline storage and Cholesky factorization.

Uploaded by

Rija Hossain
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 10

Numerical Methods for Engineers 13.

002 Lecture 5
Introduction to Numerical Analysis for Engineers
Systems of Linear Equations Mathews
Cramers Rule
Gaussian Elimination 3.3-3.5
Numerical implementation
3.3-3.4
Numerical stability
Partial Pivoting
Equilibration
Full Pivoting
Multiple right hand sides
Computation count
LU factorization 3.5
Error Analysis for Linear Systems 3.4
Condition Number
Special Matrices
Iterative Methods 3.6
Jacobis method
Gauss-Seidel iteration
Convergence
Numerical Methods for Engineers 13.002 Lecture 5
Linear Systems of Equations
Tri-diagonal Systems
y(x,t)
x
i
Forced Vibration of a String
f(x,t)
Harmonic excitation
f(x,t) = f(x) cos(t)
Differential Equation
Boundary Conditions
Discrete Difference Equations
Matrix Form
Tridiagonal Matrix
Finite Difference
Symmetric, positive definite: No pivoting needed
Numerical Methods for Engineers 13.002 Lecture 5
Linear Systems of Equations
Tri-diagonal Systems
General Tri-diagonal Systems
LU Factorization
Numerical Methods for Engineers 13.002 Lecture 5
Linear Systems of Equations
Tri-diagonal Systems
LU Factorization
Reduction
Forward Substitution
Back Substitution
LU Factorization: 2*(n-1) operations
Forward substitution: n-1 operations
Back substitution: n-1 operations
Total: 4(n-1) ~ O(n) operations
Numerical Methods for Engineers 13.002 Lecture 5
Linear Systems of Equations
Special Matrices
General, Banded Coefficient Matrix
0
0
p
q
p super-diagonals
q sub-diagonals
w = p+q+1 bandwidth
b is half-bandwidth
Banded Symmetric Matrix
Numerical Methods for Engineers 13.002 Lecture 5
0
0
q
0
0
p
Linear Systems of Equations
Special Matrices
Banded Coefficient Matrix
Gaussian Elimination
No Pivoting
=
=
Numerical Methods for Engineers 13.002 Lecture 5
0
0
q
=
00
Banded Coefficient Matrix
Gaussian Elimination
With Pivoting
q
Linear Systems of Equations
Special Matrices
Numerical Methods for Engineers 13.002 Lecture 5
0
0
q
q
p
p
Banded Coefficient Matrix
Compact Storage
Linear Systems of Equations
Special Matrices
0
0
n
2
Diagonal
n(p+2q+1)
q
i
j
i
j - i
Numerical Methods for Engineers 13.002 Lecture 5
Linear Systems of Equations
Special Matrices
Sparse and Banded Coefficient Matrix
Skyline Systems
..

.
.
..
..
Skyline storage applicable when no pivoting is needed, e.g. for banded,
symmetric, and positive definite matrices: FEM and FD methods. Skyline
solvers are usually based on Choleski factorization
Skyline
Storage
Pointers
1 4 9 11 16 20
0
0
Numerical Methods for Engineers 13.002 Lecture 5
Linear Systems of Equations
Special Matrices
Symmetric, Positive Definite Coefficient Matrix
No pivoting needed
Choleski Factorization
where

You might also like