Modelos lineales
Sistemas lineales de ecuaciones como modelos computacionales
1.1 Modelo de Gauss Jordan (utiliza matrices)
1.2 Modelo de Gauss Seidel (método iteractivo)
a) Modelo con solución única
b) Modelo con infinita soluciones
c) Modelo sin solución
Gauss Jordan
(Listo)
Gauss Seidel
1. Verificar que el sistema sea diagonalmente dominante
4*X1 + X2 – X3 = 1 (4*X1= 1-X2+X3)
-X1 – 4*X2 + X3 = -1 (-4*X2= -1+X1-X3)
X1 + X2 + 4*X3 = 1 (4*X3= 1 -X1-X2)
2. Se reescriben las ecuaciones despejando cada incógnita en cada una
de las ecuaciones.
X1 = ¼ -1/4*X2 +1/4*X3 (Ecuación 1)
X2=1/4 - 1/4* X1 +1/4*X3 (Ecuación 2)
X3= ¼ - ¼*X1 – ¼*X2 (Ecuación 3)
3. Inicializo cada variable con un vector ceros para comenzar los
cálculos (X1, X2, X3) = (0, 0, 0)
i X1 X2 X3 |X1(i+1)-X1(i|)<0.001 |X2(i+1)-X2(i)|<0.001 |X3(i+1)-X3(i)|<0.001
0 0 0 0 0 0 0
1 0.2500 0.1875 0.1406 0,2500 0.1875 0.1406
2 0.2382 0.2256 0.1340 0.0118 0.0381 0.0066
3 0,2271 0.2267 0.1365 0.0111 0.0011 0.0025
4
Primera iteraccion i=1
X1 = ¼ -1/4*X2 +1/4*X3 (Ecuación 1)
X2=1/4 - 1/4* X1 +1/4*X3 (Ecuación 2)
X3= ¼ - ¼*X1 – ¼*X2 (Ecuacion 3)
X1= ¼ -1/4*(0) + ¼*(0) = 0.2500
X2= ¼ -1/4*(0.25) +1/4*(0) = 0.1875
X3= ¼ -1/4*(0.25) -1/4*(0.1875) = 0.1406
Segunda iteraccion i=2
X1 = ¼ -1/4*X2 +1/4*X3 (Ecuación 1)
X2=1/4 - 1/4* X1 +1/4*X3 (Ecuación 2)
X3= ¼ - ¼*X1 – ¼*X2 (Ecuacion 3)
X1=0.2382
X2=0.2256
X3=0.1340
Tercera iteraccion i=3
X1 = ¼ -1/4*X2 +1/4*X3 (Ecuación 1)
X2=1/4 - 1/4* X1 +1/4*X3 (Ecuación 2)
X3= ¼ - ¼*X1 – ¼*X2 (Ecuacion 3)
X1=
X2=
X3=