Facultad de Ciencias
Departamento de Matemáticas y Ciencias de la Computación
TALLER II
El desarrollo del presente taller lo puede realizar en cualquiera de las 2 opciones siguiente:
i. De forma manuscrita digitalizando posteriormente sus respuestas y pegando las imágenes
con el desarrollo y la respuesta correspondiente en cada recuadro. Los códigos MATLAB
deben ser pegados en los recuadros solicitados, junto con su respuesta.
ii. De forma digital con el desarrollo y la respuesta correspondiente en cada recuadro. Los
códigos MATLAB deben ser pegados en los recuadros solicitados, junto con su respuesta.
En cualquiera de las dos opciones debe detallar (explicar) todo su desarrollo y adjuntar gráficos y
códigos utilizados.
Es importante que cumpla con los plazos definidos en plataforma, y entregue el archivo en formato
PDF.
Tendrá disponible para consultas tanto la sesión con el profesor, como la sesión de ayudantía en la
respectiva semana.
1. Aproxime la solución del siguiente sistema lineal
6𝑥 − 2𝑦 + 𝑧 = 11
−2𝑥 + 7𝑦 + 2𝑧 = 5
𝑥 + 2𝑥 − 5 𝑧 = −1
a) Usando el método iterativo de Jacobi, con cinco iteraciones. Estudie la convergencia del
método para este sistema.
Se tiene que analizando por norma infinito este nos arroja un valor 0.6 por lo cual por este criterio
la convergencia es factible por este método.
Luego por radio espectral también es < 0 y más cercano a cero, lo que corrobora la conclusión de
convergencia por norma infinito.
Para las soluciones tenemos que:
b) Usando el método de Gauss – Seidel.
Para Gauss-Seidel tenemos que los criterios de convergencia del método son los siguientes:
Se observa los valores de norma infinito y radio espectral son más cercanos a cero que en Jacobi
Analizando los resultados se tiene que el criterio de convergencia al ser más cercano a cero en
gauss-saidel (error=0.00130) con 5 iteraciones llegamos a un error más pequeño que con Jacobi
(error=0.0367).
c) Indique el error cometido en su aproximación usando norma infinita.
El error en norma infinito se denota de la siguiente manera:
Máx {‖𝑥 𝑘+1 − 𝑥 𝑘 ‖; ‖𝑦 𝑘+1 − 𝑦 𝑘 ‖; ‖𝑧 𝑘+1 − 𝑧 𝑘 ‖}∞
Gauss-Seidel: Máx {‖𝑥 5 − 𝑥 4 ‖; ‖𝑦 5 − 𝑦 4 ‖; ‖𝑧 5 − 𝑧 4 ‖}∞ = {1.3 ∗ 10−3 }
Jacobi = {3.6 ∗ 10−2 }
En conclusión, tenemos que el error del método Gauss-Seidel es menor a Jacobi al aplicar la
misma cantidad de iteraciones.
2. un ingeniero civil que trabaja en la construcción requiere 4800, 5800 y 5700 m 3 de arena,
grava fina y grava gruesa, respectivamente, para cierto proyecto. Hay tres canteras de las
que puede obtenerse dichos materiales. La composición de dichas canteras es la siguiente
Arena (%) Grava gruesa (%) Grava fina (%)
Cantera 1 55 15 20
Cantera 2 25 30 60
Cantera 3 25 55 20
¿Cuantos metros cúbicos debe extraer de cada cantera a fin de satisfacer las necesidades
del ingeniero?
a) Platee el sistema de ecuación que define la situación del ingeniero
Ec1: 55x + 25y + 25z = 4800 (Ecuación de arena)
Ec2: 15x + 30y + 55z = 5800 (Ecuación de grava fina)
Ec3: 20x + 60y + 20z = 5700 (Ecuación de grava gruesa)
b) Aproxime la solución aplicando algún método numérico.
Utilizando Gauss-Seidel
Se observa que por norma infinito y radio espectral > 1 el método no converge, lo que se puede
observar aplicando calculando las raíces por el método. Se ve que claramente en método se
escapa de un valor determinado.
Utilizando Jacobi:
Tenemos que por norma infinito y radio espectral estamos más cerca de 0 si comparamos al
método de Gauss-Siedel, pero sigue siendo > 1 por lo cual no asegura la convergencia.
Efectivamente, el método como nos indica por norma infinito y radio espectral no converge.
Ocurriendo algo similar con el método aplicado anteriormente.
Por lo tanto, ninguno del método que manejamos hasta el momento puede satisfacer la
convergencia del sistema propuesto.
c) Estime el error cometido en norma infinito.
El error en norma infinito es:
Gauss-Seidel = {158423} no converge
Jacobi = {15734} no converge
Luego, ambos métodos a la quinta iteración tienen errores del orden de las 4 y 5 cifras,
para Jacobi tenemos que el error es menor vs G-S.