Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 40 vistas

    Roger Unocc Sihuinta

    Cargado por

    Javier Zarate Neira
    Este documento presenta dos problemas relacionados con métodos iterativos para resolver sistemas de ecuaciones lineales. El Problema 11.10 muestra el código para implementar el método de Gauss-Seidel para resolver un sistema 3x3. El Problema 11.14 pide identificar qué conjuntos de ecuaciones no podrían resolverse con Gauss-Seidel y explicar por qué no convergerían.

    Copyright:

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

    Roger Unocc Sihuinta

    Cargado por

    Javier Zarate Neira
    40 vistas6 páginas
    Este documento presenta dos problemas relacionados con métodos iterativos para resolver sistemas de ecuaciones lineales. El Problema 11.10 muestra el código para implementar el método de Gauss-Seidel para resolver un sistema 3x3. El Problema 11.14 pide identificar qué conjuntos de ecuaciones no podrían resolverse con Gauss-Seidel y explicar por qué no convergerían.

    Descripción original:

    13190167

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

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

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento presenta dos problemas relacionados con métodos iterativos para resolver sistemas de ecuaciones lineales. El Problema 11.10 muestra el código para implementar el método de Gauss-Seidel para resolver un sistema 3x3. El Problema 11.14 pide identificar qué conjuntos de ecuaciones no podrían resolverse con Gauss-Seidel y explicar por qué no convergerían.

    Copyright:

    © All Rights Reserved
    Descargue como ODT, PDF, TXT o lea en línea desde Scribd
    Descargar como odt, pdf o txt
    40 vistas6 páginas

    Roger Unocc Sihuinta

    Cargado por

    Javier Zarate Neira
    Este documento presenta dos problemas relacionados con métodos iterativos para resolver sistemas de ecuaciones lineales. El Problema 11.10 muestra el código para implementar el método de Gauss-Seidel para resolver un sistema 3x3. El Problema 11.14 pide identificar qué conjuntos de ecuaciones no podrían resolverse con Gauss-Seidel y explicar por qué no convergerían.

    Copyright:

    © All Rights Reserved
    Descargue como ODT, PDF, TXT o lea en línea desde Scribd
    Descargar como odt, pdf o txt
    de 6

    Problema 11.

    10:
    %function x=GaussSeidel(A,b)
    %% Datos
    A = [10 2 -1;-3 -6 2;1 1 5];
    b = [27 -61.5 -21.5];
    n = 3;
    % Resultados
    % x = vector solucion
    [n n] = size(A);
    x = zeros(n,1);
    y = zeros(n,1);
    error = 0.05;
    flag = 1;
    NTOL = 50;
    k = 0;
    fprintf('%5d',k);
    for m = 1:n
    fprintf('%10.5f',x(m));
    end

    % prueba de diagonalizacion
    i = 1;
    while (i <= n) & (flag == 1)
    suma = 0;
    for j = 1:n
    if i ~= j

    suma = suma + abs(A(i,j));


    end
    end
    if abs(A(i,i)) <= suma
    fprintf('\nError de diagonalizacion');
    flag = 0;
    end
    i=i+1;
    end % fin de la prueba
    if flag == 0
    break
    end
    while 1
    flag = 1;
    for i = 1:n
    suma = 0;
    for j = 1:n
    if i ~= j
    suma = suma+A(i,j)*x(j)/A(i,i);
    end
    end
    y(i) = b(i)/A(i,i)-suma;
    if abs(y(i) - x(i)) > error
    flag = 0;
    end
    x(i) = y(i);
    end
    k=k+1;

    fprintf('\n%5d',k);
    for i=1:n
    fprintf('%10.5f',x(i));
    end
    if (NTOL == k) | (flag == 1)
    break
    end
    end

    los resultados son:

    0 0.00000 0.00000 0.00000


    1 2.70000 8.90000 -6.62000
    2 0.25800 7.91433 -5.93447
    3 0.52369 8.01000 -6.00674
    4 0.49733 7.99909 -5.99928

    Problema 11.14:
    De los tres conjuntos siguientes de ecuaciones
    lineales,

    identifique aquel(los) que no podra resolver con


    el uso de un
    mtodo iterativo tal como el de Gauss-Seidel.
    Demuestre que su
    solucin no converge, utilizando cualquier
    nmero de iteraciones
    que sea necesario. Enuncie con claridad su
    criterio de convergencia (es decir, cmo se sabe
    que no est convergiendo).
    Conjunto uno

    Conjunto dos

    9x+ 3y+ z= 13
    5z= 6
    6x+ 8z= 2
    3

    Conjunto tres

    x+ y+ 6z= 8
    x+ 5y z= 5

    2x+ 5y z= 6

    4x+ 2y 2z= 4

    3x+ 4y+
    2x+ 2y 3z=
    2y z = 1

    Al ejecutar el primer conjunto en el matlab se debe


    recordar que debe ser diagonal positiva.
    Por lo tanto para que converga la ecuacin se debe
    cambiar la fila 2 por la fila 3.
    9x+ 3y+ z= 13
    2x+ 5y z= 6
    6x+ 8z= 2

    La solucin es:
    0

    0.00000

    0.00000

    0.00000

    1.44444

    0.62222

    1.33333

    1.08889

    1.03111

    1.06667

    0.98222

    1.02044

    0.98667

    0.99467

    0.99947

    0.99600

    Para hallar los conjuntos 2 y 3 por el mtodo de


    Gauss-Seidel no se puede debido a que no
    converge ya que la suma de las diagonales
    principales debe ser mayor a la suma de sus
    componentes de la misma fila.

    Conjunto dos
    x+ y+ 6z= 8
    x+ 5y z= 5
    4x+ 2y 2z= 4

    Conjunto tres
    3x+ 4y+ 5z= 6
    2x+ 2y 3z= 3
    2y z = 1

    conjunto 2:
    4 debe ser mayor que (2+2) para que sea convergente.
    En el conjunto 3 :
    La fila 1 debe haber un numero mayor que las demas
    de la fila para que converga.
    La fila 2 tambien sucede lo mismo.
    Por lo tanto el conjunto 2 y 3 no convergen .

    Problema11.15:
    Emplee la librera o paquete de software de su
    preferencia
    para obtener una solucin, calcular la inversa y
    determinar el
    nmero de condicin (sin dar escala) con base en la
    norma de
    suma de renglones, para los sistemas
    a)

    En ambos casos, las respuestas para todas las xdeben


    ser 1.

    También podría gustarte