REPUBLICA BOLIVARIANA DE VENEZUELA
MINISTERIO DEL PODER POPULAR PARA LA
EDUCACION UNIVERSITARIA, CIENCIA Y
TECNOLOGIA.
INSTITUTO UNIVERSITARIO POLITECNICO
“SANTIAGO MARIÑO”
EXTENSION MERIDA
METODOS NUMERICOS (METODO
DE GAUSS-SEIDEL)
13 DE MARZO DE 2021
ANDRY BETHANIA MONTILLA MENDEZ CI 21.284.307
INGENIERIA INDUSTRIAL
1
Método de Gauss-Seidel
Este es uno de los métodos más interesantes del análisis numérico y particularmente útil ya
que nos permite encontrar la solución de un sistema de “n” ecuaciones con “n” incógnitas.
Para comenzar es preciso mencionar que es un método iterativo, es decir que debe aplicarse
recursivamente hasta encontrar una solución adecuada o con un error considerablemente pequeño.
En cada iteración obtenemos una solución posible del sistema con un error determinado, a
medida que aplicamos nuevamente el método, la solución puede ser más precisa, entonces se dice
que el sistema converge, pero si al aplicar el método reiteradas veces la solución tiene un error (ya
explicaremos como se calcula este error) cada vez mayor se dice que el sistema no converge y no
se puede resolver el sistema de ecuaciones por este método.
Teniendo el siguiente sistema de ecuaciones:
Despejamos x1 de la ecuación 1, x2 de la ecuación 2,…, xn de la ecuación n, quedando:
Desde la fórmula anterior resultan las fórmulas que se deberán ir aplicando en las diferentes
iteraciones. Para comenzar a aplicar el método debemos asignar un valor arbitrario a las variables
x2,…xn con el fin de obtener x1. Lo más conveniente en este caso es que los valores comiencen en
cero, lo cual nos facilitaría el trabajo ya que se reduce el cálculo de las primeras soluciones, entonces
de esto resulta que:
2
Ahora despejamos x2 de la ecuación 2 y reemplazamos a x1 por el valor
obtenido en la ecuación anterior. De esto nos queda:
Una vez que tenemos x2, despejamos x3 de la ecuación 3 y así
sucesivamente con las n ecuaciones, cada vez asignando el valor de las x1, x2,…
xn-1 obtenido en el paso anterior.
Cuando hemos despejado las xn, tenemos lo que se conoce como primera
solución o solución de la primera iteración:
Con los nuevos valores de x1, x2,…, xn aplicamos los mismos pasos
anteriores, pero con los nuevos valores de las xn, de esta manera conseguimos una
segunda solución:
Al tener esta segunda solución estamos en condiciones de calcular el error
que se calcula como sigue:
3
Así, repetimos el método tantas veces hasta que el error sea muy pequeño o
los suficientemente aceptable.
Ahora solo queda mencionar que para que un sistema sea convergente se
debe cumplir que la matriz de coeficientes sea diagonalmente dominante, y para
ello se debe verificar la siguiente expresión:
Si no se cumple esa condición, se puede permutar las filas de la matriz, con
el fin de poder convertirla en una diagonalmente dominante.
Este método es prácticamente idéntico al de Jacob, la única diferencia estriba
en que el método de Gauss-Seidel se acerca más rápido a la solución cuando el
método converge, debido a que una vez que se calcula la componente Xi (k 1)
i la utiliza inmediatamente en la misma iteración, esto es:
El criterio de convergencia de este método es el mismo que el de Jacobi.
4
EJERCICIOS PROPUESTOS.
1.- RESOLVER EL SIGUIENTE SISTEMA DE ECUACIONES, APLICANDO EL
METODO DE GAUSS-SEIDEL
X -7*Y + 2*Z = 19 (A)
2*X + Y + 9*Z = 13 (B)
8*X + Y – Z = 21 (C)
PASO 1. DEBEMOS ORDENAR NUESTRO SISTEMA DE MANERA QUE
OBTENGAMOS LA DIAGONAL DOMINANTE, CON ESTO ASEGURAMOS LA
CONVERGENCIA DEL SISTEMA. ESTO LO HACEMOS DE LA SIGUIENTE
MANERA:
MATRIZ DOMINANTE
8*X + Y – Z = 21
X - 7*Y + 2*Z = 19
2*X + Y + 9*Z = 13
COMO PODEMOS VER TENEMOS LA DIAGONAL DOMINANTE.
PASO 2. DESPEJAMOS LAS VARIABLES X, Y, Z DE CADA UNA DE LAS
ECUACIONES:
X= (21 - Y + Z)/8
Y= (-19 + X + 2*Z)/7
Z= (13 - 2*X - Y)/9
PASO 3. DAMOS VALORES INICIALES A LAS VARIABLES, CON EL FIN DE
ENCONTRAR ALGUN VALOR INCIAL PARA COMENZAR LA ITERACION.
PRIMERA ITERACION
DANDO VALORES INICIALES PARA
Y=0
Z=0
ENCUENTRO VALOR INICIAL PARA X
X= (21 - (0) + (0))/8
X= 2,62500
CON VALORES PARA
X= 2,62500
Z= 0
ENCUENTRO VALOR INICIAL PARA Y
5
Y= (-19 + (2,62500) + 2*(0))/7
Y= -2,33929
CON VALORES PARA
X= 2,62500
Y= -2,33929
ENCUENTRO VALOR INICIAL PARA Z
Z= (13 - 2*(2,62500) - (-2,33929))/9
Z= 1,12103
PASO 4. CON LOS PRIMEROS VALORES OBTENIDOS DE LAS VARIABLES,
HALLO NUEVOS VALORES PARA LAS VARIABLES.
SEGUNDA ITERACION
CON MIS NUEVOS VALORES PARA
Y= -2,33929
Z= 1,12103
ENCUENTRO UN NUEVO VALOR PARA X
X= (21 - (-2,33929) + (1,12103))/8
X= 3,05754
CON MIS NUEVOS VALORES PARA
X= 3,05754
Z= 1,12103
ENCUENTRO UN NUEVO VALOR PARA Y
Y= (-19 + (3,05754) + 2*(1,12103))/7
Y= -1,95720
CON MIS NUEVOS VALORES PARA
X= 3,05754
Y= -1,95720
ENCUENTRO VALOR INICIAL PARA Z
Z= (13 - 2*(3,05754) - (-1,95720))/9
Z= 0,98246
PASO 5. CON LOS VALORES DE X, Y, Z DE LA SEGUNDA ITERACION, BUSCO EL
ERROR CORRESPONDIENTE A CADA UNO DE ELLOS.
6
PARA EL ERROR DE X, USAMOS:
E Xn= ((|Xi+1 - Xi)/Xi+1) |*100
E Xn= ((|(3,05754) - (2,62500))/3,05754)
|*100%
E Xn= 14,14667%
PARA EL ERROR DE Y, USAMOS:
E Yn= ((|Yi+1 - Yi|)/Yi+1)*100
E Yn= ((|(-1,95720) - (-2,33929))/(-1,95720))
|*100%
E Yn= 19,52206
PARA EL ERROR DE Z, USAMOS:
E Zn= ((|Zi+1 - Zi)/Zi+1)|*100
E Zn= ((|( 0,98246) - (1,12103))/
0,98246)|*100%
E Zn= 14,10439%
PASO 6. COMO VEMOS NUESTRO ERROR SOBREPASA EL ERROR DE
PRECISION, EL CUAL DEBE SER MENOR O IGUAL A 0,00001%.
CONTINUAREMOS REALIZANDO ITERACIONES DE LA MISMA MANERA
EXPUESTA ANTERIORMENTE HASTA LLEGAR AL ERROR DESEADO. LOS
RESULTADOS SE MOSTRARAN EN LA SIGUIENTE TABLA.
ITERACION VALOR X VALOR Y VALOR Z % EXn % EYn % EZn
1 2,62500 -2,33929 1,12103 100,00000 100,00000 100,00000
2 3,05754 -1,95720 0,98246 14,14667 19,52206 14,10439
3 2,99246 -2,00609 1,00235 2,17489 2,43706 1,98483
4 3,00106 -1,99918 0,99967 0,28650 0,34577 0,26796
5 2,99986 -2,00011 1,00004 0,03997 0,04683 0,03705
6 3,00002 -1,99998 0,99999 0,00545 0,00646 0,00507
7 3,00000 -2,00000 1,00000 0,00075 0,00088 0,00070
8 3,00000 -2,00000 1,00000 0,00010 0,00012 0,00010
9 3,00000 -2,00000 1,00000 0,00001 0,00002 0,00001
10 3,00000 -2,00000 1,00000 0,00000 0,00000 0,00000
CON UN ERROR DEL 0,00000 PARA CADA UNA DE LAS VARIABLES
OBTENEMOS LOS SIGUIENTES VALORES PARA CADA UNA DE ELLAS:
7
X= 3
Y= -2
Z= 1
2.- RESOLVER EL SIGUIENTE SISTEMA DE ECUACIONES, APLICANDO EL
METODO DE GAUSS-SEIDEL
- 2*X + 10*Y + 3*Z + 2*W = 35 (A)
- X + 3*Y - 2*Z +9*W = 35 (B)
- 8*X + Y + Z + 5*W = 17 (C)
3*X – 4*Y + 7*Z + W = 20 (D)
PASO 1. DEBEMOS ORDENAR NUESTRO SISTEMA DE MANERA QUE
OBTENGAMOS LA DIAGONAL DOMINANTE, CON ESTO ASEGURAMOS LA
CONVERGENCIA DEL SISTEMA. ESTO LO HACEMOS DE LA SIGUIENTE
MANERA:
MATRIZ DOMINANTE
- 8*X + Y + Z + 5*W = 17
- 2*X + 10*Y + 3*Z + 2*W = 35
3*X – 4*Y + 7*Z + W = 20
- X + 3*Y - 2*Z +9*W = 35
COMO PODEMOS VER TENEMOS LA DIAGONAL DOMINANTE.
PASO 2. DESPEJAMOS LAS VARIABLES X, Y, Z DE CADA UNA DE LAS
ECUACIONES:
X= (-17 + Y + Z + 5*W)/8
Y= (35 + 2*X - 3*Z - 2*W)/10
Z= (20 - 3*X + 4*Y - W)/7
W= (35 + X - 3*Y + 2*Z)/9
PASO 3. DAMOS VALORES INICIALES A LAS VARIABLES, CON EL FIN DE
ENCONTRAR ALGUN VALOR INCIAL PARA COMENZAR LA ITERACION.
8
PRIMERA ITERACION
DANDO VALORES INICIALES PARA
Y=0
Z=0
W= 0
ENCUENTRO VALOR INICIAL PARA X
X= (-17 + (0) + (0) + 5*(0))/8
X= -2,12500
CON VALORES PARA
X= -2,12500
Z=0
W= 0
ENCUENTRO VALOR INICIAL PARA Y
Y= (35 + 2*(-2,12500) - 3*(0) - 2*(0))/10
Y= 3,07500
CON VALORES PARA
X= -2,12500
Y= 3,07500
W= 0
ENCUENTRO VALOR INICIAL PARA Z
Z= (20 - 3*(-2,12500) + 4*(3,07500) - (0))/7
Z= 5,52500
CON VALORES PARA
X= -2,12500
Y= 3,07500
Z= 5,52500
ENCUENTRO VALOR INICIAL PARA W
W= (35 + (-2,12500) - 3*(3,07500) +
2*(5,52500))/9
W= 3,85556
PASO 4. CON LOS PRIMEROS VALORES OBTENIDOS DE LAS VARIABLES,
HALLO NUEVOS VALORES PARA LAS VARIABLES.
SEGUNDA ITERACION
CON MIS NUEVOS VALORES PARA
Y= 3,07500
9
Z= 5,52500
W= 3,85556
ENCUENTRO UN NUEVO VALOR PARA X
X= (-17 + (3,07500) + (5,52500) +
5*(3,85556))/8
X= 1,35973
CON MIS NUEVOS VALORES PARA
X= 1,35973
Z= 5,52500
W= 3,85556
ENCUENTRO UN NUEVO VALOR PARA Y
Y= (35 + 2*(1,35973) - 3*(5,52500) -
2*(3,85556))/10
Y= 1,34333
CON MIS NUEVOS VALORES PARA
X= 1,35973
Y= 1,34333
W= 3,85556
ENCUENTRO UN NUEVO VALOR PARA Z
Z= (20 - 3*(1,35973) + 4*(1,343330) -
(3,85556))/7
Z= 2,49122
CON MIS NUEVOS VALORES PARA
X= 1,35973
Y= 1,34333
Z= 2,49122
ENCUENTRO UN NUEVO VALOR PARA W
W= (35 + (1,35973) - 3*1,34333) + 2*(
2,49122))/9
W= 4,14580
PASO 5. CON LOS VALORES DE X, Y, Z DE LA SEGUNDA ITERACION, BUSCO EL
ERROR CORRESPONDIENTE A CADA UNO DE ELLOS.
PARA EL ERROR DE X, USAMOS:
E Xn= (|(Xi+1 - Xi)/(Xi+1)|)*100
E Xn= ((|(1,35973) - ( -
2,12500)/(1,35973)|*100%
1
0
E Xn= 256,28103%
PARA EL ERROR DE Y, USAMOS:
E Yn= (|(Yi+1 - Yi)/(Yi+1)|)*100
E Yn= ((|(1,34333) - (
3,07500)/(1,34333)|*100%
E Yn= 128,90819%
PARA EL ERROR DE Z, USAMOS:
E Zn= (|(Zi+1 - Zi)/(Zi+1)|)*100
E Zn= ((|(2,49122) - (
5,52500)/(2,49122)|*100%
E Zn= 121,77798%
PARA EL ERROR DE W, USAMOS:
E Wn= (|(Wi+1 - Wi)/(Wi+1)|)*100
E Wn= ((|(4,14580) - (
5,52500)/(3,85556)|*100%
E Wn= 7,00088%
PASO 6. COMO VEMOS NUESTRO ERROR SOBREPASA EL ERROR DE
PRECISION, EL CUAL DEBE SER MENOR O IGUAL A 0,00001%.
CONTINUAREMOS REALIZANDO ITERACIONES DE LA MISMA MANERA
EXPUESTA ANTERIORMENTE HASTA LLEGAR AL ERROR DESEADO. LOS
RESULTADOS SE MOSTRARAN EN LA SIGUIENTE TABLA.
ITERACION VALOR X VALOR Y VALOR Z VALOR W % EXn % EYn % EZn % Ewn
1 -2,12500 3,07500 5,52500 3,85556 100,00000 100,00000 100,00000 100,00000
2 1,35973 1,34333 2,49122 4,14580 256,28103 128,90876 121,77889 7,00082
3 0,94544 2,11256 3,06687 3,97128 43,81924 36,41230 18,77006 4,39460
4 1,00448 1,98658 2,99451 4,00375 5,87716 6,34181 2,41643 0,81110
5 0,99998 2,00089 2,99998 3,99970 0,44968 0,71538 0,18225 0,10141
6 0,99992 2,00005 3,00011 4,00000 0,00623 0,04207 0,00418 0,00754
7 1,00002 1,99997 2,99998 4,00001 0,00989 0,00390 0,00434 0,00020
8 1,00000 2,00001 3,00000 4,00000 0,00209 0,00166 0,00089 0,00019
9 1,00000 2,00000 3,00000 4,00000 0,00028 0,00030 0,00012 0,00004
10 1,00000 2,00000 3,00000 4,00000 0,00002 0,00004 0,00001 0,00001
1,00000 2,00000 3,00000 4,00000 0,00000 0,00000 0,00000 0,00000
CON UN ERROR DEL 0,00000 PARA CADA UNA DE LAS VARIABLES
OBTENEMOS LOS SIGUIENTES VALORES PARA CADA UNA DE ELLAS:
1
1
X= 1
Y= 2
Z= 3
W= 4
1
2
1
3
