Cálculo Numérico
Universidad Nacional del Comahue –Facultad de Informática- Depto de Computación Aplicada-
TRABAJO PRACTICO 3: Sistemas de Ecuaciones Lineales
1. Resolver A x = b por el método de eliminación de Gauss, utilizando calculadora y
trabajando con 5 cifras significativas, siendo A | b la matriz ampliada del sistema de 3
ecuaciones con 3 incógnitas.
A|b
3 1 0 -2
-2 0 2 1
1 3 4 0
a) ¿Es aconsejable transformar en 1 los elementos de la diagonal? Justifique
b) Resuelva el sistema utilizando Scilab y compare los resultados.
2. Implementar el algoritmo de eliminación de Gauss en Scilab.
3. Resolver el siguiente sistema utilizando:
a. Gauss
b. Gauss con Pivoteo parcial
c. Comparar los resultados
Obs: Trabajar con 3, 4 y 5 cifras significativas:
0.0003 𝑥1 + 3.0000 𝑥2 = 2.0001
1.0000 𝑥1 + 1.0000 𝑥2 = 1.0000
4. Resolver el siguiente sistema utilizando Gauss Jordan:
3 𝑥1 + 𝑥2 = −2
−2 𝑥1 + 2 𝑥3 = 1
𝑥1 + 3𝑥2 + 4𝑥3 = 0
5. Dado el siguiente sistema:
−𝑥1 + 2 𝑥2 + 42 𝑥3 = 83
72 𝑥1 − 41 𝑥2 − 14 𝑥3 = 44
35 𝑥1 + 10 𝑥2 − 5 𝑥3 = 25
a) Resolver utilizando Gauss Jordan
b) Utilizando Gauss Jordan con pivoteo parcial (trabajar con tres cifras significativas)
c) Comparar los resultados
6. Descomponer cada una de las siguientes matrices, en una matriz superior U y una
matriz inferior L. Comprobar resultados con Scilab.
2 −1 1 2 −1.5 3
[3 3 9] [−1 0 2]
3 3 5 4 −4.5 5
Cálculo Numérico
Universidad Nacional del Comahue –Facultad de Informática- Depto de Computación Aplicada-
7. Implementar los algoritmos en Scilab para para resolver los sistemas triangulares:
L y = b (substituciones hacia delante) ; U x = y (substituciones hacia atrás)
8. Resolver el siguiente sistema utilizando descomposición LU.
4 𝑥1 − 2𝑥2 − 𝑥3 = 9
5𝑥1 + 𝑥2 − 𝑥3 = 7
𝑥1 + 2𝑥2 − 𝑥3 = 12
9. Implementar los algoritmos en Scilab para generar las matrices de la descomposición L
y U. Luego verifica los resultados obtenidos en el ejercicio anterior.
10. El siguiente ejemplo muestra un fenómeno llamado “mal condicionamiento” de un
sistema, por el cual, pequeños cambios en los elementos de A o del vector b producen
grandes cambios en el vector solución, a pesar de que los cálculos se realizan con
exactitud.
Considerando el siguiente sistema A | b:
𝑥1 − 11 𝑥2 = −1
9 𝑥1 + 100 𝑥2 = 9.1
a) Compruebe que la solución es x = (0.1, 0.1)T
b) ¿Cómo cambia la solución del sistema si en lugar de considerar a b= (-1, 9.1)T,
consideramos b’=(-1, 9)T
c) Analice la “sensibilidad” de la solución cuando se cambia a12 por -11.1.
d) Usando Scilab calcular el número de condición en cada matriz de los sistemas.
Comentar los resultados.
11. Resolver el siguiente sistema por Jacobi y por Gauss Seidel, verificando previamente el
criterio de convergencia.
𝑥1 − 4𝑥2 − 𝑥3 = 11
2 𝑥2 + 3𝑥3 = 0
8 𝑥1 + 2𝑥2 = 2
12. Implementar en Scilab los métodos de Jacobi y Gauss Seidel. Luego compruebe los
resultados del ejercicio anterior.
13. Con ayuda de Scilab, hallar las matrices de Jacobi del sistema A x= f, donde A y f son:
2 −1 0 1
[−1 2 −1] ; [1]
0 −1 2 5
Recordatorio: Bj =-(L + U) en el proceso iterativo X(k+1)=B X(k)+c = -(L+U) x(k)+ D-1f
Cálculo Numérico
Universidad Nacional del Comahue –Facultad de Informática- Depto de Computación Aplicada-
14. Dado el sistema lineal:
𝑥1 − 𝑥2 = −1
𝑥1 − 6 𝑥2 = -24
Y el valor inicial x(0)=(0,0)
a) Calcular mediante la forma algebraica de Gauss Seidel las siguientes 3
iteraciones.
b) Sabiendo que la matriz de iteración del método es : BGS= - (I+L)-1 * U, encontrar
una condición de convergencia sobre dicha matriz. (utilizar la ayuda de un
software matemático para realizar los cálculos)
c) ¿Qué conclusión podría realizar considerando los incisos anteriores?
15. Una aplicación sencilla de sistemas de ecuaciones lineales a las leyes de Kirchooff para
circuitos eléctricos. Se aplica un potencial de 15 Volts entre los puntos A y G de un
circuito, formado por los nodos A,B,C,D,E, y F ; sabe que los potenciales en dichos
puntos son vb, vc, vd,, ve y vf , respectivamente. Entonces dichos potenciales satisfacen
el siguiente sistema de ecuaciones:
31 Vb - 10 𝑉c -6 𝑉f = 15
2 𝑉𝑏 − 8 𝑉𝑐 + 3 𝑉e = 0
𝑉𝑐 − 3 𝑉d + 2𝑉𝑒 = 0
2𝑉𝑐 + 4 𝑉𝑑 - 7 𝑉𝑒 + 𝑉𝑓 = 0
12𝑉b + 15 𝑉𝑑 − 47 𝑉𝑓 = 0
Determinar los potenciales en los nodos A, B, C, D, E y F del circuito, usando los
métodos de Jacobi y Gauss-Seidel. Comparar la rapidez de convergencia.
