Cálculo de la inversa (Método de Gauss-Jordan)

Cálculo de la Matriz Inversa. El método del Espejo Mágico de Gauss-Jordan

En las lecciones anteriores exploramos el concepto teórico de la Matriz Inversa (A^-1) y comprendimos geométricamente por qué sólo las matrices de rango máximo (no singulares) poseen esa "máquina del tiempo" capaz de revertir las transformaciones del espacio.

Entender la teoría está bien, pero en un examen universitario o en el día a día de la física computacional necesitamos un procedimiento práctico, ágil y sistemático para calcular la matriz inversa sin perdernos en un mar de números. Aunque existen varios caminos (como el uso de determinantes y adjuntos que estudiaremos más adelante), el método más elegante, directo y preferido por los algoritmos informáticos es el **Método de Gauss-Jordan con la Matriz Ampliada**.

1. El Principio del Espejo: La Matriz Ampliada (A | I_n)

¿En qué consiste el truco de Gauss-Jordan? La idea es de una belleza aplastante. Construimos una super-matriz "ampliada" pegando dos bloques lado a lado:

• En el **bloque izquierdo** colocamos la matriz cuadrada A que queremos invertir (de tamaño n × n).
• En el **bloque derecho** colocamos la Matriz Identidad I_n del mismo tamaño (llena de unos en la diagonal y ceros en el resto).

Matriz Ampliada Inicial = ( A | I_n )

A partir de este momento, aplicamos las **operaciones elementales por filas** (Tipo I, Tipo II y Tipo III) sobre la matriz entera, afectando tanto al lado izquierdo como al lado derecho. Nuestro objetivo es transformar el bloque izquierdo hasta convertirlo en la Matriz Identidad I_n (alcanzando la Forma Canónica de Hermite por Filas).

2. El Chivato de Matrices Singulares

¿Qué ocurre si intentas aplicar este método a una matriz que **no tiene inversa** (una matriz singular)?

El método de Gauss-Jordan incluye un sistema de autoprotección infalible: si durante el proceso de escalonamiento en el bloque izquierdo **una fila entera se convierte en ceros (0 0 ... 0)**, el algoritmo se detiene inmediatamente. Esa fila de ceros es la prueba irrefutable de que rg(A) < n. No es posible formar la Matriz Identidad en el lado izquierdo y concluyes, sin lugar a dudas, que **la matriz A no es invertible**.

3. Ejemplo Desarrollado Paso a Paso (Matriz 3×3)

Vamos a calcular la inversa de una matriz de orden 3 con números limpios, mostrando cada operación elemental para que veas el baile de filas en directo.

Queremos calcular la inversa de la siguiente matriz A:

Paso 1: Construimos la Matriz Ampliada ( A | I₃ )

Paso 2: Hacemos ceros en la primera columna por debajo del primer pivote

Usamos el 1 de la esquina superior izquierda para limpiar su columna en las filas 2 y 3:

• A la Fila 2 le restamos el doble de la Fila 1 (F₂ → F₂ - 2F₁)
• A la Fila 3 le restamos el cuádruple de la Fila 1 (F₃ → F₃ - 4F₁)

Paso 3: Cambiamos el signo de la Fila 2 para tener un pivote igual a 1

Multiplicamos toda la Fila 2 por -1 (F₂ → -1 · F₂):

Paso 4: Eliminamos el elemento debajo del segundo pivote

Restamos la Fila 2 a la Fila 3 (F₃ → F₃ - F₂):

Paso 5: Preparamos el tercer pivote y limpiamos hacia arriba (Jordan)

Multiplicamos la Fila 3 por -1 para fijar su pivote como +1 (F₃ → -1 · F₃):

Ahora limpiamos la tercera columna de las filas superiores usando el pivote de la Fila 3:

• A la Fila 2 le restamos la Fila 3 (F₂ → F₂ - F₃)
• A la Fila 1 le restamos el doble de la Fila 3 (F₁ → F₁ - 2F₃)

¡Objetivo cumplido! El bloque izquierdo se ha transformado en la Matriz Identidad I₃. Por consiguiente, el bloque derecho contiene, sin lugar a dudas, la **Matriz Inversa A^-1**:

El método de Gauss-Jordan nos da el control absoluto sobre la inversión de matrices. En nuestro próximo capítulo entraremos en la **FASE 6: El Determinante**, donde descubriremos ese "número mágico" que mide cuánto se estira el universo y nos dirá, en un solo segundo, si una matriz tiene inversa o es singular.