lunes, 21 de septiembre de 2009

Unidad ll : Sistema de ecuaciones lineales.

 Definicion.

Ecuaciones lineales con más de dos variables.
Para sistemas de ecuaciones lineales con más de dos variables, podemos usar el método de eliminación por sustitución o el método de eliminación por suma o resta (por adición o sustracción).
El método de eliminación por suma o resta es la técnica más breve y fácil de hallar soluciones. Además, lleva la técnica de matrices que se estudia en esta sección.
Cualquier sistema de ecuaciones lineales con tres variables tiene una solución única, un número infinito de soluciones o no tiene solución.
Método de eliminación para resolver un 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 cuerpoanillo conmutativo. Un ejemplo de sistema lineal de ecuaciones sería el siguiente: o un
    \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 x1, x2 y x3 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.

 Clasificacion.

Sistema general

La forma genérica de un sistema de m\, ecuaciones y n\, incógnitas es la siguiente:
\left\{\begin{matrix}F_1(x_1,...,x_n)=0 \\
\vdots \\
F_m(x_1,...,x_n)=0\end{matrix}\right.
donde F_1, \ldots, F_m son funciones de las incógnitas. La solución, perteneciente al espacio euclídeo  \R^n , será tal que el resultado de evaluar cualquier expresión F_i\, con los valores de dicha solución, verifique la ecuación.

Representación gráfica.

Los sistemas de 2 o 3 incógnitas admiten representaciones gráficas cuando las funciones F_i\, en (1) son continuas a tramos. En cada ecuación se representa como una curva o una superficie curva. La existencia de soluciones en ese caso puede deducirse a partir de la existencia de intersecciones comunes a dichas curvas o superficies curvas.

Clasificación de los sistemas.

Un sistema de ecuaciones sobre \R^n puede clasificarse de acuerdo con el número de soluciones en:
  • Sistema incompatible cuando no admite ninguna solución.
  • Sistema compatible cuando admite alguna solución que a su vez pueden dividirse en:
    • Sistemas compatibles indeterminados cuando existe un número infinito de soluciones que forman una variedad continua.
    • Sistemas compatibles determinados cuando admiten un conjunto finito de soluciones, o un conjunto infinito de soluciones aisladas con a lo sumo un número finito de puntos de acumulación.

Sistema lineal.

Un sistema como el anterior en que las anteriores ecuaciones son funciones afines. A diferencia del caso general, la solución de los sistemas de ecuaciones lineales son fáciles de encontrar cuando los coeficientes de las ecuaciones son números reales o complejos. También existen medios generales cuando los coeficientes pertenecen a un anillo, aunque la búsqueda de las soluciones en ese caso puede ser un poco más complicada.
Una característica importante de los sistemas lineales de ecuaciones es que admiten la llamada forma matricial. Esa forma permite representar el sistema usando tres matrices, de la siguiente forma:
\begin{pmatrix}
a_{11} & a_{12} & \cdots & a_{1Y} \\
a_{21} & a_{22} & \cdots & a_{2Y} \\
\vdots & \vdots & \ddots & \vdots \\
a_{X1} & a_{X2} & \cdots & a_{XY} \end{pmatrix} 

\begin{pmatrix}
x_1 \\
x_2 \\
\vdots \\
x_Y
\end{pmatrix} 
=
\begin{pmatrix}
b_1 \\
b_2 \\
\vdots \\
b_X
\end{pmatrix}
La primera es la matriz de coeficientes, donde el término a_{ij}\, representa al coeficiente que acompaña a la j-ésima incógnita de la ecuación i-ésima. La segunda es la matriz de incógnitas, donde cada término se corresponde con una de las Y\, incógnitas que queremos averiguar. Y la tercera matriz es la de términos independientes, donde el cada b_i\, representa al término independiente de la ecuación i-ésima.
Esta representación matricial facilita el uso de algunos métodos de resolución, como el método de Gauss, en el que, partiendo de la matriz aumentada (matriz de coeficientes a la que se le ha acoplado la matriz de términos independientes), y aplicando transformaciones lineales sobre las ecuaciones, pretendemos llegar a una matriz de este tipo:
\begin{pmatrix}
1 & 0 & \cdots & 0 & b_1 \\
0 & 1 & \cdots & 0 & b_2 \\
\vdots & \vdots & \ddots & \vdots & \vdots \\
0 & 0 & \cdots & 1 & b_X \end{pmatrix}
Una vez la matriz se ha triangulado, el valor de cada término b_i\, se corresponderá con el de la incógnita x_i\,. Si nos encontramos alguna fila del tipo \begin{pmatrix}0 & 0 & 0 & b_X \\\end{pmatrix}, con b_X\ne0\,, el sistema no tendrá solución.


Tipos de soluciones.

Existencia de soluciones.

El teorema de la función inversa proporciona condiciones suficientes de existencia de solución, de un sistema como con m = n\,. Si sucede que la función vectorial:
\mathbf{F}:\R^n \longrightarrow \R^n, \qquad
(x_1,\dots,x_n) \mapsto (F_1(x_1,\dots,x_n),\dots,F_n(x_1,\dots,x_n))
Es diferenciable con continuidad, es decir, es de clase C^1(\R^n) y su jacobiano no se anula en ningún punto entonces existe una única solución del sistema. Ya que en ese caso existirá una función inversa, y podremos escribir la solución buscada simplemente como:
(x_1,\dots,x_n)=\mathbf{F}^{-1}(\mathbf{0})
Sin embargo, la condición de diferenciabilidad anterior aún siendo condición suficiente, no es una condición necesaria, por lo que existen sistemas de ecuaciones en que las funciones F_i\, no son diferenciables y sin embargo, existen soluciones. Más aún, en casos en que existe más de una solución si la función es diferenciable entonces el jacobiano se anula en algún punto, pero eso no impide que existan varias soluciones.
En casos de un menor número de ecuaciones que de incógnitas, cuando m < n\,, entonces el sistema es compatible indeterminado o carece de soluciones. En esos casos, el teorema de la función implícita proporciona condiciones suficientes, aunque no necesarias, para la existencia de soluciones de un modo similar a como el teorema de la función inversa las proporciona en el caso m = n\,.

Métodos de resolución.

Si bien para los sistemas de ecuaciones lineales existen multitud de técnicas del álgebra lineal, para los sistemas de ecuaciones no-lineales el problema es técnicamente bastante más difícil.

Métodos analíticos.

Los métodos analíticos se restringen casi exclusivamente a sistemas de ecuaciones lineales. Ni siquiera se conoce una solución analítica para el sistema de ecuaciones de segundo grado general:
\begin{matrix}
a_{1,11}x^2_1 + a_{1,12}x_1x_2 + \dots + a_{1,nn}x_n^2 + b_{1,1}x_1 + \dots + b_{1,n}x_n + c_1= 0 \\
a_{2,11}x^2_1 + a_{2,12}x_1x_2 + \dots + a_{2,nn}x_n^2 + b_{2,1}x_1 + \dots + b_{2,n}x_n + c_2= 0 \\
\dots \\
a_{n,11}x^2_1 + a_{n,12}x_1x_2 + \dots + a_{n,nn}x_n^2 + b_{n,1}x_1 + \dots + b_{n,n}x_n + c_n= 0
\end{matrix}

Métodos numéricos

Las aplicaciones técnicas generalmente recurren a algoritmos numéricos que permiten calcular aproximaciones numéricas a las solciones de un sistema de ecuaciones.
Uno de los métodos numéricos que puede generalizarse a sistemas no-lineales es el método de Newton-Raphson. En el caso multidimensional la resolución numérica del sistema de n ecuaciones \scriptstyle \mathbf{f}(\mathbf{x}) = \mathbf{f}(x_1, \dots, x_n)=0 puede hacerse a partir del conocimiento de una solución aproximada \scriptstyle \mathbf{x}^{(0)} = (x_1^{(0)}, \dots, x_n^{(0)}), siempre y cuando la aplicación anterior sea difernciable, mediante el esquema iterativo:
\mathbf{x}^{(m+1)} = \mathbf{x}^{(m)} - [D\mathbf{f}(\mathbf{x}^{(m)})]^{-1}(\mathbf{f}(\mathbf{x}^{(m))}), \qquad \mathbf{f}:\R^n \to \R^n,\ \mathbf{f}\in C^{(1)}(\R^n;\R^n)
O más explícitamente:
\begin{bmatrix} x^{(m+1)}_1\\ \dots \\ x^{(m+1)}_n \end{bmatrix} =
\begin{bmatrix} x^{(m)}_1\\ \dots \\ x^{(m)}_n \end{bmatrix} - 
\begin{bmatrix} D_1f_1 & \dots & D_nf_1 \\ \dots \\ D_1f_n & \dots & D_nf_n \end{bmatrix}^{-1} 
\begin{bmatrix} f_1(\mathbf{x}^{(m)})\\ \dots \\ f_n(\mathbf{x}^{(m)}) \end{bmatrix}
Lamentablemente la convergencia del esquema iterativo anterior no está garantizada y en casos de soluciones múltiples la convergencia puede darse hacia la solución no deseada.

Métodos gráficos

Los métodos gráficos son didácticos e elustrativos, aunque en general carecen de interés práctico en las aplicaciones técnicas de importancia. Además están restringidos generalmente a sistemas de dos o tres ecuaciones reales.
Dos sistemas de ecuaciones con dos incognitas de valor real, suelen aparecer como uno de los cinco tipos diferentes mencionados a continuación. Tienen una relación con el número de soluciones:
  1. Aquellos sitemas de ecuaciones que representan gráficamente rectas y curvas que se intersectan entre si. Este tipo de sistema de ecuación es considerado como el normal. Suele tener un numero de soluciones finito cada uno formado por las coordenadas de los punto de intersección.
  2. Sistemas que tienen simplificaciones falsas. Por ejemplo: 1 = 0. Graficamente se representan como un conjunto de lineas que nunca se intersectan entre si, como lineas paralelas.
  3. Sistemas de ecuaciones en las que ambos simplificar a una identidad (por ejemplo, x = 2x - y 0 y x = 0). Cualquier asignación de valores a las variables desconocidas satisface las ecuaciones. Por lo tanto, hay un número infinito de soluciones, que gráficamente, se representa como todos los puntos del plano que representa la solución.
  4. Sistemas en los que las dos ecuaciones representan el mismo conjunto de puntos: son matemáticamente equivalentes (una ecuación general puede ser transformada en otra a través de la manipulación algebraica). Estos sistemas representan completamente la superposición de líneas o curvas, etc Una de las dos ecuaciones es redundante y puede ser desechada. Cada punto de la serie de puntos corresponde a una solución. Generalmente, esto significa que hay un número infinito de soluciones.
  5. Sistemas en los que una (y sólo una) de las dos ecuaciones se simplifica a una identidad. Por lo tanto, es redundante y puede ser descartada, según el tipo anterior. Cada punto de la serie de puntos representados por los demás es una solución de la ecuación de los que hay a continuación, por lo general un número infinito.
La ecuación x2 + y2 = 0 puede ser pensada como la ecuación de un círculo cuyo radio se ha reducido a cero, por lo que representa un único punto: (x = 0, y = 0), a diferencia de una normal de un círculo que contiene infinito número de puntos. Este y otros casos similares muestran la razón por la cual los dos últimos tipos anteriormente descritos necesitan la calificación de "normalmente". Un ejemplo de un sistema de ecuaciones del primer tipo descrito anteriormente, con un número infinito de soluciones viene dada por x = | x |, y = | y | (donde la notación | • | indica el valor absoluto de la función), cuyas soluciones de forma un cuadrante de la x - y plano. Otro ejemplo es x = | y |, y = | x |, cuya solución representa un rayo.


Interpretacion Geometrica.


Resolución y representación gráfica de sistemas de ecuaciones lineales

Resolvemos gráficamente el sistema x + y = 6; x - y = 2}


  1. o Despejamos y en las dos ecuaciones. x + y = 6 → y = 6 - x
    x - y = 2 → y = x - 2
  2. o Dando valores a x, formamos una tabla de valores para cada una de las dos ecuaciones.

y = 6 - x

x 0 1 2 3 4
y 6 5 4 3 2
y = x - 2

x 0 1 2 3 4
y -2 -1 0 1 2

  1. o Representamos estos puntos sobre un sistema de ejes.


Uniendo los puntos de cada ecuación, obtenemos dos rectas que representan todas las soluciones de cada una de las ecuaciones Uniendo los puntos de cada ecuación, obtenemos dos rectas que representan todas las soluciones de cada una de las ecuaciones

  1. o Puede ocurrir uno de los siguientes casos:
    • Si las rectas no se cortan, es decir, son paralelas, el sistema es incompatible, no tiene solución.
    • Si las rectas se cortan en un punto, el sistema tiene solución única. Decimos que es compatible determinado.
    • Si las dos rectas coinciden, esto es, son la misma, el sistema tiene infinitas soluciones. Es un sistema compatible indeterminado.
    En nuestro caso, las rectas se cortan en el punto (4, 2). La solución del sistema es x = 4 e y = 2.

Clasificamos los siguientes sistemas de ecuaciones lineales


    a) 2 x + y = 6 2 x - y = 2 } b) x + y = 3 2 x + 2 y = 6 } c) x + y = 3 x + y = - 1 }


    a) Dibujamos las rectas que representan las soluciones de cada ecuación: Dos soluciones de la primera ecuación son: x = 1, y = 4; x = 2, y = 2 Dos soluciones de la segunda ecuación son: x = 1, y= 0; x = 2, y = 2 Las rectas se cortan en un punto que será la solución:x = 2, y = 2. Por tanto, el sistema será compatible determinado. Vemos la representación en el margen. b) Dibujamos las rectas que representan las soluciones de cada ecuación: Dos soluciones de la primera ecuación son: x = 0, y = 3; x = 3, y = 0 Dos soluciones de la segunda ecuación son: x = 1, y = 2; x = 2, y = 1 Las rectas coinciden, toda la recta es solución del sistema (infinitas soluciones). Por tanto, el sistema será compatible indeterminado. Vemos la representación en el margen. c) Dibujamos las rectas que representan las soluciones de cada ecuación: Dos soluciones de la primera ecuación son: x = 0,y = 3; x = 3,y = 0 Dos soluciones de la segunda ecuación son: x = 0, y =-1; x = -2, y = 1 Las rectas son paralelas, no tienen ningún punto en común, luego el sistema no tiene solución. Por tanto, el sistema será incompatible. Vemos la representación siguiente:


Representación gráfica de los tres sistemas Representación gráfica de los tres sistemas



 Metodos de solucion  de sistemas de ecuaciones lineales.

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í:
   -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 \, :
   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 \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 Aj 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.

Artículo principal: ecuación diofántica
Los métodos para resolver el sistema 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 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 \ne \varnothing.

 Obtenido de:
http://es.wikipedia.org/wiki/Sistema_de_ecuaciones

http://mx.kalipedia.com/matematicas-algebra/tema/resolucion-representacion-grafica-sistemas.html?x1=20070926klpmatalg_145.Kes&x=20070926klpmatalg_146.Kes


http://es.wikipedia.org/wiki/Sistema_de_ecuaciones_lineales