[justify]El problema central del álgebra lineal es la solución de ecuaciones y es más sencillo cuando el numero de incógnitas es igual al número de ecuaciones[br]empecemos con n = 2:[br][/justify][center][br][math][br]\large [br]1x + 2y=3 \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \:ecuación 1\\[br]\large [br]4x+5y=6 \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: \: ecuación 2[br][/math][br][/center][br]las incógnitas son x y y. Para resolver estas ecuaciones podemos usar el método de eliminación o el de determinantes[br]1. [b]Eliminación :[/b]De la segunda ecuación, reste 4 veces la primera ecuación. Así se elimina x de la segunda ecuación, y queda una ecuación para y:[br][center][br][math][br]\large (ecuación 2) - 4(ecuación 1) \hspace{2cm} -3y=-6 \\[br][br][/math][br][/center] [br]De inmediato se sabe que [math]\large y=2[/math], Luego se averigua [math]\large x[/math] a partir de la primera ecuación, Con sustitución hacia atrás:[br][br][center][br][math][br]\large \text{sustitución hacia atrás} \hspace{1.5cm} 1x+2(2) = 3 \hspace{1.5cm} x+4=1\hspace{1.5cm} x=-1\\[br][br][/math][br][/center][br][br]Procediendo cuidadosamente, se comprueba que [math]\large x[/math] e [math]\large y[/math], cumplen también la segunda ecuación, esto debe funcionar, como es el caso: 4 veces ([math]\large x=-1 [/math] ) mas cinco veces ([math]\large y=2[/math],) es igual a 6.[br][br][br]2. [b]Determinantes[/b] : La solución [math]\large y=2[/math], depende completamente de los seis números en las ecuaciones. Debe existir una fórmula para [math]\large y[/math], (y también para [math]\large x [/math],). Se trata de una razón de determinantes, que espero, el lector me permita escribir directamente:[br][center][br][math][br]\large [br]y=\frac{\begin{vmatrix}[br] 1 & 3\\ [br] 4 & 6[br]\end{vmatrix}}{\begin{vmatrix}[br] 1 & 2\\ [br]4 & 5[br]\end{vmatrix}}= \frac{1.6-3.4}{1.5-2.4}=\frac{-6}{-3}=2[br][/math][br][/center][br][br]Lo anterior puede parecer algo misterioso, a menos que el lector ya conozca algo sobre el cálculo de determinantes de 2 por 2. Estos determinantes proporcionan la misma respuesta y = 2, proveniente de la misma razón de -6 a -3. Si nos quedamos con los determinantes (lo cual no pensamos hacer), hay una fórmula semejante para calcular la otra incógnita, [math]\large x[/math],:[br][br][center][br][math][br]\large[br]x=\frac{\begin{vmatrix}[br] 3 & 2\\ [br] 6 & 5[br]\end{vmatrix}}{\begin{vmatrix}[br] 1 & 2\\ [br]4 & 5[br]\end{vmatrix}}= \frac{3.5-2.6}{1.5-2.4}=\frac{3}{-3}=-1[br][/math][br][/center][justify][/justify][justify][/justify]
[justify]A continuación se compararán ambos métodos, pensando en futuros problemas reales en los que n es mucho más grande (n = 1000 es un tamaño bastante moderado en cálculos científicos). Lo cierto es que el uso directo de la fórmula de los determinantes para 1000 ecuaciones puede ser un desastre total, ya que el millón de números a la izquierda se utilizaría correcta pero ineficazmente. Esta fórmula se encontrará en el capítulo 4 (regla de Cramer), aunque en el capítulo 1 se presenta un método aceptable para resolver 1000 ecuaciones. Este método aceptable es la eliminación gaussiana. Se trata del algoritmo que suele aplicarse de manera constante para resolver grandes sistemas de ecuaciones. A partir de los ejemplos en un libro de texto (n = 3 se aproxima al limite superior de la paciencia del autor y del lector), quizá el lector no puede apreciar mucha diferencia. En las ecuaciones (2) y ( 4) se siguieron esencialmente los mismos pasos para encontrar y = 2. Ciertamente x se conoció más rápido por la sustitución hacia atrás en la ecuación (3) que por la razón en (5). Para sistemas con más incógnitas y ecuaciones no hay caso, Gana la eliminación (a propósito el mejor método para calcular determinantes ¡¡¡es eliminación!!!). [br]La idea de eliminación es engañosamente simple, el lector la dominará luego de unos cuantos ejemplos. Constituye la base de la mitad de este libro, su objetivo es simplificar (maximizar el número de ceros) de una matriz de modo que sea posible resolver un sistema de ecuaciones, pero esta matriz nos ayudará a entender también las nociones de espacio vectorial, dimensión, independencia etc... Junto con la mecánica del algoritmo, en este capítulo es necesario explicar cuatro aspectos más profundos.[/justify]
[justify]1. Las ecuaciones lineales llevan a la geometría de planos. No es fácil visualizar un plano nueve-dimensional en un espacio de diez dimensiones. Es más difícil ver diez de estos planos, que se cortan en la solución de diez ecuaciones, aunque de alguna manera esto es casi posible. Nuestro ejemplo tiene dos rectas en la figura , que se encuentran en el punto (x, y) = (-1, 2). El álgebra lineal mueve esta imagen hacia diez dimensiones, donde la intuición debe imaginar la geometría.[br][br]2. Pasamos a la notación matricial al escribir las n incógnitas como un vector [math]\large \vec{x}[/math] y las n ecuaciones como [math]\large A\vec{x}=\vec{b}[/math]. Transformamos A con la finalidad de obtener una matriz con ceros bajo la diagonal U. Con lo anterior, resolver el sistema de ecuaciones es inmediato.[br][br]3. En la mayor parte de los casos, la eliminación se realiza sin dificultades. La matriz tiene una inversa, y el sistema [math]\large A\vec{x}=\vec{b}[/math] tiene una solución. En otros casos el método falla, en este caso las ecuaciones no tienen una solución única. El caso singular aparece si 8 se sustituye por 5 en nuestro ejemplo:[br][br][math]\large \text{ Caso Singular \: \: \: \: \: \: \: \: \: \: 1x+2y=3 }[/math][br][math]\large \text{ Dos rectas paralelas \: \: \: \: 4x+8y=6 }[/math][br][br]La eliminación resta simplemente 4 veces la primera ecuación de la segunda. Sin embargo, ¡observe el resultado!:[br][br](ecuación 2) - 4(ecuación 1) 0=-6.[br][br][br]Este caso singular no tiene solución. Otros casos singulares tienen una cantidad infinita de soluciones. (Cambie 6 por 12 en el ejemplo anterior y la eliminación producirá O = O) . Entonces la variable y puede asumir [b]cualquier valor.[/b]) Cuando la eliminación falla, se quiere encontrar toda solución posible.[br][br]4. Se requiere una estimación aproximada del número de pasos de eliminación necesarios para resolver un sistema de tamaño n. El costo de cómputo a menudo determina la precisión del modelo. Cien ecuaciones requieren alrededor de 300 000 pasos (multiplicaciones y restas). La computadora es capaz de hacer estos pasos rápidamente, pero no es así para el caso de varios billones de pasos. Y después de un millón de pasos, el error por redondeo puede ser significativo. (Algunos problemas son sensibles; otros no.) Sin entrar en todos los detalles, pretendemos considerar grandes sistemas que se presentan en la práctica, así como la manera en que se resuelven realmente.[br][br]El resultado final de este capítulo, es un algoritmo de eliminación que es lo más eficaz posible. Se trata del algoritmo que suele usarse en una numerosa variedad de aplicaciones. Y al mismo tiempo, comprenderlo en términos de la matriz de coeficientes A es un fundamento esencial de la teoría. Pasamos entonces a explicar de manera general el método de eliminación.[/justify]