Sistemas Escalonados y el Efecto Dominó: El arte de la Sustitución Hacia Atrás
En las lecciones anteriores aprendimos a analizar el destino de un sistema (Discusión) y a manipular sus ecuaciones mediante las operaciones de Gauss. Vimos que el objetivo de Gauss es limpiar la matriz, haciendo ceros por debajo de la diagonal principal para construir una **estructura en escalera**.
Pero, ¿qué ocurre una vez que la matriz ha alcanzado esa codiciada forma escalonada? Ocurre la fase más gratificante de todo el álgebra lineal: la **Sustitución Hacia Atrás** (o Back-Substitution).
Si la eliminación de Gauss es el duro trabajo de construcción, la sustitución hacia atrás es el **efecto dominó**: la última ecuación te entrega el valor de la primera incógnita sin esfuerzo, y ese dato cae hacia arriba resolviendo, una a una y en cascada, todas las demás incógnitas del problema.
1. La Anatomía de la Matriz Escalonada
Recordemos qué aspecto tiene una matriz ampliada (A|B) una vez que hemos completado el escalonamiento de Gauss. Todos los elementos por debajo de los pivotes principales se han convertido en ceros:
( 0 a22 a23 | b2 )
( 0 0 a33 | b3 )
Si traducimos esta matriz de nuevo al lenguaje de ecuaciones explícitas, la magia salta a la vista inmediatamente:
(Fila 1) a11x + a12y + a13z = b1
(Fila 2) a22y + a23z = b2
(Fila 3) a33z = b3
Observa la Fila 3. Como los coeficientes de x e y son ceros, la variable z ha quedado **completamente aislada**. Ya no depende de las demás variables. La solución de la última incógnita está a un solo paso aritmético de distancia.
2. El Algoritmo del Dominó Paso a Paso
El proceso de sustitución hacia atrás sigue un orden estricto de abajo hacia arriba (desde la fila m hasta la fila 1):
Fase 1: Despeje del Pivote Inferior
Tomamos la última fila no nula y despejamos la última incógnita principal. Como a33 ≠ 0 (es un pivote válido), simplemente dividimos entre él:
z = b3 / a33
Fase 2: Propagación Ascendente
Tomamos el número concreto que acabamos de obtener para z y lo **inyectamos** en la Fila 2. Al sustituir z por su valor numérico, la Fila 2 pasa de tener dos incógnitas a tener **una sola incógnita (y)**. Despejamos y:
a22y + a23(z_conocido) = b2 ➔ y = [b2 - a23(z_conocido)] / a22
Fase 3: Cúspide del Sistema
Con los valores numéricos de y y z ya calculados, subimos a la Fila 1 y sustituimos ambas variables. La primera ecuación queda reducida a una simple aritmética de primer grado para despejar x.
¿Por qué la informática ama la Sustitución Hacia Atrás?
En la física computacional, la eficiencia lo es todo. Eliminar números con Gauss requiere un esfuerzo de cálculo proporcional a n3 operaciones (donde n es el número de ecuaciones). Sin embargo, una vez que el sistema está escalonado, la sustitución hacia atrás sólo requiere n2 operaciones.
Es un proceso **ultra-rápido**: requiere poquísima memoria RAM y un tiempo de cómputo despreciable para las supercomputadoras. Por eso, casi todos los solvers numéricos del mundo triangulan primero y sustituyen después.
3. Ejemplo Numérico Desarrollado en Directo
Vamos a resolver paso a paso un sistema que ya ha sido escalonado por Gauss para experimentar el efecto dominó en directo.
Supongamos que tras aplicar Gauss a un problema de ingeniería, obtenemos la siguiente matriz ampliada (A|B):
( 0 -2 4 | 6 )
( 0 0 3 | 9 )
Paso 1: Fila 3 (La base del dominó)
Escribimos la tercera ecuación: 3z = 9.
z = 9 / 3 ➔ z = 3
Paso 2: Fila 2 (Primera caída)
Escribimos la segunda ecuación: -2y + 4z = 6. Inyectamos z = 3:
-2y + 4(3) = 6 ➔ -2y + 12 = 6 ➔ -2y = 6 - 12 ➔ -2y = -6
y = -6 / -2 ➔ y = 3
Paso 3: Fila 1 (Cúspide final)
Escribimos la primera ecuación: 2x + 3y - z = 5. Inyectamos y = 3 y z = 3:
2x + 3(3) - (3) = 5 ➔ 2x + 9 - 3 = 5 ➔ 2x + 6 = 5 ➔ 2x = 5 - 6 ➔ 2x = -1
x = -1 / 2 ➔ x = -0.5
Vector Solución Único (S.C.D.):
(x, y, z) = (-0.5, 3, 3)
4. La analogía física: Cadenas de Tensión en Equilibrio
En la física del mundo real, la sustitución hacia atrás se vive constantemente. Imagina una cadena suspendida del techo con 3 eslabones pesados sujetos unos a otros.
- Para saber la tensión que soporta el **último eslabón (el de abajo del todo)**, solo necesitas medir la fuerza del peso colgado de él (Fila 3, solución directa).
- Para calcular la tensión del **eslabón intermedio**, necesitas sumar su propio peso más la tensión que ya calculaste del eslabón inferior (Fila 2, inyección de datos).
- Para calcular la tensión del **primer eslabón (el sujeto al techo)**, debes sumar su peso a la tensión acumulada de todos los eslabones inferiores (Fila 1, cúspide final).
La naturaleza resuelve la distribución de fuerzas de abajo hacia arriba de la misma manera que el método de Gauss. En la siguiente lección profundizaremos en la gestión de variables cuando el sistema no está completamente cerrado: **Jerarquía de variables: Incógnitas Principales vs. Incógnitas Secundarias (Parámetros)**.
Todas las Etiquetas:
Comentarios
Todavía no hay comentarios. ¡Sé el primero en comentar!