Un m�todo iterativo con el cual se resuleve el sistema lineal Ax = b comienza con una aproximaci�n inicial x⁽⁰⁾a la solucion 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 Ax = b en otro equivalente de la forma x = Tx + c para alguna matriz fija T y un vector c.

    Luego de seleccionar el vector inicial x⁽⁰⁾ 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) = Tx^(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 notacion A = D-L-U, entonces transformamos la ecuaci�n Ax = b, o (D-L-U)x = b, en

Dx = (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 T_j = D^-1(L+U) y c_j, esta t�cnica tiene la forma

x^(k) = Tx^(k-1) + c

Seudoc�digo

    Es de mencionar el siguiente teorema: " Si A es estrictamente diagonal dominante, entonces con cualquier eleccion de la aproximaci�n inicial, el m�todo de Jacobi da una sucesion que converge a la soluci�n �nica de Ax = b"

ANTERIOR   SIGUIENTE