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:

  1. 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.
  2. 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.