Método de jacobi Métodos Númericos

25
Clase 10 12-Marzo-2015

Transcript of Método de jacobi Métodos Númericos

Page 1: Método de jacobi Métodos Númericos

Clase 10

12-Marzo-2015

Page 2: Método de jacobi Métodos Númericos

Esta técnica muestra cierta similitud con el método de iteración de punto fijo, ya

que consiste en despejar una de las incógnitas de una ecuación dejándola en

función de las otras. La manera mas sencilla es despejar 𝑥1 de la primera

ecuación; 𝑥2 de la segunda ecuación; 𝑥𝑖 de la i-esima ecuación, hasta 𝑥𝑛 de la n-

esima ecuación. Es necesario, por razones obvias que todos los elementos de la

diagonal principal de la matriz de coeficientes del sistema lineal, sean diferentes

de cero.

Page 3: Método de jacobi Métodos Númericos

Sea el sistema lineal:

11 1 12 2 13 3 1n n 1

21 1 22 2 23 3 2n n 2

31 1 32 2 33 3 3n n 3

n1 1 n2 2 n3 3 nm n n

a x a x a x ... a x C

a x a x a x ... a x C

a x a x a x ... a x C

.

.

.

a x a x a x ... a x C

Page 4: Método de jacobi Métodos Númericos

Al realizar los despejes

propuestos se obtiene 𝑥1 de la

primera ecuación, 𝑥2 de la

segunda ecuación, etc., se

obtiene:

131 12 1n1 2 3 n

11 11 11 11

232 21 2n2 1 3 n

22 22 22 22

3 31 32 3n3 1 2 n

33 33 33 33

n 1,nn n1 n2n 1 2 n 1

mn mn mn mn

aC a ax x x ... x

a a a a

aC a ax x x ... x

a a a a

C a a ax x x ... x

a a a a

.

.

.

aC a ax x x ... x

a a a a

Page 5: Método de jacobi Métodos Númericos

Para estimar la primera aproximación a la solución se debe partir de un vector

inicial, el cual puede ser un vector 𝒙𝟎 = 𝟎, o algún otro que se encuentre próximo

al vector solución 𝒙.

Page 6: Método de jacobi Métodos Númericos

Resolver el sistema lineal por medio del método de Jacobi. Emplear el vector

inicial 𝑥0 = 0

1 2 3 4

1 2 3 4

1 2 3 4

1 2 3 4

6x x x 4x 17

x 10x 2x x 17

3x 2x 8x x 19

x x x 5x 14

Page 7: Método de jacobi Métodos Númericos

Al despejar las incógnitas correspondientes al esquema se tiene

1 2 3 4

2 1 3 4

3 1 2 4

4 1 2 3

x 17 x x x / 6

x 17 x 2x x / 10

x 19 3x 2x x / 8

x 14 x x x / 5

Page 8: Método de jacobi Métodos Númericos

Si se inicia el proceso iterativo con el vector cero se obtiene:

1

1

1

2

1

3

1

4

x 2.833333

x 1.7

x 2.375

x 2.8

Page 9: Método de jacobi Métodos Númericos

Los resultados del vector 𝑥 1 se utilizan para estimar el vector 𝑥(2), los del vector

𝑥 3 y así sucesivamente. Los resultados del proceso iterativo se muestran en la

tabla 1

Tabla 1. Resultado

de las iteraciones

Page 10: Método de jacobi Métodos Númericos

En general, el vector aproximación a la solución después de las iteraciones se

puede calcular de la siguiente manera:

Page 11: Método de jacobi Métodos Númericos

k 1 k k k131 12 1n1 2 3 n

11 11 11 11

k 1 k k k232 21 2n2 1 3 n

22 22 22 22

k 1 k k k3 31 32 3n3 1 2 n

33 33 33 33

k 1 k k kn 1,nn n1 n2n 1 2 n 1

mn mn mn mn

aC a ax x x ... x

a a a a

aC a ax x x ... x

a a a a

C a a ax x x ... x

a a a a

.

.

.

aC a ax x x ... x

a a a a

Page 12: Método de jacobi Métodos Númericos

O bien escrito en forma compacta:

𝑥𝑖𝑘+1 =

1

𝑎𝑖𝑖𝐶𝑖 − 𝑗=1

𝑗≠1

𝑛 𝑎𝑖𝑗𝑥𝑗𝑘 𝑝𝑎𝑟𝑎 1 ≤ 𝑖 ≤ 𝑛……………………(A)

Page 13: Método de jacobi Métodos Númericos

1. Utilizar la herramienta de Excel para generar la tabla de la figura 1 que

contiene a la matriz aumentada, el vector inicial y la programación de los

despejes que se generen al utilizar la ecuación (A)

Page 14: Método de jacobi Métodos Númericos

Para el método de Jacobi, considere un sistema 𝐴𝑥 = 𝑏

Sea 𝐴 = 𝐷 − 𝐸 − 𝐹, donde 𝐷 es la diagonal de 𝐴, −𝐸 la triangular inferior y −𝐹 la triangular

superior.

Así, la sucesión que se construye con este método iterativo será:

𝐴𝑥 = 𝑏

𝐷 − 𝐸 − 𝐹 𝑥 = 𝑏

𝐷𝑥 = 𝐸 + 𝐹 𝑥 + 𝑏

𝑥𝑘 = 𝐷−1 𝐸 + 𝐹 𝑥𝑘−1 + 𝐷−1𝑏

Page 15: Método de jacobi Métodos Númericos

El siguiente programa resuelve mediante el método de Jacobi un sistema de

ecuaciones 𝐴𝑥 = 𝑏 con un error menor que una tolerancia dada tol.

Note que el programa necesita un dato inicial 𝑥0.

Además, el programa se detiene si se alcanza un número máximo de iteraciones

maxit sin que se satisfaga el criterio de convergencia.

Page 16: Método de jacobi Métodos Númericos
Page 17: Método de jacobi Métodos Númericos

Ejemplo:

Resuelva 𝐴𝑥(0) = 𝑏, con una aproximación inicial 𝑥(0) = [0 0 0]′

Page 18: Método de jacobi Métodos Númericos
Page 19: Método de jacobi Métodos Númericos
Page 20: Método de jacobi Métodos Númericos
Page 21: Método de jacobi Métodos Númericos
Page 22: Método de jacobi Métodos Númericos
Page 23: Método de jacobi Métodos Númericos
Page 24: Método de jacobi Métodos Númericos
Page 25: Método de jacobi Métodos Númericos