{"id":225,"date":"2025-10-09T09:15:41","date_gmt":"2025-10-09T07:15:41","guid":{"rendered":"https:\/\/clases.jesussoto.es\/?p=225"},"modified":"2025-10-05T21:07:56","modified_gmt":"2025-10-05T19:07:56","slug":"biomath-sistemas-de-ecuaciones-y-minimos-cuadrados","status":"publish","type":"post","link":"https:\/\/clases.jesussoto.es\/?p=225","title":{"rendered":"MathBio: Sistemas de ecuaciones y m\u00ednimos cuadrados"},"content":{"rendered":"<p>Recordad que aprendimos c\u00f3mo deducir las ecuaciones impl\u00edcitas de una variedad af\u00edn, en concreto de una recta y de un plano. Determinar los puntos que pertenece a una variedad es equivalente a resolver el sistema de ecuaciones lineales que definen la variedad. Para resolverlas utilizamos las matrices. As\u00ed todo sistema de ecuaciones lineales lo podemos plantear como un sistema matricial de la forma <strong>Ax<\/strong>=<strong>b<\/strong>, donde <strong>A<\/strong> es la matriz de coeficientes del sistema, <strong>x<\/strong> es la matriz columna de inc\u00f3gnitas y <strong>b<\/strong> es la matriz columna de t\u00e9rminos independientes:<\/p>\n<p>\\[{\\displaystyle {\\begin{bmatrix}a_{11}&#038;a_{12}&#038;\\cdots &#038;a_{1n}\\\\a_{21}&#038;a_{22}&#038;\\cdots &#038;a_{2n}\\\\\\vdots &#038;\\vdots &#038;\\ddots &#038;\\vdots \\\\a_{m1}&#038;a_{m2}&#038;\\cdots &#038;a_{mn}\\end{bmatrix}}{\\begin{bmatrix}x_{1}\\\\x_{2}\\\\\\vdots \\\\x_{n}\\end{bmatrix}}={\\begin{bmatrix}b_{1}\\\\b_{2}\\\\\\vdots \\\\b_{m}\\end{bmatrix}}}\\]<\/p>\n<p>Para tratarlos mejor podemos intentar transformalos en sistemas escalonados, que es resultado de transformar la matriz ampliada <strong>[A b]<\/strong>, mediante operaciones elementales de fila, en una matriz escalonada. Es te el m\u00e9todo que conocemos como <a href=\"http:\/\/es.wikipedia.org\/wiki\/Eliminaci%C3%B3n_de_Gauss-Jordan\" target=\"_blank\" rel=\"noopener noreferrer\">m\u00e9todo de Gauss<\/a>.<\/p>\n<p>Los sistemas de ecuaciones m\u00e1s sencillos resultan aquellos que podemos emplear la<a href=\"http:\/\/es.wikipedia.org\/wiki\/Regla_de_Cramer\" target=\"_blank\" rel=\"noopener noreferrer\"> regla de Cramer<\/a>.<\/p>\n<p>La importancia de <a href=\"http:\/\/es.wikipedia.org\/wiki\/Teorema_de_Rouch%C3%A9%E2%80%93Frobenius\" target=\"_blank\" rel=\"noopener noreferrer\">Teorema de Rouch\u00e9-Frobenius <\/a>estriba en que determina cuando un sistema tiene soluci\u00f3n o no. Este resultado junto con el anterior nos permiten resolver con facilidad los sistemas de ecuaciones como los ejercicios que hemos realizado.<\/p>\n<p>El Teorema de Rouch\u00e9-Fr\u00f6benius nos afirma que Existen soluciones para el sistema si y solo si el rango de la matriz de coeficientes es igual al rango de la matriz ampliada.<\/p>\n<p>As\u00ed un sistema ser\u00e1:<\/p>\n<p>\\[\\left\\{\\begin{array}{l} \\begin{array}{c} Compatible \\\\ rang(A)=rang(A|\\textbf{b}) \\end{array}\\left\\{\\begin{array}{l} \\begin{array}{c} Determinado \\\\ rang(A)=\\mbox{N\u00famero de inc\u00f3gnitas}<br \/>\\end{array} \\\\ \\begin{array}{c} Indeterminado \\\\ rang(A)<\\mbox{N\u00famero de inc\u00f3gnitas}\\end{array} \\\\\\end{array}\\right.\\\\\\begin{array}{c} Incompatible \\\\ rang(A)\\neq rang(A|\\textbf{b}) \\end{array}\\\\ \\end{array}\\right.\\]<\/p>\n<p>Para resolver un sistema compatible s\u00f3lo tenemos que encontrar un menor de \\(A\\) distinto de cero y del mismo orden que en rango de \\(A\\). Supongamos que  \\(\\bar{A}\\) es la submatriz de \\(A\\) cuyo menor es el que buscamos. Entonces \\([A|\\textbf{b}]\\) se puede transformar mediante operaciones elementales por filas en una matriz<\/p>\n<p>\\[[A|\\textbf{b}]\\sim\\left[\\begin{array}{c}<br \/>\n\\bar{A}\\,\\bar{P}\\\\<br \/>\n0\\end{array}\\left|\\begin{array}{c}<br \/>\n\\bar{\\textbf{b}}\\\\<br \/>\n0\\end{array}\\right.\\right]\\]<br \/>\nDonde \\(\\bar{P}\\) son o \\(0\\) o las columnas de la martiz \\(A\\) tales que \\[rang(A)+\\mbox{n\u00bacolumnas}(\\bar{P})=\\mbox{N\u00famero de inc\u00f3gnitas}.\\]<br \/>\nDe este modo el sistema tendr\u00e1 por soluci\u00f3n<br \/>\n\\[\\bar{X}=\\bar{A}^{-1}\\cdot [\\bar{\\textbf{b}}-\\bar{P}K],\\]<br \/>\ndonde \\(K\\) son las variables, en forma de par\u00e1metros, que faltan en el menor de \\(\\bar{A}\\), y tales que \\(X^t=\\bar{X}^t K^t\\).<\/p>\n<blockquote>\n<p><strong>Ejercicio:<\/strong> Determinar las ecuaciones param\u00e9tricas de la recta dada por la intersecci\u00f3n de los planos \\(\\pi_1:2x+3y-z=0\\) y \\(\\pi_2:-x-2y+z=0.\\)<\/p>\n<\/blockquote>\n<p><script>\nfunction showHtmlDiv20() {\n  var htmlShow20 = document.getElementById(\"html-show20\");\n  if (htmlShow20.style.display === \"none\") {\n    htmlShow20.style.display = \"block\";\n  } else {\n    htmlShow20.style.display = \"none\";\n  }\n}\n<\/script><\/p>\n<p><button onclick=\"showHtmlDiv20()\">Soluci\u00f3n:<\/button><\/p>\n<div id=\"html-show20\" style=\"display: none;\">\n<iframe loading=\"lazy\" title=\"\u00c1lgebra Lineal - Ecuaci\u00f3n impl\u00edcita: EJ. 3 - Jes\u00fas Soto\" width=\"640\" height=\"360\" src=\"https:\/\/www.youtube.com\/embed\/WOz2oJpCDQU?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> Determinar las ecuaciones param\u00e9tricas del plano \\(3x-6y+z=23.\\)<\/p>\n<\/blockquote>\n<p><script>\nfunction showHtmlDiv2() {\n  var htmlShow2 = document.getElementById(\"html-show2\");\n  if (htmlShow2.style.display === \"none\") {\n    htmlShow2.style.display = \"block\";\n  } else {\n    htmlShow2.style.display = \"none\";\n  }\n}\n<\/script><\/p>\n<p><button onclick=\"showHtmlDiv2()\">Soluci\u00f3n:<\/button><\/p>\n<div id=\"html-show2\" style=\"display: none;\">\n<iframe loading=\"lazy\" title=\"\u00c1lgebra Lineal - Ecuaci\u00f3n Impl\u00edcita Ejemplo 1 - Jes\u00fas Soto\" width=\"640\" height=\"360\" src=\"https:\/\/www.youtube.com\/embed\/PJtL7pah4cI?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> Sea la aplicaci\u00f3n lineal \\(f:\\mathbb{R}^3\\to\\mathbb{R}^3\\), dada por \\[f(x,y,z)=(2x-y+z,3x-2y+z,2x+2y-2z),\\] determinar la imagen rec\u00edproca del vector (-5,-9,-8).<\/p>\n<\/blockquote>\n<p><script>\nfunction showHtmlDiv2b() {\n  var htmlShow2b = document.getElementById(\"html-show2b\");\n  if (htmlShow2b.style.display === \"none\") {\n    htmlShow2b.style.display = \"block\";\n  } else {\n    htmlShow2b.style.display = \"none\";\n  }\n}\n<\/script><\/p>\n<p><button onclick=\"showHtmlDiv2b()\">Soluci\u00f3n:<\/button><\/p>\n<div id=\"html-show2b\" style=\"display: none;\">\n<iframe loading=\"lazy\" title=\"\u00c1lgebra Lineal - Imagen Rec\u00edproca. Ej. 2 - Jes\u00fas Soto\" width=\"640\" height=\"360\" src=\"https:\/\/www.youtube.com\/embed\/y3mT4R1a4DA?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> Sea la aplicaci\u00f3n lineal \\(f:\\mathbb{R}^4\\to\\mathbb{R}^3\\), dada por \\[f(x,y,z,t)=(x-y+z+t,x+2z-t,x+y+3z-3t),\\] determinar su n\u00facleo.<\/p>\n<\/blockquote>\n<p><script>\nfunction showHtmlDiv2a() {\n  var htmlShow2a = document.getElementById(\"html-show2a\");\n  if (htmlShow2a.style.display === \"none\") {\n    htmlShow2a.style.display = \"block\";\n  } else {\n    htmlShow2a.style.display = \"none\";\n  }\n}\n<\/script><\/p>\n<p><button onclick=\"showHtmlDiv2a()\">Soluci\u00f3n:<\/button><\/p>\n<div id=\"html-show2a\" style=\"display: none;\">\n<iframe loading=\"lazy\" title=\"\u00c1lgebra Lineal - N\u00facleo de una aplicaci\u00f3n lineal. Ejemplo 2 - Jes\u00fas Soto\" width=\"640\" height=\"360\" src=\"https:\/\/www.youtube.com\/embed\/N1UD5XE0Fis?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>M\u00ednimos cuadrados<\/h2>\n<p>Hemos visto c\u00f3mo solucionar sistemas de ecuaciones lineales. Sin embargo, hay en ocasiones que los sistemas no tiene soluci\u00f3n. En estos casos podemos buscar el punto m\u00e1s cercano a la soluci\u00f3n. Recordemos que todo sistema podemos plantearlo en su forma matricial como<br \/>\n\\[    A\\ x=\\ \\textbf{b},\\]<br \/>\ndonde \\(A\\in\\mathcal{M}_{m\\times n}(\\mathbb{R})\\), \\(x\\in\\mathcal{M}_{n\\times 1}(\\mathbb{R})\\) y \\(\\textbf{b}\\in\\mathcal{M}_{m\\times 1}(\\mathbb{R})\\). <\/p>\n<p>Cuando se necesita una soluci\u00f3n pero no hay ninguna, lo mejor que puede hacerse es encontrar una soluci\u00f3n \\(x\\) que deje a \\(A\\ x\\) tan cercana a \\(\\textbf{b}\\) como sea posible.<\/p>\n<blockquote>\n<p>\nSea \\( A\\in\\mathcal{M}_{m\\times n}(\\mathbb{R})\\), \\( x^t\\in\\mathbb{R}^n\\) y \\(\\textbf{b}^t\\in\\mathbb{R}^m\\), llamamos soluci\u00f3n por m\u00ednimos cuadrados de la ecuaci\u00f3n, a una aproximaci\u00f3n \\( \\hat{x}^t\\in\\mathbb{R}^n\\), tal que<br \/>\n\\[\\parallel\\textbf{b}-A\\ \\hat{x}\\parallel\\leq \\parallel \\textbf{b}-A\\ x\\parallel\\ \\forall\\ x^t\\in\\mathbb{R}^n.\\]\n<\/p>\n<\/blockquote>\n<p>El conjunto de soluciones por m\u00ednimos cuadrados de \\( A\\ x=\\ \\textbf{b}\\) coincide con el conjunto no vac\u00edo de soluciones de<br \/>\n\\[    A^t A\\ x=\\ A^t\\textbf{b},\\]<\/p>\n<p>De la propiedad anterior se deduce un resultado concluyente:<\/p>\n<blockquote>\n<p>\nSi las columnas de \\(A\\) son linealmente independientes, entonces \\(A^t A\\) es invertible y la ecuaci\u00f3n \\(A\\ x=\\ \\textbf{b}\\) tiene solamente una soluci\u00f3n por m\u00ednimos cuadrados dada por<br \/>\n\\[ \\hat{x}=\\ (A^t A)^{-1} A^t\\textbf{b}.\\]\n<\/p>\n<\/blockquote>\n<p>Esta forma de calcular la soluci\u00f3n por m\u00ednimos cuadrados ser\u00eda equivalente a considerar \\(\\bar{A}\\) el subespacio vectorial de \\(\\mathbb{R}^m\\) generado por los vectores columna de \\(A\\) y determinar<br \/>\n\\[   \\hat{\\textbf{b}}=proy_{\\bar{A}}(\\textbf{b}).\\]<br \/>\nEntonces<br \/>\n\\[ A\\ \\hat{x}=\\ \\hat{\\textbf{b}}.\\]<\/p>\n<hr>\n<h3>Bibliograf\u00eda<\/h3>\n<ul>\n<li>Cap\u00edtulo 6 del libro \u00c1lgebra lineal y sus aplicaciones. 4\u00ba edici\u00f3n, David C. Lay. Pearson<\/li>\n<\/ul>\n<hr>\n<table id=\"yzpi\" border=\"0\" width=\"100%\" cellspacing=\"0\" cellpadding=\"3\" bgcolor=\"#999999\">\n<tbody>\n<tr>\n<td width=\"100%\">\n<p><strong>Ejercicio:<\/strong> Cu\u00e1l es el error de una soluci\u00f3n por m\u00ednimos cuadrados del sistema incompatible<br \/>\n\\[\\begin{array}{r}<br \/>\n    -x +2y=4, \\\\<br \/>\n    2x-3y=1, \\\\<br \/>\n    -x+3y=2.<br \/>\n\\end{array}\\]<\/p>\n<div id=\"menu-a\">\n<ul>\n<li>\\(\\sqrt{11}\\)<\/li>\n<li>\\(\\sqrt{5}\\)<\/li>\n<li>\\(\\sqrt{21}\\)<\/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>A.)<\/strong><\/p>\n<p><iframe loading=\"lazy\" src=\"https:\/\/uploads.jesussoto.es\/maxima\/Ejer_minimos_cuadrados.html\" width=\"650\" height=\"300\" allow=\"fullscreen\"><\/iframe>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>Recordad que aprendimos c\u00f3mo deducir las ecuaciones impl\u00edcitas de una variedad af\u00edn, en concreto de una recta y de un plano. Determinar los puntos que pertenece a una variedad es equivalente a&#8230;<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[4],"tags":[],"class_list":["post-225","post","type-post","status-publish","format-standard","hentry","category-mathbio"],"_links":{"self":[{"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=\/wp\/v2\/posts\/225","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=225"}],"version-history":[{"count":7,"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=\/wp\/v2\/posts\/225\/revisions"}],"predecessor-version":[{"id":238,"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=\/wp\/v2\/posts\/225\/revisions\/238"}],"wp:attachment":[{"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=225"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=225"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/clases.jesussoto.es\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=225"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}