Gauss Jordan elimination

Algebra ->  Algebra  -> College  -> Linear Algebra -> Gauss Jordan elimination      Log On


   

Gauss–Jordan elimination

  (Redirected from Gauss Jordan elimination)
Jump to: navigation, search

In linear algebra, Gauss–Jordan elimination is a version of Gaussian elimination that puts zeros both above and below each pivot element as it goes from the top row of the given matrix to the bottom. In other words, Gauss–Jordan elimination brings a matrix to reduced row echelon form, whereas Gaussian elimination takes it only as far as row echelon form. Every matrix has a reduced row echelon form, and this algorithm is guaranteed to produce it.

Gauss–Jordan elimination is considerably less efficient than Gaussian elimination with backsubstitution when solving a system of linear equations. However, it is well suited for calculating the matrix inverse.

It is named in after Carl Friedrich Gauss and Wilhelm Jordan, because it is a modification of Gaussian elimination as described by Jordan in 1887. However, the method also appears in an article by Clasen published in the same year. Jordan and Clasen probably discovered Gauss–Jordan elimination independently (Althoen & McLaughlin 1987).

In computer science, Gauss-Jordan elimination as an algorithm has a time complexity of O(n3).

[ Application to finding inverses

If Gauss–Jordan elimination is applied on a square matrix, it can be used to calculate the matrix's inverse. This can be done by augmenting the square matrix with the identity matrix of the same dimensions, and through the following matrix operations:

[ A I ] \Rightarrow
A^{-1} [ A I ] \Rightarrow
[ I A^{-1} ].

If the original square matrix, A, is given by the following expression:

 A =
\begin{bmatrix}
2 & -1 & 0 \\
-1 & 2 & -1 \\
0 & -1 & 2
\end{bmatrix}.

Then, after augmenting by the identity, the following is obtained:

 [ A I ] = 
\begin{bmatrix}
2 & -1 & 0 & 1 & 0 & 0\\
-1 & 2 & -1 & 0 & 1 & 0\\
0 & -1 & 2 & 0 & 0 & 1
\end{bmatrix}.

By performing elementary row operations on the [AI] matrix until it reaches reduced row echelon form, the following is the final result:

 [ I A^{-1} ] = 
\begin{bmatrix}
1 & 0 & 0 & \frac{3}{4} & \frac{1}{2} & \frac{1}{4}\\
0 & 1 & 0 & \frac{1}{2} & 1 & \frac{1}{2}\\
0 & 0 & 1 & \frac{1}{4} & \frac{1}{2} & \frac{3}{4}
\end{bmatrix}.

The matrix augmentation can now be undone, which gives the following:

 I =
\begin{bmatrix}
1 & 0 & 0 \\
0 & 1 & 0 \\
0 & 0 & 1
\end{bmatrix}\qquad
 A^{-1} =
\begin{bmatrix}
\frac{3}{4} & \frac{1}{2} & \frac{1}{4}\\
\frac{1}{2} & 1 & \frac{1}{2}\\
\frac{1}{4} & \frac{1}{2} & \frac{3}{4}
\end{bmatrix}.

or

 
 A^{-1} =\frac{1}{4}
\begin{bmatrix}
3 & 2 & 1\\
2 & 4 & 2\\
1 & 2 & 3
\end{bmatrix}=\frac{1}{det(A)}
\begin{bmatrix}
3 & 2 & 1\\
2 & 4 & 2\\
1 & 2 & 3
\end{bmatrix}.

A matrix is non-singular (meaning that it has an inverse matrix) if and only if thee identity matrix can be obtained using only elementary row operations.

[ References

[ External links

Source: this wikipedia article, under CC-BY-SA.

Tutors Answer Your Questions about Linear Algebra (FREE)


Older solutions: 1..45, 46..90, 91..135, 136..180, 181..225, 226..270, 271..315, 316..360, 361..405, 406..450, 451..495, 496..540, 541..585, 586..630, 631..675, 676..720, 721..765, 766..810, 811..855, 856..900, 901..945, 946..990, 991..1035, 1036..1080, 1081..1125, 1126..1170, 1171..1215, 1216..1260, 1261..1305, 1306..1350, 1351..1395, 1396..1440, 1441..1485, 1486..1530, 1531..1575, 1576..1620, 1621..1665, 1666..1710, 1711..1755, 1756..1800, 1801..1845, 1846..1890