Descomposición lu
Descargar como DOCX, PDF2 recomendaciones8,336 vistas
La descomposición LU involucra descomponer una matriz A en dos matrices triangulares, una inferior L y una superior U, mediante operaciones sobre los coeficientes de A. Esto proporciona una forma eficiente de resolver sistemas de ecuaciones lineales. El proceso implica hacer ceros los valores por debajo y por encima de los pivotes en L y U respectivamente, a través de factores multiplicativos.
![DESCOMPOSICIÓN LU<br />Su nombre se deriva de las palabras inglesas “Lower\"
y “Upper”, que en español se traducen como “Inferior” y “Superior”. Estudiando el proceso que se sigue en la descomposición LU es posible comprender el porqué de este nombre, analizando cómo una matriz original se descompone en dos matrices triangulares, una superior y otra inferior.<br />La descomposición LU involucra solo operaciones sobre los coeficientes de la matriz [A], proporcionando un medio eficiente para calcular la matriz inversa o resolver sistemas de álgebra lineal.<br />Primeramente se debe obtener la matriz [L] y la matriz [U].<br />[L] es una matriz diagonal inferior con números 1 sobre la diagonal. [U] es una matriz diagonal superior en la que sobre la diagonal no necesariamente tiene que haber números 1.<br />El primer paso es descomponer o transformar [A] en [L] y [U], es decir obtener la matriz triangular inferior [L] y la matriz triangular superior [U].<br />PASOS PARA ENCONTRAR LA MATRIZ TRIANGULAR SUPERIOR (MATRIZ [U])<br />Hacer cero todos los valores abajo del pivote sin convertir este en 1.<br />Para lograr lo anterior se requiere obtener un factor el cual es necesario para convertir a cero los valores abajo del pivote.<br />Dicho factor es igual al número que se desea convertir en cero entre el número pivote.<br />Este factor multiplicado por -1 se multiplica luego por el pivote y a ese resultado se le suma el valor que se encuentra en la posición a cambiar (el valor en la posición que se convertirá en cero). <br />PASOS PARA ENCONTRAR LA MATRIZ TRIANGULAR INFERIOR (MATRIZ [L])<br />Para encontrar la matriz triangular inferior se busca hacer ceros los valores de arriba de cada pivote, así como también convertir en 1 cada pivote. Se utiliza el mismo concepto de “factor” explicado anteriormente y se ubican todos los “factores” debajo de la diagonal según corresponda en cada uno.<br />Esquemáticamente se busca lo siguiente:<br /> <br />Originalmente se tenía:<br />Debido a que [A] = [L][U], al encontrar [L] y [U] a partir de [A] no se altera en nada la ecuación y se tiene lo siguiente:<br />Por lo tanto, si Ax = b, entonces LUx = b, de manera que Ax = LUx = b.<br />PASOS PARA RESOLVER UN SISTEMA DE ECUACIONES POR EL MÉTODO DE DESCOMPOSICIÓN LU<br />Obtener la matriz triangular inferior L y la matriz triangular superior U.<br />Resolver Ly = b (para encontrar y).<br />El resultado del paso anterior se guarda en una matriz nueva de nombre “y”.<br />Realizar Ux = y (para encontrar x).<br />El resultado del paso anterior se almacena en una matriz nueva llamada “x”, la cual brinda los valores correspondientes a las incógnitas de la ecuación.<br />EJEMPLO 1 DE DESCOMPOSICIÓN LU<br />PROBLEMA: Encontrar los valores de x1, x2 y x3 para el siguiente sistema de ecuaciones:<br />NOTA: Recuérdese que si la matriz es 2x2 se hará 1 iteración; si es 3x3, 2 iteraciones; si es 4x4, 3 iteraciones; y así sucesivamente.<br />SOLUCIÓN:<br />4- 2- 19[A] =51- 1[B] =712- 412<br />ITERACIÓN 1<br />factor 1 = (a21 / a11) = 5 / 4 = 1.25<br />factor 2 = (a31 / a11) = 1 / 4 = 0.25<br />Encontrando [U]<br />fila 2 = - (factor 1) * (fila 1) + (fila 2)<br />fila 3 = - (factor 2) * (fila 1) + (fila 3)<br />a11 = a11<br />a12 = a12<br />a13 = a13<br />a21 = - (1.25) * (4) + (5) = 0<br />a22 = - (1.25) * (- 2) + (1) = 3.5<br />a23 = - (1.25) + (- 1) + (- 1) = 0.25<br />a31 = - (0.25) * (4) + (1) = 0<br />a32 = - (0.25) * (- 2) + (2) = 2.5<br />a33 = - (0.25) * (- 1) + (- 1) = - 0.75<br /> <br />4- 2- 1[U] =03.50.2502.5- 0.75<br />Encontrando [L]<br />100[L] =1.25000.2500<br />ITERACIÓN 2<br />factor 3 = (u32 / u22) = 2.5 / 3.5 = 0.7142857143<br />Encontrando [U]<br />fila 3 = - (factor 3) * (fila 2) + (fila 3)<br />a31 = - (2.5 / 3.5) * (0) + (0) = 0<br />a32 = - (2.5 / 3.5) * (3.5) + (2.5) = 0<br />a33 = - (2.5 / 3.5) * (0.25) + (- 0.75) = - 0.9285714286<br />4- 2- 1[U] =03.50.2500- 0.9285714286<br />Encontrando [L]<br />100[L] =1.25100.250.71428571431<br />Ahora ya se tiene la matriz [U] y la matriz [L]. El siguiente paso es resolver <br />Ly = b para encontrar la matriz y. En pocas palabras es como que se pidiera resolver el siguiente sistema de ecuaciones, encontrando los valores de y1, y2 y y3:<br />Al resolver el sistema anterior, se obtienen los siguientes valores para y1, y2 y y3:<br />El último paso es resolver Ux = y para encontrar la matriz x. En otras palabras es como que se pidiera resolver el siguiente sistema de ecuaciones, encontrando los valores de x1, x2 y x3:<br />La solución del sistema es:<br />Este es finalmente el valor de x1, x2 y x3; es decir, la respuesta del ejercicio utilizando la descomposición LU.<br />EJEMPLO 2 DE DESCOMPOSICIÓN LU<br />PROBLEMA: Encontrar los valores de x1, x2 y x3 para el siguiente sistema de ecuaciones:<br />SOLUCIÓN:<br />11- 3- 218[A] =5- 2- 8[B] =134- 722<br />ITERACIÓN 1<br />factor 1 = (a21 / a11) = 5/11 = 0.4545454545<br />factor 2 = (a31 / a11) = 4/11 = 0.3636363636<br />Encontrando [U]<br />fila 2 = - (factor 1) * (fila 1) + (fila 2)<br />fila 3 = - (factor 2) * (fila 1) + (fila 3)<br />a11 = a11<br />a12 = a12<br />a13 = a13<br />a21 = - (0.4545454545) * (11) + (5) = 0<br />a22 = - (0.4545454545) * (- 3) + (- 2) = - 0.6363636365<br />a23 = - (0.4545454545) + (- 2) + (- 8) = - 7.0909090919<br />a31 = - (0.3636363636) * (11) + (4) = 0<br />a32 = - (0.3636363636) * (- 3) + (- 7) = - 5.909090909<br />a33 = - (0.3636363636) * (- 2) + (2) = 2.7272727272<br />11-3-2[U] =0- 0.6363636365- 7.09090909190- 5.9090909092.7272727272<br />Encontrando [L]<br />100[L] =0.45454545000.3636363600<br />ITERACIÓN 2<br />factor 3 = (u32/u22) = - 5.909090909 / - 0.6363636365 = 9.285714284<br />Encontrando [U]<br />fila 3 = - (factor 3) * (fila 2) + (fila 3)<br />a31 = - (9.285714284) * (0) + (0) = 0<br />a32 = - (9.285714284) * (- 0.6363636365) + (- 5.909090909) = 0<br />a33 = - (9.285714284) * (- 7.0909090919) + (2.7272727272) = 68.57142857<br />11- 3- 2[U] =0- 0.6363636365- 7.09090909190068.57142857<br />Encontrando [L]<br />100[L] =0.4545454545100.36363636369.2857142841<br />Ahora ya se tiene la matriz [U] y la matriz [L]. El siguiente paso es resolver <br />Ly = b para encontrar la matriz y. En pocas palabras es como que se pidiera resolver el siguiente sistema de ecuaciones, encontrando los valores de y1, y2 y y3:<br />Al resolver el sistema anterior, se obtienen los siguientes valores para y1, y2 y y3:<br />El último paso es resolver Ux = y para encontrar la matriz x. En otras palabras es como que se pidiera resolver el siguiente sistema de ecuaciones, encontrando los valores de x1, x2 y x3:<br />La solución del sistema es:<br />Este es finalmente el valor de x1, x2 y x3; es decir, la respuesta del ejercicio utilizando la descomposición LU.<br />](https://image.slidesharecdn.com/descomposicinlu-101221215534-phpapp01/85/Descomposicion-lu-1-320.jpg)
Descomposición lu
- 1. DESCOMPOSICIÓN LU<br />Su nombre se deriva de las palabras inglesas “Lower\" y “Upper”, que en español se traducen como “Inferior” y “Superior”. Estudiando el proceso que se sigue en la descomposición LU es posible comprender el porqué de este nombre, analizando cómo una matriz original se descompone en dos matrices triangulares, una superior y otra inferior.<br />La descomposición LU involucra solo operaciones sobre los coeficientes de la matriz [A], proporcionando un medio eficiente para calcular la matriz inversa o resolver sistemas de álgebra lineal.<br />Primeramente se debe obtener la matriz [L] y la matriz [U].<br />[L] es una matriz diagonal inferior con números 1 sobre la diagonal. [U] es una matriz diagonal superior en la que sobre la diagonal no necesariamente tiene que haber números 1.<br />El primer paso es descomponer o transformar [A] en [L] y [U], es decir obtener la matriz triangular inferior [L] y la matriz triangular superior [U].<br />PASOS PARA ENCONTRAR LA MATRIZ TRIANGULAR SUPERIOR (MATRIZ [U])<br />Hacer cero todos los valores abajo del pivote sin convertir este en 1.<br />Para lograr lo anterior se requiere obtener un factor el cual es necesario para convertir a cero los valores abajo del pivote.<br />Dicho factor es igual al número que se desea convertir en cero entre el número pivote.<br />Este factor multiplicado por -1 se multiplica luego por el pivote y a ese resultado se le suma el valor que se encuentra en la posición a cambiar (el valor en la posición que se convertirá en cero). <br />PASOS PARA ENCONTRAR LA MATRIZ TRIANGULAR INFERIOR (MATRIZ [L])<br />Para encontrar la matriz triangular inferior se busca hacer ceros los valores de arriba de cada pivote, así como también convertir en 1 cada pivote. Se utiliza el mismo concepto de “factor” explicado anteriormente y se ubican todos los “factores” debajo de la diagonal según corresponda en cada uno.<br />Esquemáticamente se busca lo siguiente:<br /> <br />Originalmente se tenía:<br />Debido a que [A] = [L][U], al encontrar [L] y [U] a partir de [A] no se altera en nada la ecuación y se tiene lo siguiente:<br />Por lo tanto, si Ax = b, entonces LUx = b, de manera que Ax = LUx = b.<br />PASOS PARA RESOLVER UN SISTEMA DE ECUACIONES POR EL MÉTODO DE DESCOMPOSICIÓN LU<br />Obtener la matriz triangular inferior L y la matriz triangular superior U.<br />Resolver Ly = b (para encontrar y).<br />El resultado del paso anterior se guarda en una matriz nueva de nombre “y”.<br />Realizar Ux = y (para encontrar x).<br />El resultado del paso anterior se almacena en una matriz nueva llamada “x”, la cual brinda los valores correspondientes a las incógnitas de la ecuación.<br />EJEMPLO 1 DE DESCOMPOSICIÓN LU<br />PROBLEMA: Encontrar los valores de x1, x2 y x3 para el siguiente sistema de ecuaciones:<br />NOTA: Recuérdese que si la matriz es 2x2 se hará 1 iteración; si es 3x3, 2 iteraciones; si es 4x4, 3 iteraciones; y así sucesivamente.<br />SOLUCIÓN:<br />4- 2- 19[A] =51- 1[B] =712- 412<br />ITERACIÓN 1<br />factor 1 = (a21 / a11) = 5 / 4 = 1.25<br />factor 2 = (a31 / a11) = 1 / 4 = 0.25<br />Encontrando [U]<br />fila 2 = - (factor 1) * (fila 1) + (fila 2)<br />fila 3 = - (factor 2) * (fila 1) + (fila 3)<br />a11 = a11<br />a12 = a12<br />a13 = a13<br />a21 = - (1.25) * (4) + (5) = 0<br />a22 = - (1.25) * (- 2) + (1) = 3.5<br />a23 = - (1.25) + (- 1) + (- 1) = 0.25<br />a31 = - (0.25) * (4) + (1) = 0<br />a32 = - (0.25) * (- 2) + (2) = 2.5<br />a33 = - (0.25) * (- 1) + (- 1) = - 0.75<br /> <br />4- 2- 1[U] =03.50.2502.5- 0.75<br />Encontrando [L]<br />100[L] =1.25000.2500<br />ITERACIÓN 2<br />factor 3 = (u32 / u22) = 2.5 / 3.5 = 0.7142857143<br />Encontrando [U]<br />fila 3 = - (factor 3) * (fila 2) + (fila 3)<br />a31 = - (2.5 / 3.5) * (0) + (0) = 0<br />a32 = - (2.5 / 3.5) * (3.5) + (2.5) = 0<br />a33 = - (2.5 / 3.5) * (0.25) + (- 0.75) = - 0.9285714286<br />4- 2- 1[U] =03.50.2500- 0.9285714286<br />Encontrando [L]<br />100[L] =1.25100.250.71428571431<br />Ahora ya se tiene la matriz [U] y la matriz [L]. El siguiente paso es resolver <br />Ly = b para encontrar la matriz y. En pocas palabras es como que se pidiera resolver el siguiente sistema de ecuaciones, encontrando los valores de y1, y2 y y3:<br />Al resolver el sistema anterior, se obtienen los siguientes valores para y1, y2 y y3:<br />El último paso es resolver Ux = y para encontrar la matriz x. En otras palabras es como que se pidiera resolver el siguiente sistema de ecuaciones, encontrando los valores de x1, x2 y x3:<br />La solución del sistema es:<br />Este es finalmente el valor de x1, x2 y x3; es decir, la respuesta del ejercicio utilizando la descomposición LU.<br />EJEMPLO 2 DE DESCOMPOSICIÓN LU<br />PROBLEMA: Encontrar los valores de x1, x2 y x3 para el siguiente sistema de ecuaciones:<br />SOLUCIÓN:<br />11- 3- 218[A] =5- 2- 8[B] =134- 722<br />ITERACIÓN 1<br />factor 1 = (a21 / a11) = 5/11 = 0.4545454545<br />factor 2 = (a31 / a11) = 4/11 = 0.3636363636<br />Encontrando [U]<br />fila 2 = - (factor 1) * (fila 1) + (fila 2)<br />fila 3 = - (factor 2) * (fila 1) + (fila 3)<br />a11 = a11<br />a12 = a12<br />a13 = a13<br />a21 = - (0.4545454545) * (11) + (5) = 0<br />a22 = - (0.4545454545) * (- 3) + (- 2) = - 0.6363636365<br />a23 = - (0.4545454545) + (- 2) + (- 8) = - 7.0909090919<br />a31 = - (0.3636363636) * (11) + (4) = 0<br />a32 = - (0.3636363636) * (- 3) + (- 7) = - 5.909090909<br />a33 = - (0.3636363636) * (- 2) + (2) = 2.7272727272<br />11-3-2[U] =0- 0.6363636365- 7.09090909190- 5.9090909092.7272727272<br />Encontrando [L]<br />100[L] =0.45454545000.3636363600<br />ITERACIÓN 2<br />factor 3 = (u32/u22) = - 5.909090909 / - 0.6363636365 = 9.285714284<br />Encontrando [U]<br />fila 3 = - (factor 3) * (fila 2) + (fila 3)<br />a31 = - (9.285714284) * (0) + (0) = 0<br />a32 = - (9.285714284) * (- 0.6363636365) + (- 5.909090909) = 0<br />a33 = - (9.285714284) * (- 7.0909090919) + (2.7272727272) = 68.57142857<br />11- 3- 2[U] =0- 0.6363636365- 7.09090909190068.57142857<br />Encontrando [L]<br />100[L] =0.4545454545100.36363636369.2857142841<br />Ahora ya se tiene la matriz [U] y la matriz [L]. El siguiente paso es resolver <br />Ly = b para encontrar la matriz y. En pocas palabras es como que se pidiera resolver el siguiente sistema de ecuaciones, encontrando los valores de y1, y2 y y3:<br />Al resolver el sistema anterior, se obtienen los siguientes valores para y1, y2 y y3:<br />El último paso es resolver Ux = y para encontrar la matriz x. En otras palabras es como que se pidiera resolver el siguiente sistema de ecuaciones, encontrando los valores de x1, x2 y x3:<br />La solución del sistema es:<br />Este es finalmente el valor de x1, x2 y x3; es decir, la respuesta del ejercicio utilizando la descomposición LU.<br />