jueves, 20 de noviembre de 2014

4.1.5 Eliminación de Gauss y Gauss-Jordan

En matemáticas, la eliminación de Gauss-Jordan, llamada así debido a Carl Friedrich Gauss y Wilhelm Jordan, es un algoritmo del álgebra lineal para determinar las soluciones de un sistema de ecuaciones lineales, encontrar matrices e inversas. Un sistema de ecuaciones se resuelve por el método de Gauss cuando se obtienen sus soluciones mediante la reducción del sistema dado a otro equivalente en el que cada ecuación tiene una incógnita menos que la anterior. El método de Gauss transforma la matriz de coeficientes en una matriz triangular superior. El método de Gauss-Jordan continúa el proceso de transformación hasta obtener una matriz diagonal.


El método de Gauss consiste en transformar un sistema de ecuaciones en otro equivalente de forma que este sea escalonado.
Para facilitar el cálculo vamos a transformar el sistema en una matriz, en la que pondremos los coeficientes de las variables y los términos independientes (separados por una recta).
matriz matriz

Ejemplos

 1 
3x+ 2y+ z=1cierre
5x+ 3y+ 4z=2
x+ y− z=1
Primer paso
solución
solución
solución
solución
solución

Método de Gauss-Jordan

Como hemos visto, el método de Gauss transforma la matriz de coeficientes en una matriz triangular superior. El método de Gauss-Jordan continúa el proceso de transformación hasta obtener una matriz diagonal unitaria (aij=0 para cualquier $i \neq j$).
Veamos el método de Gauss-Jordan siguiendo con el ejemplo empleado en el apartado anterior. Aplicando el método de Gauss habíamos llegado a la siguiente ecuación: 
\begin{displaymath}\begin{array}{llll}
\left(
\begin{array}{rrrr}
6 & -2 & 2 ...
...ray}{r}
12 \\ 10 \\ -9 \\ -3
\end{array} \right)
\end{array}\end{displaymath}

Ahora seguiremos un procedimiento similar al empleado en el método de Gauss. Tomaremos como pivote el elemento a44=-3; multiplicamos la cuarta ecuación por $\frac{-3}{4}$ y la restamos a la primera: 
\begin{displaymath}\begin{array}{llll}
\left(
\begin{array}{rrrr}
6 & -2 & 2 ...
...rray}{r}
8 \\ 10 \\ -9 \\ -3
\end{array} \right)
\end{array}\end{displaymath}

Realizamos la misma operación con la segunda y tercera fila, obteniendo: 
\begin{displaymath}\begin{array}{llll}
\left(
\begin{array}{rrrr}
6 & -2 & 2 ...
...array}{r}
8 \\ 8 \\ -4 \\ -3
\end{array} \right)
\end{array}\end{displaymath}

Ahora tomamos como pivote el elemento a33=2, multiplicamos la tercera ecuación por $\frac{2}{2}=1$ y la restamos a la primera: 
\begin{displaymath}\begin{array}{llll}
\left(
\begin{array}{rrrr}
6 & -2 & 0 ...
...rray}{r}
12 \\ 8 \\ -4 \\ -3
\end{array} \right)
\end{array}\end{displaymath}

Repetimos la operación con la segunda fila: 
\begin{displaymath}\begin{array}{llll}
\left(
\begin{array}{rrrr}
6 & -2 & 0 ...
...ray}{r}
12 \\ 12 \\ -4 \\ -3
\end{array} \right)
\end{array}\end{displaymath}

Finalmente, tomamos como pivote a22=-4, multiplicamos la segunda ecuación por $\frac{-2}{-4}$ y la sumamos a la primera: 
\begin{displaymath}\begin{array}{llll}
\left(
\begin{array}{rrrr}
6 & 0 & 0 &...
...rray}{r}
6 \\ 12 \\ -4 \\ -3
\end{array} \right)
\end{array}\end{displaymath}

El sistema de ecuaciones anterior es, como hemos visto, fácil de resolver. Empleando la ecuación (46) obtenemos las soluciones: 
\begin{displaymath}x = \left( \begin{array}{r} 1 \\ -3 \\ -2 \\ 1 \end{array} \right)
\end{displaymath}






No hay comentarios:

Publicar un comentario