Calculadora de Descomposición LU de Matriz

Bienvenido a la Calculadora de Descomposición LU de Matriz, una herramienta integral de álgebra lineal que factoriza cualquier matriz cuadrada en el producto de una matriz triangular inferior (L) y una matriz triangular superior (U) mediante la eliminación gaussiana con pivoteo parcial. Obtenga una eliminación detallada paso a paso, animación interactiva de la descomposición y verificación automática. Ideal para estudiantes, ingenieros y cualquier persona que trabaje con sistemas de ecuaciones lineales.

¿Qué es la descomposición LU?

La descomposición LU (también llamada factorización LU) expresa una matriz cuadrada \(A\) como el producto de dos matrices triangulares:

Donde:

• L (Triangular inferior): tiene unos en la diagonal y entradas distintas de cero solo debajo de la diagonal. Estas entradas son los multiplicadores utilizados durante la eliminación gaussiana.
• U (Triangular superior): tiene entradas distintas de cero solo en y por encima de la diagonal. Esta es la forma escalonada por filas de la matriz.

Cuando se utiliza pivoteo parcial (para evitar pivotes cero y mejorar la estabilidad numérica), la factorización se convierte en:

Donde \(P\) es una matriz de permutación que registra los intercambios de filas realizados durante la eliminación.

Cómo usar esta calculadora

1. Ingrese su matriz: Introduzca una matriz cuadrada con las filas en líneas separadas o separadas por puntos y comas. Los elementos pueden separarse por espacios, comas o pestañas. Admite hasta 8×8.
2. Establezca la precisión decimal: Elija cuántos decimales (2-10) desea mostrar en los resultados.
3. Haga clic en Descomponer: La calculadora realiza la factorización LU con pivoteo parcial y muestra los resultados.
4. Revise los resultados: Examine las matrices L, U y P, la descomposición animada y los detalles de la eliminación paso a paso.

Resolución de sistemas lineales con descomposición LU

La descomposición LU es particularmente potente para resolver sistemas de ecuaciones lineales \(Ax = b\). Una vez que se tiene \(PA = LU\), la resolución se convierte en un proceso de dos pasos:

Paso 1: Sustitución hacia adelante

Resuelva \(Ly = Pb\) para \(y\). Dado que \(L\) es triangular inferior, esto es directo: comience desde la ecuación superior y avance hacia abajo:

Paso 2: Sustitución hacia atrás

Resuelva \(Ux = y\) para \(x\). Dado que \(U\) es triangular superior, comience desde la ecuación inferior y avance hacia arriba:

Cálculo del determinante

El determinante de \(A\) se puede calcular eficientemente a partir de los factores LU:

Donde \(s\) es el número de intercambios de filas (pivotes) y \(U_{ii}\) son las entradas de la diagonal de \(U\). Dado que \(\det(L) = 1\) (todas las entradas de la diagonal son 1) y \(\det(P) = (-1)^s\), la fórmula se deduce de \(\det(P)\det(A) = \det(L)\det(U)\).

¿Por qué el pivoteo parcial?

Sin pivoteo, la descomposición LU falla si algún elemento pivote es cero. Incluso cuando los pivotes son distintos de cero pero pequeños, el resultado calculado puede sufrir errores numéricos severos. El pivoteo parcial selecciona el pivote más grande disponible en cada columna, lo que:

• Evita la división por cero
• Minimiza el crecimiento de errores de redondeo
• Garantiza que los multiplicadores en L satisfagan \(|L_{ij}| \leq 1\)
• Asegura que toda matriz no singular pueda ser descompuesta

Aplicaciones de la descomposición LU

LU frente a otras descomposiciones

• LU vs QR: LU es más rápido (\(O(\frac{2}{3}n^3)\) frente a \(O(\frac{4}{3}n^3)\)) pero menos estable numéricamente. Se prefiere QR para problemas de mínimos cuadrados.
• LU vs Cholesky: Cholesky (\(A = LL^T\)) funciona solo para matrices simétricas definidas positivas, pero es el doble de rápido y más estable que la LU general.
• LU vs Eliminación gaussiana: LU es eliminación gaussiana, pero la forma factorizada L y U se puede reutilizar para resolver múltiples términos independientes de manera eficiente.

Preguntas frecuentes

¿Qué es la descomposición LU?

La descomposición LU (también llamada factorización LU) es un método que factoriza una matriz cuadrada A en el producto de una matriz triangular inferior L y una matriz triangular superior U, de modo que A = LU (o PA = LU con pivoteo parcial). La matriz L tiene unos en la diagonal y almacena los multiplicadores de eliminación, mientras que U es el resultado de la eliminación gaussiana.

¿Por qué es necesario el pivoteo parcial en la descomposición LU?

El pivoteo parcial intercambia filas para colocar el valor absoluto más grande en la posición del pivote. Esto evita la división por cero cuando un elemento pivote es cero y reduce los errores numéricos causados por la división por números pequeños. Con el pivoteo parcial, la factorización se convierte en PA = LU, donde P es una matriz de permutación que registra los intercambios de filas.

¿Cuáles son las aplicaciones de la descomposición LU?

La descomposición LU se utiliza para resolver sistemas de ecuaciones lineales (Ax = b) de manera eficiente, calcular determinantes de matrices, encontrar matrices inversas y analizar la estabilidad numérica. Es especialmente eficiente al resolver múltiples sistemas con la misma matriz de coeficientes pero diferentes términos independientes, ya que la factorización solo debe realizarse una vez.

¿Cómo se resuelve Ax = b usando la descomposición LU?

Después de calcular PA = LU, resolver Ax = b consiste en: primero resolver Ly = Pb mediante sustitución hacia adelante (fácil porque L es triangular inferior), luego resolver Ux = y mediante sustitución hacia atrás (fácil porque U es triangular superior). Este proceso de dos pasos es mucho más rápido que la eliminación gaussiana cuando se resuelven múltiples sistemas.

¿Se puede descomponer por LU cualquier matriz cuadrada?

No todas las matrices cuadradas tienen una descomposición LU sin pivoteo. Una matriz tiene una factorización LU si y solo si todos sus menores principales líderes son distintos de cero. Sin embargo, con pivoteo parcial (PA = LU), toda matriz cuadrada no singular puede descomponerse. Las matrices singulares pueden fallar si se encuentra un pivote cero.

Recursos adicionales

• Descomposición LU - Wikipedia
• Matriz triangular - Wikipedia
• Eliminación gaussiana - Wikipedia