LH.rar_The Method Method

clc;clear all;close all; fprintf(' Universidad Autónoma de Nuevo León \n'); fprintf(' Facultad de Ingeniería Mecánica y Eléctrica\n'); fprintf('M.C. de la Ing. Electrica con Orientacion en S.E.P. \n'); fprintf(' ANALISIS DE SISTEMAS ELECTRICOS DE POTENCIA \n'); fprintf(' Salvador Alvarez Caballero \n'); fprintf(' Solucion de Matriz Cuadrada por Factorizacion LH \n\n') fprintf('Nota: A=LH en el sistema:\n') fprintf('AV(Volts)=I(Corriente); entonces\n') fprintf('LX=I y HV=X\n\n') N=input('Indique el Orden de la Matriz='); %Matriz de Trabajo for m=1:N for n=1:N fprintf('Valor del Elemento A(%d,%d)=' , m,n ) %Matriz Primaria A(m,n)=input(''); end end fprintf('\n\nMatriz Capturada\n\n'); disp(A) [X,Y]=lu(A); fprintf('Matriz L:\n\n'); L=Y'; disp(L) fprintf('Matriz H:\n\n'); H=X'; disp(H) n=length(L); fprintf(' Capture el Vector de Corrientes\n\n'); %Vector de Corrientes for i=1:N fprintf('Valor del Elemento I(1,%d)=', i); %Valores de la Matriz I(i,1)=input(''); end L=[L,I]; L(1,n+1)=L(1,n+1)/L(1,1); for i=2:n suma=0; for j=1:i-1 suma=suma+L(i,j)*L(j,n+1); end L(i,n+1)=(L(i,n+1)-suma)/L(i,i); end fprintf('\n\nResultados de Solucion Progresiva (L)\n\n'); for i=1:N fprintf('Valor %d = %d\n', i,L(i,n+1)); %Valores de la Matriz Z(i,1)=L(i,n+1); end H=[H,Z]; fprintf('\n\nResultados de Solucion Regresiva(H)\n\n'); V(N)=H(N,N+1)/H(N,N); % Cálculo de la última incógnita for i=N-1:-1:1 Suma=0; for j=i+1:N Suma=Suma+H(i,j)*V(j); end V(i)=(H(i,N+1)-Suma)/H(i,i); % Cálculo de las incógnitas end for i=1:N fprintf('V(%d) = %d \n',i,V(i)); % Presentación de resultados end