Análisis numéricos
2do examen 20%
Prof.: Ing. Víctor Gómez
1. Usa el método de eliminación gaussiana con pivoteo para resolver el
siguiente sistema:
Solución. Escribimos la matriz aumentada del sistema:
( )
2 1 1 5
4 −6 0 −2
−2 7 2 9
Buscamos el pivote, para ello observamos la columna 1 y buscamos el número
mayor sin tomar en cuenta el signo:
Pivote = 4, hacemos el intercambio de la fila 1 por la fila 2:
( )
4 −6 0 −2
2 1 1 5
−2 7 2 9
Ahora empezamos a hacer cero debajo de nuestro pivote que es 4:
Queremos hacer cero el 2 de la fila 2 debemos multiplicar toda esa fila por -2 y el
resultado lo sumamos con la fila 1:
−2 F2 + F 1
−2∗2+ 4=0
−2∗1−6=−8
−2∗1+ 0=−2
−2∗5−2=−12
Ahora hacemos cero el -2 de la fila 3, para ello multiplicamos por 2 y le sumamos
el valor de la fila 3, y así hacemos con toda esa fila:
2 F3 + F 1
2∗(−2)+4=0
2∗7−6=8
2∗2+ 0=4
2∗9−2=16
Ahora mostramos la nueva matriz que nos va quedando con estos nuevos valores:
( )
4 −6 0 −2
0 −8 −2 −12
0 8 4 16
Ahora debemos buscar un nuevo pivote de la columna 2 sin tomar en cuenta la fila
1, el valor mayor es 8 (tomamos el de la fila 2 para evitar intercambiar filas), ese
será nuestro nuevo pivote.
Hacemos cero debajo del pivote para ello la fila 2 con la fila 3:
F 2+ F 3
−8+ 8=0
−2+ 4=2
−12+16=4
Podemos observar que nos queda una matriz escalonada
( )
4 −6 0 −2
0 −8 −2 −12
0 0 2 4
Por lo tanto, podemos despejar los valores x1, x2, x3
4 x1 −6 x2 =−2
−8 x 2−2 x3 =−12
2 x3 =4
Despejamos x3 de la última ecuación. Luego sustituimos x3 en la ecuación 2 y
despejamos x2 luego con x2 sustituimos en ecuación 1 y despejamos x1
4
2 x3 =4 → x 3= =2
2
−12+2∗2
−8 x 2−2∗2=−12 → x 2= =1
−8
−2+6∗1
4 x1 −6∗1=−2 → x 1= =1
4
Por lo tanto, las soluciones del sistema son:
x 1=1
x 2=1
x 3=2
2. Usa el método de Gauss-Jordan para resolver el siguiente sistema:
Solución. Escribimos la matriz aumentada del sistema:
( )
1 0 2 6
−3 4 6 30
−1 −2 3 8
Buscamos el pivote, para ello observamos la columna 1 y buscamos el número
mayor sin tomar en cuenta el signo:
Pivote = 3, hacemos el intercambio de la fila 1 por la fila 2:
( )
−3 4 6 30
1 0 2 6
−1 −2 3 8
Ahora hacemos ceros debajo del pivote que es 3:
Queremos hacer cero el 1 de la fila 2, aplicamos la fórmula:
3 F2 + F 1
3∗1−3=0
3∗0+ 4=4
3∗2+6=12
3∗6+30=4 8
Ahora hacemos cero el -1 de la fila 3, aplicamos la fórmula:
−3 F3 + F 1
−3∗(−1 )−3=0
−3∗(−2 ) +4=10
−3∗3+6=−3
−3∗8+30=6
Ahora mostramos la nueva matriz que nos va quedando con estos nuevos valores:
( )
−3 4 6 30
0 4 12 4 8
0 10 −3 6
Buscamos un nuevo pivote de la columna 2 sin tomar en cuenta la fila 1, el valor
mayor es 10, ese será nuestro nuevo pivote.
Intercambiamos la fila 2 por la 3
( )
−3 4 6 30
0 10 −3 6
0 4 12 4 8
Pivote 10
Hacemos cero debajo del pivote, aplicando la fórmula:
−2,5 F3 + F 2
−2,5∗4 +10=0
−2,5∗12−3=−33
−2,5∗4 8+6=−114
Hecho esto obtenemos la matriz nos queda así:
( )
−3 4 6 30
0 10 −3 6
0 0 −33 −114
Debemos hacer ahora cero los valores por encima de la diagonal principal
Aprovechando el pivote 10 vamos hacer cero el valor encima de el para ello
aplicamos la siguiente fórmula:
−2,5 F1 + F 2
−2,5∗(−3)+0=7,5
−2,5∗4 +10=0
−2,5∗6−3=−18
−2,5∗30+6=−6 9
( )
7,5 0 −18 −69
0 10 −3 6
0 0 −33 −114
Ahora hacemos cero lo que está encima del -33, aplicando la fórmula (para hacer
cero el -3):
−3 F 3
+F2
33
−3∗(−33)
−3=0
33
−3∗(−114)
+6=16,36364
33
Aplicamos la fórmula (para hacer cero el -18):
−18 F 3
+ F1
33
−18∗(−33)
−18=0
33
−18∗(−114)
−69=−6,81818
33
La nueva matriz nos queda así:
( )
7,5 0 0 −6,81818
0 10 0 16,36364
0 0 −33 −114
Ahora aplicamos las siguientes fórmulas para obtener unos en la diagonal
principal:
F1
7,5
7,5
=1
7,5
−6,81818
=−0,90909
7,5
F2
10
10
=1
10
16,36364
=1 , 63636
10
F3
=¿
−33
−33
=1
−33
−114
=3,45455
−33
Finalmente, la matriz nos queda así:
( )
1 0 0 −0,90909
0 1 0 1,63636
0 0 1 3,45455
Luego concluimos que las soluciones del sistema son:
x 1=−0,90909
x 2=1,63636
x 3=3,45455
Solución. Escribimos la matriz aumentada del sistema:
( )
0,7 2,7 −6 0,7 1, 6487
2 −0,8 3 −1 −2,342
−1 −1 , 5 1 , 4 3 −4,189
0 7 −1 , 56 5 15,792
Buscamos el pivote, para ello observamos la columna 1 y buscamos el número
mayor sin tomar en cuenta el signo:
Pivote = 2, hacemos el intercambio de la fila 1 por la fila 2:
( )
2 −0,8 3 −1 −2,342
0,7 2,7 −6 0,7 1,6487
−1 −1,5 1,4 3 −4,189
0 7 −1,56 5 15,792
Ahora hacemos ceros debajo del pivote que es 2:
Para hacer cero el primer elemento de la fila 1 aplicamos la fórmula siguiente:
2 F2−0,7 F 1
2∗0,7−0,7∗2=0
2∗2,7−0,7∗(−0,8)=5,96
2∗(−6)−0,7∗3=14,1
2∗0,7−0,7∗(−1)=2,1
2∗1,6487−0,7∗(−2,342)=4,9368
Para hacer cero el primer elemento de la fila 1 aplicamos la fórmula siguiente:
2 F3 + F 1
2∗ (−1 )+ 2=0
2∗(−1,5 )−0,8=−3,8
2∗ (1,4 ) +3=5,8
2∗ ( 3 )−1=5
2∗(−4,189 ) −2,3 42=−10,72
Como el primer elemento de la fila 3 ya vale cero, lo dejamos como está.
La nueva matriz es:
( )
2 −0,8 3 −1 −2,342
0 5,96 −14,1 2,1 4,9368
0 −3,8 5,8 5 −10,72
0 7 −1,56 5 15,792
El nuevo pivote es 7 de la columna 2.
Intercambiamos la fila 2 por la fila 4.
( )
2 −0,8 3 −1 −2,342
0 7 −1,56 5 15,792
0 −3,8 5,8 5 −10,72
0 5,96 −14,1 2,1 4,9368
Ahora hacemos ceros debajo del pivote que es 7:
Para hacer cero el segundo elemento de la fila 3 aplicamos la fórmula siguiente:
7 F 3+3,8 F 2
7∗(−3,8 ) +3,8∗7=0
7∗( 5,8 ) +3,8∗(−1,56)=34,672
7∗( 5 ) +3,8∗5=54
7∗(−10,72 ) +3,8∗15,792=−15,0304
Para hacer cero el segundo elemento de la fila 4 aplicamos la fórmula siguiente:
7 F 4−5,96 F 2
7∗5,96−5,96∗7=0
7∗(−14,1 )−5,96∗(−1,56 ) =−89,4024
7∗2,1−5,96∗5=−15,1
7∗4,9368−5,96∗15,792=−59,56272
La nueva matriz es:
( )
2 −0,8 3 −1 −2,342
0 7 −1,56 5 15,792
0 0 34,672 54 −15,0304
0 0 −89,4024 −15,1 −59,56272
El nuevo pivote es 89,4024 de la columna 3.
Se intercambia la columna 3 con la 4.
( )
2 −0,8 3 −1 −2,342
0 7 −1,56 5 15,792
0 0 −89,4024 −15,1 −59,56272
0 0 34,672 54 −15,0304
Ahora hacemos ceros debajo del pivote que es 89,4024:
Para hacer cero el tercer elemento de la fila 4 aplicamos la fórmula siguiente:
−89,4024 F 4 −34,672 F 3
−89,4024∗34,672−34,672∗(−89,4024 )=0
−89,4024∗54−34,672∗(−15,1 )=−4304,1824
−89,4024∗(−15,0304)−34,672∗(−59,56272 )=3408,91246
La nueva matriz es:
( )
2 −0,8 3 −1 −2,342
0 7 −1,56 5 15,792
0 0 −89,4024 −15,1 −59,56272
0 0 0 −4304,1824 3408,91246
Si dividimos la fila 4 entre -4304,1824 los cálculos se harán más sencillos en
los pasos siguientes y no habrá ningún problema con el algoritmo porque
los elemento que se usarán para hacer ceros son diferentes de cero y
además ningún elemento es cero cero aproximadamente igual a cero:
( )
2 −0,8 3 −1 −2,342
0 7 −1,56 5 15,792
0 0 −89,4024 −15,1 −59,56272
0 0 0 1 −0,79200
El 1 de la fila 4 actuará como elemento pivote:
Para hacer cero el cuarto elemento en la fila 1 aplicamos la fórmula:
F 4 + F1
0+2=2
0−0,8=−0,8
0+3=3
1−1=0
−0,792 00−2,342=−3,134 00
Para hacer cero el cuarto elemento en la fila 2 aplicamos la fórmula:
F 2−5 F 4
7−5∗0=7
−1,56−5∗0=−1,56
5−5∗1=0
15,792−5∗(−0,792 )=19,752
Para hacer cero el cuarto elemento en la fila 3 aplicamos la fórmula:
F 3+15,1 F 4
−89,4024+ 15,1∗0=−89,4024
−15,1+15,1∗1=0
−59,56272+15,1∗(−0,79200 )=−71,52192
La nueva matriz es:
( )
2 −0,8 3 0 −3,13400
0 7 −1,56 0 19,75200
0 0 −89,4024 0 −71,52192
0 0 0 1 −0,79200
Si dividimos la fila 3 entre -89,4024 la matriz nos queda así::
( )
2 −0,8 3 0 −3,13400
0 7 −1,56 0 19,75200
0 0 1 0 0,80000
0 0 0 1 −0,79200
Ahora el 1 en la columna 2 será el elemento pivote, por lo cual haremos cero los
elementos por encima de él.
Para hacer cero el tercer elemento en la fila 1 aplicamos la fórmula:
F 1−3 F 3
2−3∗0=2
−0,8−3∗0=−0,8
3−3∗1=0
0−3∗0=0
−3,13400−3∗0,80000=−5,53400
Para hacer cero el tercer elemento en la fila 2 aplicamos la fórmula:
F 2+1,56 F 3
7+1,56∗0=7
−1,56+1,56∗1=0
0+1,56∗0=0
19,75200+1,56∗0,80000=21,00000
La nueva matriz es:
( )
2 −0,8 0 0 −5,53400
0 7 0 0 21,00000
0 0 1 0 0,80000
0 0 0 1 −0,79200
Ahora dividimos la fila 2 entre 7 y obtenemos:
( )
2 −0,8 0 0 −5,53400
0 1 0 0 3
0 0 1 0 0,80000
0 0 0 1 −0,79200
Para hacer cero el segundo elemento en la fila 1 aplicamos la fórmula:
F 1+ 0,8 F 2
2+0,8∗0=2
−0,8+ 0,8∗1=0
−5,53400+0,8∗3=−3,134
La nueva matriz es:
( )
2 0 0 0 −3,134
0 1 0 0 3 , 00000
0 0 1 0 0,8 0000
0 0 0 1 −0,79 200
Ahora dividimos la fila 1 entre 2 y obtemenos:
( )
1 0 0 0 −1,56700
0 1 0 0 3,00000
0 0 1 0 0,80000
0 0 0 1 −0,792000
Por lo tanto, podemos concluimos que la solución del sistema de ecuaciones es:
x 1=−1,56700=−1,567
x 2=3,00000=3
x 3=0,80000=0,8
x 3=−0,79200=−0,792
3. Usa el método de Gauss-Seidel hasta que el error aproximado sea igual a 0
para aproximar la solución del siguiente sistema de ecuaciones:
10x1 – x2 + 2x3 = 6
-1x1 + 11x2 – 1x3 + 3x4 = 25
2x1 – 1x2 + 10x3 – 1x4 = -11
3x2 – 1x3 + 8x4 = 15
De la ecuación 1 despejemos x1, de la ecuación 2 despejemos x2, de la
ecuación 3 despejemos x3. Esto nos da el siguiente conjunto de
ecuaciones:
6+ x 2−2 x3
x 1=
10
25+ x 1+ x 3−3 x 4
x 2=
11
−11−2 x 1+ x 2 + x 4
x 3=
10
15−3 x2 + x 3
x4=
8
Iteración 1:
x 2=x 3=0
6+0−2∗0
x 1= =0,6
10
x 3=x 4 =0
25+0,6+0−3∗0
x 2= =2,32727
11
x 4 =0
−11−2∗0,6+2,32727+ 0
x 3= =−0,98727
10
15−3∗2,32727 +(−0,98727)
x4= =0,87886
8
Iteración 2:
x 2=2,32727 ; x 3=−0,98727 ;
6+2,32727−2∗(−0,98727)
x 1= =1,03018
10
x 3=−0,98727 ; x 4 =0,87886
25+1,03018+(−0,98727)−3∗(0,87886)
x 2= =2,03694
11
x 4 =0,87886
−11−2∗1,03018+2,03694+0,87886
x 3= =−1,01446
10
15−3∗1,03018+(−1,01446)
x4= =0,98434
8
Calculamos el error:
|1,03018−0,6
|∈a ,1|=
1,03018
∗100|=41,75775 %
|∈ |=| ∗100|=14 , 25374 %
2,03694−2,32727
a,2
2,03694
|∈ |=| ∗100|=2 , 68024 %
−1,01446−(−0,98727)
a ,3
−1,01446
|∈ |=| ∗100|=10 ,71581 %
0,98434−0,87886
a,4
0,98434
Iteración 3:
x 2=2,03694 ; x 3=−1,01446 ;
6+2,03694−2∗(−1,01446)
x 1= =1,00659
10
x 3=−1,01446 ; x 4 =0,98434
25+1,00659+(−1,01446)−3∗(0,98434)
x 2= =2 , 0 0356
11
x 4 =0 , 98434
−11−2∗1,00659+ 2 ,00356+ 0 , 98434
x 3= =−1 , 0 0253
10
15−3∗2,00356 +(−1 , 00253)
x4= =0 , 99835
8
Calculamos el error:
|1,00659−1,03018
|∈a ,1|= 1,00659
∗100|=2 , 34356 %
|∈ |=| ∗100|=1 ,66603 %
2,00356−2,03694
a,2
2,00356
|∈ |=| ∗100|=1, 18999 %
−1,00253−(−1,01446)
a ,3
−1,00253
|∈a , 4|= |0,99835−0,98434
0,99835
∗100|=1 , 40332 %
Iteración 4:
x 2=2,00356 ; x 3=−1 , 00253;
6+2 , 00356−2∗(−1 , 00253)
x 1= =1 ,00 086
10
x 3=−1 , 00253 ; x 4=0 , 99835
25+1 ,00086+(−1 , 00253)−3∗(0 , 99835)
x 2= =2 , 00030
11
x 4 =0 , 99835
−11−2∗1, 00086+ 2, 00030+0 , 99835
x 3= =−1,000 3 1
10
15−3∗2 , 00030+(−1 ,00031)
x4= =0,99 9 85
8
Calculamos el error:
|1 ,00086−1
|∈a ,1|=
, 00659
1 , 00086
∗100|=0,57251%
|∈ |=| ∗100|=0,16298 %
2 ,00030−2 , 00356
a, 2
2 , 00030
|∈ |=| ∗100|= 0,22193 %
−1 , 00031−(−1 , 00253)
a ,3
−1 , 00031
|∈ |=| ∗100|=0,15002 %
0 ,99985−0 , 99835
a,4
0 , 99985
Continuamos iterando hasta que el error sea cero para todas las variables.
Iteración 5:
x 2=2 , 00030; x 3=−1, 00031 ;
6+2 , 00030−2∗(−1 , 00031)
x 1= =1 , 000 09
10
x 3=−1 , 00031; x 4 =0 , 99985
25+1 ,00009+(−1 , 00031)−3∗(0 , 99985)
x 2= =2, 000 02
11
x 4 =0 , 99985
−11−2∗1, 00009+2 , 00002+0 , 99985
x 3= =−1 , 000 03
10
15−3∗2 , 00002+(−1, 00003)
x4= =0 , 999 99
8
Calculamos el error:
|1 ,00009−1
|∈a ,1|=
,00086
1 , 00009
∗100|=0,07699%
|∈ |=| ∗100|=0,01400 %
2 ,00002−2 ,00030
a, 2
2 , 00002
|∈ |=| ∗100|= 0,02800 %
−1 , 00003−(−1, 00031)
a ,3
−1 , 00003
|∈ |=| ∗100|=0,01400 %
0 ,99999−0 , 99985
a,4
0 , 99999
Iteración 6:
x 2=2 , 00002; x 3=−1 , 00003 ;
6+2 , 00002−2∗(−1 , 00003)
x 1= =1 , 0000 08
10
x 3=−1 , 00003 ; x 4=0 , 99999
25+1 ,00001+(−1 ,00003)−3∗(0 , 99999)
x 2= =2, 0000 0
11
x 4 =0 , 99999
−11−2∗1, 00001+2,0000 0+0 , 99999
x 3= =−1 , 0000 0
10
15−3∗2,0000 0+(−1,0000 0)
x4= =1,00000
8
x 1=1,000008 ≅ 1
x 2=2,00000=2
x 3=−1,00000=−1
x 4 =1,00000=1
Calculamos el error:
|1 ,00001−1,
|∈a ,1|=
00009
1 , 00001
∗100|=0,00799 % ≅ 0 %
|∈ |=| ∗100|=0,00 100 % ≅ 0 %
2 ,00000−2 , 00002
a,2
2 , 00000
|∈ |=| ∗100|=0,00 300 % ≅ 0 %
−1 , 00000−(−1 ,00003)
a ,3
−1 , 00000
|∈ |=| ∗100|=0,00 100 % ≅ 0 %
1,00000−0 , 99999
a,4
1,00000
En la iteración 6 podemos notar que el error es aproximadamente igual a
cero, por lo tanto, las soluciones del sistema por el método de Gauss-Seidel
son:
4. Calcula la matriz inversa de la siguiente matriz usando el método de Gauss-
Jordan
Solución. En una matriz, colocamos la matriz A y a su derecha agregamos la
matriz identidad
( )
1 3 3 1 0 0
1 4 3 0 1 0
1 3 4 0 0 1
Puesto que ningún elemento de la matriz dada es cero ni aproximadamente cero,
para hallar la matriz inversa por el método de Gauss-Jordan haremos el
procedimiento un poco diferente y más sencillo sin intercambiar filas:
Pivote 1 (primer elemento de la fila 1):
Para hacer cero el primer elemento de la fila 2, aplicamos la fórmula:
F 2−F1
1−1=0
4−3=1
3−3=0
0−1=−1
1−0=1
0−0=0
Para hacer cero el primer elemento de la fila 3, aplicamos la fórmula:
F 2−F1
1−1=0
3−3=0
4−3=1
0−1=−1
0−0=0
1−0=1
La nueva matriz es:
( )
1 3 3 1 0 0
0 1 0 −1 1 0
0 0 1 −1 0 1
Pivote 1 de la columna 2 (segundo elemento de la fila 2)
Para hacer cero el segundo elemento de la fila 1, aplicamos la fórmula:
F 1−3 F 2
1−3∗0=1
3−3∗1=0
3−3∗0=3
1−3∗(−1)=4
0−3∗1=−3
0−3∗0=0
La nueva matriz es:
( )
1 0 3 4 −3 0
0 1 0 −1 1 0
0 0 1 −1 0 1
Pivote 1 de la columna 3 (tercer elemento de la fila 3)
Para hacer cero el tercer elemento de la fila 1, aplicamos la fórmula:
F 1−3 F 3
1−3∗0=1
0−3∗0=0
3−3∗1=0
4−3∗(−1)=7
−3−3∗0=−3
0−3∗1=−3
La nueva matriz es:
( )
1 0 0 7 −3 −3
0 1 0 −1 1 0
0 0 1 −1 0 1
Hemos terminado porque a la izquierda está la matriz identidad y a la derecha
obtuvimos la matriz inversa.
Por lo tanto, la matriz inversa es:
( )
7 −3 −3
A−1= −1 1 0
−1 0 1
En todos los ejercicios, redondea tus resultados a cinco decimales.
Recuerden dejar evidencia con todas las muestras de cálculo
Enviar en pdf hasta la fecha y hora señalada en plataforma. Recuerden que les había
comentado que esta misma semana debo entregar notas por tanto, procuren no perder
el día del examen
