MN Guia2

Universidad de O’Higgins Métodos numéricos
Escuela de Ingeniería Profesor Daniel Casagrande

Semestre de primavera 2018
GUÍA DE EJERCICIOS N°2
Resolución de sistemas de ecuaciones
1. Clasifique las siguientes matrices según su condicionamiento de acuerdo al valor

del determinante:
1 2 3 2.11 −0.8 1.72
a) [2 3 4] b) [−1.84 3.03 1.29]
3 4 5 −1.57 5.25 4.30
2 −1 0 4 3 −1
c) [−1 2 −1] d) [7 −2 3]
0 −1 2 5 −18 13
2. Use el método de Gauss-Jordan sin pivoteo para resolver el siguiente sistema:

2𝑥1 + 𝑥2 − 𝑥3 = 1
5𝑥1 + 2𝑥2 + 2𝑥3 = −4
3𝑥1 + 𝑥2 + 2𝑥3 = 5
Sustituir los resultados en las ecuaciones originales para comprobar las soluciones.
3. Dadas las siguientes ecuaciones:

2𝑥1 − 6𝑥2 − 𝑥3 = −38
−3𝑥1 − 𝑥2 + 7𝑥3 = −34
−8𝑥1 + 𝑥2 + 2𝑥3 = −20
Resolver empleando el método de Gauss con pivoteo parcial. Sustituir los
resultados en las ecuaciones originales para comprobar las soluciones.
4. Resuelva el siguiente sistema empleando el método de Gauss sin y con pivoteo

y el método de Gauss-Jordan:
𝑥1 + 𝑥2 − 𝑥3 = −3
6𝑥1 + 2𝑥2 + 2𝑥3 = 2
−3𝑥1 + 4𝑥2 + 𝑥3 = 1
5. Determine los coeficientes del polinomio 𝑦 = 𝑎0 + 𝑎1 𝑥1 + 𝑎2 𝑥 2 + 𝑎3 𝑥 3 que pase

por los puntos (0,10), (1,35), (3,31) y (4,2).
6. Determine los coeficientes del polinomio de cuarto orden que pasa por los puntos
(0,-1), (1,1), (3,3), (5,2) y (6,-2).
7. Desarrolle un algoritmo que multiplique dos matrices. Compruebe el

funcionamiento con las funciones predefinidas de MATLAB o GNU Octave.
8. Desarrolle un algoritmo que calcule la traspuesta de una matriz. Compruebe el

funcionamiento con las funciones predefinidas de MATLAB o GNU Octave.
9. Calcule el número de condición de la siguiente matriz empleando (a) la norma
euclídea y (b) la norma infinito. Puede emplear MATLAB o GNU Octave para
calcular la matriz inversa.
1 −1 −1
[0 1 −2]
0 0 1
10. Realice la descomposición LU del siguiente sistema:

10𝑥1 + 2𝑥2 − 𝑥3 = 27
−3𝑥1 − 6𝑥2 + 2𝑥3 = −61.5
𝑥1 + 𝑥2 + 5𝑥3 = −21.5
Multiplique las matrices L y U y compruebe que el resultado sea la matriz de
coeficientes. Resuelva el sistema empleando las matrices L y U. Resuelva también
para el vector 𝐛 = [12 18 − 6]𝑇 .
11. Resuelva el siguiente sistema de ecuaciones empleando el método de

descomposición en LU sin pivoteo y con el resultado calcule el determinante.
8 2 1
[3 7 2]
2 3 9
12. Realice la descomposición LU del siguiente sistema sin emplear pivoteo:

8𝑥1 + 4𝑥2 − 𝑥3 = 11
−2𝑥1 + 5𝑥2 + 𝑥3 = 4
2𝑥1 − 𝑥2 + 6𝑥3 = 7
Determine la inversa de A y compruebe que A por su inversa es la matriz identidad.
13. Emplee descomposición LU sin pivoteo para determinar la matriz inversa del
sistema, compruebe el resultado.
10𝑥1 + 2𝑥2 − 𝑥3 = 27
−3𝑥1 − 6𝑥2 + 2𝑥3 = −61.5
𝑥1 + 𝑥2 + 5𝑥3 = −21.5
14. El siguiente sistema de ecuaciones se emplea para determinar las

concentraciones (c en [g/m3]) en una serie de reactores acoplados en función de
la cantidad de masa ingresada a cada reactor (vector b en [g/día]).
15𝑐1 − 3𝑐2 − 𝑐3 = 3800
−3𝑐1 + 18𝑐2 − 6𝑐3 = 1200
−4𝑐1 − 𝑐2 + 12𝑐3 = 2350
a) Determine la inversa de A.
b) Emplee la matriz inversa para calcular la solución.
c) Determine cuánto se debe incrementar el ingreso de masa al reactor 3 para
inducir un incremento de 10 g/m3 en la concentración del reactor 1.
d) Determine cuánto se reducirá la concentración en el reactor 3 si el ingreso
de masa en los reactores 1 y 2 se reduce por 500 y 250 g/día
respectivamente.
15. Resuelva el siguiente sistema de ecuaciones empleando el método de
descomposición en LU:
3𝑥1 − 2𝑥2 + 𝑥3 = −10
2𝑥1 + 6𝑥2 − 4𝑥3 = 44
−𝑥1 − 2𝑥2 + 5𝑥3 = −26
16. Los cuatro bloques de la figura están conectados por cuerdas de masa
despreciable. Tres de los bloques se encuentran en un plano inclinado, cuyo
coeficiente de fricción está dado por el vector 𝝁𝑖 . Las ecuaciones de movimiento
para los bloques son las siguientes:
𝑇1 + 𝑚1 𝑎 = 𝑚1 𝑔(sin 𝜃 − 𝜇1 cos 𝜃)
−𝑇1 + 𝑇2 + 𝑚2 𝑎 = 𝑚2 𝑔(sin 𝜃 − 𝜇2 cos 𝜃)
−𝑇2 + 𝑇3 + 𝑚3 𝑎 = 𝑚3 𝑔(sin 𝜃 − 𝜇3 cos 𝜃)
−𝑇3 + 𝑚4 𝑎 = −𝑚4 𝑔
Donde 𝐓𝑖 son las fuerzas en las cuerdas y 𝑎 es la aceleración. Determine 𝐓𝑖 y 𝑎

para 𝜃 = 45°, 𝑔 = 9.81 𝑚/𝑠 y
𝐦 = [10 4 5 6]𝑇 𝑘𝑔
𝛍 = [0.25 0.3 0.2]𝑇
17. Resuelva el siguiente sistema tridiagonal empleando el método de

descomposición en LU y el método iterativo de Gauss-Seidel:
0.8𝑥1 − 0.4𝑥2 = 41
−0.4𝑥1 + 0.8𝑥2 − 0.4𝑥3 = 25
−0.4𝑥2 + 0.8𝑥3 = 105
Calcule la matriz inversa empleando la descomposición LU.
18. El siguiente es un sistema tridiagonal que debe ser resuelto como parte de
un algoritmo más grande (método de Crank-Nicholson) empleado para resolver
ecuaciones diferenciales parciales:
2.01475 −0.020875 0 0 𝑇1 4.175
−0.020875 2.01475 −0.020875 0 𝑇 0
[ ] × [ 2] [ ]
0 −0.020875 2.01475 −0.020875 𝑇3 0
0 0 −0.020875 2.01475 𝑇4 2.0875
Encuentre el vector de soluciones 𝐓𝑖 .
19. Realice la descomposición de Cholesky de la siguiente matriz:

8𝑥1 + 20𝑥2 + 15𝑥3 = 50
20𝑥1 + 80𝑥2 + 50𝑥3 = 250
15𝑥1 + 50𝑥2 + 60𝑥3 = 100
20. Realice la descomposición de Cholesky de la siguiente matriz simétrica y
compruebe que el producto de L y su traspuesta dan como resultado la matriz de
coeficientes A:
16 15 55
[15 55 225]
55 225 979
Resuelva el sistema considerando un vector 𝐛 = [152.6 585.6 2488.8]𝑇 .
21. El siguiente circuito eléctrico puede ser analizado como tres mallas, y
aplicando las leyes de Kirchoff a cada una de ellas se pueden obtener las
siguientes ecuaciones:
(50 + 𝑅)𝑖1 − 𝑅𝑖2 − 30𝑖3 = 0

−𝑅𝑖1 + (65 + 𝑅)𝑖2 − 15𝑖3 = 0
−30𝑖1 − 15𝑖2 + 45𝑖3 = 120
Resuelva el sistema para 𝑅 = 5 Ω, 𝑅 = 10 Ω y 𝑅 = 20 Ω.
22. Resuelva el sistema del problema 17 empleando el método de Gauss-Seidel

con sobrerelajación. Emplee 𝜆 = 1.2 y un error 𝜀 = 5%.
23. Resuelva el sistema del problema 14 empleando el método de Gauss-Seidel

con 𝜀 = 5%.
24. Resuelva el siguiente sistema empleando el método de Gauss-Seidel hasta

que el error relativo sea menor al 5%.
10𝑥1 + 2𝑥2 − 1𝑥3 = 27
−3𝑥1 − 6𝑥2 + 2𝑥3 = −61.5
𝑥1 + 𝑥2 + 5𝑥3 = −21.5
25. Emplee el método de Gauss-Seidel sin y con relajación (𝜆 = 0.95) para

resolver el siguiente sistema hasta que el error relativo sea menor al 5%. Si es
necesario reacomode las ecuaciones para lograr la convergencia.
10𝑥1 + 2𝑥2 − 𝑥3 = 27
−3𝑥1 − 6𝑥2 + 2𝑥3 = −61.5
𝑥1 + 𝑥2 + 5𝑥3 = −21.5
26. Emplee el método de Gauss-Seidel sin y con relajación (𝜆 = 1.2) para

resolver el siguiente sistema hasta que el error relativo sea menor al 5%. Si es
necesario reacomode las ecuaciones para lograr la convergencia.
2𝑥1 − 6𝑥2 − 𝑥3 = −38
−3𝑥1 − 𝑥2 + 7𝑥3 = −34
−8𝑥1 + 𝑥2 − 2𝑥3 = −20
27. Considerando una matriz cuadrada A, desarrolle un comando que devuelva
la matriz de coeficientes aumentada con una matriz identidad empleando sólo
una línea de código para MATLAB o GNU Octave.
28. Escriba el siguiente set de ecuaciones en forma matricial y emplee Matlab o

GNU Octave para encontrar las soluciones, calcular la traspuesta y la inversa de
la matriz de coeficientes.
50 = 5𝑥3 − 7𝑥2
4𝑥2 + 7𝑥3 + 30 = 0
𝑥1 − 7𝑥3 = 40 − 3𝑥2 + 5𝑥1
29. El diagrama muestra 5 vasijas de mezclado conectadas por tuberías. Se

suministras agua con flujo constante a través de los tubos, tal como se muestra
en la figura. El agua que entra al sistema contiene un químico cuya concentración
se especifica a través de c [mg/m3].
Aplicando el principio de conservación de masa (la masa entrante es igual a la masa

que sale) a cada vasija se obtiene el siguiente sistema de ecuaciones para la
concentración 𝑐𝑖 en cada vasija:
−8𝑐1 + 4𝑐2 = −80
8𝑐1 − 10𝑐2 + 2𝑐3 = 0
6𝑐2 −11𝑐3 + 5𝑐4 = 0
3𝑐3 −7𝑐4 + 4𝑐5 = 0
2𝑐4 − 4𝑐5 = −30
Tenga en cuenta que el flujo de materia entrante con el químico está dado por el
producto del caudal de agua y concentración del químico. Verifique las ecuaciones
y determine las concentraciones en cada vasija.
30. Los bordes del panel mostrado en la figura se mantienen a las temperaturas
indicadas.
Asumiendo que la transferencia de calor se encuentra en estado estable, la
ecuación diferencial que gobierna la temperatura 𝑇 en el interior está dada por:
𝜕 2𝑇 𝜕 2𝑇
+ =0
𝜕𝑥 2 𝜕𝑦 2
Usando la malla indicada y aproximando esta ecuación por diferencias finitas, se
obtiene el siguiente sistema de ecuaciones algebraicas para la temperatura en
los nodos de la malla:
Resuelva el sistema empleando el método del gradiente conjugado.

MN Guia2

Cargado por

Copyright:

Formatos disponibles

MN Guia2

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

MN Guia2

Cargado por

Copyright:

Formatos disponibles

Universidad de O’Higgins Métodos numéricos

Escuela de Ingeniería Profesor Daniel Casagrande

Resolución de sistemas de ecuaciones

1. Clasifique las siguientes matrices según su condicionamiento de acuerdo al valor

2. Use el método de Gauss-Jordan sin pivoteo para resolver el siguiente sistema:

3. Dadas las siguientes ecuaciones:

4. Resuelva el siguiente sistema empleando el método de Gauss sin y con pivoteo

5. Determine los coeficientes del polinomio 𝑦 = 𝑎0 + 𝑎1 𝑥1 + 𝑎2 𝑥 2 + 𝑎3 𝑥 3 que pase

7. Desarrolle un algoritmo que multiplique dos matrices. Compruebe el

8. Desarrolle un algoritmo que calcule la traspuesta de una matriz. Compruebe el

10. Realice la descomposición LU del siguiente sistema:

11. Resuelva el siguiente sistema de ecuaciones empleando el método de

12. Realice la descomposición LU del siguiente sistema sin emplear pivoteo:

14. El siguiente sistema de ecuaciones se emplea para determinar las

Donde 𝐓𝑖 son las fuerzas en las cuerdas y 𝑎 es la aceleración. Determine 𝐓𝑖 y 𝑎

17. Resuelva el siguiente sistema tridiagonal empleando el método de

19. Realice la descomposición de Cholesky de la siguiente matriz:

(50 + 𝑅)𝑖1 − 𝑅𝑖2 − 30𝑖3 = 0

Resuelva el sistema para 𝑅 = 5 Ω, 𝑅 = 10 Ω y 𝑅 = 20 Ω.

22. Resuelva el sistema del problema 17 empleando el método de Gauss-Seidel

23. Resuelva el sistema del problema 14 empleando el método de Gauss-Seidel

24. Resuelva el siguiente sistema empleando el método de Gauss-Seidel hasta

25. Emplee el método de Gauss-Seidel sin y con relajación (𝜆 = 0.95) para

26. Emplee el método de Gauss-Seidel sin y con relajación (𝜆 = 1.2) para

28. Escriba el siguiente set de ecuaciones en forma matricial y emplee Matlab o

29. El diagrama muestra 5 vasijas de mezclado conectadas por tuberías. Se

Aplicando el principio de conservación de masa (la masa entrante es igual a la masa

Resuelva el sistema empleando el método del gradiente conjugado.

También podría gustarte