martes, 20 de julio de 2010

ITERATIVE METHODS FOR SOLUTION OF SYSTEMS OF LINEAR EQUATIONS


ITERATIVE METHODS FOR SOLUTION OF SYSTEMS OF LINEAR EQUATIONS


In general, these methods are based on fixed point method and the process is iterated so substitute in a formula.

Iterative methods compared with direct, we do not guarantee a better approach, however, are more efficient when working with large matrices.

In computational mathematics, is an iterative method to solve a problem (as an equation or a system of equations) by successive approximations to the solution, starting from an initial estimate. This approach contrasts with the direct methods, which attempt to solve the problem once (like solving a system of equations Ax = b by finding the inverse of the matrix A). Iterative methods are useful for solving problems involving a large number of variables (sometimes in the millions), where direct methods would be prohibitively expensive even with better computer power available.


  • JACOBI METHOD

The basis of the method is to construct a convergent sequence defined iteratively. The limit of this sequence is precisely the solution of the system. For practical purposes if the algorithm stops after a finite number of steps leads to an approximation of the value of x in the solution of the system.


  • GAUSS SEIDEL METHOD


The methods of Gauss and Cholesky methods are part of direct or finite. After a number of operations nito, in the absence of errors rounding, we get x solution of the system Ax = b. The Gauss-Seidel method is part of the so-called indirect methods or iterative. They start with x0 = (x01, X02,:::; x0 n), an approximation initial solution. Since x0 is building a new approach of the solution, x1 = (x11, x12;:::; x1n).

From built x1 x2 (here the superscript indicates the iteration and does not indicate a power). So on construye fxkg a sequence of vectors, with the aim, not always guaranteed to quelimk! 1xk = x:
Generally, indirect thods are a good option when the matrix is very large and dispersed or sparse (sparse), ie when the number of onzero elements is small compared to n2, total number of elements A.

In these cases you must use an appropriate data structure lets you store only the nonzero elements. In each iteration of the Gauss-Seidel method, there are n subiteraciones. In ca first subiteracion be amended only x1. The other coordinates x2, x3, ..., xn are not modified can.



fuente

2 comentarios:

  1. I like your blog, it's very amazing,, but i can not find where is the gatget to add follower

    ResponderEliminar
  2. ah!! congratulations , you has done a gret job!! =)

    ResponderEliminar