Introducción Método de Jacobi
Introducción Método de Jacobi
Introducción Método de Jacobi
Métodos Directos
Son aquellos que no involucran procesos infinitos y que, en consecuencia, serían exactos, si no
existieran errores de redondeo.
Parten de una aproximación inicial y que, por aplicación de un determinado algoritmo, conducen a
aproximaciones sucesivamente mejores en cada iteración.
Método de Jacobi
En análisis numérico el método de Jacobi es un método iterativo, usado para resolver sistemas de
ecuaciones lineales. El método de Jacobi consiste en usar fórmulas como iteración de punto fijo,
fórmulas donde se tiene despejada cada una de las variables que conforman el sistema de
ecuaciones.
𝑥0 = 0, 𝑦0 = 0, 𝑧0 = 0
𝑎1 𝑥 𝑏1 𝑦 𝑐1 𝑧 = 𝐷1
𝑎2 𝑥 𝑏2 𝑦 𝑐2 𝑧 = 𝐷2
𝑎3 𝑥 𝑏3 𝑦 𝑐3 𝑧 = 𝐷3
𝐷1 + 𝑏1 𝑦0 + 𝑐1 𝑧0 𝐷2 + 𝑎2 𝑥0 + 𝑐2 𝑧0 𝐷3 + 𝑎3 𝑥0 + 𝑏3 𝑦0
𝑥1 = 𝑦1 = 𝑧1 =
𝑎1 𝑏2 𝑐3
El método de Jacobi siempre converge si la matriz es dominante y puede converger incluso si esta
condición no se cumple, esto debido a que la iteración es de punto fijo.
El método de Jacobi, cuentan con los siguientes requisitos para poder resolverse:
La matriz no debe tener elementos nulos. (Si una matriz cuadrada tiene dos filas (o dos
columnas) iguales, su determinante es nulo)
La matriz puede o no ser dominante, aunque es mejor trabajar con la matriz dominante.
El sistema de ecuaciones debe de ser cuadrado, es decir, el número de ecuaciones es igual
al número de incógnitas.
Ejemplo:
Resolver el siguiente sistema de tres ecuaciones por el método de Jacobi, tome en cuenta los
siguientes valores iniciales: 𝑥0 = 0, 𝑦0 = 0, 𝑧0 = 0
Despejamos las variables x, y z de la 1ra, 2da y 3er ecuaciones respectivamente y para comenzar a
iterar sustituimos los valores iniciales.
500 + 2𝑦 + 3𝑧 200 + 5𝑥 + 2𝑧 30 + 5𝑥 + 5𝑦
𝑥= 𝑦= 𝑧=
17 −21 22
1ra iteración 𝑥0 = 0, 𝑦0 = 0, 𝑧0 = 0
30 + 5(0) + 5(0)
𝑧1 = = 1.3636
22
𝑥𝑛𝑢𝑒𝑣𝑎 − 𝑥𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟
𝐸𝑥 = | | ∗ 100
𝑥𝑛𝑢𝑒𝑣𝑎
28.1620 − 28.1579
𝐸𝑥 = | | ∗ 100 = 0.01455%
28.1620
𝑦𝑛𝑢𝑒𝑣𝑎 − 𝑦𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟
𝐸𝑦 = | | ∗ 100
𝑦𝑛𝑢𝑒𝑣𝑎
−16.6079 − (−16.6048)
𝐸𝑦 = | | ∗ 100 = 0.0186%
−16.6079
𝑧𝑛𝑢𝑒𝑣𝑎 − 𝑧𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟
𝐸𝑧 = | | ∗ 100
𝑧𝑛𝑢𝑒𝑣𝑎
3.9893 − 3.9882
𝐸𝑧 = | | ∗ 100 = 0.0275%
3.9893
El cálculo del error nos indica que los valores obtenidos en la 7a iteración son una muy buena
aproximación a la solución del sistema.
Valores iniciales: 𝑥0 = 0, 𝑦0 = 0, 𝑧0 = 0
9𝑥 +2𝑦 −𝑧 = −2
7𝑥 +8𝑦 +5𝑧 = 3
3𝑥 +4𝑦 −10𝑧 = 6