Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 6 vistas

    Tarea 2 2do Parcial

    Cargado por

    mafreire8

    Copyright:

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

    Tarea 2 2do Parcial

    Cargado por

    mafreire8
    6 vistas5 páginas

    Título original

    Tarea_2_2do_parcial

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    6 vistas5 páginas

    Tarea 2 2do Parcial

    Cargado por

    mafreire8

    Copyright:

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

    UNIVERSIDAD DE LAS FUERZAS ARMADAS - ESPE

    ÁREA DE ANÁLISIS FUNCIONAL


    MÉTODOS NUMÉRICOS

    Nombre: Marcelo Agustı́n Freire Poveda NRC:14020


    Fecha: 8 de febrero de 2024

    TAREA 2-P2

    Instrucciones
    Construir un script en el cual se pueda realizar la resolución de un sistema de ecuaciones de Mediante el
    método de Descomposición LU ingresando una matriz nxn

    Desarollo del script


    (1) Definimos las va variables a ser ingresadas, las cuales serán: la función que deseamos analizar
    descomplu(A,B); generaremos la Matriz A de tamaño nxn y una matriz B de tamaño 1xn con las que
    vamos a resolver el sistemas de Ecuaciones.
    En esta parte del código vamos a ubicar de manera estética visualmente la matriz A y B para ve se pueda
    observar la resolución.
    function x1 = descomplu (A , B )
    clc
    A2 = A ;
    B2 = B ;
    fprintf ( ' METODO DESCOMPOSICION LU ')
    disp ( newline ) ;
    n2 = size ( A2 ,1) ;
    (2)Para empezar a resolver el sistema de Ecuaciones vamos a generar un código con el cual se reorde-
    nara la matriz A y B .
    En donde se va a tratar de colocar los valores mas altos en la diagonal de la Matriz A antes de empezar
    la resolución de Descomposición LU.

    %Pivoteo parcial
    disp ( ' Paso 0: Reordenar ecuaciones ') ;
    for k = 1: length ( B )

    if A (k , k ) ~= max ( abs ( A (: , k ) ) )
    [ filapivote , z ] = find ( abs ( A ) == max ( abs ( A (: , k ) ) ) ) ;
    m =1;
    if length ( filapivote ) >1
    m = length ( filapivote ) ;
    end
    A ([ k , filapivote ( m ) ] ,:) = A ([ filapivote ( m ) ,k ] ,:) ;
    B ([ k , filapivote ( m ) ]) = B ([ filapivote ( m ) ,k ]) ;

    1
    end
    end
    (3)Ahora se realizara el código para resolver por el método de Descomposición LU.

    % D e s c o m p o s i c i n LU
    %~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ E l i m i n a c i n hacia adelante
    ~~~~~~~~~~~~~~~~~~~~~~~~~ %
    U = A2
    D = B2
    L = eye ( n2 ) ;
    for k2 = 1: n2 - 1
    for i2 = k2 + 1: n2
    factor = U ( i2 , k2 ) / U ( k2 , k2 ) ;
    L ( i2 , k2 ) = factor ;
    U ( i2 ,:) = U ( i2 ,:) - factor * U ( k2 ,:) ;
    D ( i2 ) = D ( i2 ) - factor * D ( k2 ) ;
    disp ( newline ) ;
    disp ( ' Matriz U ') ;
    disp ( U ) ;
    disp ( newline ) ;
    disp ( ' Matriz L ') ;
    disp ( L ) ;
    disp ( ' Matriz D ') ;
    disp ( D ) ;
    end
    end

    %~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ S u s t i t u c i n hacia a t r s
    ~~~~~~~~~~~~~~~~~~~~~~~~~ %

    x1 ( n2 ) = D ( n2 ) / U ( n2 , n2 ) ;
    for i2 = n2 - 1: -1:1
    sum = D ( i2 ) ;
    for j2 = i2 + 1: n2
    sum = sum - U ( i2 , j2 ) * x1 ( j2 ) ;
    end
    x1 ( i2 ) = sum / U ( i2 , i2 ) ;
    end
    x1 = x1 ';

    for i2 = 1: n2
    fprintf ( 'x %d = %f ' ,i2 , x1 ( i2 ) ) ;
    fprintf ( '\ n ') ;
    end
    (4)A continuación ejecutaremos el script en Matlab para verificar que se cumplan con todas la condi-
    ciones declaradas.

    Para ello utilizaremos un ejemplo de Sistema de Ecuaciones para resolver:

    En la siguiente Figura 1 se captura el sistema de ecuaciones a resolver el cual seria A y b

    2
    Figura 1: Ingresamos el sistema de ecuaciones

    A continuación se puede observar como se reordena las ecuaciones.

    Figura 2: Reordenar el Sistema de Ecuaciones

    Después que se a reordenado comienza a ejecutarse el código para resolver el sistema de ecuaciones
    mediante Descomposición LU, en este se puede observar como e revuelve paso por paso .

    3
    Figura 3: Proceso de resolución

    Al finalizar el programa entrega las respectivas soluciones .

    4
    Figura 4: Resutados Finales

    También podría gustarte