Universidad Nacional de Salta Programación Numérica
Facultad de Ciencias Exactas Cálculo Numérico
Departamento de Informática
Trabajo Práctico N° 5
Resolución de Sistemas de Ecuaciones Lineales
Ejercicio N° 1. Sea la resolución de A· x = b, con |A|0. Expresar el modelo matemático y
determinar el número de operaciones para los siguientes casos, en donde A es:
• Diagonal
• Triangular Superior
• Triangular Inferior
Ejercicio N° 2. Investigar la demostración del teorema de la descomposición LU para una
matriz cuadrada de orden n.
Ejercicio Nº 3. Descomponer la siguiente matriz en L·U, tal que A = L·U.
6 −2 2 4
12 −8 6 10
A=
3 −13 9 3
−6 4 1 −18
Ejercicio N° 4. En la resolución de A· x = b, con |A|0. Utilice el método de Gauss y
Gauss-Jordan para encontrar la solución. Utilícese la matriz del ejercicio anterior y
bT = (12, 34, 27, −38).
Ejercicio Nº 5. Determine el modelo matemático y el número de operaciones de los
métodos de Crout y Cholesky
Ejercicio N° 6. Resolver los siguientes sistemas lineales:
2 x1 + x2 − 5 x3 + x4 = 8 0.05 x1 + 0.07 x2 + 0.06 x3 + 0.05 x4 = 0.23
x − 3x − x4 = 9 0.07 x + 0.1 x + 0.08 x + 0.07 x = 0.32
a) 1 2
b) 1 2 3 4
2 x2 − x3 + 2 x4 = −5 0.06 x1 + 0.08 x 2 + 0.1 x3 + 0.09 x4 = 0.33
x1 + 4 x2 − 7 x3 + 6 x4 = 0 0.05 x1 + 0.07 x2 + 0.09 x3 + 0.1 x4 = 0.31
Mediante los siguientes métodos:
1. Gauss con pivoteo parcial. 4. Crout
2. Gauss con pivoteo completo. 5. Cholesky
3. Gauss - Jordan.
Ejercicio Nº 7. Sean los siguientes problemas:
10 7 8 7 x1 32 10 7 8 7 x1 32.1
a)
7 5 6 5 x2 = 23 b)
7 5 6 5 x2 = 22.9
8 6 10 9 x3 33 8 6 10 9 x3 33.1
7 5 9 10 x4 31 7 5 9 10 x4 30.9
¿Son problemas equivalentes?
¿Las soluciones son aproximadamente iguales? Justifique.
1
Universidad Nacional de Salta Programación Numérica
Facultad de Ciencias Exactas Cálculo Numérico
Departamento de Informática
Ejercicio Nº 8. Dado el siguiente sistema lineal, obtener la solución aproximada mediante
los métodos que se indica, usando x 0 = (0, 0, 0, 0) y xn+1 − xn 0.001:
• Jacobi • Gauss-Seidel • Relajamiento
10 x1 − x2 + 2 x3 = 6
− x1 + 11x2 − x3 + 3x4 = 25
2 x1 − x2 + 10 x3 − x4 = −11
3x2 − x3 + 8 x4 = 15
⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯⎯
Bloque 2 – Programación
— Diseñar programas para:
• Determinar normas de VECTORES.
• Determinar normas de MATRICES.
• Resolver SISTEMAS LINEALES mediante los métodos de:
Gauss Gauss – Jordan
Crout Cholesky
Gauss-Seidel Jacobi
Relajamiento Mejoramiento Iterativo
SOR
Nota: Para los métodos directos programar estrategias de pivoteo:
• Pivoteo Simple
• Pivoteo Parcial
• Pivoteo Total