Anuncio

Colapsar
No hay ningún anuncio todavía.

Factorizacion LU (sin pivotaje)

Colapsar
X
 
  • Filtro
  • Hora
  • Mostrar
Borrar todo
nuevos mensajes

  • Matlab Factorizacion LU (sin pivotaje)

    Buenas a todos. Estoy intentando programar la Factorización LU en Matlab. Intenté el siguiente código pero sorprendentemente no me da el resultado esperado:

    Código:
    function B=milu_sept(A)
    dim=size(A); n=dim(1); %n filas=columnas
    
    for col=1:(n-1)
        for fil=(col+1):n
            A(fil,col:n)=A(fil,col:n)-(A(fil,col)/A(col,col))*A(fil,col:n);
        end
    end
    B=A;
    end
    Según lo que yo creo, este código debería efectivamente devolverme la factorización LU de una matriz A cualquiera... ¿alguien puede ofrecerme ayuda?

    Un saludo y muchas gracias!
    Última edición por skinner; 06/07/2012, 00:27:10.

  • #2
    Re: Factorizacion LU (sin pivotaje)

    ¡Ya lo he solucionado! Por si a alguien le interesa, el código correcto es (simplemente cambiar un "fil" por un "col"):

    Código:
    function B=milu_sept(A)
    dim=size(A); n=dim(1); %n filas=columnas
    
    for col=1:(n-1)
        for fil=(col+1):n
            A(fil,col:n)=A(fil,col:n)-(A(fil,col)/A(col,col))*A(col,col:n);
        end
    end
    B=A;
    end
    Un saludo!
    Última edición por skinner; 06/07/2012, 08:59:35.

    Comentario

    Contenido relacionado

    Colapsar

    Trabajando...
    X