El Método Iterativo Jacobi

El método iterativo Jacobi.
Es un método iterativo, usado para resolver sistemas de ecuaciones lineales del tipo. El
algoritmo toma su nombre del matemático alemán Carl Gustav Jakob Jacobi. El método de
Jacobi consiste en usar fórmulas como iteración de punto fijo.
Un método iterativo con el cual se resuelve el sistema lineal A x = b comienza con una
aproximación inicial x (0) a la solución x y genera una sucesión de vectores x (k) que
converge a x. Los métodos iterativos traen consigo un proceso que convierte el sistema A x
= b en otro equivalente de la forma x = T x + c para alguna matriz fija T y un vector c.
Luego de seleccionar el vector inicial x (0) la sucesión de los vectores de la solución

aproximada se genera calculando:
X (k) = Tx(k-1) + c
Para cada k = 1, 2,3,....
El método se escribe en la forma x (k) = T x (k-1) + c separando A en sus partes diagonal

D y fuera de la diagonal. Sea D la matriz diagonal cuya diagonal es la misma que A, sea -L
la parte estrictamente triangular inferior de la parte A y sea -U la parte estrictamente
triangular superior de A.
Con esta notación A = D-L-U, entonces transformamos la ecuación A x = b, o (D-L-U) x =

b, en
D x = (L+U) x + b
Y, si D-1 existe, es decir, si a i,i es distinto de cero para cada i, entonces
x = D-1(L+U) x + D-1b.
Esto da origen a la forma matricial del método iterativo de Jacobi:
X k) = D-1(L+U) x (k-1) + D-1b, k = 1,2,...
Al introducir la notación Tj = D-1(L+U) y c esta técnica tiene la forma
X (k) = T x (k-1) + c
Es de mencionar el siguiente teorema: " Si A es estrictamente diagonal dominante, entonces

con cualquier elección de la aproximación inicial, el método de Jacobi da una sucesión que
converge a la solución única de A x = b"
Sucesión
Se construye descomponiendo la matriz del sistema A en la forma siguiente: A= D+L+U
Donde
D es una matriz diagonal. Partiendo de Ax=b, podemos reescribir dicha ecuación como:
Dx(L+U) x=b
Luego
Si aii ≠0 para cada i. Por la regla iterativa, la definición del Método de Jacobi puede ser
expresado de la forma.
Donde k es el contador de iteración, finalmente tenemos:
Algoritmo
EJEMPLO
Con un vector inicial:
x1=0 x2=0 x3=0
Resolver por el método de Jacobi, el siguiente sistema de ecuaciones.
6x1 + 2x2 + x3 = 22 -x1 + 8x2 + 2x3 = 30 x1 - x2 + 6x3 = 23 Solución:
2. Para un vector inicial (0; 0; 0) hallar los valores de x1, x2, x3.
x1 = (22 - 2x2 - x3)/6 x2 = (30 + x 1 - 2x3)/8 x3 = (23 - x1 + x2)/6
2. Para un vector inicial (0; 0; 0) hallar los valores de x1, x2, x3.
x1 = (22 - 2(0) - (0)) /6 x2 = (30 + (0) - 2(0)) /8 x3 = (23 - (0) + (0)) /6
3. Teniendo para nuestra primera iteración los siguientes valores:
x1 = 3.66 x2 = 3.75 x3 = 3.83
4. Así sucesivamente sustituimos los valores de la primera iteración en las
ecuaciones para seguir a los valores de las variables.
A continuación se muestra un ejemplo de cómo se resuelve un ejercicio
Jacobi
3x0-x1-x2=1
-x0+3x1+x2=3
2x0+x1+4x2=7
Una vez que se tiene la ecuación se hacen las iteraciones hasta encontrar la raíz.
MÉTODO DE GAUSS SEIDEL
El Método de Gauss-Seidel consiste en hacer iteraciones, a partir de un vector inicial, para
encontrar los valores de las incógnitas hasta llegar a una tolerancia deseada, la diferencia
radica en que cada vez que se desee encontrar un nuevo valor de una xi, además de usar los
valores anteriores de las x, también utiliza valores actuales de las x encontradas antes
(desde x0 hasta xi-1). La ecuación es la siguiente:
El método de Gauss-Seidel surgio como una modificación del método de Jacobi que acelera
la convergencia de éste.
El método de Gauss-Seidel recorta sustancialmente el número de iteraciones a realizar para

obtener una cierta precisión en la solución. Evidentemente los criterios de convergencia son
similares a los de Jacobi.
Este criterio no solo se aplica a las ecuaciones lineales que se resuelven con el método de
Gauss-Seidel sino también para el método iterativo del punto fijo y el método de jacobi. Por
tanto, al aplicar este criterio sobre las ecuaciones de Gauss-Seidel y evaluando con respecto
a cada una de las incógnitas, obtenemos la expresión siguiente:
El valor absoluto de las pendientes en la ecuación, deben ser menor que la unidad para asegurar
la convergencia.
Es decir, el elemento diagonal debe ser mayor que el

elemento fuera de la diagonal para cada reglón de ecuaciones. La generalización del criterio
anterior para un sistema de n ecuaciones es:
El método de Gauss-Seidel está basado en el concepto de punto fijo, es decir ( xi = gi (x), i = 1..
n), para resolver sistemas de ecuaciones lineales.Para garantizar la convergencia se debe de
cumplir que el sistema tenga una diagonal dominante, es decir que se cumpla la desigualdad
siguiente, si se cambió el orden de las ecuaciones esta puede divergir.
Resolver el siguiente sistema de ecuación por el método Gauss-Seidel utilizando un
ξ= 0.001.
0.1 X1 + 7.0 X2 – 0.3 X3 = -19.30
3.0 X1 – 0.1 X2 – 0.2 X3 = 7.85
0.3 X1 – 0.2 X2 – 10.0 X3 = 71.40
SOLUCIÓN:
Primero ordenamos las ecuaciones, de modo que en la diagonal principal esten los
coeficientes mayores para asegurar la convergencia.
3.0 X1 – 0.1 X2 – 0.2 X3 = 7.85

0.1 X1 + 7.0 X2 – 0.3 X3 = -19.30
0.3 X1 – 0.2 X2 – 10.0 X3 = 71.40
Despejamos cada una de las variables sobre la diagonal:
Suponemos los valores
iniciales X2 = 0 y X3 = 0 y calculamos X1
Este valor junto con el de X3 se puede utilizar para obtener X2
La primera iteración se
completa sustituyendo los valores de X1 y X2 calculados obteniendo:
En la segunda iteración, se repite el mismo procedimiento:

Comparando los valores calculados entre la primera y la segunda iteración
Como podemos observar, no se cumple la condición
Entonces tomamos los valores

calculados en la última iteración y se toman como supuestos para la siguiente iteración. Se
repite entonces el proceso:
Comparando de nuevo los valores obtenidos

Como se observa todavía no se cumple la condición
Así que hacemos otra iteración
Comparando los valores obtenidos
Dado que se cumple la condición, el resultado es:
X1 = 3.0
X2 = -2.5
X3 = 7.0
Anuncios

El Método Iterativo Jacobi

Cargado por

Copyright:

Formatos disponibles

El Método Iterativo Jacobi

Cargado por

Información del documento

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

El Método Iterativo Jacobi

Cargado por

Copyright:

Formatos disponibles

El método iterativo Jacobi.

Luego de seleccionar el vector inicial x (0) la sucesión de los vectores de la solución

El método se escribe en la forma x (k) = T x (k-1) + c separando A en sus partes diagonal

Con esta notación A = D-L-U, entonces transformamos la ecuación A x = b, o (D-L-U) x =

Es de mencionar el siguiente teorema: " Si A es estrictamente diagonal dominante, entonces

Donde k es el contador de iteración, finalmente tenemos:

A continuación se muestra un ejemplo de cómo se resuelve un ejercicio

El método de Gauss-Seidel recorta sustancialmente el número de iteraciones a realizar para

Es decir, el elemento diagonal debe ser mayor que el

Resolver el siguiente sistema de ecuación por el método Gauss-Seidel utilizando un

3.0 X1 – 0.1 X2 – 0.2 X3 = 7.85

Este valor junto con el de X3 se puede utilizar para obtener X2

En la segunda iteración, se repite el mismo procedimiento:

Como podemos observar, no se cumple la condición

Entonces tomamos los valores

Comparando de nuevo los valores obtenidos

Así que hacemos otra iteración

Comparando los valores obtenidos

Dado que se cumple la condición, el resultado es:

También podría gustarte