Lab 4
Lab 4
Lab 4
EJERCICIO 1
55 I 1 +0 I 3−25 I 4=200
−20 I 1−25 I 3−0 I 4 =−450
0 I 1+37 I 3−4 I 4=250
Datos de entrada
clc
clear all
A=load('MATRIZ1.txt')
[f c]=size(A);
for k=1:c-1
A(k,:)=A(k,:)/A(k,k)
for i=k+1:f
A(i,:)=A(i,:)-A(k,:)*A(i,k)
end
end
% gauss-jordan
for k=f:-1:2
for i=k-1:-1:1
A(i,:)=A(i,:)-A(k,:)*A(i,k)
i=i-1
end
Resultados
A=
55 0 -25 200
-20 -25 0 -450
0 37 -4 250
A=
1.0000 0 -0.4545 3.6364
-20.0000 -25.0000 0 -450.0000
0 37.0000 -4.0000 250.0000
A=
1.0000 0 -0.4545 3.6364
0 -25.0000 -9.0909 -377.2727
0 37.0000 -4.0000 250.0000
A=
1.0000 0 -0.4545 3.6364
0 -25.0000 -9.0909 -377.2727
0 37.0000 -4.0000 250.0000
A=
1.0000 0 -0.4545 3.6364
0 1.0000 0.3636 15.0909
0 37.0000 -4.0000 250.0000
A=
1.0000 0 -0.4545 3.6364
0 1.0000 0.3636 15.0909
0 0 -17.4545 -308.3636
A=
1.0000 0 -0.4545 3.6364
0 1.0000 0.3636 15.0909
0 0 1.0000 17.6667
x=
0
0
0
x3=11.67
x2=8.67
xl=17.67
EJERCICIO 2
clc
clear all
A=load('MATRIZ2.txt')
[f c]=size(A);
for k=1:c-1
A(k,:)=A(k,:)/A(k,k)
for i=k+1:f
A(i,:)=A(i,:)-A(k,:)*A(i,k)
end
end
% gauss-simple
b=A(:,c);
x=zeros(f,1)
for i=f:-1:1
x(i)=(b(i)-A(i,i+1:f)*x(i+1:f))/A(i,i);
end
fprintf('x3=%2.2f\n x2=%2.2f\n xl=%2.2f\n ',x)
Resultados
A=
30 -7 0 0
-7 47 -10 0
0 -10 50 140
A=
1.0000 -0.2333 0 0
-7.0000 47.0000 -10.0000 0
0 -10.0000 50.0000 140.0000
A=
1.0000 -0.2333 0 0
0 45.3667 -10.0000 0
0 -10.0000 50.0000 140.0000
A=
1.0000 -0.2333 0 0
0 45.3667 -10.0000 0
0 -10.0000 50.0000 140.0000
A=
1.0000 -0.2333 0 0
0 1.0000 -0.2204 0
0 -10.0000 50.0000 140.0000
A=
1.0000 -0.2333 0 0
0 1.0000 -0.2204 0
0 0 47.7957 140.0000
A=
1.0000 -0.2333 0 0
0 1.0000 -0.2204 0
0 0 1.0000 2.9291
x=
0
0
0
x3=0.15
x2=0.65
xl=2.93
EJERCICIO 3
Para la producción de tres tipos de componentes eléctricos se requiere tres clases de material:
(metal, plástico y caucho), de lo cual se obtuvo las siguientes ecuaciones:
15 I x + 17 I y +19 I z=3.89
0.3 I x +0.4 I y + 0.55 I z=0 .095
1 I x +1.2 I y + 1.5 I z =0 .282
Datos de entrada
clc
clear all
A=load('MATRIZ3.txt')
[f c]=size(A);
for k=1:c-1
A(k,:)=A(k,:)/A(k,k)
for i=k+1:f
A(i,:)=A(i,:)-A(k,:)*A(i,k)
end
end
% gauss-simple
b=A(:,c);
x=zeros(f,1)
for i=f:-1:1
x(i)=(b(i)-A(i,i+1:f)*x(i+1:f))/A(i,i);
end
fprintf('x3=%2.2f\n x2=%2.2f\n xl=%2.2f\n ',x)
Resultados
A=
15.0000 17.0000 19.0000 3.8900
0.3000 0.4000 0.5500 0.0950
1.0000 1.2000 1.5000 0.2820
A=
1.0000 1.1333 1.2667 0.2593
0.3000 0.4000 0.5500 0.0950
1.0000 1.2000 1.5000 0.2820
A=
1.0000 1.1333 1.2667 0.2593
0 0.0600 0.1700 0.0172
1.0000 1.2000 1.5000 0.2820
A=
1.0000 1.1333 1.2667 0.2593
0 0.0600 0.1700 0.0172
0 0.0667 0.2333 0.0227
A=
1.0000 1.1333 1.2667 0.2593
0 1.0000 2.8333 0.2867
0 0.0667 0.2333 0.0227
A=
1.0000 1.1333 1.2667 0.2593
0 1.0000 2.8333 0.2867
0 0 0.0444 0.0036
A=
1.0000 1.1333 1.2667 0.2593
0 1.0000 2.8333 0.2867
0 0 1.0000 0.0800
x=
0
0
0
x3=0.09
x2=0.06
xl=0.08
EJERCICIO 4
Para el circuito mediante la ley de malla de corrientes se obtuvo las siguientes ecuaciones:
Datos de entrada
clc
clear all
A=load('MATRIZ4.txt')
[f c]=size(A);
for k=1:c-1
A(k,:)=A(k,:)/A(k,k)
for i=k+1:f
A(i,:)=A(i,:)-A(k,:)*A(i,k)
end
end
% gauss-jordan
for k=f:-1:2
for i=k-1:-1:1
A(i,:)=A(i,:)-A(k,:)*A(i,k)
i=i-1
end
Resultados
A=
60 -40 0 200
-40 150 -100 0
0 -100 130 -220
A=
1.0000 -0.6667 0 3.3333
-40.0000 150.0000 -100.0000 0
0 -100.0000 130.0000 -220.0000
A=
1.0000 -0.6667 0 3.3333
0 123.3333 -100.0000 133.3333
0 -100.0000 130.0000 -220.0000
A=
1.0000 -0.6667 0 3.3333
0 123.3333 -100.0000 133.3333
0 -100.0000 130.0000 -220.0000
A=
1.0000 -0.6667 0 3.3333
0 1.0000 -0.8108 1.0811
0 -100.0000 130.0000 -220.0000
A=
1.0000 -0.6667 0 3.3333
0 1.0000 -0.8108 1.0811
0 0 48.9189 -111.8919
A=
1.0000 -0.6667 0 3.3333
0 1.0000 -0.8108 1.0811
0 0 1.0000 -2.2873
x=
0
0
0
x3=2.82
x2=-0.77
xl=-2.29