Método de Gauss-Jordán para Matrizes

Vista previa del material en texto

[INGENIERIA AMBIENTAL]
[BALANCE DE ENERGIA]
[14/12/2021] 
Código
%método de Gauss-Jordán (con procedimiento) 
%en el siguiente programa se puede determinar la solución de una matriz nxn 
%si en la matriz 2 se especifican los coeficientes o se puede determinar la
%inversa de una matriz si a la matriz 2 se le asigna el valor de la matriz
%identidad. La idea de este programa es ir realizando el procedimiento 
%paso a paso para finalmente llegar a la respuesta 
%ejemplo1: A= [1 2 3; 4 5 6; 7 8 9] y B= [3 4 5]'
A=input('ingrese la matriz 1 '); %% entrada de
B=input('ingrese la matriz 2 '); % datos %%
C=[A B]; %%unión de los datos en una solo matriz
for i=1:length(C(:,1)) %%para i desde la primera fila hasta el número de filas existentes 
if C(i,i)~=1 %%si el elemento i,i de la diagonal es diferente de 1 
 C(i,:)= C(i,:)./C(i,i); %entonces se convierte a 1 dividiendo toda la fila por dicho elemento
 disp(C) %salida de datos
end
%además el resto de elementos de la columna deben convertirse a 0 : 
%es decir si n es diferente de i ya que si i y n son iguales entonces el
%elemento se encuentra en la diagonal 
 for n=1:length(C(:,1)) %para n desde la primera fila hasta el número de filas existentes
 if n~=i % si n en la columna i no está en la diagonal es decir si i no es igual a n
 C(n,:)=-C(n,i).*C(i,:)+C(n,:); %entonces se convierte a 0 
 disp(C) 
 end 
 end 
end
manual
 1 0 -1 3
 0 1 1 2
 0 2 1 5
 1 0 -1 3
 0 1 1 2
 0 0 -1 1
 1 0 -1 3
 0 1 1 2
 0 0 1 -1
 1 0 0 2
 0 1 1 2
 0 0 1 -1
 1 0 0 2
 0 1 0 3
 0 0 1 -1
12