Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Introducción Método de Jacobi

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 5

METODOS NUMERICOS

Solución de sistemas de ecuaciones lineales

Métodos Directos

Son aquellos que no involucran procesos infinitos y que, en consecuencia, serían exactos, si no
existieran errores de redondeo.

Métodos Indirectos o Iterativos

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.

Ing. Rafael González Betancourt 29-sep-2021


METODOS NUMERICOS

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

17𝑥 −2𝑦 −3𝑧 = 500


−5𝑥 −21𝑦 −2𝑧 = 200
−5𝑥 −5𝑦 22𝑧 = 30

En este caso la matriz es dominante

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

500 + 2𝑦0 + 3𝑧0 200 + 5𝑥0 + 2𝑧0 30 + 5𝑥0 + 5𝑦0


𝑥1 = 𝑦1 = 𝑧1 =
17 −21 22

500 + 2(0) + 3(0)


𝑥1 = = 29.4117
17

200 + 5(0) + 2(0)


𝑦1 = = −9.5238
−21

30 + 5(0) + 5(0)
𝑧1 = = 1.3636
22

Ing. Rafael González Betancourt 29-sep-2021


METODOS NUMERICOS

2da iteración 𝑥1 = 29.4117, 𝑦1 = −9.5238, 𝑧1 = 1.3636

500 + 2𝑦1 + 3𝑧1 200 + 5𝑥1 + 2𝑧1 30 + 5𝑥1 + 5𝑦1


𝑥2 = 𝑦2 = 𝑧2 =
17 −21 22
500 + 2(−9.5238) + 3(1.3636)
𝑥2 = = 28.5319
17
200 + 5(29.4117) + 2(1.3636)
𝑦2 = = −16.6564
−21
30 + 5(29.4117) + 5(−9.5238)
𝑧2 = = 5.8836
22

3ra iteración 𝑥2 = 28.5319, 𝑦2 = −16.6564, 𝑧2 = 5.8836

500 + 2𝑦2 + 3𝑧2 200 + 5𝑥2 + 2𝑧2 30 + 5𝑥2 + 5𝑦2


𝑥3 = 𝑦3 = 𝑧3 =
17 −21 22
500 + 2(−16.6564) + 3(5.8836)
𝑥3 = = 28.4904
17
200 + 5(28.5319) + 2(5.8836)
𝑦3 = = −16.8774
−21
30 + 5(28.5319) + 5(−16.6564)
𝑧3 = = 4.0626
22

4ta iteración 𝑥3 = 28.4904, 𝑦3 = −16.8774, 𝑧3 = 4.0626

500 + 2𝑦3 + 3𝑧3 200 + 5𝑥3 + 2𝑧3 30 + 5𝑥3 + 5𝑦3


𝑥4 = 𝑦4 = 𝑧4 =
17 −21 22
500 + 2(−16.8774) + 3(4.0626)
𝑥4 = = 28.1431
17
200 + 5(28.4904) + 2(4.0626)
𝑦4 = = −16.6941
−21
30 + 5(28.4904) + 5(−16.8774)
𝑧4 = = 4.0029
22

Ing. Rafael González Betancourt 29-sep-2021


METODOS NUMERICOS

5ta iteración 𝑥4 = 28.1431, 𝑦4 = −16.6941, 𝑧4 = 4.0029

500 + 2𝑦4 + 3𝑧4 200 + 5𝑥4 + 2𝑧4 30 + 5𝑥4 + 5𝑦4


𝑥5 = 𝑦5 = 𝑧5 =
17 −21 22
500 + 2(−16.6941) + 3(4.0029)
𝑥5 = = 28.1541
17
200 + 5(28.1431) + 2(4.0029)
𝑦5 = = −16.6057
−21
30 + 5(28.1431) + 5(−16.6941)
𝑧5 = = 3.9656
22

6ta iteración 𝑥5 = 28.1541, 𝑦5 = −16.6057, 𝑧5 = 3.9656

500 + 2𝑦5 + 3𝑧5 200 + 5𝑥5 + 2𝑧5 30 + 5𝑥5 + 5𝑦5


𝑥6 = 𝑦6 = 𝑧6 =
17 −21 22
500 + 2(−16.6057) + 3(3.9656)
𝑥6 = = 28.1579
17
200 + 5(28.1541) + 2(3.9656)
𝑦6 = = −16.6048
−21
30 + 5(28.1541) + 5(−16.6057)
𝑧6 = = 3.9882
22

7a iteración 𝑥6 = 28.1579, 𝑦6 = −16.6048, 𝑧6 = 3.9882

500 + 2𝑦6 + 3𝑧6 200 + 5𝑥6 + 2𝑧6 30 + 5𝑥6 + 5𝑦6


𝑥7 = 𝑦7 = 𝑧7 =
17 −21 22
500 + 2(−16.6048) + 3(3.9882)
𝑥7 = = 28.1620
17
200 + 5(28.1579) + 2(3.9882)
𝑦7 = = −16.6079
−21
30 + 5(28.1579) + 5(−16.6048)
𝑧7 = = 3.9893
22

Ing. Rafael González Betancourt 29-sep-2021


METODOS NUMERICOS

Calculemos el error hasta la 7a iteración

𝑥6 = 28.1579, 𝑦6 = −16.6048 𝑧6 = 3.9882

𝑥7 = 28.1620, 𝑦7 = −16.6079 𝑧7 = 3.9893

𝑥𝑛𝑢𝑒𝑣𝑎 − 𝑥𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟
𝐸𝑥 = | | ∗ 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.

𝑥7 = 28.1620, 𝑦7 = −16.6079 𝑧7 = 3.9893

Ejercicio para practicar

Resuelva el siguiente sistema de ecuaciones por el método de Jacobi.

Valores iniciales: 𝑥0 = 0, 𝑦0 = 0, 𝑧0 = 0

9𝑥 +2𝑦 −𝑧 = −2
7𝑥 +8𝑦 +5𝑧 = 3
3𝑥 +4𝑦 −10𝑧 = 6

Ing. Rafael González Betancourt 29-sep-2021

También podría gustarte