Metodos Numericos ICC-1027
Metodos Numericos ICC-1027
Metodos Numericos ICC-1027
DE ACAYUCAN
MANUAL DE PRCTICAS DE
METODOS NUMERICOS
CLAVE: ICC 1027
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 2
MANUAL DE PRCTICAS DE
METODOS NUMERICOS
CLAVE: ICC 1027
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 3
Contenido
Objetivo .......................................................................................................................................................................... 4
Introduccin ................................................................................................................................................................. 4
Competencia a desarrollar ...................................................................................................................................... 7
Prctica 1: Tipos de errores.................................................................................................................................... 8
Prctica 2: Software de cmputo numrico .................................................................................................. 10
Prctica 3: Mtodo de biseccin y regla falsa ............................................................................................... 12
Prctica 4: Mtodo de Newton Raphson y secante ..................................................................................... 15
Prctica 5: Aplicacin de la solucin de ecuaciones no lineales ............................................................ 17
Prctica 6: Interpolacin lineal .......................................................................................................................... 19
Prctica 7: Interpolacin de Lagrange ............................................................................................................. 22
Prctica 8: Aplicaciones de interpolacin ...................................................................................................... 25
Prctica 9: Integracin numrica ...................................................................................................................... 27
Prctica 10: Aplicacin de la integracin numrica ................................................................................... 30
Prctica 11: Mtodo de Gauss Jordn .............................................................................................................. 33
Prctica 12: Mtodo de Gauss Seidel................................................................................................................ 36
Prctica 13: Mtodos de Jacobi........................................................................................................................... 38
Prctica 14: Aplicaciones de los sistemas de ecuaciones lineales ........................................................ 40
Prctica 15: Mtodo de Newton Raphson para sistemas de ecuaciones no lineales ..................... 42
Prctica 16: Aplicaciones de los sistemas de ecuaciones no lineales .................................................. 45
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 4
Objetivo
Seleccionar y aplicar mtodos (algoritmos) numricos para resolver
problemas matemticos referentes a diferentes reas de a ingeniera,
probabilidad, anlisis estadstico, entre otras, de acuerdo al tipo de
problema en particular.
Evaluar los temas vistos en cursos tradicionales de clculo, lgebra
lineal, ecuaciones diferenciales, etc. desde el punto de vista numrico,
concretando en el anlisis de una serie de mtodos o algoritmos, y su
aplicacin mediante el uso de computadoras y el software apropiado.
Introduccin
Caracterizacin
de la
asignatura
Esta
asignatura
aporta
al
perfil
del
Ingeniero
en
nuevas
referente a las
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 5
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 6
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 7
Competencia a desarrollar
Competencias especficas:
Resolver problemas relacionados con la ingeniera de procesos mediante
la aplicacin de algoritmos numricos y el uso de computadoras digitales.
Competencias genricas
Competencias
instrumentales
Solucin de problemas
Toma de decisiones.
Competencias
interpersonales
Trabajo en equipo
Habilidades interpersonales
Competencias
sistmicas
Habilidades de investigacin
Capacidad de aprender
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 8
Unidad 1
Subtemas
Introduccin
E p p*
Error por redondeo. Este error es el resultado de representar aproximadamente
nmeros exactos. Es decir, se debe a la omisin de algunas de las cifras significativas de
algn valor especfico. Los errores de redondeo se deben a que las computadoras solo
guardan un nmero finito de cifras significativas durante un clculo. Las computadoras
realizan esta funcin de maneras diferentes. Por ejemplo, si solo se guardan siete cifras
significativas, la computadora puede almacenar y usar como = 3.141592, omitiendo
los trminos restantes y generando un error de redondeo.
Errores de truncamiento. Los errores de truncamiento son aquellos que resultan al usar
una sucesin finita o infinita de pasos en el cual se realizan clculos para producir un
resultado exacto, se trunca prematuramente despus de un cierto nmero de pasos.
Truncar la siguiente cifra hasta centsimos, o hasta que sean dos las cifras significativas:
7 2 . 645751311
2.64
Como podemos ver, en este tipo de error, lo que se hace es omitir algunas de las cifras de
una cantidad, debido a que esta contiene muchos decimales, entonces se trunca o corta el
nmero, por lo que tambin cae en un error.
Los errores de truncamiento son aquellos que resultan al usar una aproximacin en lugar
de un procedimiento matemtico exacto. Estos errores se regresan a la formulacin
matemtica usada ampliamente en los mtodos numricos para expresar funciones en
forma polinomial.
Herramienta de
apoyo
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 9
Desarrollo de la practica
Instruccin I
Problema 1
Aplique el mtodo de redondeo para las siguientes cantidades segn se indique el nmero
de cifras significativas y a la vez escriba el nmero de digito que tenga el nmero inicial.
Numero
Cifras significativas
5.6723
10.406
7.3500
88.21650
1.25001
Redondeo
Dgitos significativos
Conclusin
Referencia
Chapra , S., & Canale, R. . Mtodos Nmericos para ingenieros (Quinta ed.). Mexico: Mc Graw
bibliogrfica
Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 10
Realizar una lista del software de cmputo numrico tanto comercial como libres.
Unidad 1
Subtemas
Introduccin
Herramienta de
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
Respuesta:
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 11
Instruccin II
Actividad 1
Realice una investigacin acerca de los lenguajes de programacin y con esa informacin
elabore una lnea de tiempo.
Respuesta
Conclusin
Referencia
bibliogrfica
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 12
Unidad 2
Subtemas
Introduccin
El mtodo de biseccin
El mtodo de biseccin conocido tambin como de corte binario, de particin en dos
intervalos iguales o mtodo Bolzano, es un mtodo de bsqueda incremental en el que el
intervalo se divide siempre en dos. La posicin de la raz se determina situndola en el
punto medio del subintervalos dentro del cual ocurre un cambio de signo. El proceso se
repite hasta obtener una mejor aproximacin.
xr
x1 xu
2
f ( xu ) .
xr xu
Herramienta de
f ( x1 ) y
f xu x1 xu
f x1 f xu
apoyo
Desarrollo de la practica
Instruccin I
INGENIERA
CIVIL
Problema 1
Cdigo:
Revisin: 0
Pgina 13
P 2 V b RT
V
Dnde:
He
0.03412
0.02370
Realice los clculos necesarios para resolver esta ecuacin usando los mtodos indicados:
a) Mtodo de biseccin con
V1 0.8v , Vu 1.2v ,
V1 0.8v , Vu 1.2v
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 14
Respuesta (b)
Conclusin
Referencia
Chapra , S., & Canale, R. (2007). Mtodos Nmericos para ingenieros (Quinta ed.). Mexico:
bibliogrfica
Mc Graw Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 15
Unidad 2
Subtemas
Introduccin
xi 1 xi
f xi
f xi
xi 1 xi
Herramienta de
xi xi1 f xi
f xi f xi 1
xi 1 xi
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
x0 0 para
xi 1 xi 0.001
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 16
Respuesta
Problema 2
f ( x) x 2 1 tan( x) .
EA 1%
Respuesta
Conclusin
Referencia
Chapra , S., & Canale, R. (2007). Mtodos Nmericos para ingenieros (Quinta ed.). Mexico:
bibliogrfica
Mc Graw Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 17
Unidad 2
Subtemas
Introduccin
Para hacer uso del clculo de races de ecuaciones no lineales se puede hacer uso de los
mtodos de biseccin, falsa posicin o regla falsa, Newton Raphson y mtodo de la
Secante.
Herramienta de
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
Para obtener la temperatura de burbuja de una mezcla de CCl4 y CF4 en equilibrio con su
vapor, se lleg a la ecuacin:
6.898T1221.8
6.195T376.71
227.4
241.2
f (T ) 0.75 10
0.25
10
760
105, 100
Respuesta
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 18
Conclusin
Referencia
Chapra , S., & Canale, R. . Mtodos Nmericos para ingenieros (Quinta ed.). Mexico: Mc Graw
bibliogrfica
Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 19
Unidad 3
Interpolacin
Subtemas
Introduccin
f 1( x) f ( x0 ) f ( x1 ) f ( x0 )
x x0
x1 x0
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 20
f 1( x) f ( x0 )
f ( x1 ) f ( x0 )
( x x0 )
x1 x0
f 1( x) f ( x0 ) f x0 , x1 ( x x0 )
La cual es una frmula de interpolacin lineal. La notacin f 1(x) indica que se trata de un
polinomio de interpolacin de primer orden. Ntese que adems de representar la
pendiente de la linera que conecta los dos puntos, el termino
f ( x1 ) f ( x0 )
x1 x0
Es una aproximacin de diferencias divididas finitas a la primera derivada. En general,
entre ms pequeo sea el intervalo entre dos puntos, ms exacta ser la aproximacin.
Interpolacin cuadrtica. Una estrategia que mejora la aproximacin es la introducir
cierta curvatura en la lnea que conecta a los puntos. Si se dispone de tres datos, lo anterior
se puede llevar a cabo con un polinomio de segundo orden (llamado tambin polinomio
cuadrtico o parbola). Una manera conveniente para este caso es:
f 2( x) f ( x0 ) f x0 , x1 ( x x0 ) f x0 , x1 , x2 ( x x0 )( x x1 )
Herramienta de
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
Calcular ln 9.2 a partir de ln 9.0, ln 9.5, ln 11.0 por interpolacin cuadrtica de Lagrange
y determinar el error a partir de ln 9.2 = 2.2192 (4D)
Respuesta
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 21
la siguiente tabla:
Punto
T (K)
280
650
1000
1200
1500
1700
Cp
32.7
45.4
52.15
53.7
52.9
50.3
Conclusin
Referencia
Chapra , S., & Canale, R. (2007). Mtodos Nmericos para ingenieros (Quinta ed.). Mexico:
bibliogrfica
Mc Graw Hill
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 22
Unidad 3
Interpolacin
Subtemas
Introduccin
f 1 ( x)
x x1
x x1
f x0
f x1
x0 x1
x1 x0
f 2 x
Herramienta de
x x1 x x2 f x x x0 x x2 f x x x0 x x1 f x
x0 x1 x0 x2 0 x1 x0 x1 x2 1 x2 x0 x2 x1 2
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
Calcular ln 9.2 a partir de ln 9.0, ln 9.5, ln 11.0 por interpolacin cuadrtica de Lagrange
y determinar el error a partir de ln 9.2 = 2.2192
Respuesta
Problema 2
Use la interpolacin de Lagrange, para las densidades de las soluciones acuosas del cido
sulfrico varan con la temperatura y la concentracin de acuerdo con la tabla :
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 23
T ( C)
C (%)
10
30
60
100
1.0344
1.0281
1.0140
0.9888
20
1.1453
1.1335
1.1153
1.0885
40
1.3103
1.2953
1.2732
1.2446
70
1.6923
1.6014
1.5753
1.5417
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 24
Respuesta (b)
Respuesta (c)
Conclusin
Referencia
Chapra , S., & Canale, R. (2007). Mtodos Nmericos para ingenieros (Quinta ed.). Mexico:
bibliogrfica
Mc Graw Hill
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 25
Unidad 3
Interpolacin
Subtemas
Introduccin
Herramienta de
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
Respuesta
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 26
Problema 2
Respuesta
Conclusin
Referencia
Nieves, A., & Federico, C. (2010). Metodos Numericos para ingenieros (Tercera ed.). Mxico:
bibliogrfica
Mc Graw Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 27
Calcular la integracin numrica usando el mtodo del trapecio para los problemas reales
Unidad 4
Integracin numrica
Subtemas
Introduccin
/ f ( x)dx f1 ( x)dx
b
El resultado de la integracin es :
I b a
f (a) f (b)
2
los segmentos.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 28
n 1
h
f x0 2 f xi f xn
2
i 1
I (b a)
Herramienta de
f ( x0 ) 2 f ( xi ) f ( xn )
i 1
2n
apoyo
Desarrollo de la practica
Instruccin I
Resolver los siguientes ejercicios usando el mtodo del trapecio. Utilice Excel para
resolverlos.
Problema 1
x 5Cos(t ) Cos(5t )
y 5Sen(t ) Sen(5t )
b
usando :
s
a
dy dx
dt dt dt
Hallar la distancia recorrida por el punto en una vuelta completa alrededor de crculo
mayor. Use el mtodo de trapecio de segmento mltiple con n = 10
Respuesta
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 29
Conclusin
Referencia
bibliogrfica
Chapra , S., & Canale, R. (2007). Mtodos Nmericos para ingenieros (Quinta ed.). Mexico:
Mc Graw Hill.
Larson , R., & Hostetler, R. (2000). Calculo (Sexta ed., Vol. 2). Madrid: Mc Graw Hill.
Nieves, A., & Federico, C. (2010). Metodos Numericos para ingenieros (Tercera ed.). Mxico:
Mc Graw Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 30
Calcular la integracin numrica para problemas de aplicacin usando los mtodos vistos
Unidad 4
Integracin numrica
Subtemas
Introduccin
Herramienta de
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
H (en
H mcT
Donde c esta en cal/gC, masa (g) y
podra aumentar podra aumentar con la temperatura de acuerdos con una relacin tal y
como:
c (T )
T2
T1
c(T )dT
T2 T1
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 31
T2 c(T )dT
T
T m T2 c(T )dT
H mcT m 1
T1
T2 T1
Encuentre el centro de masa de una lmina cuya forma se encuentra en la figura adjunta,
suponiendo que la densidad en un punto p(x,y) de la lmina est dada por:
p( x, y ) xsen( y )
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 32
Respuesta
Conclusin
Referencia
Nieves, A., & Federico, C. (2010). Metodos Numericos para ingenieros (Tercera ed.). Mxico:
bibliogrfica
Mc Graw Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 33
Unidad 5
Subtemas
Introduccin
Herramienta de
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
INGENIERA
CIVIL
Problema 2
Cdigo:
Revisin: 0
Pgina 34
En una columna de cinco platos, se requiere absorber benceno contenido en una corriente
de gas V, con un aceite L que circula a contracorriente del gas. Considrese que el benceno
transferido no altera sustancialmente el nmero de moles de V y L, fluyendo a
contracorriente, que la relacin de equilibrio est dada por la ley de henry (y = mx) y que
la columna opera a rgimen permanente. Calcule la composicin del benceno en cada
plato, usando el mtodo de Gauss Jordn
Datos:
V = 100 moles / min; L = 500 moles / min.
x0 0.0 Fraccin molar del benceno en L (el aceite entra por el domo sin benceno).
m = 0.12. Realice un balance de materia para el benceno en cada plato son con
E = 0.0001
Respuesta
INGENIERA
CIVIL
Problema 3
Cdigo:
Revisin: 0
Pgina 35
1
2
3
4
Totales
Horas
hombre, kg /
computadora
3
4
7
20
504
Metales kg /
computadora
Plsticos kg /
computadora
20
25
40
50
1970
10
15
20
22
970
Componentes,
unidades /
computadora
10
8
10
15
601
Respuesta
Conclusin
Referencia
Nieves, A., & Federico, C. (2010). Metodos Numericos para ingenieros (Tercera ed.). Mxico:
bibliogrfica
Mc Graw Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 36
Unidad 5
Subtemas
Introduccin
Ax B
a,i
xik xik 1
*100 s
xik
Herramienta de
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
En una columna de cinco platos, se requiere absorber benceno contenido en una corriente
de gas V, con un aceite L que circula a contracorriente del gas. Considrese que el benceno
transferido no altera sustancialmente el nmero de moles de V y L, fluyendo a
contracorriente, que la relacin de equilibrio est dada por la ley de henry (y = mx) y que
la columna opera a rgimen permanente. Calcule la composicin del benceno en cada
plato.
Datos:
V = 100 moles / min; L = 500 moles / min.
x0 0.0 Fraccin molar del benceno en L (el aceite entra por el domo sin benceno).
= 0.12. Realice un balance de materia para el benceno en cada plato son con E = 0.0001
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 37
Respuesta
Conclusin
Referencia
Chapra , S., & Canale, R. (2007). Mtodos Nmericos para ingenieros (Quinta ed.). Mexico:
bibliogrfica
Mc Graw Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 38
Unidad 5
Subtemas
Mtodo de Jacobi
Introduccin
Mtodo de Jacobi
El mtodo de Jacobi es un proceso simple de iteraciones de punto fijo en la solucin de
races de una ecuacin. La iteracin de punto fijo tiene dos problemas fundamentales:
a,i
Herramienta de
xik xik 1
*100 s
xik
apoyo
Desarrollo de la practica
Instruccin I
Use el mtodo de Jacobi para encontrar las soluciones al sistema de ecuaciones lineales,
apyese de Excel para resolver.
Actividad 1
En una columna de cinco platos, se requiere absorber benceno contenido en una corriente de
gas V, con un aceite L que circula a contracorriente del gas. Considrese que el benceno
transferido no altera sustancialmente el nmero de moles de V y L, fluyendo a contracorriente,
que la relacin de equilibrio est dada por la ley de henry (y = mx) y que la columna opera a
rgimen permanente. Calcule la composicin del benceno en cada plato.
Datos:
V = 100 moles / min; L = 500 moles / min.
y0 0.09
x0 0.0 Fraccin molar del benceno en L (el aceite entra por el domo sin benceno).
0.12. Realice un balance de materia para el benceno en cada plato son con E = 0.0001
m=
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 39
Respuesta
Conclusin
Bibliografa
Chapra , S., & Canale, R. (2007). Mtodos Nmericos para ingenieros (Quinta ed.). Mexico: Mc
Graw Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 40
Unidad 5
Subtemas
Introduccin
Para resolver los sistemas de ecuaciones lineales se pueden lograr usando el mtodo de
Gauss Jordn, o usar el mtodo de Gauss Seidel.
Herramienta de
apoyo
Desarrollo de la practica
Instruccin I
Problema
Producto 2
Producto 3
Producto 4
Maquina 1
Maquina 2
Maquina 3
Por ejemplo, en la produccin de una unidad del producto la mquina 1 se usa 1 hora, la
maquina 2 se usa 2 horas y la maquina 3 se usa 1 hora. Encuentre el nmero de unidades
que se deben producir de cada uno de los 4 productos un da de 8 horas completas.
Respuesta
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 41
Conclusin
Referencia
Nieves, A., & Federico, C. (2010). Metodos Numericos para ingenieros (Tercera ed.). Mxico:
bibliogrfica
Mc Graw Hill.
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 42
Unidad 6
Subtemas
Introduccin
f1 ( x, y) x 2 y 2 1
f 2 ( x, y) x 2 y 2 1
La grafica de
Herramienta de
f1 ( x, y) x 2 y 2 1 se muestra en la figura
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
f ( x, y) x2 y 2 4
f ( x, y) y x 2
INGENIERA
CIVIL
Respuesta
Problema 2
Respuesta
f ( x, y) 10( y x 2 )
f ( x, y) 1 x
Cdigo:
Revisin: 0
Pgina 43
INGENIERA
CIVIL
Problema 3
Cdigo:
Revisin: 0
Pgina 44
f ( x, y) sin( y)
f ( x, y) cos( x)
Respuesta
Conclusin
Referencia
Nieves, A., & Federico, C. (2010). Mtodos Numricos para ingenieros (Tercera ed.).
bibliogrfica
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 45
Unidad 6
Subtemas
Introduccin
Herramienta de
apoyo
Desarrollo de la practica
Instruccin I
Problema 1
Respuesta
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 46
Conclusin
Referencia
Nieves, A., & Federico, C. (2010). Mtodos Numricos para ingenieros (Tercera ed.).
bibliogrfica
INGENIERA
CIVIL
Cdigo:
Revisin: 0
Pgina 47