En matemáticas, la eliminación gaussiana (también llamada reducción de filas) es un método sistemático para resolver sistemas de ecuaciones lineales. Lleva el nombre de Carl Friedrich Gauss, un famoso matemático alemán que trabajó y describió el método, aunque no fue su descubridor original.
¿Qué se hace y por qué?
La idea principal es transformar el sistema en una matriz aumentada (con los coeficientes de las incógnitas y los términos independientes) y aplicar operaciones elementales sobre las filas para llevar dicha matriz a una forma más sencilla desde la que se pueda leer la solución del sistema o detectar si no existe o es indeterminado.
Operaciones elementales de fila
Las operaciones que se permiten son tres y preservan el conjunto de soluciones del sistema:
- Tipo 1: Intercambiar dos filas.
- Tipo 2: Multiplicar una fila por un escalar distinto de cero.
- Tipo 3: Sumar a una fila un múltiplo de otra fila.
Formas escalonadas
El objetivo usual es obtener la matriz en forma de fila-echelón (forma escalonada). En esa forma, cada fila no nula comienza con más ceros iniciales que la fila anterior; además los pivotes (primeros elementos no nulos de cada fila) están a la derecha de los pivotes de las filas superiores. Si, además, cada pivote es 1 y es el único elemento no nulo en su columna, la matriz está en forma reducida de fila-echelón (reduced row-echelon form). El proceso que conduce directamente a la forma reducida suele llamarse eliminación de Gauss–Jordan.
Algoritmo paso a paso (eliminación hacia adelante y sustitución hacia atrás)
Un procedimiento típico de eliminación gaussiana consta de dos fases:
- Eliminación hacia adelante: convertir la matriz aumentada en una matriz de fila-echelón usando operaciones elementales. Esto implica seleccionar un pivote (idealmente distinto de cero), usarlo para anular los elementos por debajo en su columna y continuar con la submatriz inferior.
- Sustitución hacia atrás: una vez que la matriz está en forma de fila-echelón, resolver las incógnitas empezando por la última ecuación no nula (la que tiene menos incógnitas) y retrocediendo hacia arriba.
Para mejorar la estabilidad numérica, es habitual aplicar pivoteo parcial (intercambiar filas para colocar en el pivote el elemento de mayor valor absoluto de la columna) o pivoteo total (intercambiar filas y columnas). El coste computacional de la eliminación gaussiana es aproximadamente O(n³) para un sistema de n ecuaciones con n incógnitas.
Ejemplo paso a paso
Resolver el sistema
x + 2y − z = 3
2x + 3y + z = 7
−x + y + 2z = −4
Matriz aumentada inicial:
[ 1 2 −1 | 3 ]
[ 2 3 1 | 7 ]
[−1 1 2 | −4 ]
1) Eliminar debajo del primer pivote (1):
R2 ← R2 − 2·R1 ⇒ [0 −1 3 | 1]
R3 ← R3 + 1·R1 ⇒ [0 3 1 | −1]
Ahora la matriz es:
[ 1 2 −1 | 3 ]
[ 0 −1 3 | 1 ]
[ 0 3 1 | −1 ]
2) Usamos el pivote de la segunda fila (−1) y eliminamos debajo:
R3 ← R3 + 3·R2 ⇒ [0 0 10 | 2]
La matriz escalonada queda:
[ 1 2 −1 | 3 ]
[ 0 −1 3 | 1 ]
[ 0 0 10 | 2 ]
3) Escalamos la tercera fila para tener pivote 1:
R3 ← (1/10)·R3 ⇒ [0 0 1 | 0.2]
4) Sustitución hacia atrás (o eliminación hacia arriba):
De R3: z = 0.2 = 1/5.
De R2: −y + 3z = 1 ⇒ −y + 3·0.2 = 1 ⇒ y = −0.4 = −2/5.
De R1: x + 2y − z = 3 ⇒ x + 2(−0.4) − 0.2 = 3 ⇒ x = 4.
Solución: x = 4, y = −2/5, z = 1/5.
Casos especiales y relaciones con el rango
- Si durante la eliminación aparece una fila de la forma [0 0 ... 0 | b] con b ≠ 0, el sistema es inconsistente (sin solución).
- Si el número de pivotes (el rango de la matriz de coeficientes) es menor que el número de incógnitas y no hay contradicción, el sistema tiene infinitas soluciones; quedan variables libres que se parametrizan.
- El teorema de Rouché–Capelli relaciona rango de la matriz de coeficientes y de la aumentada para determinar existencia y cantidad de soluciones.
Aplicaciones y consideraciones numéricas
- La eliminación gaussiana sirve también para calcular inversas de matrices (aplicando el proceso a la matriz aumentada [A | I]) y como paso clave para obtener la factorización LU (A = LU).
- En implementación numérica, el pivoteo parcial suele ser suficiente para garantizar estabilidad en la mayoría de los casos; sin embargo, problemas mal condicionados pueden amplificar errores y requerir técnicas más robustas o aritmética de mayor precisión.
- El método es fundamental en álgebra lineal computacional y está implementado en todas las bibliotecas numéricas estándar.
En resumen, la eliminación gaussiana es una herramienta sistemática y versátil para tratar sistemas lineales: permite encontrar soluciones únicas, detectar sistemas incompatibles o describir soluciones infinitas, y sirve como base para muchas otras técnicas en álgebra lineal y cálculo numérico.