3.3 Metodo de Gauss-Jordan..
3.3 Metodo de Gauss-Jordan..
3.3 Metodo de Gauss-Jordan..
3 Mtodo de Gauss-Jordan
Es un Mtodo exacto que nos permite obtener la solucin de un sistema de ecuaciones
lineales, llamado tambin mtodo de eliminacin completa.
PROCEDIMIENTO:
*1.- Si es posible, intercambiar las ecuaciones en el sistema tratando de que algn trmino
con valor de 1 quede en la diagonal principal
2.- Multiplicar todos los elementos de un rengln por una constante que viene a ser igual
al valor del trmino que se desea eliminar con signo opuesto.
3.- Sumar el resultado del paso anterior a la ecuacin o rengln para eliminar ese trmino.
El procedimiento anterior se realiza para todos los elementos del sistema hasta obtener al
lado izquierdo (trminos dependientes) una matriz identidad y al lado derecho un vector
de datos que corresponde al conjunto solucin.
Ejemplo.EDU A = magic(3)
A =
8
1
6
3
5
7
4
9
2
EDU b = [-4;1;3];
EDU B = [A b]
B =
8
1
6
-4
3
5
7
1
4
9
2
3
EDU
EDU B(1,:) = B(1,:)*(1/8)
B =
1.0000
0.1250
0.7500 -0.5000
3.0000
5.0000
7.0000
1.0000
4.0000
9.0000
2.0000
3.0000
EDU B(2,:) = B(1,:)*(-1)*B(2,1) + B(2,:)
B =
1.0000
0.1250
0.7500
-0.5000
0
4.6250
4.7500
2.5000
4.0000
9.0000
2.0000
3.0000
EDU B(3,:) = B(1,:)*(-1)*B(3,1) + B(3,:)
B =
1.0000
0.1250
0.7500
-0.5000
0
4.6250
4.7500
2.5000
0
8.5000
EDU B(2,:) = B(2,:) *
B =
1.0000
0.1250
0
1.0000
0
8.5000
EDU B(1,:) = B(2,:) *
B =
1.0000
0
0
1.0000
0
8.5000
EDU B(3,:) = B(2,:) *
B =
1.0000
0
0
1.0000
0
0
EDU B(3,:) = B(3,:) *
B =
1.0000
0
0
1.0000
0
0
EDU B(1,:) = B(3,:) *
B =
1.0000
0
0
1.0000
0
0
EDU B(2,:) = B(3,:) *
B =
1.0000
0
0
1.0000
0
0
EDU %ahora se hace la
EDU sol = B(:,4);
EDU A*sol
ans =
-4.00000000000000
1.00000000000000
3.00000000000000
EDU b
b =
-4
1
3
-1.0000
1/B(2,2)
5.0000
0.7500
-0.5000
1.0270
0.5405
-1.0000
5.0000
(-1)*B(1,2) + B(1,:)
0.6216
-0.5676
1.0270
0.5405
-1.0000
5.0000
(-1)*B(3,2) + B(3,:)
0.6216
1.0270
-9.7297
1/B(3,3)
-0.5676
0.5405
0.4054
0.6216
-0.5676
1.0270
0.5405
1.0000
-0.0417
(-1)*(B(1,3)) + B(1,:)
0
-0.5417
1.0270
0.5405
1.0000
-0.0417
(-1)*(B(2,3)) + B(2,:)
0
-0.5417
0
0.5833
1.0000
-0.0417
comprobacin...
x=[-5:0.1:5];
EDU y=sqrt(5-x.^2);
EDU y2=-sqrt(5-x.^2);
EDU y3=1-exp(x);
EDU plot(x,y,'g')
EDU hold on
EDU plot(x,y2,'r')
EDU plot(x,y3,'b')
EDU axis([-sqrt(5)-0.1 sqrt(5)+0.1 -sqrt(5) sqrt(5) ])
EDU axis square
x = -2.0881
EDU y = 1 - exp(x)
y = 0.8761
EDU x = - sqrt(5 - y.^2), y = 1 - exp(x)
x = -2.0573
y = 0.8722
EDU x = - sqrt(5 - y.^2), y = 1 - exp(x)
x = -2.0589
y = 0.8724
EDU
EDU x=[-2.5:0.1:2.5];
EDU y=sqrt(5-x.^2);
EDU y2=-sqrt(5-x.^2);
EDU y3=1-exp(x);
EDU plot(x,y,'g')
EDU hold on
EDU plot(x,y2,'r')
EDU plot(x,y3,'b')
EDU axis([-sqrt(5)-0.1 sqrt(5)+0.1 -sqrt(5) sqrt(5) ])
EDU axis square
EDU grid
2.- Obtenemos los valores iniciales:
1. Segundo cuadrante (-2,0.8)
2. Cuarto cuadrante (1.2,-1.9)
3.- Identificar f(x,y) y g(x,y) y sus derivadas parciales
f(x,y) = x2 + y2 5
g(x,y) = exp(x) + y 1
;
=0
g(x,y) +
=0
= f(x,y)
= g(x,y)
Regla de Cramer
*NOTA: Xn+1 quiere decir que es el siguiente valor ya que es un mtodo iterativo
En Xn+1 y yn+1
(
)
(
)
(
5.- sustituir los valores iniciales por sistema en las ecuaciones generales en forma iterativa
hasta que el criterio de evaluacin se cumpla.
a) segundo cuadrante: (-2,0.8)
b) cuarto cuadrante: (1.2,-1.9)
% Del segundo cuadrante
% con un programilla
EDU x=-2; y=0.8;
EDU for n=1:10,
XXX= x + (-(x.^2+y.^2-5)+(exp(x)+y-1)*(2*y))./(2*x-(exp(x)*2*y));
YYY= y + (-(exp(x)+y-1)*2*x +(x.^2+y.^2-5)*exp(x))./(2*x(exp(x)*2*y));
x= XXX; y= YYY;
fprintf ('%3d %13.10f %13.10f\n',n,x,y)
end
1 -2.0608405633 0.8728985916
2 -2.0588637617 0.8724013787
3 -2.0588629045 0.8724010203Este mtodo es muy rpido!!!
4 -2.0588629045 0.8724010203
5 -2.0588629045 0.8724010203
6 -2.0588629045 0.8724010203
7 -2.0588629045 0.8724010203
8 -2.0588629045 0.8724010203
9 -2.0588629045 0.8724010203
10 -2.0588629045 0.8724010203
% Del cuarto cuadrante
EDU x=1.2; y= -1.9;
EDU for n=1:7,
XXX= x + (-(x.^2+y.^2-5)+(exp(x)+y-1)*(2*y))./(2*x-(exp(x)*2*y));
YYY= y + (-(exp(x)+y-1)*2*x +(x.^2+y.^2-5)*exp(x))./(2*x(exp(x)*2*y));
x= XXX; y= YYY;
fprintf ('%3d %13.10f %13.10f\n',n,x,y)
end
1 1.0903572462 -1.9560901603
2 1.0838117343 -1.9558617079
3 1.0837904912 -1.9558625135Este mtodo es muy rpido!!!
4 1.0837904910 -1.9558625135
5 1.0837904910 -1.9558625135
6 1.0837904910 -1.9558625135
7 1.0837904910 -1.9558625135
TAREA
Encontrar las intersecciones de
f(x,y)= 3*sin(x)-(x-2).^2-y= 0
g(x,y)= (x+1).^2 + (y+1).^2-7= 0
x=[-5:0.1:5];
y1=3*sin(x)-(x-2).^2;
y2= (sqrt(7-(x+1).^2)-1);
y3= (-sqrt(7-(x+1).^2)-1);
plot(x,y1,'b',x,y2,'r',x,y3,'c'),grid, hold on
axis([-4 2 -4 2])
valores iniciales
1re. Cuadrante(0.9 0.9)
4To cuadrante(0.1 -3.5)
Con (x1,y) y (0.9 0.9)
EDU x=0.9; y=0.9;
EDU for n=1:13,
x= asin(((x-2).^2 + y)/3);
y= sqrt(7 - (x + 1).^2)-1;
fprintf('%3d %13.9f %13.9f\n',n,x,y)
end
1
0.780075841
0.957378349
2
0.953070573
0.784801204
3
0.677637355
1.045857498
4
1.198515541
0.471913522
5
0.380549039
1.257007831
6
1.570796327
0.567409702
7
0.054268723
1.509413884
8
1.029577578
1.382014039
9
0.151257885
1.627146412
10
0.698233270
1.601145136
11
0.264920156
1.682964452
12
0.535319655
1.675634274
13
0.331991019
1.706790391
************DIVERGE*******************
EDU x=0.9; y=0.9;
EDU for n=1:4,
x= sqrt(3*sin(x) - y) + 2;
y= sqrt(7 - (x + 1).^2)-1;
fprintf('%3d %13.9f %13.9f\n',n,x,y)
end
1
3.204151456 -1.000000000
2
3.445542962
0.373819629
3
2.112510470
0.987008600
4
4.122259993 -0.799617553
************DIVERGE*******************
En el 4to cuadrante (0.1 -3.5)
Con (x1,y)
EDU x=0.1; y=-3.5;
EDU for n=1:4,
x= asin(((x-2).^2 + y)/3);
y= -sqrt(7 - (x + 1).^2)-1;
fprintf('%3d %13.9f %13.9f\n',n,x,y)
end
1
0.036674888 -3.434194975
2
0.140613084 -3.387258217
3
0.023355952 -3.439824296
4
0.156402667 -3.379649737
EDU x=0.1; y=-3.5;
EDU for n=1:54,
x= asin(((x-2).^2 + y)/3);
y= -sqrt(7 - (x + 1).^2)-1;
fprintf('%3i %13.9f %13.9f\n',n,x,y)
end
1 0.036674888 -3.434194975
2 0.140613084 -3.387258217
3 0.023355952 -3.439824296
4 0.156402667 -3.379649737
5 0.006400507 -3.446866981
6 0.176776593 -3.369640659
7 -0.015166270 -3.455626707
8 0.203150592 -3.356359194
9 -0.042576662 -3.466442895
10 0.237450581 -3.338528610
11 -0.077393289 -3.479676765
12 0.282365973 -3.314203429
13 -0.121611799 -3.495683107
14 0.341801059 -3.280256547
15 -0.177813127 -3.514758188
16 0.421764113 -3.231274705
17 -0.249392600 -3.537043265
18 0.532367480 -3.156814759
19 -0.340851629 -3.562327736
20 0.693309597 -3.032904968
21 -0.457628330 -3.589562313
22 0.955827730 -2.781779416
23 -0.599013368 -3.615188276
24 1.570796327 -2.029426639
25 -0.683548966 -3.635469980
26 1.113869476 -2.950936211
27 -1.046264619 -3.667649350
28 1.154736096 -3.355504503
29 -1.271157919 -3.812300341
30 0.706044613 -3.888800105
31 -0.934031933 -4.011056893
32 0.276955734 -4.024768041
33 -0.622663738 -4.043407064
34 0.055290767 -4.018291818
35 -0.441145327 -4.020390671
36 -0.053204258 -3.993506281
37 -0.339230949 -3.993831668
38 -0.109421614 -3.973354513
39 -0.280331712 -3.973674259
40 -0.140442745 -3.959506102
41 -0.245185562 -3.959854621
42 -0.158421517 -3.950409837
43 -0.223694542 -3.950712106
44 -0.169209395 -3.944521699
45 -0.210328756 -3.944752358
46 -0.175836676 -3.940729529
47 -0.201921643 -3.940893701
48 -0.179972099 -3.938291614
49 -0.196594116 -3.938403903
50 -0.182579167 -3.936725288
51 -0.193201744 -3.936800295
52 -0.184233712 -3.935719139
53 -0.191034839 -3.935768521
54 -0.185288291 -3.935072851
************DIVERGE*******************
Con (x2,y)
EDU x=0.1; y=-3.5;
EDU for n=1:4,
x=sqrt(3*sin(x) - y) +2;
y=-sqrt(7 - (x+1).^2) -1;
fprintf('%3i %13.9f %13.9f\n',n,x,y)
end
1
3.949230682 -1.000000000
2
3.259906666 -3.006607410
3
4.578986164 -3.385683903
4
3.030928573 -4.713223747
************DIVERGE*******************
EDU x=0.9; y=0.9;
EDU for n=1:4,
y=3*sin(x)-(x-2).^2;
x= sqrt(7 - (x + 1).^2)-1;
fprintf('%3i %13.9f %13.9f\n',n,x,y)
end
1
0.841195264
1.139980729
2
0.900000000
0.893492730
3
0.841195264
1.139980729
4
0.900000000
0.893492730
************DIVERGE*******************
EDU %cuarto cuadrante (y,x2) y (-0.1,3.5)
EDU x=0.1; y=-3.5;
EDU for n=1:4,
y=3*sin(x)-(x-2).^2;
x= sqrt(7 - (y + 1).^2)-1;
fprintf('%3i %13.9f %13.9f\n',n,x,y)
end
1
0.289027116 -3.310499750
2
1.418682559 -2.072368722
3 -1.000000000
2.627429009
4 23.611170832 -18.043473358
*este mtodo no va a funcionar en todos los casos
f(x,y)= 3*sin(x)-(x-2).^2-y= 0
g(x,y)= (x+1).^2 + (y+1).^2-7= 0
Acomodando signos
(