Algebra Lineal para Estudiantes de Ingenie - Juan Carlos Del Valle Sotelo-500-1145-300-646 - 115 PDF

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

894 CAPÍTULO 8 Álgebra lineal numérica

donde T = −Q−1U y c = Q−1b. Ahora bien,

U = A−Q

y, en consecuencia,

Q−1U = Q−1 (A − Q)
= Q−1 A − In .

Luego,

−Q−1U = In − Q−1 A.

Ası́, la matriz que corresponde al formato general (8.33) (cfr. página 877) para la transformación afı́n
(8.34) en el método de Gauss-Seidel es

TQ = In − Q−1 A

y la transformación afı́n cuyos puntos fijos coinciden con las soluciones del sistema para este método es

gQ (x) = (In − Q−1 A)x +c

donde c = Q−1b. Entonces el método iterativo de Gauss-Seidel será convergente si


 
In − Q−1 A < 1 (8.53)

para alguna norma matricial natural. Esta condición es difı́cil de aplicar en la práctica, debido a la dificul-
tad que entraña calcular Q−1 . El siguiente teorema contiene condiciones suficientes para la convergencia
del método iterativo de Gauss-Seidel en términos de una caracterı́stica muy simple de identificar en la
matriz de coeficientes.

Teorema 8.7 Sean A una matriz cuadrada de orden n y b ∈ Rn . Si la matriz A es diagonalmente


dominante, entonces el método iterativo de Gauss-Seidel converge a la única solución del sistema
lineal Ax = b para cualquier vector inicial x0 ∈ Rn .

DEMOSTRACIÓN Q Sea Q la matriz (8.50), entonces el método iterativo de Gauss-Seidel en forma matricial está dado por
(8.52). Sea λ un valor propio de la matriz In −Q−1 A yx = (x1 , x2 , . . . xn ) un vector propio correspondiente
tal que x ∞ = 1. Entonces

(In − Q−1 A)x = λx

que equivale a

x − Q−1 Ax = λx;

(Q − A)x = λQx.

Claramente Q − A = −U, donde U es la matriz (8.51); entonces, para cada i = 1, 2, . . . , n,

Page (PS/TeX): 76 / 894, COMPOSITE

También podría gustarte