Codigo Matlab Factorizacion Lu Cholesky
Codigo Matlab Factorizacion Lu Cholesky
Codigo Matlab Factorizacion Lu Cholesky
SISTEMAS DE ECUACIONES
M E TO D O S
N U M E R IC O S
O2 - 5
CYNDY ARGOTE
JONATHAN CELIS
JONATHAN PEREZ
JHONATAN QUINTERO
LINA MARGARITA GOMEZ
DESCOMPOSICION LU
1. Hacer cero todos los valores abajo del pivote sin convertir
este en 1.
2. Para lograr lo anterior se requiere obtener un factor el cual
es necesario para convertir a cero los valores abajo del
pivote.
3. Dicho factor es igual al número que se desea convertir en
cero entre el número pivote.
4. Este factor multiplicado por -1 se multiplica luego por el
pivote y a ese resultado se le suma el valor que se encuentra
en la posición a cambiar (el valor en la posición que se
convertirá en cero).
PASOS PARA ENCONTRAR LA MATRIZ
TRIANGULAR INFERIOR (MATRIZ [L])
4 -2 -1 9
A=5 1 -1 b= 7
1 2 -1 12
FACTORIZACION LU
EJEMPLO N°1
• SOLUCION
1. Se halla “U”:
4 -2 -1
5 1 -1 R2 R2 – (5/4)*R1
1 2 -1 R3 R3 – (1/4)*R1
4 -2 -1
0 7/2 ¼
0 5/2 -3/4 R3 R3 – (5/2)/(7/2)*R2
FACTORIZACION LU
EJEMPLO N°1
SOLUCION
1. Se halla “U”:
4 -2 -1
U= 0 7/2 ¼
0 0 -13/14
2. Se halla “L”:
1 0 0 1 0 0
L =? 1 0 L = 5/4 1 0
? ? 1 ¼ 5/7 1
FACTORIZACION LU
EJEMPLO N°1
3. Se verifica L*U = A
1 0 0 4 -2 -1
5/4 1 0 x 0 7/2 ¼ =
¼ 5/7 1 0 0 -13/14
Y1 =9 Y1 = 9
5/4Y1 + Y2 =7 Y2 = -17/4
1/4Y1 + 5/7Y2 +Y3 = 12 Y3 = 179/14
FACTORIZACION LU
EJEMPLO N°1
11 -3 -2 18
A=5 -2 -8 b= 13
4 -7 2 2
FACTORIZACION LU
EJEMPLO N°2
SOLUCION
1. Se halla “U”:
11 -3 -2
5 -2 -8 R2 R2 – (5/11)*R1
4 -7 2 R3 R3 – (4/11)*R1
11 -3 -2
0 -7/11 -78/11
0 -65/11 30/11 R3 R3 – (-65/11)/(-7/11)*R2
FACTORIZACION LU
EJEMPLO N°2
SOLUCION
1. Se halla “U”:
11 -3 -2
U= 0 -7/11 -78/11
0 0 480/7
2. Se halla “L”:
1 0 0 1 0 0
L =? 1 0 L = 5/11 1 0
? ? 1 4/11 65/7 1
FACTORIZACION LU
EJEMPLO N°2
SOLUCION
3. Se verifica L*U = A
1 0 0 11 -3 -2
5/11 1 0 x 0 -7/11 -78/11 =
4/11 65/7 1 0 0 480/7
SOLUCION
4. Se despeja “Y” de L*Y = b
1 0 0 Y1 18
5/11 1 0 * Y2 = 13
4/11 65/7 1 Y3 2
Y1 = 18 Y1 = 18
5/11Y1 + Y2 = 13 Y2 = 53/11
4/11Y1 + 65/7Y2 +Y3 =2 Y3 = -345/7
FACTORIZACION LU
EJEMPLO N°2
SOLUCION
5. Se despeja “X” de U*X = Y
11 -3 -2 X1 18
0 -7/11 -78/11 * X2 = 53/11
0 0 480/7 X3 -345/7
j = k+1, n
Sustitución hacia
2 i = 2, n delante:
L*Y = b
sum = bi
j = 1, i-1
Se almacenan los
bi = sum nuevos “b” (Y)
3
FACTORIZACION LU
DIAGRAMA DE FLUJO
Xn = bn/An,n Sustitución hacia
3
atrás:
U*X = Y
i = n-1, 1, -1 i = 1, n
sum = 0 Xi
j = i+1, n
Fin!
Xi = (bi – sum)/Ai,j
FACTORIZACION LU
Referencias de consulta
• http://www.unizar.es/aragon_tres/unidad6/Matrices/u
6matte20.pdf
• http://www.ditutor.com/matrices/matriz_simetrica.ht
ml
• http://www.cramster.com/reference/wii.aspx?
wiki_name=Band_matrix
• Chapra, Steven; Canale, Raymond. Métodos númericos
para ingenieros. 3ra Edición. Mc Graw Hill 2000.
GRACIAS