Método de Jacobi.
Este es un método alterno e
iterativo en la búsqueda de
soluciones de sistemas de
ecuaciones y utiliza el criterio de
Cauchy como criterio de paro.
Para derivar la fórmula del método,
hágase A = D + R; esto es, la matriz
del sistema se divide en dos, de las
cuales la matriz D es una matriz
diagonal y R, por el contrario, es tal
que sólo la diagonal principal tiene
elementos cero.
a11 a12 a13 ... a1n a11 0 0 ... 0 0 a12 a13 ... a1n
a a 22 a 23 ... a 2 n 0 a 22 0 ... 0 a 21 0 a 23 ... a 2 n
21
a31 a32 a33 ... a3n 0 0 a33 ... 0 a31 a32 0 ... a3n
... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
a m1 am2 am3 ... a mn 0 0 0 ... a mn a m1 am2 a m3 ... 0
Derivando la fórmula:
AX = B forma matricial del sistema
(D + R) X = B sustituyendo A
DX+RX = B quitando los paréntesis
DX = B–RX despejando
X = B D-1 – R X D-1 multiplicando por D-1
Así, para calcular el valor de X se requiere
conocer el valor de X. Esto no es una
contradicción, sino que de aquí se manifiesta la
naturaleza iterativa del método.
Hay que dar una primera aproximación a la
matriz X para acercarse al valor real según
una ε preestablecida. Así, la fórmula
general del método es:
Xi+1 = B D-1 – R D-1 Xi
Este método converge siempre que cada
elemento sobre la diagonal principal sea
mayor, en valor absoluto, que la suma de
los otros elementos de la fila, en valor
absoluto.
De la fórmula iterativa, se tiene:
Xi+1 = D-1 (B – R Xi)
Esto equivale a despejar del sistema las xi de
la i - ésima ecuación; en forma matemática:
1 n
xi
aii
bi
aij x j con j ≠ i
j 1
Este método requiere de un valor inicial para
X (por simplicidad se toma el vector cero) y
el proceso se detiene hasta que se cumpla
que |Xi+1 – Xi| < ε
NOTA
Es importante considerar para este método, que el
valor inicial para cada Xi de la i-ésima ecuación sea
el valor cero, ya que a partir de él se puede obtener
valores positivos o negativos; el proceso repetitivo
se detiene hasta que se cumpla con |Xi+1 – Xi| < ε
(Criterio de Cauchy ) para todos los casos.
Algoritmo Jacobi:
Leer n, ε
Para i = 1 hasta n fin_si
Para j = 1 hasta n+1 fin_para
Leer aij yi = yi/aij
fin_para delta = |xi – yi|
fin_para Si delta > ε entonces
Para i = 1 hasta n
fin = falso
fin_si
xi = 0
fin_para
fin_para
Para i = 1 hasta n
Repetir
xj = yj
fin = verdadero
fin_para
Para i = 1 hasta n
hasta fin
yi = ai,n+1 Para i = 1 hasta n
Para j = 1 hasta n Imprimir xi
Si i ≠ j entonces fin_para
yi = y i – Terminar
aij*xi
