Sistema de ecuacións lineais

[miniatura|dereita|Un sistema lineal de tres variables determina unha colección de [plano (matemáticas)|planos. O punto de intersección é a solución.]]
En matemáticas, un **sistema de ecuacións lineais** (ou **sistema linear**) é unha colección de ecuacións lineais que levan o mesmo conxunto de variables. Por exemplo,
:\begin{alignat}{7}
3x &&\; + \;&& 2y             &&\; - \;&& z  &&\; = \;&& 1 & \\
2x &&\; - \;&& 2y             &&\; + \;&& 4z &&\; = \;&& -2 & \\
-x &&\; + \;&& \tfrac{1}{2} y &&\; - \;&& z  &&\; = \;&& 0 &
\end{alignat}
é un sistema de tres ecuacións das tres variables x, y e z. Unha **solución** a un sistema lineal é unha asignación de números ás variables de xeito que tódalas ecuacións sexan simultaneamente satisfeitas. Unha solución ao sistema anterior é dada por
:\begin{alignat}{2}
x & = & 1 \\
y & = & -2 \\
z & = & -2
\end{alignat}
xa que fai que as tres ecuacións sexan válidas.

En matemáticas, a teoría de sistemais lineais é unha rama da álxebra lineal, unha materia fundamental das matemáticas modernas. Os algoritmos computacionais para atopar solucións son unha parte importante da álxebra lineal numérica e, tales métodos, xogan un rol moi importante na enxeñería, a física, a química, a informática e as ciencias económicas. Adoito, un sistema de ecuacións non lineais pódese aproximar por un sistema lineal (ver linealización), unha técnica útil para facer un modelo matemático ou unha simulación computerizada dun sistema relativamente complexo.

## Exemplo elemental
Un dos tipos máis sinxelos de sistema linear implica dúas ecuacións e dúas variábeis:

: \begin{alignat}{5}
2x &&\; + \;&& 3y &&\; = \;&& 6 &. \\
4x &&\; + \;&& 9y &&\; = \;&& 15&.
\end{alignat}

Un método para resolver un sistema deste tipo é o seguinte. Primeiro, resolva a ecuación superior para x en función de y (o que usualmente se chama despexar ou resolver a x):

: x = 3 - \frac{3}{2}y.

Agora substitúa esta expresión por *x* na ecuación inferior:

: 4\left( 3 - \frac{3}{2}y \right) + 9y = 15.
: 12 - 6y  + 9y = 15.

Isto dá como resultado unha única ecuación que inclúe só a variábel y. Ao resolver temos y = 1, e substituíndo isto de novo na ecuación con x obtemos x = \frac{3}{2}.

## Forma xeral
Un sistema xeral de *m* ecuacións lineares con *n* incógnitas e os coeficientes, de cada ecuación e cada incógnita, pódese escribir como
:\begin{cases}
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 \\ 
\vdots\\
a_{m1} x_1 + a_{m2} x_2 + \dots + a_{mn} x_n = b_m,
\end{cases}
onde x_1, x_2,\dots,x_n son as incógnitas, a_{11},a_{12},\dots,a_{mn} son os coeficientes do sistema e  b_1,b_2,\dots,b_m son os termos constantes.

### Ecuación vectorial 
Unha visión moi útil é que cada incógnita é un peso para un vector columna nunha combinación linear.
:
x_1\begin{bmatrix}a_{11}\\a_{21}\\ \vdots \\a_{m1}\end{bmatrix} +
x_2\begin{bmatrix}a_{12}\\a_{22}\\ \vdots \\a_{m2}\end{bmatrix} +
\dots +
x_n\begin{bmatrix}a_{1n}\\a_{2n}\\ \vdots \\a_{mn}\end{bmatrix} = \begin{bmatrix}b_1\\b_2\\ \vdots \\b_m\end{bmatrix}

### Ecuación matricial 
A ecuación vectorial é equivalente a unha ecuación matricial da forma
: A\mathbf{x} = \mathbf{b} 
onde *A* é unha matriz de dimensión *m*×*n*, **x** é un vector columna con *n* coeficientes e **b** é un vector columna con *m* coeficientes.

:
A=
\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},\quad
\mathbf{x}=
\begin{bmatrix}
x_1 \\
x_2 \\
\vdots \\
x_n
\end{bmatrix},\quad
\mathbf{b}=
\begin{bmatrix}
b_1 \\
b_2 \\
\vdots \\
b_m
\end{bmatrix}.

Agora poderemos relacionar conceptos do sistema coas matrices como o *rango* da matriz.

## Conxunto de solucións
thumb|O conxunto de solucións para as ecuacións  e  é o punto único (2, 3).

Unha solución] dun sistema linear é unha asignación de valores ás variábeis  de forma que se cumpre para cada unha das ecuacións. O conxunto de todas as posíbeis solucións chámase *conxunto de solucións*.

Un sistema linear pode comportarse de calquera das tres formas seguintes:
# O sistema ten *infinitas solucións*.
# O sistema ten unha *solución única*.
# O sistema non ten "ningunha solución".

### Interpretación xeométrica
Para un sistema que inclúe dúas variábeis (*x* e *y*), cada ecuación linear determina unha recta no plano *xy*. Dado que unha solución dun sistema linear debe satisfacer todas as ecuacións, o conxunto solución é a intersección destas rectas e, polo tanto, é unha recta (cando as rectas son coincidentes), un só punto (cando as rectas se cruzan) ou o conxunto baleiro (cando as rectas son paralelas).

Para tres variábeis, cada ecuación linear determina un plano no espazo tridimensional, e o conxunto de solucións é a intersección destes planos. De modo parecido que en dúas dimensións, con máis combinacións de posibilidades, os planos poden coincidir, seren paralelos ou intersecar nun punto ou nunha recta.

Para *n* variábeis a interpretación sería dentro dun hiperplano.

### Comportamento xenérico
En xeral, o comportamento dun sistema linear está determinado pola relación entre o número de ecuacións e o número de incógnitas.
* Un sistema con menos ecuacións que incógnitas ten infinitas solucións, pero pode que non teña solución. Este sistema coñécese como sistema indeterminado.
* Un sistema co mesmo número de ecuacións e incógnitas ten unha solución única.
* En xeral, un sistema con máis ecuacións que incógnitas non ten solución. Este sistema tamén se coñece como sistema sobredeterminado.
No primeiro caso, a dimensión do conxunto de solucións é, en xeral, igual a , onde *n* é o número de variábeis e *m* é o número de ecuacións.

Un sistema de ecuacións lineares compórtase de forma diferente ao caso xeral se as ecuacións son *dependentes linearmente*, ou se é *inconsistente* e non ten máis ecuacións que incógnitas.

## Propiedades
### Independencia
As ecuacións dun sistema linear son **independentes** se ningunha das ecuacións se pode derivar alxébricamente das outras. Cando as ecuacións son independentes, cada ecuación contén nova información sobre as variábeis, e ao eliminar calquera das ecuacións aumenta o tamaño do conxunto de solucións.

thumb|As ecuacións ,  e  son linearmente dependentes porque a terceira é a suma das dúas primeiras.
Por exemplo, as ecuacións
: 3x+2y=6\;\;\;\;\text{e}\;\;\;\;6x+4y=12
non son independentes: son a mesma ecuación cando se escalan cun factor de dous e producirían gráficos idénticos. Este é un exemplo de equivalencia nun sistema de ecuacións lineares.

Para un exemplo máis complicado, as ecuacións
: \begin{alignat}{5}
 x &&\; - \;&& 2a &&\; = \;&& -1 & \\
 3x &&\; + \;&& 5y &&\; = \;&& 8 & \\
 4x &&\; + \;&& 3y &&\; = \;&& 7 &
\end{alignat}
non son independentes, porque a terceira ecuación é a suma das outras dúas. De feito, calquera destas ecuacións pódese derivar das outras dúas, e calquera das ecuacións pódese eliminar sen afectar ao conxunto de solucións. As gráficas destas ecuacións son tres rectas que se cortan nun só punto.

### Consistencia
thumb|As ecuacións  e  son inconsistentes.
Un sistema linear é **inconsistente** se non ten solución e no caso contrario, dise que é **consistente**. Cando o sistema é inconsistente, é posíbel derivar unha contradición das ecuacións, que sempre pode reescribirse como a declaración .

Por exemplo, as ecuacións
:3x+2y=6\;\;\;\;\text{e}\;\;\;\;3x+2y=12
son inconsistentes. De feito, restando a primeira ecuación da segunda e multiplicando os dous lados do resultado por 1/6, obtemos . As gráficas destas ecuacións no plano *xy* son un par de liñas paralelas.

É posíbel que tres ecuacións lineares sexan inconsistentes, aínda que dúas calquera delas sexan consistentes xuntas. Por exemplo, as ecuacións
: \begin{alignat}{7}
 x &&\; + \;&& y &&\; = \;&& 1 & \\
2x &&\; + \;&& y &&\; = \;&& 1 & \\
3x &&\; + \;&& 2y &&\; = \;&& 3 &
\end{alignat}
son inconsistentes. Sumando as dúas primeiras ecuacións temos , que se pode restar da terceira ecuación para obter . Dúas destas ecuacións teñen unha solución común. O mesmo fenómeno pode ocorrer para calquera número de ecuacións.

En xeral, as inconsistencias ocorren se os lados esquerdos das ecuacións dun sistema son linearmente dependentes e os termos constantes non satisfán a relación de dependencia. Un sistema de ecuacións cuxos lados esquerdos son linearmente independentes sempre é consistente.

Dito doutro xeito, segundo o teorema de Rouché-Capelli, calquera sistema de ecuacións (sobredeterminado ou non) é inconsistente se o rango da matriz aumentada é maior que o rango da matriz de coeficientes. Se, pola contra, os rangos destas dúas matrices son iguais, o sistema debe ter polo menos unha solución. A solución é única se e só se o rango é igual ao número de variábeis. En caso contrario, a solución xeral ten *k* parámetros libres onde *k* é a diferenza entre o número de variábeis e o rango; daí que en tal caso exista unha infinidade de solucións. 

Os sistemas de ecuacións pódense clasificar segundo o número de solucións que poden presentar. Segundo este caso, pódense presentar os seguintes casos:

* **Sistema compatíbel** ou **posíbel** se ten solución, neste caso tamén se pode distinguir entre:
** **Sistema compatíbel determinado** ou **posíbel determinado** cando tes unha única solución.
** **Sistema compatíbel indeterminado** ou **posíbel indeterminado** cando admite un conxunto infinito de solucións.
* **Sistema incompatíbel** ou **imposíbel** se non hai solución.

Un sistema é compatíbel determinado \Longleftrightarrow \det(\mathbf{A}) \ne 0 .

### Equivalencia
Dous sistemas lineares que usan o mesmo conxunto de variábeis son **equivalentes** se cada unha das ecuacións do segundo sistema se pode derivar alxébricamente a partir das ecuacións do primeiro sistema, e viceversa. Dous sistemas son equivalentes se ambos os dous son inconsistentes ou cada ecuación de cada un deles é unha combinación linear das ecuacións do outro. Por tano dous sistemas lineares son equivalentes se e só se teñen o mesmo conxunto de solucións.

## Resolver un sistema linear
### Eliminación de variábeis
O método máis elemental para resolver un sistema de ecuacións lineares é eliminar varias variábeis. Este método pódese describir do seguinte xeito:
# Na primeira ecuación, resolva unha das variábeis en función das outras.
# Substitúa esta expresión nas ecuacións restantes. Isto dá un sistema de ecuacións cunha ecuación e unha incógnita menos.
# Repita os pasos 1 e 2 ata que o sistema se reduza a unha única ecuación linear.
# Resolva esta ecuación 
# Percorra cara atrás as ecuacións substituíndo os valores conseguidos ata atopar a solución completa.

Por exemplo, considere o seguinte sistema:
:\begin{cases}
x+3y-2z=5\\
3x+5y+6z=7\\
2x+4y+3z=8
\end{cases}
Resolvendo a primeira ecuación para x temos x=5+2z-3y, e metendo ese resultado na segunda e terceira ecuación obtemos
:\begin{cases}
3(5+2z-3y)+5y+6z=7\\
2(5+2z-3y)+4y+3z=8
\end{cases}
que ao simplificar fica como
:\begin{cases}
y=3z+2\\
y=\tfrac{7}{2}z+1
\end{cases}
Xa que a parte esquerda destas dúas ecuacións (LHS en inglés) é igual a *y*, ao igualar as partes dereitas das ecuacións (RHS en inglés) temos unha ecuación cunha única incógnita. Agora temos:
:\begin{align}
3z+2=\tfrac{7}{2}z+1\\ \rightarrow z=2
\end{align}
 
Imos cara atrás e substituímos *z* = 2 na segunda ou terceira ecuación dá *y* = 8, e os valores de *y* e *z* na primeira ecuación producen *x* = - 15. Polo tanto, o conxunto de solucións é o terna ordenada (x,y,z)=(-15,8,2) .

### Redución de filas
Na **redución de filas** (tamén coñecida como **eliminación de Gauss**), o sistema linear represéntase como unha matriz aumentada

:\left[\begin{array}{rrr|r}
1 & 3 & -2 & 5 \\
3 & 5 & 6 & 7 \\
2 & 4 & 3 & 8
\end{array}\right]

A continuación, esta matriz modifícase usando operacións elementais de filas ata alcanzar a forma escalonada de filas reducida. Hai tres tipos de operacións elementais de filas: 
:**Tipo 1**: troca as posicións de dúas filas.
:**Tipo 2**: Multiplica unha fila por un escalar (número) distinto de cero.
:**Tipo 3**: Suma ou resta a unha fila outra fila multiplicada por un escalar.
Dado que estas operacións son reversíbeis, a matriz aumentada producida sempre representa un sistema linear que é equivalente ao orixinal.

Existen varios algoritmos específicos para reducir por filas unha matriz aumentada, o máis sinxelo é a eliminación de Gauss (ou  Gauss-Jordan). O seguinte cálculo mostra a eliminación de Gauss-Jordan aplicada á matriz anterior:
:\begin{align}\left[\begin{array}{rrr|r}
1 & 3 & -2 & 5 \\
3 & 5 & 6 & 7 \\
2 & 4 & 3 & 8
\end{array}\right]&\sim
\left[\begin{array}{rrr|r}
1 & 3 & -2 & 5 \\
0 & -4 & 12 & -8 \\
2 & 4 & 3 & 8
\end{array}\right]\sim
\left[\begin{array}{rrr|r}
1 & 3 & -2 & 5 \\
0 & -4 & 12 & -8 \\
0 & -2 & 7 & -2
\end{array}\right]\sim
\left[\begin{array}{rrr|r}
1 & 3 & -2 & 5 \\
0 & 1 & -3 & 2 \\
0 & -2 & 7 & -2
\end{array}\right]
\\
&\sim
\left[\begin{array}{rrr|r}
1 & 3 & -2 & 5 \\
0 & 1 & -3 & 2 \\
0 & 0 & 1 & 2
\end{array}\right]\sim
\left[\begin{array}{rrr|r}
1 & 3 & -2 & 5 \\
0 & 1 & 0 & 8 \\
0 & 0 & 1 & 2
\end{array}\right]\sim
\left[\begin{array}{rrr|r}
1 & 3 & 0 & 9 \\
0 & 1 & 0 & 8 \\
0 & 0 & 1 & 2
\end{array}\right]\sim
\left[\begin{array}{rrr|r}
1 & 0 & 0 & -15 \\
0 & 1 & 0 & 8 \\
0 & 0 & 1 & 2
\end{array}\right]\end{align}
A última matriz está en forma de chanzo de fila reducida e representa o sistema x = -15,  y = 8,  z = 2 . Unha comparación co exemplo da sección anterior sobre a eliminación alxébrica de variábeis mostra que estes dous métodos son de feito iguais; a diferenza radica en como se escriben os cálculos.

### Regra de Cramer
**A regra de Cramer** é unha fórmula explícita para a solución dun sistema de ecuacións lineares, con cada variábel dada por un cociente de dous determinantes.
Seguindo co exemplo deste artigo:

:
x=\frac
{\, \begin{vmatrix}5&3&-2\\7&5&6\\8&4&3\end{vmatrix} \,}
{\, \begin{vmatrix}1&3&-2\\3&5&6\\2&4&3\end{vmatrix} \,}
,\;\;\;\;
y=\frac
{\, \begin{vmatrix}1&5&-2\\3&7&6\\2&8&3\end{vmatrix} \,}
{\, \begin{vmatrix}1&3&-2\\3&5&6\\2&4&3\end{vmatrix} \,}
,\;\;\;\;
z=\frac
{\, \begin{vmatrix}1&3&5\\3&5&7\\2&4&8\end{vmatrix} \,}
{\, \begin{vmatrix}1&3&-2\\3&5&6\\2&4&3\end{vmatrix} \,}.

Para cada variábel, o denominador é o determinante da matriz de coeficientes, mentres que o numerador é o determinante dunha matriz na que unha columna foi substituída polo vector de termos constantes.

Aínda que a regra de Cramer é importante teoricamente, ten pouco valor práctico para matrices grandes, xa que o cálculo de grandes determinantes é un bocadiño engorroso.

### Solución matricial
Se o sistema de ecuacións se expresa na forma matricial A\mathbf{x}=\mathbf{b}, o conxunto de solucións enteiro tamén se pode expresar en forma matricial. Se a matriz *A* é cadrada (ten *m* filas e *n*=*m* columnas) e ten rango completo (todas as filas *m* son independentes), entón o sistema ten unha solución única dada por

:\mathbf{x}=A^{-1}\mathbf{b}

onde A^{-1} é a inversa de *A*. De forma máis xeral, independentemente de se *m*=*n* ou non e independentemente do rango de *A*, todas as solucións (se as existen) danse usando a inversa de Moore-Penrose de *A*, denotada como A^+, como segue:

:\mathbf{x}=A^+ \mathbf{b} + \left(I - A^+ A\right)\mathbf{w}

onde \mathbf{w} é un vector de parámetros libres que abrangue todos os posíbeis vectores de dimensión n \times 1 .

### Outros métodos
Aínda que os sistemas de tres ou catro ecuacións pódense resolver facilmente a man, as computadoras úsanse a miúdo para sistemas máis grandes. O algoritmo estándar para resolver un sistema de ecuacións lineares baséase na eliminación de Gauss con algunhas modificacións. En primeiro lugar, é esencial evitar a división por números pequenos, o que pode levar a resultados inexactos. Isto pódese facer reordenando as ecuacións se é necesario, un proceso coñecido como *pivotación*. En segundo lugar, o algoritmo non fai exactamente a eliminación de Gauss, mais calcula a descomposición LU da matriz *A*.

Se a matriz *A* ten algunha estrutura especial, esta pode ser usada para obter algoritmos máis rápidos ou precisos. Por exemplo, os sistemas cunha matriz simétrica definida positiva pódense resolver dúas veces máis rápido coa descomposición de Cholesky. A recursividade de Levinson é un método rápido para matrices de Toeplitz. Tamén existen métodos especiais para matrices con moitos elementos cero (as chamadas matrices dispersas), que aparecen a miúdo nas aplicacións.

Adoita adoptarse un enfoque completamente diferente para sistemas moi grandes, que doutro xeito levarían demasiado tempo ou memoria. A idea é comezar cunha aproximación inicial á solución (que non ten por que ser precisa en absoluto), e mudar esta aproximación en varios pasos para achegala á verdadeira solución. Unha vez que a aproximación é suficientemente precisa, esta é a solución do sistema. Isto son os métodos iterativos. Para algunhas matrices dispersas, a introdución da aleatoriedade mellora a velocidade dos métodos iterativos. Un exemplo de método iterativo é o método de Jacobi, onde a matriz A divídese na súa compoñente diagonal D e a súa compoñente non diagonal L+U. Unha suposición inicial {\bold x}^{(0)} úsase ao comezo do algoritmo. Cada iteración posterior calcúlase mediante a ecuación iterativa:

:{\bold x}^{(k+1)} = D^{-1}({\bold b} - (L+U){\bold x}^{(k)}).

Cando a diferenza entre {\bold x}^{(k)} e {\bold x}^{(k+1)} é suficientemente pequena, o algoritmo dise que *converxe* na solución.

## Sistemas homoxéneos
Un sistema de ecuacións lineares é **homoxéneo** se todos os termos constantes son cero:
:\begin{alignat}{7}
a_{11} x_1 &&\; + \;&& a_{12} x_2 &&\; + \cdots + \;&& a_{1n} x_n &&\; = \;&&& 0 \\
a_{21} x_1 &&\; + \;&& a_{22} x_2 &&\; + \cdots + \;&& a_{2n} x_n &&\; = \;&&& 0 \\
 && && && && && \vdots\;\ &&& \\
a_{m1} x_1 &&\; + \;&& a_{m2} x_2 &&\; + \cdots + \;&& a_{mn} x_n &&\; = \;&&& 0. \\
\end{alignat}
Un sistema homoxéneo equivale a unha ecuación matricial da forma
:A\mathbf{x}=\mathbf{0}
onde *A* é unha matriz , **x** é un vector columna con *n* coeficientes e ** 0** é o vector cero con *m* coeficientes.

### Conxunto de solucións homoxéneas
Todo sistema homoxéneo ten polo menos unha solución, coñecida como solución *cero* (ou *trivial*), que se obtén asignando o valor de cero a cada unha das variábeis. Se o sistema ten unha matriz non singular (), entón tamén é a única solución. Se o sistema ten unha matriz singular entón hai un conxunto de solucións cun número infinito de solucións. Este conxunto de solucións ten as seguintes propiedades adicionais:
# Se **u** e **v** son dous vectores que representan solucións a un sistema homoxéneo, entón a suma vectorial  tamén é unha solución para o sistema.
# Se **u** é un vector que representa unha solución a un sistema homoxéneo, e *r* é calquera escalar, entón *r* **u** tamén é unha solución do sistema.
Estas son exactamente as propiedades necesarias para que o conxunto de solucións sexa un subespazo linear de **R***n*. En particular, o conxunto de solucións para un sistema homoxéneo é o mesmo que o espazo nulo da matriz *A* correspondente.

### Relación con sistemas non homoxéneos
Existe unha estreita relación entre as solucións dun sistema linear e as solucións do sistema homoxéneo correspondente:
:A\mathbf{x}=\mathbf{b}\qquad \text{e}\qquad A\mathbf{x}=\mathbf{0}.
Especificamente, se **p** é algunha solución específica para o sistema linear A\mathbf{x}=\mathbf{b}, entón a totalidade do conxunto de solucións pódese describir como
:\left\{ \mathbf{p}+\mathbf{v} : \mathbf{v}\text{ é calquera solución a }A\mathbf{x}=\mathbf{0} \right\}. 
Xeometricamente, isto di que o conxunto de solucións de  é unha translación do conxunto de solucións para . En concreto, o plano para o primeiro sistema pódese obter trasladando o subespazo linear do sistema homoxéneo polo vector **p**.

Este razoamento só se aplica se o sistema  ten polo menos unha solución. Isto ocorre se e só se o vector **b** reside na imaxe da transformación linear *A*.

## Notas 
## Véxase tamén 
### Bibliografía 
* 
* 
* 
* 
* 
* 
* 
* 
* 

### Outros artigos 
* Eliminación de Gauss

### Ligazóns externas 
* Sistemas_de_Equacoes_Lineares .