Modelos lineales(sistemas de ecuaciones lineales)
a) Método de Gauss Jordan (Matrices)
b) Método de Gauss Seidel (método iteractivo)
1. Sistema con solución única
2. Sistema con infinitas soluciones
3. Sistema sin solución
Metodo de Gauss Seidel
1. El sistema debe ser diagonalmente dominante
4*X1 –X2 +X3 = 1
-X1 +4*X2 +X3 = -1
X1 +X2 -4*X3 = 1
2. Se despejan cada una de las incógnitas en cada una de las ecuaciones
X1= 1/4 +1/4*X2 -1/4*X3 (Ecuación 1)
X2= -1/4 +1/4*X1 -1/4*X3 (Ecuación 2)
X3= -1/4 +1/4*X1 +1/4*X2 (Ecuación 3)
3. Inicializar las incógnitas con un vector de ceros
(X1=0 ; X2=0 ; X3=0)
Primera iteración i=1
i X1 X2 X3 |X1(K+1) – X1(k)|<0.001 |X2(k+1) – X2(k)|<0.001 |X3(k+1) – X3(k)|<0.001
0 0 0 0 0 0 0
1 0.25 -0.1875 -0.2344 0.25 0.1875 0.2344
2
X1= 1/4 +1/4*X2 -1/4*X3 (Ecuación 1)
X2= -1/4 +1/4*X1 -1/4*X3 (Ecuación 2)
X3= -1/4 +1/4*X1 +1/4*X2 (Ecuación 3)
X1= ¼ + ¼*(0) -1/4*(0) = 0.2500
X2= -0.1875
X3= -0.2344
Segunda iteración i=2
i X1 X2 X3 |X1(K+1) – X1(k)|<0.001 |X2(k+1) – X2(k)|<0.001 |X3(k+1) – X3(k)|<0.001
0 0 0 0 0 0 0
1 0.25 -0.1875 -0.2344 0.25 0.1875 0.2344
2 0.2617 -0.1260 -0.2161 0.0117 0.0615 0.0183
3
X1= 1/4 +1/4*X2 -1/4*X3 (Ecuación 1)
X2= -1/4 +1/4*X1 -1/4*X3 (Ecuación 2)
X3= -1/4 +1/4*X1 +1/4*X2 (Ecuación 3)
tercera iteración i=3
i X1 X2 X3 |X1(K+1) – X1(k)|<0.001 |X2(k+1) – X2(k)|<0.001 |X3(k+1) – X3(k)|<0.001
0 0 0 0 0 0 0
1 0.25 -0.1875 -0.2344 0.25 0.1875 0.2344
2 0.2617 -0.1260 -0.2161 0.0117 0.0615 0.0183
3 0.2725 -0.1278 -0.2138 0.0108 0.0018 0.0023
4 0.2715 -0.1287 -0.2143 0.001 0.0009 0,0005
5 0.2714 -0.1286 -0.2143 0.0001 0.0001 0.0000
X1= 1/4 +1/4*X2 -1/4*X3 (Ecuación 1)
X2= -1/4 +1/4*X1 -1/4*X3 (Ecuación 2)
X3= -1/4 +1/4*X1 +1/4*X2 (Ecuación 3)