Metodo de jacobi
Páginas: 2 (366 palabras)
Publicado: 12 de enero de 2010
Método de Jacobi
En la iteración de Jacobi, se escoge una matriz Q que es diagonal y cuyos elementos diagonales son los mismos que los de la matriz A. La matriz Q toma la forma:
y la ecuacióngeneral (63) se puede escribir como
Qx(k) = (Q-A)x(k-1) + b | (65) |
Si denominamos R a la matriz A-Q:
la ecuación (65) se puede reescribir como:
Qx(k) = -Rx(k-1) + b
El producto de la matrizQ por el vector columna x(k) será un vector columna. De modo análogo, el producto de la matriz R por el vector columna x(k-1) será también un vector columna. La expresión anterior, que es una ecuaciónvectorial, se puede expresar por necuaciones escalares (una para cada componente del vector). De este modo, podemos escribir, para un elemento i cualquiera y teniendo en cuenta que se trata de unproducto matriz-vector:
Si tenemos en cuenta que en la matriz Q todos los elementos fuera de la diagonal son cero, en el primer miembro el único término no nulo del sumatorio es el que contiene elelemento diagonal qii, que es precisamente aii. Más aún, los elementos de la diagonal de Rson cero, por lo que podemos eliminar el término i=j en el sumatorio del segundo miembro. De acuerdo con lo dicho,la expresión anterior se puede reescribir como:
de donde despejando xi(k) obtenemos:
que es la expresión que nos proporciona las nuevas componentes del vector x(k) en función de vector anteriorx(k-1) en la iteración de Jacobi. En la figura (14) se presenta un algoritmo para el método de Jacobi.
Figure: Implementación del método de Jacobi. |
|
El método de Jacobi se basa en escribirel sistema de ecuaciones en la forma:
| (66) |
Partimos de una aproximación inicial para las soluciones al sistema de ecuaciones y sustituimos estos valores en la ecuación (66). De esta forma,se genera una nueva aproximación a la solución del sistema, que en determinadas condiciones, es mejor que la aproximación inicial. Esta nueva aproximación se puede sustituir de nuevo en la parte...
En la iteración de Jacobi, se escoge una matriz Q que es diagonal y cuyos elementos diagonales son los mismos que los de la matriz A. La matriz Q toma la forma:
y la ecuacióngeneral (63) se puede escribir como
Qx(k) = (Q-A)x(k-1) + b | (65) |
Si denominamos R a la matriz A-Q:
la ecuación (65) se puede reescribir como:
Qx(k) = -Rx(k-1) + b
El producto de la matrizQ por el vector columna x(k) será un vector columna. De modo análogo, el producto de la matriz R por el vector columna x(k-1) será también un vector columna. La expresión anterior, que es una ecuaciónvectorial, se puede expresar por necuaciones escalares (una para cada componente del vector). De este modo, podemos escribir, para un elemento i cualquiera y teniendo en cuenta que se trata de unproducto matriz-vector:
Si tenemos en cuenta que en la matriz Q todos los elementos fuera de la diagonal son cero, en el primer miembro el único término no nulo del sumatorio es el que contiene elelemento diagonal qii, que es precisamente aii. Más aún, los elementos de la diagonal de Rson cero, por lo que podemos eliminar el término i=j en el sumatorio del segundo miembro. De acuerdo con lo dicho,la expresión anterior se puede reescribir como:
de donde despejando xi(k) obtenemos:
que es la expresión que nos proporciona las nuevas componentes del vector x(k) en función de vector anteriorx(k-1) en la iteración de Jacobi. En la figura (14) se presenta un algoritmo para el método de Jacobi.
Figure: Implementación del método de Jacobi. |
|
El método de Jacobi se basa en escribirel sistema de ecuaciones en la forma:
| (66) |
Partimos de una aproximación inicial para las soluciones al sistema de ecuaciones y sustituimos estos valores en la ecuación (66). De esta forma,se genera una nueva aproximación a la solución del sistema, que en determinadas condiciones, es mejor que la aproximación inicial. Esta nueva aproximación se puede sustituir de nuevo en la parte...
Leer documento completo
Regístrate para leer el documento completo.