MÉTODOS ITERATIVOS PARA RESOLVER
SISTEMAS LINEALES
Un método iterativo para resolver el sistema lineal AX = b comienza con una aproximación
inicial x(0) a la solución x, y se genera una sucesión de vectores {x(k)} que cuando converge
lo hace al vector solución x.
Los métodos iterativos convierten el sistema Ax = b en otro sistma equivalente dela forma x
= Tx + C, para alguna matriz fija T y C un vector.
Se selecciona el vector inicial x(0),la sucesión de vectores se genera así:
x(k) = Tx(k-1) + C, k = 1, 2, 3, …
Y se itera hasta que se satisfacen alguno de los siguientes criterios de paro:
i) El proceso se continúa hasta que se cumpla alguno de los siguientes criterios:
a) < TOL, TOL = 10-s, s IN
b) < TOL, TOL = 10-s, s IN
ii) Controlar el número máximo de iteraciones.
Se usará la norma euclidiana (norma l2) y la norma l, las cuales se definen de la
siguiente manera:
Definición:
Sea x = (x1, x2, …, xn) y las normas se definen:
L2 = = L = =
Método de Jacobi
Sea el sistema de ecuaciones lineales Ax = b:
E1: a11x1 + a12x2 + … + a1nxn = b1
E1: a21x1 + a22x2 + … + a2nxn = b2
……………………………………………………..
Ei: ai1x1 + ai2x2 + … + ainxn = bi
………………………………………………………
En: an1x1 + an2x2 + … + annxn = bn
De la ecuación Ei se despeja xi:
xi = , aii 0, i = 1, …, n
Y para generar a partir de , la iteración es:
=, i = 1, …, n
Ejemplo:
Sea el sistema :
E1: 3x1 – 0.1x2 – 0.2x3 = 7.85
E2: 0.1x1 + 7x2 – 0.3x3 = -19.3
E3: 0.3x1 – 0.2x2 + 10x3 = 71.4
a) Hallar la matriz T y el vector C
b) Resuelva el sistema usando el método de Jacobi, tomando como aproximación inicial
x(0) = (0,0,0)t, e itere hasta que se cumpla que < 10-4.
Solución:
a) De E1: despejamos x1 =
De E2: despejamos x2 =
De E3: despejamos x3 =
Luego:
T = y C=
b) Esquema iterativo de Jacobi:
Hallando la primera aproximación (k=0):
=
= = 2.616667
=
= = -2,757143
=, k0
= = 7.140000
Así
x(1) = (2.616667, -2,757143, 7.140000)t
y = 1 < 10-4, no se cumple
Hallando x(2), (k=1):
= = = 3.000762
= = = -2,488524
= = = 7.006357
Así x(2) = (3.000762, -2, 488524, 7. 006357)t
y = = 0.05420929 < 10-4, no se cumple
Hallando x(3), (k=2):
= = = 3.000806
= = = -2,499738
= = = 7.000207
Así x(3) = (3.000806, -2, 499738, 7. 000207)t
y = = 0.001602 < 10-4, no se cumple
Hallando x(4), (k=3):
= = = 3.000022
= = = -2,500003
= = = 6.999981
Así x(4) = (3.000806, -2, 499738, 7. 000207)t
y = = 0.000112 < 10-4, no se cumple
Hallando x(5), (k=4):
= = = 2.999999
= = = -2,500001
= = = 6.999999
Así x(5) = (2.999999, -2,500001, 6.999999)t
y = = 0.000112 < 10-4, se cumple
x = (2.999999, -2,500001, 6.999999)t
Método de Gauss-Seidel
Sea el sistema de ecuaciones lineales Ax = b:
E1: a11x1 + a12x2 + … + a1nxn = b1
E1: a21x1 + a22x2 + … + a2nxn = b2
……………………………………………………..
Ei: ai1x1 + ai2x2 + … + ainxn = bi
………………………………………………………
En: an1x1 + an2x2 + … + annxn = bn
De la ecuación Ei se despeja xi:
xi = , aii 0, i = 1, …, n
Y para generar a partir de , la iteración es:
=, i = 1, …, n
Luego la aproximación es: x(k) = (x1(k), x2(k),…., xk(k), xk+1(k-1),…. xn(k-1),)
Ejemplo:
Sea el sistema :
E1: 3x1 – 0.1x2 – 0.2x3 = 7.85
E2: 0.1x1 + 7x2 – 0.3x3 = -19.3
E3: 0.3x1 – 0.2x2 + 10x3 = 71.4
Resuelva el sistema usando el método de Gauss - Seidel, tomando como aproximación
inicial x(0) = (0,0,0)t, e itere hasta que se cumpla que < 10-4.
Solución:
a) De E1: despejamos x1 =
De E2: despejamos x2 =
De E3: despejamos x3 =
Luego:
Esquema iterativo de Gauss - Seidel:
=
=
=, k0
Hallando la primera aproximación (k=1):
= = 2.616667
= = = -2,794524
= = = = 7.005610
Así x(1) = (2.616667, -2,794524, 7.005610)t y = 1 < 10-4, no se cumple
Hallando x(2), (k=2):
= = = 2.990557
= = = -2,499625
= = = = 7.000291
Así x(2) = (2.990557, -2,499625, 7.0000291)t y = < 10-4, no se cumple
Hallando x(3), (k=3):
= = = 3.000032
= = = -2,499988
= = = = 6.999999
Así x(3) = (3.000032, -2,499988, 6.999999)t y = = 0.00135 < 10-4, no se cumple
Hallando x(4), (k=4):
= = = 3.000000
= = = -2,500000
= = = = 7.000000
Así x(4) = (3.000000, -2,500000, 7.000000)t y = = 0.000005 < 10-4, se cumple
Ejemplo:
Resuelva el sistema :
3.8x1 + 1.6 x2 + 0.9x3 = 3.72
-0.7x1 + 5.4x2 + 1.6x3 = 3.16
1.5x1 + 1.1x2 – 3.2x3 = 43.78
Resuelva el sistema usando el método de Jacobi y el de Gauss - Seidel, tomando como
aproximación inicial x(0) = (0,0,0)t, Realice 5 iteraciones redondeando los cálculos hasta
cinco cifras significativos.
