0% found this document useful (0 votes)
50 views2 pages

Inverse Matriz

Gaussian elimination is a method to solve systems of linear equations by transforming the original system into an upper triangular system. This is done by applying row operations like row interchanging and replacing a row with a linear combination of rows to the augmented matrix. The original system is then equivalent to the upper triangular system, which can be solved using back substitution. The matrix inverse is the inverse of a square matrix S that satisfies the property SX = I and XS = I. It can be computed by solving n linear systems with S and the unit vectors as the right hand sides to obtain the columns of the inverse matrix X.

Uploaded by

Julian Gutierrez
Copyright
© Attribution Non-Commercial (BY-NC)
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
0% found this document useful (0 votes)
50 views2 pages

Inverse Matriz

Gaussian elimination is a method to solve systems of linear equations by transforming the original system into an upper triangular system. This is done by applying row operations like row interchanging and replacing a row with a linear combination of rows to the augmented matrix. The original system is then equivalent to the upper triangular system, which can be solved using back substitution. The matrix inverse is the inverse of a square matrix S that satisfies the property SX = I and XS = I. It can be computed by solving n linear systems with S and the unit vectors as the right hand sides to obtain the columns of the inverse matrix X.

Uploaded by

Julian Gutierrez
Copyright
© Attribution Non-Commercial (BY-NC)
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/ 2

Gaussian Elimination Example 1

Gaussian elimination for the solution of a linear system transforms the Suppose that we want to solve

system Sx = f into an equivalent system U x = c with upper triangular


    
2 4 −2 x1 2
matrix U (that means all entries in U below the diagonal are zero). 4 9 −3 x2  =  8  . (1)
−2 −3 7 x3 10
This transformation is done by applying three types of transformations to
the augmented matrix (S | f ). We apply Gaussian elimination. To keep track of the operations, we use,
e.g., R2 = R2 − 2 ∗ R1 , which means that the new row 2 is computed by
Type 1: Interchange two equations; and subtracting 2 times row 1 from row 2.

Type 2: Replace an equation with the sum of the same equation and a
   
2 4 −2 2 2 4 −2 2
 4 9 −3 8  →  0 1 1 4  R2 = R2 − 2 ∗ R1
multiple of another equation. −2 −3 7 10 0 1 5 12 R3 = R3 + R1
Once the augmented matrix (U | f ) is transformed into (U | c), where U 
2 4 −2 2

is an upper triangular matrix, we can solve this transformed system → 0 1 1 4 


0 0 4 8 R3 = R3 − R2
U x = c using backsubstitution.

The original system (1) is equivalent to


    
2 4 −2 x1 2
0 1 1  x2  = 4 . (2)
0 0 4 x3 8

The system (2) can be solved by backsubstitution. We get the solution

x3 = 8/4 = 2, x2 = (4−1∗2)/1 = 2, x1 = (2−4∗2−(−2)∗2)/2 = −1.

M. Heinkenschloss - CAAM335 Matrix Analysis Gaussian Elimination and Matrix Inverse – 1 M. Heinkenschloss - CAAM335 Matrix Analysis Gaussian Elimination and Matrix Inverse – 2

Example 2 Example 2 (cont.)


Suppose that we want to solve If f3 − 2f2 − f1 = 0, then x3 can be chosen arbitrarily and x2 , x1 can be

2 3 −2
   
x1 f1 determined by backsubstitution.
1 −2 3  x2  = f2  . (3) If f3 − 2f2 − f1 = 0, then
4 −1 4 x3 f3
x3 = any scalar, x2 = (f2 −f1 /2−4x3 )∗(−2/7), x1 = (f1 −3x2 +2x3 )/2.
We apply Gaussian elimination.

For example if f1 = f2 = f3 = 1, and if we choose x3 = 0, then


   
2 3 −2 f1 2 3 −2 f1
 1 −2 3 f2  → 0 −7/2 4 f2 − f1 /2  R2 = R2 − 0.5 ∗ R1
4 −1 4 f3 0 −7 8 f3 − 2f1 R3 = R3 − 2 ∗ R1 x3 = 0, x2 = −1/7, x1 = 5/7.
 
2 3 −2 f1
→ 0 −7/2 4 f2 − f1 /2 
0 0 0 f3 − 2f2 − f1 R3 = R3 − 2R2

The original system (3) is equivalent to


    
2 3 −2 x1 f1
0 −7/2 4  x2  =  f2 − f1 /2  . (4)
0 0 0 x3 f3 − 2f2 − f1

The last equation in system (4) reads 0x1 + 0x2 + 0x3 = f3 − 2f2 − f1 .
This can only be satisfied of the right hand side satisfies
f3 − 2f2 − f1 = 0, for example if f1 = f2 = f3 = 1.
M. Heinkenschloss - CAAM335 Matrix Analysis Gaussian Elimination and Matrix Inverse – 3 M. Heinkenschloss - CAAM335 Matrix Analysis Gaussian Elimination and Matrix Inverse – 4
The Matrix Inverse Computation of the Matrix Inverse
A square matrix S ∈ Rn×n is invertible if there exists a matrix We want to find the inverse of S ∈ Rn×n , that is we want to find a
X ∈ Rn×n such that matrix X ∈ Rn×n such that SX = I.
XS = I and SX = I. I Let X:,j denote the jth column of X, i.e., X = (X:,1 , . . . , X:,n ).
Consider the matrix-matrix product SX. The jth column of SX is
The matrix X is called the inverse of S and is denoted by S −1 .
the matrix-vector product SX:,j , i.e., SX = (SX:,1 , . . . , SX:,n ).
I An invertible matrix is also called non-singular. The jth column of the identity I is the jth unit vector
A matrix is called non-invertible or singular if it is not invertible. ej = (0, . . . , 0, 1, 0, . . . , 0)T .
I A matrix S ∈ Rn×n cannot have two different inverses. Hence SX = (SX:,1 , . . . , SX:,n ) = (e1 , . . . , en ) = I implies that we
In fact, if X, Y ∈ Rn×n are two matrices with XS = I and SY = I,
can compute the columns X:,1 , . . . , X:,n of the inverse of S by
then
X = XI = X(SY ) = (XS)Y = IY = Y. solving n systems of linear equations
I If S ∈ Rn×n is invertible, then Sx = f implies x = S −1 Sx = S −1 f , SX:,1 = e1 ,
i.e., for every f the linear system Sx = f has a solution x = S −1 f .
The linear system Sx = f cannot have more than one solution
..
.
because Sx = f and Sy = f imply
S(x − y) = Sx − Sy = f − f = 0 and x − y = S −1 0 = 0. SX:,n = en .
Hence if S is invertible, then for every f the linear system Sx = f
has the unique solution x = S −1 f . Note that if for every f the linear system Sx = f has a unique
I We will see later that if for every f the linear system Sx = f has a solution x, then there exists a unique X = (X:,1 , . . . , X:,n ) with
unique solution x, then S is invertible. SX = I.
M. Heinkenschloss - CAAM335 Matrix Analysis Gaussian Elimination and Matrix Inverse – 5 M. Heinkenschloss - CAAM335 Matrix Analysis Gaussian Elimination and Matrix Inverse – 6

Example 3 Example 4
Suppose that we want the inverse of Suppose that we want the inverse of
 
2 4 −2  
2 3 −2
S= 4 9 −3 .
−2 −3 7 S = 1 −2 3 .
4 −1 4
We can use Gaussian Elimination to solve the systems
SX:,1 = e1 , SX:,2 = e2 , SX:,3 = e3 for the three columns of X = S −1 We can use Gaussian Elimination to solve the systems
    SX:,1 = e1 , SX:,2 = e2 , SX:,3 = e3 for the three columns of X = S −1
2 4 −2 1 0 0 2 4 −2 1 0 0
 4 9 −3 0 1 0  →  0 1 1 −2 1 0     
2 3 −2 1 0 0 2 3 −2 1 0 0
−2 −3 7 0 0 1 0 1 5 1 0 1
   1 −2 3 0 1 0  →  0 −7/2 4 −1/2 1 0 
2 4 −2 1 0 0
4 −1 4 0 0 1 0 −7 8 −2 0 1
→ 0 1 1 −2 1 0   
0 0 4 3 −1 1 2 3 −2 1 0 0
→  0 −7/2 4 −1/2 1 0 .
   
2 4 0 5/2 −1/2 1/2 1 0 0 27/4 −11/4 3/4
→  0 1 0 −11/4 5/4 −1/4  →  0 1 0 −11/4 5/4 −1/4  . 0 0 0 −1 −2 1
0 0 1 3/4 −1/4 1/4 0 0 1 3/4 −1/4 1/4
  None of the linear systems SX:,1 = e1 , SX:,2 = e2 , SX:,3 = e3 has a
27 −11 3 solution. Therefore, S is not invertible.
1
S −1 = −11 5 −1  .
4
3 −1 1

M. Heinkenschloss - CAAM335 Matrix Analysis Gaussian Elimination and Matrix Inverse – 7 M. Heinkenschloss - CAAM335 Matrix Analysis Gaussian Elimination and Matrix Inverse – 8

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy