Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Sistema Dee Cuac I Ones

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 8

1)Resuelva el sistema por el método de iteración de gauss seidel continuando el proceso hasta

que la próxima diferencia entre valores sucesivos de x y o z sea menor que 0.02.

10𝑥 + 2𝑦 + 6𝑧 = 28
𝑥 + 10𝑦 + 9𝑧 = 7
−2𝑥 + 7𝑦 + 10𝑧 = 17
Termine los cálculos si no se verifica dicho creiterio de detención luego de haber efectuado 5
iteraciones.

Iteracion 1
X=2.8 ea=100
Y=0.42 ea=100
Z=1.9660 ea=100
--------------------------------------------
Iteracion 2
X=1.5364 ea=82.2442%
Y=-1.223 ea=134.3407%
Z=2.8634 ea=31.3406%
--------------------------------------------
Iteracion 3
X=1.3266 ea=15.8181
Y=-2.0097 ea=39.1439
Z=3.3721 ea=15.0858
--------------------------------------------
Iteracion 4
X=1.1787 ea=12.5472%
Y=-2.4528 ea=18.0633%
Z=3.6527 ea=7.6809%
--------------------------------------------
Iteracion 5
X=1.0989 ea=7.2546%
Y=-2.6973 ea=9.0657%
Z=3.8079 ea=4.0764%
--------------------------------------------
3)un ingeniero supervisa la producción de tres tipos de automóviles. Se requieren tres clases
de materiales (metal plástico y caucho) para la producción. La cantidad necesaria para producir
cada modelo de auto es de:
Automóvil Metal plastico Caucho
Modelo1 1500 25 100
Modelo2 1700 33 120
Modelo3 1900 42 160

Si se dispone de un total de 106 toneladas de metal, 2.17 toneladas de plástico y 8.2 toneladas
de caucho ¿Cuántos automóviles se pueden producir por dia de cada modelo?.

1500𝑥 + 1700𝑦 + 1900𝑧 = 106000


25𝑥 + 33𝑦 + 42𝑧 = 2170
100𝑥 + 120𝑦 + 160𝑧 = 8200

1500 1700 1900


A= 25 33 42 .
100 120 160

0.0018462 -0.3384615 0.6692308

A-1= 0.0015385 0.3846154 -1.1923077

-0.0023077 -0.0769231 0.5384615

106000
B= 2170
8200

10
A-1xB =20
30

X=10 y=20 z=30


4)Considere el método de gauss-seidel para obtener la solución al siguiente sistema de
ecuaciones lineales. Asegure 3 cifras significativas osea un error porcentual de 0.05%

3𝑥 − 0.1𝑦 − 0.2𝑧 = 7.85


0.1𝑥 + 7𝑦 − 0.3𝑧 = −19.3
0.3𝑥 − 0.2𝑦 + 10𝑧 = 71.4
Iteracion 1
X=2.6167 ea=100
Y=-2.7945 ea=100
Z=7.0056 ea=100
--------------------------------------------
Iteracion 2
X= 2.9906 ea=12.5023
Y=-2.4996 ea=11.7977
Z= 7.0003 ea=0.076

--------------------------------------------
Iteracion 3
X=3 ea=0.31580
Y=-2.5 ea=0.0145
Z=7 ea=0.0042

--------------------------------------------
Iteracion 4
X=3 ea=0.0011
Y=-2.5 ea=0.0005
Z=7 ea=0.0000

--------------------------------------------
5)Use el método de Gauss jordan para obtener la solución al siguiente sistema de ecuaciones
lineales.

3𝑥 − 0.1𝑦 − 0.2𝑧 = 7.85


0.1𝑥 + 7𝑦 − 0.3𝑧 = −19.3
0.3𝑥 − 0.2𝑦 + 10𝑧 = 71.4
A =

3. -0.1 -0.2 7.85


0.1 7. -0.3 -19.3
0.3 -0.2 10. 71.4

--> rref(A)
ans =

1. 0. 0. 3.
0. 1. 0. -2.5
0. 0. 1. 7.

Se observa que la matriz de coeficientes se ha transformado en la matriz identidad y la


solución se obtiene en el lado derecho de la matriz.
2.Se desea hallar los voltajes en los nodos a, b y c de una red eléctrica conectada en tres
terminales con voltajes conocidos. Se obtuvieron las tres ecuaciones siguientes:

(ea-20) /2+(ea-eb) /4+(ea-ec) /3=0

(eb-ea) /4+(eb) /3+(eb-ec) /3=0

(ec-5) /3+(ec-ea) /3+(ec-eb) /3=0

Emplee el método de Gauss-Jordán y muestre todas las etapas del cálculo necesario ea, eb ,ec
voltaje en voltios.

Luego de despejar, obtenemos la siguiente matriz dominante:

13 | -3 |-4

15 | 47 |-12

-3 | -3 | 9

A partir de operaciones en filas y columnas obtenemos la matriz identidad:

13 | -3 | -4 :120

15 | 47 | -12 :0

-3 | -3 | 9 :15

Cambiar columna 1 con 2:

15 | 47 | -12 :0

13 | -3 | -4:120

-3 | -3 | 9 : 15

- 13/15 × (fila 1) con fila 2:

15 | 47 | -12 :0

0 | -656/15 | 17/5 :120

-3 | -3 |9 :15

+ 1/5 × (fila 1) con fila 3:

15 | 47 | -12 :0

0 | -656/15 | 17/5 :120

0 | 32/5 | 33/5: 15

+6/41 × (fila 2) con fila 3:

15 | 47 | -12 :0
0 | -656/15 | 17/5 :120

0 |0 | 291/41 :1335/41

Por fila 3 entre 41/291:

15 | 47 | -12 :0

0 | -656/15 | 17/5 :120

0 |0 |1 :445/97

-17/5 × (fila 3) con fila 2:

15 | 47 | -12 :0

0 | -656/15 | 0 :10127/97

0 |0 | 1 :445/97

+ 12 × fila 3 con fila 1:

15 | 47 | 0 :5340/97

0 | -656/15 | 0 :10127/97

0 |0 | 1 :445/97

Por fila 2entre -15/656:

15 | 47 | 0 :5340/97

0 | 1 | 0: -3705/1552

0 | 0 | 1 :445/97

-47 × (fila 2) con fila 1:

15 | 0 | 0 :167.251933

0 | 1 | 0: -3705/1552

0 | 0 | 1 :445/97

divido fila 1entre 15:

1 | 0 | 0 :10.072816

0 | 1 |0: -3705/1552

0 | 0 | 1 :445/97

X=10.072816;

Y=-2.1116505;

Z=4.3203883
function gauss_seidel
A=input('ingrese la matriz de coeficientes:');
B=input('ingrese las constantes:');
es=input('ingrese el %de error:');
x=[0 0 0];
xs=[0 0 0];
i=input('ingrese las iteraciones maximas: ');
for k=1:i
//calculando las nuevas aproximaciones
x(1)=(B(1)-A(1,2)*x(2)-A(1,3)*x(3))/A(1,1);
x(2)=(B(2)-A(2,1)*x(1)-A(2,3)*x(3))/A(2,2);
x(3)=(B(3)-A(3,1)*x(1)-A(3,2)*x(2))/A(3,3);
//calculando los porcentajes de error
error1=abs((xs(1)-x(1))/x(1))*100;
error2=abs((xs(2)-x(2))/x(2))*100;
error3=abs((xs(3)-x(3))/x(3))*100;
//muestra los resultados en la pantalla
printf('Iteracion%2d\n',k);
printf('\t X\t Y\t Z\n');
printf('\t%1.4f\t%1.4f\t%1.4f\n',x(1),x(2),x(3));
printf('\terror1\terror2\terror3\n');
printf('\t%1.4f\t%1.4f\t%1.4f\n',error1,error2,error3);
disp('--------------------------------------------');
//termina si se llega al error
if error1<=es & error2<=es & error3<=es
break
end
xs=x;
end
endfunction

También podría gustarte