{"id":217,"date":"2025-10-08T08:15:35","date_gmt":"2025-10-08T06:15:35","guid":{"rendered":"https:\/\/clases.jesussoto.es\/?p=217"},"modified":"2025-10-03T20:01:00","modified_gmt":"2025-10-03T18:01:00","slug":"alg-factorizacion-lu","status":"publish","type":"post","link":"https:\/\/clases.jesussoto.es\/?p=217","title":{"rendered":"ALG: Factorizaci\u00f3n LU"},"content":{"rendered":"<p>La factorizaci\u00f3n <strong>LU <\/strong>es una forma de factorizaci\u00f3n de una matriz como el producto de una matriz triangular inferior y una superior.<\/p>\n<p>El prop\u00f3sito es dada una matriz \\(A\\) conseguir descomponer esta en un producto \\[\\textbf{A}=\\textbf{L}\\textbf{U}\\] de manera \\(\\textbf{L}\\) sea triangular inferior y \\(\\textbf{U}\\) triangular superior.<\/p>\n<p>Recordad que una operaci\u00f3n elemental entre filas se puede considerar como una matriz. De esta manera Podemos realizar una serie de operaciones elementales entre filas para transformar la matriz de partida \\(A\\), en una matriz escalonada (triangular superior). Es decir;<br \/>\n\\[E_k\\,E_{k-1}\\cdots E_1\\, \\textbf{A}= \\textbf{U}.\\]<br \/>\nComo cada matriz \\(E_i\\) es regular (por sus propiedades), entonces:<br \/>\n\\[\\textbf{A}= E_1^{-1}\\, E_2^{-1}\\cdots E_k^{-1}\\, \\textbf{U}\\]<br \/>\nAs\u00ed ser\u00e1 \\[\\textbf{L}=E_1^{-1}\\, E_2^{-1}\\cdots E_k^{-1}.\\]<br \/>\nEs f\u00e1cil comprobar que \\(L\\) es triangular inferior. El prop\u00f3sito es que la diagonal principal de \\(L\\) ser\u00e1n todo unos; as\u00ed<br \/>\n\\[|\\textbf{A}|=|\\textbf{L}\\textbf{U}|=|\\textbf{L}|\\cdot|\\textbf{U}|=\\prod_{i=1}^nu_{ii}\\]<\/p>\n<blockquote>\n<p><strong>Ejercicio:<\/strong> Dada la matriz \\[A=\\begin{bmatrix}<br \/>\n1 &#038; 2 &#038; 3\\\\<br \/>\n-3 &#038; -4 &#038; 13 \\\\<br \/>\n 2 &#038; 1 &#038; -5<br \/>\n\\end{bmatrix}\\in\\mathcal{M}_3(\\mathbb{R}),\\]  \u00bfcu\u00e1l es la traza de la matriz \\(\\textbf{U}\\) resultado de su factorizaci\u00f3n \\(\\textbf{LU}\\)? <\/p>\n<\/blockquote>\n<p><script>\nfunction showHtmlDiv4() {\n  var htmlShow4 = document.getElementById(\"html-show4\");\n  if (htmlShow4.style.display === \"none\") {\n    htmlShow4.style.display = \"block\";\n  } else {\n    htmlShow4.style.display = \"none\";\n  }\n}\n<\/script> <\/p>\n<p><button onclick=\"showHtmlDiv4()\">Soluci\u00f3n:<\/button><\/p>\n<div id=\"html-show4\" style=\"display: none;\">\n<iframe loading=\"lazy\" title=\"\u00c1lgebra Lineal - Factorizaci\u00f3n de matrices (LU) EJ.1 - Jes\u00fas Soto\" width=\"640\" height=\"360\" src=\"https:\/\/www.youtube.com\/embed\/9oJsBnttdDE?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe>\n<\/div>\n<hr \/>\n<blockquote>\n<p><strong>Ejercicio:<\/strong> Dada la matriz \\[A=\\begin{bmatrix}<br \/>\n2 &#038; -4 &#038; 3\\\\<br \/>\n6 &#038; -8 &#038; 5 \\\\<br \/>\n6 &#038; 1 &#038; 7<br \/>\n\\end{bmatrix}\\in\\mathcal{M}_3(\\mathbb{R}),\\]  \u00bfcu\u00e1nto suman los elementos de la matriz \\(\\textbf{L}\\) resultado de su factorizaci\u00f3n \\(\\textbf{LU}\\)? <\/p>\n<\/blockquote>\n<p><script>\nfunction showHtmlDiv41() {\n  var htmlShow41 = document.getElementById(\"html-show41\");\n  if (htmlShow41.style.display === \"none\") {\n    htmlShow41.style.display = \"block\";\n  } else {\n    htmlShow41.style.display = \"none\";\n  }\n}\n<\/script> <\/p>\n<p><button onclick=\"showHtmlDiv41()\">Soluci\u00f3n:<\/button><\/p>\n<div id=\"html-show41\" style=\"display: none;\">\n<iframe loading=\"lazy\" title=\"\u00c1lgebra lineal - Factorizaci\u00f3n (LU) EJ. 2 - Jes\u00fas Soto\" width=\"640\" height=\"360\" src=\"https:\/\/www.youtube.com\/embed\/mbevLrcbHns?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe>\n<\/div>\n<hr \/>\n<blockquote>\n<p><strong>Teorema<\/strong>: Para que la factorizaci\u00f3n LU de una matriz \\(A\\) sea \u00fanica, se deben cumplir las siguientes condiciones<\/p>\n<ol>\n<li><strong>Matriz no singular<\/strong>: La matriz \\(A\\) debe ser invertible, es decir, su determinante debe ser distinto de cero.<\/li>\n<li><strong>Pivotes no nulos<\/strong>: En el proceso de eliminaci\u00f3n gaussiana, todos los pivotes (elementos diagonales de la matriz \\(U\\)) deben ser distintos de cero.<\/li>\n<\/ol>\n<\/blockquote>\n<h2>Sistemas<\/h2>\n<p>Uno de los usos est\u00e1 en la posibilidad de resolver sistemas de ecuaciones. Consideremos queremos resolver el sistema de ecuaciones \\[\\textbf{A}x=\\textbf{b},\\] donde \\(\\textbf{A}\\in\\mathcal{M}_{n\\times m}(\\mathbb{K})\\). Si conseguimos una factorizaci\u00f3n \\[\\textbf{A}=\\textbf{L}\\textbf{U},\\] donde \\(\\textbf{L}\\in\\mathcal{M}_{n\\times n}(\\mathbb{K})\\), y, \\(\\textbf{U}\\in\\mathcal{M}_{n\\times m}(\\mathbb{K})\\), resultar\u00e1<\/p>\n<p>\\[\\textbf{A}x=(\\textbf{L}\\textbf{U})x=\\textbf{L}(\\textbf{U}x)=\\textbf{b}.\\]<br \/>\nPara resolver el problema podemos afrontar la estrategia de resolver primero:<br \/>\n\\[\\textbf{L}y=\\textbf{b},\\] para despu\u00e9s<br \/>\n\\[\\textbf{U}x=\\textbf{y}.\\]<\/p>\n<p>Como ambas matrices \\(\\textbf{L}\\) y \\(\\textbf{U}\\) son triangulares su soluci\u00f3n es f\u00e1cil mediante sustituci\u00f3n.<\/p>\n<blockquote>\n<p><strong>Ejercicio:<\/strong> Encontrar la factorizaci\u00f3n \\(\\textbf{LU}\\) que permita resolver el sistema \\(Ax=b\\) donde \\[A=\\begin{bmatrix}<br \/>\n3 &#038; -7 &#038; -2 &#038; 2\\\\<br \/>\n-3 &#038; 5 &#038; 1 &#038; 0 \\\\<br \/>\n 6 &#038; -4 &#038; 0 &#038; -5\\\\<br \/>\n-9 &#038; 5 &#038; -5 &#038; 12<br \/>\n\\end{bmatrix},\\ b=\\begin{bmatrix} -9\\\\ 5\\\\ 7\\\\ 11\\end{bmatrix}\\]  <\/p>\n<\/blockquote>\n<p><script>\nfunction showHtmlDiv43() {\n  var htmlShow43 = document.getElementById(\"html-show43\");\n  if (htmlShow43.style.display === \"none\") {\n    htmlShow43.style.display = \"block\";\n  } else {\n    htmlShow43.style.display = \"none\";\n  }\n}\n<\/script> <\/p>\n<p><button onclick=\"showHtmlDiv43()\">Soluci\u00f3n:<\/button><\/p>\n<div id=\"html-show43\" style=\"display: none;\">\n<iframe loading=\"lazy\" title=\"\u00c1lgebra Lineal - Factorizaci\u00f3n LU  Ej. 6 - Jes\u00fas Antonio Soto Espinosa\" width=\"640\" height=\"360\" src=\"https:\/\/www.youtube.com\/embed\/0J9RNjAbdxs?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe>\n<\/div>\n<hr \/>\n<h2>Matriz Permutaci\u00f3n<\/h2>\n<p>En ocasiones no es posible encontrar una factorizaci\u00f3n LU as\u00ed; por ejemplo si nos aparece un cero en la diagonal principal de la matriz U. En tal caso debemos permutar las filas o columnas de la matriz \\(A\\) para que no ocurra. Pero si lo hacemos debemos observar que ahora buscaremos una factorizaci\u00f3n de \\(PA\\) no de \\(A\\). Es decir, \\[PA=LU.\\]<\/p>\n<blockquote>\n<p><strong>Ejercicio:<\/strong> Dada la matriz \\[A=\\begin{bmatrix}<br \/>\n1 &#038; 4 &#038; -3\\\\<br \/>\n2 &#038; 8 &#038; 1 \\\\<br \/>\n-5 &#038; -9 &#038; 7<br \/>\n\\end{bmatrix}\\in\\mathcal{M}_3(\\mathbb{R}),\\]  \u00bfcu\u00e1nto suman los elementos de la matriz \\(\\textbf{L}\\) resultado de su factorizaci\u00f3n \\(\\textbf{LU}\\)? <\/p>\n<\/blockquote>\n<p><script>\nfunction showHtmlDiv42() {\n  var htmlShow42 = document.getElementById(\"html-show42\");\n  if (htmlShow42.style.display === \"none\") {\n    htmlShow42.style.display = \"block\";\n  } else {\n    htmlShow42.style.display = \"none\";\n  }\n}\n<\/script> <\/p>\n<p><button onclick=\"showHtmlDiv42()\">Soluci\u00f3n:<\/button><\/p>\n<div id=\"html-show42\" style=\"display: none;\">\n<iframe loading=\"lazy\" title=\"\u00c1lgebra Lineal - Factorizaci\u00f3n LU EJ.4 - Jes\u00fas Soto\" width=\"640\" height=\"360\" src=\"https:\/\/www.youtube.com\/embed\/sj8SbYWPd8w?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe>\n<\/div>\n<hr \/>\n<h2>Matrices rectangulares<\/h2>\n<p>Un enfoque curioso es abordar una pseudo-factorizaci\u00f3n <strong>LU<\/strong> cuando la matriz \\(A\\) no es cuadrada. Supongamos que \\(A\\in\\mathcal{M}_{n\\times m}(\\mathbb{K})\\). En este caso podemos buscar una factorizaci\u00f3n \\[\\textbf{A}=\\textbf{L}\\textbf{U}\\] de manera \\(\\textbf{L}\\) sea triangular inferior, de orden \\(n\\) y \\(\\textbf{U}\\) escalonada, de manera que, la submatriz cuadrada cuya diagonal principal es la diagonal principal de \\(\\textbf{U}\\), es una matriz triangular superior. Por ejemplo,<\/p>\n<p>\\[\\begin{bmatrix}<br \/>\n*&#038;*&#038;*&#038;*\\\\<br \/>\n*&#038;*&#038;*&#038;*\\\\<br \/>\n*&#038;*&#038;*&#038;*<br \/>\n\\end{bmatrix}<br \/>\n=<br \/>\n\\begin{bmatrix}<br \/>\n1&#038;0&#038;0\\\\<br \/>\n*&#038;1&#038;0\\\\<br \/>\n*&#038;*&#038;1<br \/>\n\\end{bmatrix}<br \/>\n\\begin{bmatrix}<br \/>\n*&#038;*&#038;*&#038;*\\\\<br \/>\n0&#038;*&#038;*&#038;*\\\\<br \/>\n0&#038;0&#038;*&#038;*<br \/>\n\\end{bmatrix}<br \/>\n\\]<br \/>\no<br \/>\n\\[\\begin{bmatrix}<br \/>\n*&#038;*&#038;*&#038;*\\\\<br \/>\n*&#038;*&#038;*&#038;*\\\\<br \/>\n*&#038;*&#038;*&#038;*\\\\<br \/>\n*&#038;*&#038;*&#038;*\\\\<br \/>\n*&#038;*&#038;*&#038;*<br \/>\n\\end{bmatrix}<br \/>\n=<br \/>\n\\begin{bmatrix}<br \/>\n1&#038;0&#038;0&#038;0&#038;0\\\\<br \/>\n*&#038;1&#038;0&#038;0&#038;0\\\\<br \/>\n*&#038;*&#038;1&#038;0&#038;0\\\\<br \/>\n*&#038;*&#038;*&#038;1&#038;0\\\\<br \/>\n*&#038;*&#038;*&#038;*&#038;1<br \/>\n\\end{bmatrix}<br \/>\n\\begin{bmatrix}<br \/>\n*&#038;*&#038;*&#038;*\\\\<br \/>\n0&#038;*&#038;*&#038;*\\\\<br \/>\n0&#038;0&#038;*&#038;*\\\\<br \/>\n0&#038;0&#038;0&#038;*\\\\<br \/>\n0&#038;0&#038;0&#038;0\\\\<br \/>\n\\end{bmatrix}<br \/>\n\\]<\/p>\n<p>Podremos hacer esta factorizaci\u00f3n, si los menores principales de \\(A\\) son todos distintos de cero. <\/p>\n<p>Un resultado equivalente ser\u00eda:<\/p>\n<blockquote>\n<p><strong>Teorema<\/strong>: Sea \\(A\\in\\mathcal{M}_{n\\times m}(\\mathbb{R})\\) una matriz que se puede reducir a una forma escalonada efectuando \u00fanicamente operaciones elementales de eliminaci\u00f3n (operaciones del tipo \\(\\alpha f_i+f_j\\) con \\(i &lt;j\\)). Entonces existe una matriz \\(n\\times n\\) triangular inferior \\(L\\) con unos en la diagonal principal y  una matriz \\(n\\times m\\), \\(U\\) con \\(u_{ij} = 0\\) si \\(i &gt;j\\) tales que \\(A=LU\\).\n<\/p>\n<\/blockquote>\n<p>Lo dicho anteriormente no garantiza la existencia o unicidad de la factorizaci\u00f3n; sin embargo, <\/p>\n<blockquote>\n<p><strong>Teorema<\/strong>: Sea \\(A\\in\\mathcal{M}_{n\\times m}(\\mathbb{K})\\), con una factorizaci\u00f3n \\(A=LU\\). Si \\(\\mathbf{rank}(A)=n\\)(\\(U\\) no tiene una fila de ceros) la factorizaci\u00f3n es \u00fanica<\/p>\n<\/blockquote>\n<blockquote>\n<p><strong>Ejercicio:<\/strong> Dada la matriz \\[A=\\begin{bmatrix}<br \/>\n3 &#038; -1 &#038; 4 &#038; 2\\\\<br \/>\n1 &#038; 2 &#038; -3 &#038; 5 \\\\<br \/>\n2 &#038; 4 &#038; 1&#038; 5<br \/>\n\\end{bmatrix}\\in\\mathcal{M}_{3\\times 4}(\\mathbb{R}),\\]  \u00bfcu\u00e1nto suman los elementos de la primera columna de la matriz \\(\\textbf{L}\\) resultado de su factorizaci\u00f3n \\(\\textbf{LU}\\)? <\/p>\n<\/blockquote>\n<p><script>\nfunction showHtmlDiv542() {\n  var htmlShow542 = document.getElementById(\"html-show542\");\n  if (htmlShow542.style.display === \"none\") {\n    htmlShow542.style.display = \"block\";\n  } else {\n    htmlShow542.style.display = \"none\";\n  }\n}\n<\/script> <\/p>\n<p><button onclick=\"showHtmlDiv542()\">Soluci\u00f3n:<\/button><\/p>\n<div id=\"html-show542\" style=\"display: none;\">\nSoluci\u00f3n: 2.\n<\/div>\n<hr \/>\n<hr \/>\n<h2>Aplicaci\u00f3n<\/h2>\n<p>El Boeing X-48 es un veh\u00edculo a\u00e9reo no tripulado (UAV) experimental de investigaci\u00f3n de las caracter\u00edsticas de aviones de ala integrada (BWB), un tipo de ala volante. Para su dise\u00f1o, los ingenieros de Boeing\u2019s Phantom Works usan el modelado tridimensional (3D) y la din\u00e1mica de fluidos computacional (DFC). Con la teor\u00eda DFC utiliza sistemas de ecuaciones que describen el flujo del aire sobre la superficie de la aeronave. El proceso para encontrar el flujo de aire alrededor de la aeronave implica la soluci\u00f3n repetida de un sistema de ecuaciones lineales \\(\\textbf{A}x=\\textbf{b},\\) que puede implicar hasta dos millones de ecuaciones y variables, adem\u00e1s de que el vector \\(\\textbf{b}\\) cambia a cada momento. Los sistemas resultantes son muy complicados. El programa computacional DFC aplicado en el Boeing utiliza la factorizaci\u00f3n LU de la matriz de coeficientes de estos sistemas.<\/p>\n<p>As\u00ed, para poder analizar una soluci\u00f3n de un sistema de flujo de aire, los ingenieros tienen que visualizar el flujo de aire sobre la superficie de la aeronave; para ello, utilizan gr\u00e1ficos generados por computadora, y el \u00e1lgebra lineal proporciona las herramientas para trazarlas.<\/p>\n<p>&nbsp;<\/p>\n<div style=\"text-align:center\">\n<a href=\"https:\/\/commons.wikimedia.org\/wiki\/File:X-48B_from_above.jpg#\/media\/Archivo:X-48B_from_above.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter\" src=\"https:\/\/upload.wikimedia.org\/wikipedia\/commons\/thumb\/d\/d4\/X-48B_from_above.jpg\/1200px-X-48B_from_above.jpg\" alt=\"X-48B from above.jpg\" width=\"542\" height=\"369\" \/><\/a><\/p>\n<p style=\"margin-left:100px;text-align:left;width: 40vw; font-size: x-small; \">\nDe Tony Landis for NASA &#8211; Este archivo fue catalogado por Armstrong Flight Research Center de la &lt;a href=\u00bb\/\/commons.wikimedia.org\/wiki\/National_Aeronautics_and_Space_Administration\u00bb class=\u00bbmw-redirect\u00bb title=\u00bbNational Aeronautics and Space Administration\u00bb&gt;Administraci\u00f3n Nacional de Aeron\u00e1utica y del Espacio&lt;\/a&gt; (&lt;a href=\u00bb\/\/commons.wikimedia.org\/wiki\/NASA\u00bb title=\u00bbNASA\u00bb&gt;NASA&lt;\/a&gt;) de los Estados Unidos de Am\u00e9ricabajo el identificador de foto:&amp;nbsp;&lt;a rel=\u00bbnofollow\u00bb class=\u00bbexternal text\u00bb href=\u00bbhttps:\/\/images.nasa.gov\/search-results?q=ED06-0198-62&#8243;&gt;ED06-0198-62&lt;\/a&gt;., Dominio p\u00fablico, <a href=\"https:\/\/commons.wikimedia.org\/w\/index.php?curid=1478722\">Enlace<\/a><\/p>\n<\/div>\n<hr \/>\n<table id=\"yzpi\" border=\"0\" width=\"100%\" cellspacing=\"0\" cellpadding=\"3\" bgcolor=\"#999999\">\n<tbody>\n<tr>\n<td width=\"100%\"><strong>Ejercicio:<\/strong> Sea \\(A\\)=[[-1,-7,-3], [2,15,6], [1,3,2]], \u00bfcu\u00e1l es el valor de la traza de su matriz adjunta? <\/td>\n<\/tr>\n<tr>\n<td>\n<div id=\"menu-a\">\n<ul>\n<li>-1<\/li>\n<li>1<\/li>\n<li>12<\/li>\n<\/ul>\n<\/div>\n<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p><script>\nfunction showHtmlDiv() {\n  var htmlShow = document.getElementById(\"html-show\");\n  if (htmlShow.style.display === \"none\") {\n    htmlShow.style.display = \"block\";\n  } else {\n    htmlShow.style.display = \"none\";\n  }\n}\n<\/script> <\/p>\n<p><button onclick=\"showHtmlDiv()\">Soluci\u00f3n:<\/button><\/p>\n<div id=\"html-show\" style=\"display: none;\">\n<p><strong>C.)<\/strong><\/p>\n<p><iframe loading=\"lazy\" src=\"https:\/\/uploads.jesussoto.es\/maxima\/ejrALGadj01.html\" width=\"650\" height=\"300\" allow=\"fullscreen\"><\/iframe>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>La factorizaci\u00f3n LU es una forma de factorizaci\u00f3n de una matriz como el producto de una matriz triangular inferior y una superior. El prop\u00f3sito es dada una matriz \\(A\\) conseguir descomponer esta&#8230;<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[6],"tags":[],"class_list":["post-217","post","type-post","status-publish","format-standard","hentry","category-algebra"],"_links":{"self":[{"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=\/wp\/v2\/posts\/217","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=217"}],"version-history":[{"count":2,"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=\/wp\/v2\/posts\/217\/revisions"}],"predecessor-version":[{"id":223,"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=\/wp\/v2\/posts\/217\/revisions\/223"}],"wp:attachment":[{"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=217"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=217"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=217"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}