Álgebra Lineal con Scilab Scilab herramienta informática para el álgebra lineal

Con la tecnología de Blogger.

Uso de la factorización LU para resolver un sistema de ecuaciones


Suponga que se quiere resolver el sistema Ax = b, donde A es invertible. Si A satisface la hipótesis del teorema se puede escribir
Como L es invertible, existe un vector único y tal que Ly = b. Como U también es invertible, existe un vector único x tal que Ux = y. Entonces Ax = L(Ux) = Ly = b y nuestro sistema está resuelto. Observe que Ly = b se puede resolver directamente mediante la sustitución hacia adelante, mientras que el sistema Ux= y se puede resolver por sustitución hacia atrás. Esto se ilustra en el siguiente ejemplo.

Ejemplo: Uso de la factorización LU para resolver un sistema
Resuelva el sistema Ax = b, donde



Solución: Del ejemplo se puede escribir A = LU, donde




El sistema Ly=b conduce a las ecuaciones
Se acaba de realizar la sustitución hacia delante. Ahora, de Ux= y se obtiene
La solución es:
En scilab:  Vamos a comprobor los resultados con Scilab

En la anterior entrada usamos una función fact_LU(A) para obtener las Matrices L y U. Con este cálculo obtenemos y con la siguiente formula:

en la cual se obtienen los valores de y con la inversa de L.
A=[
2,3,2,4;
4,10,-4,0;
-3,-2,-5,-2;
-2,4,4,-7;
]
B=[4,-8,-4,-1]';
m=fact_LU(A);

y=m.L\B
disp(y) 

Obteniendose:


Y de la misma manera obtemos x:
Código:
 
A=[
2,3,2,4;
4,10,-4,0;
-3,-2,-5,-2;
-2,4,4,-7;
]
B=[4,-8,-4,-1]';
m=fact_LU(A);

y=m.L\B
x=m.U\y
disp(x)


Resultado:


NOTA: el segundo valor nos da un valor muy próximo a cero




Share on Google Plus

Acerca del Autor Matemática Positiva

Sitio Web Dedicado al Universo Matemático

0 comments :

Publicar un comentario