Scribd
Cargar
88 vistas

3 4

Este documento presenta la resolución de un sistema de ecuaciones lineales de 3 ecuaciones y 3 incógnitas utilizando el método de eliminación de Gauss. Se obtiene una matriz triangular superior y se resuelve para encontrar que x3=0.5, x2=-6.5 y x1=1.5. También se resuelve otro sistema de 3 ecuaciones utilizando eliminación de Gauss simple, con pivoteo parcial y el método de Gauss-Jordan, obteniendo en los tres casos que x3=9/4, x2=-1/2 y x1=-1

Cargado por

Ximena Ac
Derechos de autor
© © All Rights Reserved
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
88 vistas

3 4

Este documento presenta la resolución de un sistema de ecuaciones lineales de 3 ecuaciones y 3 incógnitas utilizando el método de eliminación de Gauss. Se obtiene una matriz triangular superior y se resuelve para encontrar que x3=0.5, x2=-6.5 y x1=1.5. También se resuelve otro sistema de 3 ecuaciones utilizando eliminación de Gauss simple, con pivoteo parcial y el método de Gauss-Jordan, obteniendo en los tres casos que x3=9/4, x2=-1/2 y x1=-1

Cargado por

Ximena Ac
Derechos de autor
© © All Rights Reserved
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
Está en la página 1/ 5

3.

Use la eliminación de Gauss para resolver:


8X1 + 2X2 - 2X3 = -2
10X1 + 2X2 + 4X3 = 4
12X1 + 2X2 + 2X3 = 6
12 2 26
8 2 −2−2
10 2 44
Para hacer uso de la eliminación de Gauss, tomaré la fila 2 y se sumará con -8/12 de toda la fila
uno, donde la matriz nos quedará de la sgte. forma:
12 2 26
0 −2/3 10 /3 6
10 2 44
Luego, tomaré la fila 3 y la sumaré con 10/12 de la fila uno, y la matriz quedará:
12 2 26
0 −2/3 10/3 6
0 −1 /3 7 /3 1
Finalmente, para obtener la matriz triangular superior, tomaremos la fila tres y la sumaremos
con -1/3 de la fila dos.
12 2 26
0 −2/3 10/3 6
0 0 42
Para encontrar los valores de las variables, pasaremos la matriz a un sistema de ecuaciones:
12X1 + 2X2 - 2X3 = 6
-2/3X2 + 10/3X3 = 6
4X3 = 2 Entonces: x3=0.5 , x2=-6.5 y x1=1.5
Finalmente, reemplazando los valores en el primer sistema de ecuaciones:
8(1.5) + 2(6.5) – 2(0.5) = -2
10(1.5) + 2(6.5) + 4(0.5) = 4
12(1.5) + 2(6.5) + 2(0.5) = 6
Código en Scilab:
function x=gauss(a, b)
ab = [a,b];
[m, n] = size(ab);
for j = 1:m-1
for i = j+1:m
ab(i,j:n) = ab(i,j:n)-ab(i,j)/ab(j,j)*ab(j,j:n);
end
end
x = zeros(m,1);
x(m) = ab(m,n)/ab(m,m);
for i = m-1:-1:1
x(i)=(ab(i,n)-ab(i,i+1:m)*x(i+1:m))/ab(i,i);
end
endfunction

4.- Resuelva el sistema:


X1 + X2 - X3 = -3
6 X1 + 2X2 + 2X3 = 4
-3 X1 + 4X2 + X3 = 1
Por medio de a) eliminación de Gauss simple, b) eliminación de Gauss con pivoteo parcial
y c) método de Gauss-Jordan sin pivoteo parcial.
a) Eliminación de Gauss.
1 1 −1−3
6 2 24
−3 4 11
Para obtener la matriz triangular superior y resolver por eliminación de Gauss, se tomará la
segunda fila y se restará con 6 veces la fila uno:
12 2 26
0 4 −8−20
−3 4 11
Siguiendo con el problema, se tomará la tercera fila y se restará con -3 la primera fila:
12 2 26
0 4 −8−20
0 −7 28
Finalmente, tomaré la segunda fila y la restaré con -7/4 de la segunda fila:
12 2 26
0 4 −8−20
0 0 1 227
Entonces pasando a un sistema de ecuaciones, se tendrá:
X1 + X2 - X3 = -3
4X2 - 8X3 = -20
12X3 = 27 donde: x3=9/4 , x2=-1/2 y x1=-1/4
b) Eliminación de Gauss con pivoteo parcial.
6 2 22
1 1 −1−3
−3 4 11

Teniendo esta matriz, se tomará la fila dos y se sumará con 1/6 de la fila uno:
6 2 22
4
0 −2/3 10/ 3
3
−3 4 11
Luego tomaré la tercera fila y la sumaré con -1/2 de la primera fila:
6 2 22
4
0 −2/3 10 /3
3
0 −5 −2−2
Y finalmente se tomará la tercera fila y se sumará con 15/2 de la segunda fila:
6 2 22
4
0 −2/3 10 /3
3
0 0 1 2 27
Finalmente se tomará el sistema de ecuaciones con la matriz :
6X1 + 2X2 +2 X3 = 2
-2/3X2 +4/3X3 = 10/3
12X3 = 27 donde: x3=9/4 , x2=-0.5000000001 y x1=-1/4

c) Gauss Jordan sin pivoteo parcial.


Reordenaremos la matriz:
1 1 −1 3
6 2 22
−3 4 1 1
Luego, se tomará la segunda fila y se sumará con 6 veces la fila uno:
1 1 −1 3
0 4 −8−2 0
−3 4 11
Tomaré la tercera fila y se sumará con -3 la primera fila:
1 1 −13
0 4 −8−20
0 −7 28
Finalmente tomaré la tercera fila y la sumaré con -7 veces la fila dos:
1 1 −1 3
0 4 −8−20
0 0 1 2 27
Se pasa esta matriz a un sistema de ecuaciones:
X1 +X2 - X3 = 3
4X2 -8/3X3 = -20
12X3 = 27 donde: x3=9/4 , x2=-1/2 y x1=-1/4
Códigos en Scilab:
a)
function x=gauss(a, b)
ab = [a,b];
[m, n] = size(ab);
for j = 1:m-1
for i = j+1:m
ab(i,j:n) = ab(i,j:n)-ab(i,j)/ab(j,j)*ab(j,j:n);
end
end
x = zeros(m,1);
x(m) = ab(m,n)/ab(m,m);
for i = m-1:-1:1
x(i)=(ab(i,n)-ab(i,i+1:m)*x(i+1:m))/ab(i,i);
end
endfunction
b)
AUM = [A b]
X = AUM
[rX cX] = size(X)
//Eliminación hacia adelante
for p = 1 : 1 : rX-1 //Fila del pivote
pivote = X(p,p);
for i = p + 1 : 1 : rX //Fila a procesar
prim_fila = X(i,p)
X(i,:) = X(i,:) - (X(p, :)/pivote)*prim_fila
end
end
x(rX) = X(rX, cX)/X(rX, cX - 1)
for i = rX : -1 : 1
s=0
for j = i + 1 : rX
s = s + X(i, j) * x(j)
end
x(i) =( X(i,cX) - s ) / X(i,i)
end

c)
function vector=EliminacionGaussiana(A,B)

%A=input('Ingrese Los Coeficientes de las ecuaciones:\n ');


%B=input('Ingrese Los Terminos del lado derecho de las ecuaciones:\n ');
n=length(B); %Muestra el largo del vector de soluciones
mat=[A,B]; %Matriz aumentada
for k=1: n-1
for i=k+1: n
multi=mat(i,k)/mat(k,k);
for j=1: n+1
mat(i,j)=mat(i,j)-(multi*mat(k,j));
end
end
end
fprintf('\n La matriz final es:\n');
disp(' ');
disp(mat)
fprintf('\n Las soluciones son :\n');
for i = n : -1 : 1 % Va solucionando el sistema de ecuaciones
suma = 0;
for s = i + 1 : n
suma = suma + mat(i,s) * x(s,1);
end
x(i,1) = (mat(i,n+1) - suma) / mat(i,i); % despeje final
end
for i=1:n
xi=x(i,1);
fprintf('\nX%g=',i)
disp(xi)
end
vector=1;
end

También podría gustarte