Código Matlab Factorizacion Lu Cholesky

Para indicar que se trata de una matriz transpuesta usamos un apóstrofe precedido del nombre de la matriz. Ejemplo >> M=[2 3 1;1 4 6; 6 7 1] M= 2 3 1 1 4 6 6 7 1 >> M' ans = 2 1 3 4 1 6

6 7 1

Para calcular el determinante de una matriz se utiliza el comando det(A) >> det(M) ans = 12 Para calcular la matriz inversa de de M usamos inv(A) >> inv(M) ans = -3.1667 0.3333 1.1667 2.9167 -0.3333 -0.9167 -1.4167 0.3333 0.4167

CÓDIGO MATLAB FACTORIZACION LU CHOLESKY

% FACTORIZACION lU CHOLESKY clc %permite borrar el area de trabajo clear %permite borrar las variables almacenadas almacenadas format long %permite utilizar la maxima capacidad de la maquina fprintf(' FACTORIZACION LU CHOLESKY\n\n\n'); %fprintf me permite ingresar comentarios de manera textual que pueden %orientar al usuario en el uso del programa %input es un comando de solicitud de entrada de datos del usuario. A=input('Ingrese la matriz A = \n');

% Las matrices A y b deben ser ingresadas entre corchetes separando las %columnas mediante coma ',' y las filas mediante punto y coma ';'. [n,m]=size(A); if n==m for k=1:n %La instrucción iterativa for permite repetir estamentos a un %numero específico de veces suma1=0; for p=1:k-1 suma1=suma1+L(k,p)*u(p,k); end L(k,k)=sqrt(A(k,k)-suma1); u(k,k)=L(k,k); u(k,k)=L(k,k); %princio del metodo for i=k+1:n suma2=0; for q=1:k-1 suma2=suma2+L(i,q)*u(q,k); end L(i,k)=(A(i,k)-suma2)/L(k,k); %obtencion de la matriz L

endfprintf('\n Matriz L:\n') disp(L)