Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 68 vistas

    Metodos Iterativos

    Cargado por

    LuisFloresRamirez
    Este documento describe los métodos iterativos de Jacobi y Gauss-Seidel para resolver sistemas de ecuaciones lineales. El método de Jacobi genera sucesivas aproximaciones a la solución mediante la iteración x(k)=Tx(k-1)+C. El método de Gauss-Seidel actualiza las aproximaciones de forma parcial, usando los valores más recientes disponibles en cada paso. Se presentan ejemplos numéricos ilustrativos para ambos métodos.

    Copyright:

    © All Rights Reserved
    Descargue como PPTX, PDF, TXT o lea en línea desde Scribd

    Metodos Iterativos

    Cargado por

    LuisFloresRamirez
    68 vistas16 páginas
    Este documento describe los métodos iterativos de Jacobi y Gauss-Seidel para resolver sistemas de ecuaciones lineales. El método de Jacobi genera sucesivas aproximaciones a la solución mediante la iteración x(k)=Tx(k-1)+C. El método de Gauss-Seidel actualiza las aproximaciones de forma parcial, usando los valores más recientes disponibles en cada paso. Se presentan ejemplos numéricos ilustrativos para ambos métodos.

    Título original

    METODOS ITERATIVOS

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PPTX, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento describe los métodos iterativos de Jacobi y Gauss-Seidel para resolver sistemas de ecuaciones lineales. El método de Jacobi genera sucesivas aproximaciones a la solución mediante la iteración x(k)=Tx(k-1)+C. El método de Gauss-Seidel actualiza las aproximaciones de forma parcial, usando los valores más recientes disponibles en cada paso. Se presentan ejemplos numéricos ilustrativos para ambos métodos.

    Copyright:

    © All Rights Reserved
    Descargue como PPTX, PDF, TXT o lea en línea desde Scribd
    Descargar como pptx, pdf o txt
    68 vistas16 páginas

    Metodos Iterativos

    Cargado por

    LuisFloresRamirez
    Este documento describe los métodos iterativos de Jacobi y Gauss-Seidel para resolver sistemas de ecuaciones lineales. El método de Jacobi genera sucesivas aproximaciones a la solución mediante la iteración x(k)=Tx(k-1)+C. El método de Gauss-Seidel actualiza las aproximaciones de forma parcial, usando los valores más recientes disponibles en cada paso. Se presentan ejemplos numéricos ilustrativos para ambos métodos.

    Copyright:

    © All Rights Reserved
    Descargue como PPTX, PDF, TXT o lea en línea desde Scribd
    Descargar como pptx, pdf o txt
    de 16

    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
    =, k0
    = = 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:
    =
    =
    =, k0

     
    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.

    También podría gustarte