Taller 2

Ejercicio 1.
Solución de Sistemas de Ecuaciones Lineales
1.19x1 + 2.11x2 − 100x3 + x4 = 1.12
14.2x1 − 0.122x2 + 12.2x3 − x4 = 3.44
100x2 − 99.9x3 + x4 = 2.15
15.3x1 + 0.110x2 − 13.1x3 − x4 = 4.16
1. Determinar si el SEL tiene, o no, solución. ¿Si la tiene, es única? Realice una breve explicación basándose
en la teoría revisada.
2. Resolver el SEL por cada uno de los esquemas de:
a. Eliminación Gaussiana simple
b. Eliminación de Gauss - Jordan.
c. Gauss – Seidel
d. Jacobi
e. S.O.R.
3. Debe establecer los criterios de convergencia en cada método. Itere la solución hasta que ésta converja a
una tolerancia de ‖A‖ ≤ 10−6 , donde ‖A‖ es una norma matricial para la matriz A (escoja la que considere más
adecuada). Valide que el vector x es solución del SEL (mediante la operación Ax = b).
4. Realizar una gráfica de la forma como va convergiendo la solución (Número de iteraciones vs norma del
error). Realice una breve explicación, sustentándose en la teoría revisada, acerca de los resultados. ¿Cuál
considera que es el mejor método para este SEL en particular? ¿por qué?
Ejercicio 2. Interpolación Numérica
La siguiente tabla presenta la velocidad de caída de un paracaidista en función del tiempo:
Tiempo (s) 1 3 5 7 20
Velocidad (cm/s) 800 2310 3090 3949 800
1. Estime el valor de la velocidad en el instante t = 10 s utilizando los métodos solicitados en el Ejercicio 1 de este tema por
medio de un polinomio de interpolación de grado adecuado.
El ajuste polinomial se realizará usando un polinomio cuadrático, y para esto se debe formular el sistema
Cuya solución serán los valores a 0 , a1 y a2, que son los coeficientes del polinomio v=f ( s ) =a2 s 2+ a1 s+ a0
Ahora, para armar el sistema, usamos la tabla:
s v s^2 s^3 s^4 s*v s^2*v

1 800 1 1 1 800 800
3 2310 9 27 81 6930 20790
5 3090 25 125 625 15450 77250
7 3949 49 343 2401 27643 193501
20 800 400 8000 160000 16000 320000
Total 36 10949 484 8496 163108 66823 612341
Como s es la variable independiente tomará los valores de x , mientras que v depende de s, por eso se le asigna y , y con
esto se arma el sistema:
5 a0 + 36 a1 + 484 a2= 10949

36 a 0+ 484 a1+ 8496 a 2= 66823
484 a0+ 8496 a 1+ 163108 a2= 612341
Para solucionar el sistema, se usarán los métodos del punto 1.
a. Eliminación Gaussiana simple
5 a0 + 36 a1 + 484 a2= 10949

36 a 0+ 484 a1+ 8496 a 2= 66823 F2 se cambia F2-F1*36/5
484 a0+ 8496 a 1+ 163108 a2= 612341 F3 se cambia F3-F1*484/5
5 a0 + 36 a1 + 484 a2= 10949
224,8 a1+ 5011,2 a2 = -12009,8
5011,2 a1+ 116256,8 a 2= -447522,2 F3 se cambia F3-F2*5011,2/224,8
5 a0 + 36 a1 + 484 a2= 10949

224,8 a1+ 5011,2 a2 = -12009,8
4548,0569395017 a 2= -179801,960854
Así
−179801,960854
a 2= =−39,55097800945,
4548,0569395017
−12009,8−5011,2∗(−39,55097800945)
a 1= =828,238 4341637
224,8
10949−36∗828,2384341637−484∗(−39,55097800945 )
a 0= =55,03466
5
Por lo tanto, el polinomio cuadrático es: v=f ( s ) =a2 s 2+ a1 s+ a0

v=f ( s ) =−39,55097800945 s2 +828,2384341637 s+55,03466
y como piden la velocidad en s=10, entonces
v=f ( 10 ) =−39,55097800945¿ 102 +828,2384341637∗10+55,03466=4382,32866
la velocidad al cabo de 10 segundos es 4382,32866 cm/ seg
b. Eliminación de Gauss - Jordan.
Retomamos el último paso de eliminación gaussiana simple
5 a0 + 36 a1 + 484 a2= 1094 9 F1 se cambia F1-F3*484

224,8 a1+ 5011,2 a2 = −12009,8 F2 se cambia F2-F3*5011,2
a 2= −39,5509780094 5
5 a0 + 36 a1 + 3 0 0 9 1,673356573 F1 se cambia F1-
36*F2/224,8
224,8 a1+ 186188,06100095
a 2= −39,5509780094 5
5 a0 + 275,0897266806
224,8 a1+ 186188,06100095
a 2= −39,5509780094 5
186188,06100095 275,0897266806
Luego a 1= =828,2384341637 , a 0= =55,03466
224,8 5
c. Gauss – Seidel
Tomamos el sistema
5 a0 + 36 a1 + 484 a2= 10949

36 a 0+ 484 a1+ 8496 a 2= 66823
484 a0+ 8496 a 1+ 163108 a2= 612341
Y despejamos a 0 en la primera ecuación, a 1 en la segunda y a 2 en la tercera del sistema obtenido en el primer punto
10949−36 a1−484 a 2 66823−36 a0 −8496 a2 612341−484 a0−8 4 96 a2
a 0= , a 1= , a 2=
5 484 163108
Primera iteración
Paso 1. Para a 0, toma a 1=a2=0
Paso 2. Para a 1tomar a 0 del paso 1 y a 2=0
Paso 3. Para a 2 tomar a 0 del paso 1 y a 1 del paso 2
Segunda iteración
Paso 4. Para a 0, toma a 1 del paso 2 y a 2 del paso 3

Paso 5. Para a 1tomar a 0 del paso 4 y a 2 del paso 3
Paso 6. Para a 2 tomar a 0 del paso 4 y a 1 del paso 5
Y asi sucesivamente, los resultados se resumen en la tabla:
iteración a_0 a_1 a_2 a_1 a_0 a_2 a_2 a_0 a_1
1 2189,8 0 0 -24,813636 2189,8 0 -1,4512197 2189,8 -24,813636
2 2508,93625 -24,813636 -1,4512197 -23,07674 2508,93625 -1,4512197 -2,4886833 2508,93625 -23,07674
3 2596,85708 -23,07674 -2,4886833 -11,404961 2596,85708 -2,4886833 -3,3575378 2596,85708 -11,404961
4 2596,92538 -11,404961 -3,3575378 3,84158565 2596,92538 -3,3575378 -4,1519055 2596,92538 3,84158565
5 2564,04503 3,84158565 -4,1519055 20,2313381 2564,04503 -4,1519055 -4,9080501 2564,04503 20,2313381
6 2519,23362 20,2313381 -4,9080501 36,8375697 2519,23362 -4,9080501 -5,6400671 2519,23362 36,8375697
7 2470,52799 36,8375697 -5,6400671 53,3099221 2470,52799 -5,6400671 -6,353555 2470,52799 53,3099221
8 2420,99268 53,3099221 -6,353555 69,5187326 2420,99268 -6,353555 -7,0508535 2420,99268 69,5187326
9 2371,78774 69,5187326 -7,0508535 85,4187861 2371,78774 -7,0508535 -7,7330497 2371,78774 85,4187861
10 2323,34395 85,4187861 -7,7330497 100,997124 2323,34395 -7,7330497 -8,400747 2323,34395 100,997124
Continuando con las iteraciones
380 56,9941031 827,566351 -39,521403 827,572426 56,9941031 -39,521403 -39,521663 56,9941031 827,572426
d. Jacobi
e. S.O.R.
2. Realice los cálculos solicitados por cada uno de los siguientes métodos:
• Interpolación de Lagrange
• Interpolación de Diferencias Divididas de Newton
• Polinomio de interpolación de Diferencias Finitas de Newton
3. Cada estudiante deberá realizar un análisis comparativo de los resultados obtenidos apoyándose en la teoría revisada.
Emplee una gráfica como ayuda para el análisis de resultados. Ejercicio sin un análisis adecuado, apoyado en la teoría y
debidamente referenciado, no será tenido en cuenta y la nota de esta etapa será de cero.cero (0.0)
Observaciones Generales (todos los estudiantes):
No se aceptarán resultados presentados en forma fraccionaria (por ejemplo, 3/7), debe usarse la notación decimal de cada
resultado (por ejemplo, 0.428571). Deberá emplear al menos seis (6) cifras decimales en cada cálculo, si emplea menos, se
considerará malo el ejercicio.
Haga uso de alguna herramienta de cálculo (por ejemplo, Excel, Geogebra o lenguaje de programación) anexando los
archivos en donde se muestre el desarrollo completo del ejercicio y la forma como se desarrollaron cada una de las
operaciones (los cálculos). No se aceptarán ejercicios resueltos empleando calculadores on-line o que sea solo colocar los
parámetros y obtener los resultados sin tener idea de lo que se hizo. Deben implementar cada uno de los métodos de
forma individual y basándose en la teoría vista. Si no anexa los archivos de apoyo, o emplea alguna calculadora on-line que
resuelva los diferentes ejercicios, la actividad no será tenido en cuenta y su nota será de cero.cero (0.0)

Taller 2

Cargado por

Copyright:

Formatos disponibles

Taller 2

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Taller 2

Cargado por

Copyright:

Formatos disponibles

Ejercicio 1.

Solución de Sistemas de Ecuaciones Lineales

1.19x1 + 2.11x2 − 100x3 + x4 = 1.12

14.2x1 − 0.122x2 + 12.2x3 − x4 = 3.44

100x2 − 99.9x3 + x4 = 2.15

15.3x1 + 0.110x2 − 13.1x3 − x4 = 4.16

2. Resolver el SEL por cada uno de los esquemas de:

a. Eliminación Gaussiana simple

b. Eliminación de Gauss - Jordan.

Ejercicio 2. Interpolación Numérica

La siguiente tabla presenta la velocidad de caída de un paracaidista en función del tiempo:

Ahora, para armar el sistema, usamos la tabla:

s v s^2 s^3 s^4 s*v s^2*v

5 a0 + 36 a1 + 484 a2= 10949

a. Eliminación Gaussiana simple

5 a0 + 36 a1 + 484 a2= 10949

5 a0 + 36 a1 + 484 a2= 10949

Por lo tanto, el polinomio cuadrático es: v=f ( s ) =a2 s 2+ a1 s+ a0

y como piden la velocidad en s=10, entonces

v=f ( 10 ) =−39,55097800945¿ 102 +828,2384341637∗10+55,03466=4382,32866

la velocidad al cabo de 10 segundos es 4382,32866 cm/ seg

b. Eliminación de Gauss - Jordan.

Retomamos el último paso de eliminación gaussiana simple

5 a0 + 36 a1 + 484 a2= 1094 9 F1 se cambia F1-F3*484

5 a0 + 36 a1 + 484 a2= 10949

Paso 1. Para a 0, toma a 1=a2=0

Paso 2. Para a 1tomar a 0 del paso 1 y a 2=0

Paso 3. Para a 2 tomar a 0 del paso 1 y a 1 del paso 2

Paso 4. Para a 0, toma a 1 del paso 2 y a 2 del paso 3

Paso 6. Para a 2 tomar a 0 del paso 4 y a 1 del paso 5

Y asi sucesivamente, los resultados se resumen en la tabla:

• Interpolación de Diferencias Divididas de Newton

• Polinomio de interpolación de Diferencias Finitas de Newton

Observaciones Generales (todos los estudiantes):

También podría gustarte

s v s^2 s^3 s^4 sv s^2v