Sistema de ecuaciones lineales

En matemática y álgebra lineal, un sistema de ecuaciones lineales, también conocido como sistema lineal de ecuaciones o simplemente sistema lineal, es un conjunto de ecuaciones lineales sobre un cuerpo o un anillo conmutativo. Un ejemplo de sistema lineal de ecuaciones sería el siguiente:

$\left\{{\begin{array}{rcrcrcr}3\,x_{1}&+&2\,x_{2}&+&\,x_{3}&=&1\\2\,x_{1}&+&2\,x_{2}&+&4\,x_{3}&=&-2\\-\,x_{1}&+&{\frac {1}{2}}\,x_{2}&-&\,x_{3}&=&0\end{array}}\right.$

El problema consiste en encontrar los valores desconocidos de las variables x₁, x₂ y x₃ que satisfacen las tres ecuaciones.

El problema de los sistemas lineales de ecuaciones es uno de los más antiguos de la matemática y tiene una infinidad de aplicaciones, como en procesamiento digital de señales, estimación, predicción y más generalmente en programación lineal así como en la aproximación de problemas no lineales de análisis numérico.

Introducción

En general, un sistema con m ecuaciones lineales n incógnitas puede ser escrito en forma ordinaria como:

${\begin{matrix}a_{11}x_{1}&+a_{12}x_{2}&+\dots &+a_{1n}x_{n}&=b_{1}\\a_{21}x_{1}&+a_{22}x_{2}&+\dots &+a_{2n}x_{n}&=b_{2}\\\dots &\dots &\dots &\dots &\dots \\a_{m1}x_{1}&+a_{m2}x_{2}&+\dots &+a_{mn}x_{n}&=b_{m}\end{matrix}}$

Donde $x_{1},\dots ,x_{n}\,$ son las incógnitas y los números $a_{ij}\in \mathbb {K}$ son los coeficientes del sistema sobre el cuerpo $\mathbb {K} \ [=\mathbb {R} ,\mathbb {C} ,\dots ]$. Es posible reescribir el sistema separando con coeficientes con notación matricial:

(1)${\begin{bmatrix}a_{11}&a_{12}&\cdots &a_{1n}\\a_{21}&a_{22}&\cdots &a_{2n}\\\vdots &\vdots &\ddots &\vdots \\a_{m1}&a_{m2}&\cdots &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}}$

Si representamos cada matriz con una única letra obtenemos:

$\mathbf {Ax} =\mathbf {b}$

Donde A es una matriz m por n, x es un vector columna de longitud n y b es otro vector columna de longitud m. El sistema de eliminación de Gauss-Jordan se aplica a este tipo de sistemas, sea cual sea el cuerpo del que provengan los coeficientes.

Sistemas lineales reales

En esta sección se analizan las propiedades de los sistemas de ecuaciones lineales sobre el cuerpo $\mathbb {R}$, es decir, los sistemas lineales en los coeficientes de las ecuaciones son números reales.

Representación gráfica

La intersección de dos planos no paralelos es una recta

Un sistema con $n\,$ incógnitas se puede representar en el n-espacio correspondiente.

En los sistemas con 2 incógnitas, el universo de nuestro sistema será el plano bidimensional, mientras que cada una de las ecuaciones será representada por una recta, si es lineal, o por una curva, si no lo es. La solución será el punto (o línea) donde intersecten todas las rectas y curvas que representan a las ecuaciones. Si no existe ningún punto en el que intersecten al mismo tiempo todas las líneas, el sistema es incompatible, o lo que es lo mismo, no tiene solución.

En el caso de un sistema con 3 incógnitas, el universo será el espacio tridimensional, siendo cada ecuación un plano dentro del mismo. Si todos los planos intersectan en un único punto, las coordenadas de éste serán la solución al sistema. Si, por el contrario, la intersección de todos ellos es una recta o incluso un plano, el sistema tendrá infinitas soluciones, que serán las coordenadas de los puntos que forman dicha línea o superficie.

Para sistemas de 4 ó más incógnitas, la representación gráfica no es intuitiva para el ser humano, por lo que dichos problemas no suelen enfocarse desde esta óptica.

Tipos de sistemas

Los sistemas de ecuaciones se pueden clasificar según el número de soluciones que pueden presentar. De acuerdo con ese caso se pueden presentar los siguientes casos:

• Sistema incompatible si no tiene ninguna solución.
• Sistema compatible si tiene alguna solución, en este caso además puede distinguirse entre:
  • Sistema compatible determinado cuando tiene un número finito de soluciones.
  • Sistema compatible indeterminado cuando admite un conjunto infinito de soluciones.

Quedando así la clasificación:

${\mbox{Tipos de sistemas}}{\begin{cases}{\mbox{Compatible}}&{\begin{cases}{\mbox{Determinado}}\\{\mbox{Indeterminado}}\end{cases}}\\{\mbox{Incompatible}}\end{cases}}$

Los sistemas incompatibles geométricamente se caracterizan por (hiper)planos o rectas que se cruzan sin cortarse. Los sistemas compatibles determinados se caracterizan por un conjunto de (hiper)planos o rectas que se cortan en un único punto. Los sistemas compatibles indeterminados se caracterizan por (hiper)planos que se cortan a lo largo de una recta [o más generalmente un hiperplano de dimensión menor]. Desde un punto de vista algebraico los sistemas compatibles determinados se caracterizan porque el determinante de la matriz es diferente de cero:

${\mbox{Sistema compatible determinado}}\Longleftrightarrow det(\mathbf {A} )\neq 0$

Sistemas compatibles indeterminados

Un sistema sobre un cuerpo K es compatible indeterminado cuando posee un número infinito de soluciones. Por ejemplo, el siguiente sistema:

$\left\{{\begin{matrix}x&+2y&=1\\2x&+4y&=2\end{matrix}}\right.$

Tanto la primera como la segunda ecuación se corresponden con la recta cuya pendiente es $-0,5\,$ y que pasa por el punto $(-1,1)\,$, por lo que ambas intersectan en todos los puntos de dicha recta. El sistema es compatible por haber solución o intersección entre las rectas, pero es indeterminado al ocurrir esto en infinitos puntos.

• En este tipo de sistemas, la solución genérica consiste en expresar una o más variables como función matemática del resto. En los sistemas lineales compatibles indeterminados, al menos una de sus ecuaciones se puede hallar como combinación lineal del resto, es decir, es linealmente dependiente.

• Una condición necesaria para que un sistema sea compatible indeterminado es que el determinante de la matriz del sistema sea cero (y por tanto uno de sus autovalores será 0):

${\mbox{sistema compatible indeterminado}}\Rightarrow \det \mathbf {A} =0$

• De hecho, de las dos condiciones anteriores se desprende, que el conjunto de soluciones de un sistema compatible indeterminado es un subespacio vectorial. Y la dimensión de ese espacio vectorial coincidirá con la multiplicidad geométrica del autovalor cero.

Texto en negrita==== Sistemas incompatibles ==== De un sistema se dice que es incompatible cuando no presenta ninguna solución. Por ejemplo, supongamos el siguiente sistema:

$\left\{{\begin{matrix}x&+2y&=4\\2x&+4y&=7\end{matrix}}\right.$

Las ecuaciones se corresponden gráficamente con dos rectas, ambas con la misma pendiente, Al ser paralelas, no se cortan en ningún punto, es decir, no existe ningún valor que satisfaga a la vez ambas ecuaciones.

Matemáticamente un sistema de estos es incompatible cuando el rango de la matriz del sistema es inferior al rango de la matriz ampliada. Una condición necesaria para que esto suceda es que el determinante de la matriz del sistema sea cero:

${\mbox{sistema incompatible}}\Rightarrow \det \mathbf {A} =0$

Los sistemas de ecuaciones lineales nos sirven en la vida diaria para

Métodos de resolución

Sustitución

El metodo de sustitución consiste en despejar en una de las ecuaciones cualquier incógnita, preferiblemente la que tenga menor coeficiente, para, a continuación, sustituirla en otra ecuación por su valor.

En caso de sistemas con más de dos incógnitas, la seleccionada debe ser sustituida por su valor equivalente en todas las ecuaciones excepto en la que la hemos despejado. En ese instante, tendremos un sistema con una ecuación y una incógnita menos que el inicial, en el que podemos seguir aplicando este método reiteradamente. Por ejemplo, supongamos que queremos resolver por sustitución este sistema:

$\left\{{\begin{matrix}3x&+y&=&22\\4x&-3y&=&-1\end{matrix}}\right.$

En la primera ecuación, seleccionamos la incógnita $y\,$ por ser la de menor coeficiente y que posiblemente nos facilite más las operaciones, y la despejamos, obteniendo la siguiente ecuación.

$y=22-3x\,$

El siguiente paso será sustituir cada ocurrencia de la incógnita $y\,$ en la otra ecuación, para así obtener una ecuación donde la única incógnita sea la $x\,$.

$4x-3(22-3x)=-1\qquad \Rightarrow 13x-66=-1,\qquad \Rightarrow 13x=65\,$

Al resolver la ecuación obtenemos el resultado $x=5\,$, y si ahora sustituimos esta incógnita por su valor en alguna de las ecuaciones originales obtendremos $y=7\,$, con lo que el sistema queda ya resuelto.

Igualación

El método de igualación se puede entender como un caso particular del método de sustitución en el que se despeja la misma incógnita en dos ecuaciones y a continuación se igualan entre sí la parte derecha de ambas ecuaciones.

Tomando el mismo sistema utilizado como ejemplo para el método de sustitución, si despejamos la incógnita $y\,$ en ambas ecuaciones nos queda de la siguiente manera:

$\left\{{\begin{matrix}y=&22-3x\\y=&{\cfrac {4x+1}{3}}\end{matrix}}\right.$

Como se puede observar, ambas ecuaciones comparten la misma parte izquierda, por lo que podemos afirmar que las partes derechas también son iguales entre sí.

$22-3x={\frac {4x+1}{3}}\,$

Llegados a este punto, la ecuación resultante es resoluble y podemos obtener el valor de la incógnita $x\,$, y a partir de aquí, sustituyendo dicho valor en una de las ecuaciones originales, obtener el valor de la $y\,$, que además ya se encuentra despejada.

La forma más fácil de tener el método de sustitución es realizando un cambio para despejar x después de averiguar el valor de la y.

Reducción

Este método suele emplearse mayoritariamente en los sistemas lineales, siendo pocos los casos en que se utiliza para resolver sistemas no lineales. El procedimiento, diseñado para sistemas con dos ecuaciones e incógnitas, consiste en transformar una de las ecuaciones (generalmente, mediante productos), de manera que obtengamos dos ecuaciones en la que una misma incógnita aparezca con el mismo coeficiente y distinto signo. A continuación, se suman ambas ecuaciones produciéndose así la reducción o cancelación de dicha incógnita, obteniendo así una ecuación con una sola incógnita, donde el método de resolución es simple.

Por ejemplo, en el sistema:

$\left\{{\begin{matrix}2x&+3y&=5\\5x&+6y&=4\end{matrix}}\right.$

no tenemos más que multiplicar la primera ecuación por $-2\,$ para poder cancelar la incógnita $y\,$. Al multiplicar, dicha ecuación nos queda así:

$-2(2x+3y=5)\quad \longrightarrow \quad -4x-6y=-10$

Si sumamos esta ecuación a la segunda del sistema original, obtenemos una nueva ecuación donde la incógnita $y\,$ ha sido reducida y que, en este caso, nos da directamente el valor de la incógnita $x\,$:

${\begin{array}{rrcr}-4x&-6y&=&-10\\5x&+6y&=&4\\\hline x&&=&-6\end{array}}$

$x=-6\,$

El siguiente paso consiste únicamente en sustituir el valor de la incógnita $x\,$ en cualquiera de las ecuaciones donde aparecían ambas incógnitas, y obtener así que el valor de $y\,$ es igual a:

$y={\frac {17}{3}}$

Método de Gauss

La eliminación de Gauss-Jordan, más conocida como método de Gauss, es un método aplicable únicamente a los sistemas lineales de ecuaciones, y consistente en triangular la matriz aumentada del sistema mediante transformaciones elementales, hasta obtener ecuaciones de una sola incógnita, cuyo valor será igual al coeficiente situado en la misma fila de la matriz. Este procedimiento es similar al anterior de reducción, pero ejecutado de manera reiterada y siguiendo un cierto orden algorítmico.

Tomemos como ejemplo el siguiente sistema:

$\left\{{\begin{array}{rrrrr}2x&+y&-z&=&8\\-3x&-y&+2z&=&-11\\-2x&+y&+2z&=&-3\end{array}}\right.$

Su matriz aumentada será esta:

$\left({\begin{array}{rrrr}2&1&-1&8\\-3&-1&2&-11\\-2&1&2&-3\end{array}}\right)$

En primer lugar, reducimos la incógnita $x\,$, sumando a la segunda fila, la primera multiplicada por ${\frac {3}{2}}$, y a la tercera, la primera fila. La matriz queda así:

$\left({\begin{array}{rrrr}2&1&-1&8\\0&{\frac {1}{2}}&{\frac {1}{2}}&1\\0&2&1&5\end{array}}\right)$

El siguiente paso consiste en eliminar la incógnita $y\,$ en la primera y tercera fila, para lo cual les sumamos la segunda multiplicada por $-2\,$ y por $-4\,$, respectivamente.

$\left({\begin{array}{rrrr}2&0&-2&6\\0&{\frac {1}{2}}&{\frac {1}{2}}&1\\0&0&-1&1\end{array}}\right)$

Por último, eliminamos la $z\,$, tanto de la primera como de la segunda fila, sumándoles la tercera multiplicada por $-2\,$ y por ${\frac {1}{2}}$, respectivamente:

$\left({\begin{array}{rrrr}2&0&0&4\\0&{\frac {1}{2}}&0&{\frac {3}{2}}\\0&0&-1&1\end{array}}\right)$

Llegados a este punto podemos resolver directamente las ecuaciones que se nos plantean:

$\left\{{\begin{matrix}2x=4\\{\cfrac {y}{2}}={\cfrac {3}{2}}\\-z=1\end{matrix}}\right.$

O, si lo preferimos, podemos multiplicar las tres filas de la matriz por: ${\frac {1}{2}}$, $2\,$ y $-1\,$ respectivamente, y obtener así automáticamente los valores de las incógnitas en la última columna.

$\left\{{\begin{matrix}x=&2\\y=&3\\z=&-1\end{matrix}}\right.$

Regla de Cramer

La regla de Cramer da una solución para sistemas compatibles determinados en términos de determinantes y adjuntos dada por:

$x_{j}={\cfrac {\det(A_{j})}{\det(\mathbf {A} )}}$

Donde A_j es la matriz resultante de remplazar la j-ésima columna de A por el vector columna b. Para un sistema de dos ecuaciones y dos incónitas:

$\left\{{\begin{matrix}a\,x&+&b\,y&=e\\c\,x&+&d\,y&=f\end{matrix}}\right.$

La regla de Cramer da la siguiente solución:

$x={\frac {\begin{vmatrix}e&b\\f&d\end{vmatrix}}{\begin{vmatrix}a&b\\c&d\end{vmatrix}}}={ed-bf \over ad-bc},\qquad y={\frac {\begin{vmatrix}a&e\\c&f\end{vmatrix}}{\begin{vmatrix}a&b\\c&d\end{vmatrix}}}={af-ec \over ad-bc}$

Nota: Cuando en la determinante original det(A) el resultado es 0, el sistema indica múltiples o sin coincidencia.

Sistemas lineales en un cuerpo arbitrario

Cuando consideramos ecuaciones lineales cuyas soluciones son números racionales, reales o complejos o más generalmente un cuerpo $\mathbb {K}$, la solución puede encontrarse mediante Regla de Cramer. Para sistemas de muchas ecuaciones la regla de Cramer puede ser computacionalmente más costosa y suelen usarse otros métodos más "económicos" en número de operaciones como la eliminación de Gauss-Jordan y la descomposición de Cholesky. Existen también métodos indirectos (basados en iteraciones) como el método de Gauss-Seidel.

Si el cuerpo es infinito (como es el caso de los números reales o complejos), entonces solo puede darse una de las tres siguientes situaciones:

• el sistema no tiene solución (en dicho caso decimos que el sistema está sobredeterminado o que es incompatible)
• el sistema tiene una única solución (el sistema es compatible determinado)
• el sistema tiene un número infinito de soluciones (el sistema es compatible indeterminado).

Un sistema de la forma

Ax = 0

se le llama sistema homogéneo. El conjunto de todas las soluciones de este tipo de sistema se le llama núcleo de la matriz y se escribe como Nuc A.

Se han diseñado algoritmos alternativos mucho más eficientes a la eliminación de Gauss-Jordan para una gran cantidad de casos específicos. La mayoría de estos algoritmos mejorados tienen una complejidad computacional de O(n²). Algunos de los métodos más usados son:

• Para los problemas de la forma Ax = b, donde A es una matriz de Toeplitz simétrica, se puede utilizar la recursión de Levinson o alguno de los métodos derivados de éste. Un método derivado de la recursión de Levinson es la recursión de Schur, que es ampliamente usado en el campo del procesamiento digital de señales.
• Para los problemas de la forma Ax = b, donde A es una matriz singular o casi singular, la matriz A se descompone en el producto de tres matrices en un proceso llamado descomposición de valores singulares.

Solución de sistemas lineales en un anillo

Los métodos para resolver el sistema (1) sobre un anillo son muy diferentes a los considerados anteriormente. De hecho la mayoría de métodos usados en cuerpos, como la regla de Cramer, son inaplicables en anillos debido a que no existen inversos multiplicativos.

La existencia de solución del sistema (1) sobre los enteros requiere varias condiciones:

1. Para cada i ${\mbox{mcd}}(a_{i1},a_{i2},...,a_{in})\;$ es divisor de $b_{i}\;$.
2. Si la condición anterior se cumple para un determinado i existe un conjunto de enteros ${\mathcal {S}}_{i}\;$ formado por el conjunto de enteros que satisface la i-ésima ecuación, y existirá solución si la intersección ${\mathcal {S}}_{1}\cap ...\cap {\mathcal {S}}_{n}\neq \varnothing$.

Véase también

• Sistema de ecuaciones
• Sistema lineal de dos ecuaciones con dos incógnitas

Enlaces externos

• Resolusor de ecuaciones lineales en línea