Unidad IV Métodos
Unidad IV Métodos
Unidad IV Métodos
AGUASCALIENTES
MÉTODOS NUMÉRICOS
04 DE JUNIO DE 2021
INTRODUCCIÓN
Para entrar un poco en contexto la integración numérica constituye una herramienta para
aproximar numéricamente el valor de integrales definidas de difícil cálculo o cuya solución
no se puede lograr mediante medios analíticos. Su importancia resulta en sus aplicaciones,
tales como, el cálculo de la capacidad de un lago a partir de datos topográficos, la fuerza
ejercida por el aire sobre las alas de un avión, la distancia recorrida a partir de un tacómetro,
por mencionar algunos de ellas.
2
UNIDAD IV “INTEGRACIÓN NÚMERICA”
Este método consiste en la aproximación del cálculo del área plana bajo una curva utilizando
trapecios curvilíneos a partir una interpolación con una función cuadrática:
𝑏
𝑏−𝑎 𝑎+𝑏
∫ 𝑓(𝑥)𝑑𝑥 ≈ [𝑓(𝑎) + 4𝑓 ( ) + 𝑓(𝑏)]
𝑎 6 2
La regla de Simpson tiene mayor aproximación que regla de los trapecios. En la regla de los
trapecios los puntos sucesivos de la gráfica y = f(x) se unen mediante líneas que forman los
trapecios, en la regla de Simpson los puntos se unen mediante segmentos de parábolas.
Esta aproximación es denominada “simple” debido a que utiliza tan solo un polinomio.
Requiere el conocimiento de tres puntos equiespaciados: los extremos y un punto central.
Aplicando esta expresión utilizando mayor cantidad de puntos intermedios (es decir,
realizando más de un Simpson 1/3 dentro del intervalo) puede definirse la variante
“compuesta” del método para el cual se utilizan N puntos que corresponden a n = N − 1
subintervalos. Este caso requiere que la cantidad de subintervalos sean pares (el caso simple
utiliza dos, por lo tanto ´este debe ser un múltiplo). Por lo tanto, se define un valor h que
corresponde a el ancho del subintervalo o el paso que hay entre puntos. Se calcula como h =
b−a n. Finalmente, la aproximación del área se puede calcular como:
𝑛−2 𝑛−2
𝑏 2 2
ℎ
∫ 𝑓(𝑥)𝑑𝑥 ≈ 𝑓(𝑎) + 𝑓(𝑏) + 4 ∑ 𝑓(𝑥2𝑘+1 ) + 2 ∑ 𝑓(𝑥2𝑘 )
𝑎 3
𝑘=0 𝑘=1
[ ]
A modo de interpretación, las sumatorias corresponden a los valores de la función en los
puntos impares y pares respectivamente sin contar los extremos.
3
EJEMPLO
2
Calcular el área de la siguiente función 𝑓(𝑥) = 𝑒 −𝑥 desde x= -1 hasta x=1.
Para calcular el área debe obtenerse la integral de la función, y por lo tanto la primitiva de la
función. En este caso es posible encontrarla en términos de funciones elementales; por lo
tanto, se realizará una aproximación simple, para ello se utiliza la siguiente fórmula:
1
2 1 − (−1) 1−1
∫ 𝑒 −𝑥 𝑑𝑥 ≈ [𝑓(−1) + 4𝑓 ( ) + 𝑓(1)]
−1 6 2
1
≈ [0.3678 + 4 ∙ 1 + 0.3678] ≈ 1.5786
3
A medida que aumente la cantidad de puntos, la separación entre ellos disminuye y se acerca
al valor de la integral definida. Por lo tanto, se plantea para el caso compuesto utilizando la
fórmula de aproximación de área con 4 puntos:
1 − (−1)
𝑛=4⟹ℎ= = 0.5
4
Luego se realiza una tabla para identificar qué valores es necesario calcular y en que parte de
la ecuación se utiliza. En este caso, se construye utilizando los valores de x espaciados por h.
Posteriormente, se calculan los valores correspondientes de y y se observan los puntos pares
e impares. En este caso, los rojos son los impares y los azules impares.
4
▪ Regla de Simpson 3/8
El esquema de integración de Simpson 3/8 considera pasar polinomios de tercer grado entre
cada cuatro puntos de la función y(x) y está dado por la integral siguiente
Para obtener:
Donde h = (b – a) /3. Esta ecuación se llama regla de Simpson 3/8 debido a que h se multiplica
por 3/8.
Para que la fórmula de Simpson 3/8 sea aplicable n debe ser múltiplo de 3. Los esquemas de
Simpson 1/3 y Simpson 3/8 proporcionan resultados similares, aunque la función
corresponda a un polinomio de grado superior a tres; estas dos expresiones son las que
producen menores errores en la estimación de la integral, pero hay que tomar en cuenta que
tienen restricciones en su aplicación. También se puede calcular el valor de la integral
combinando las dos expresiones anteriores, utilizándolas por tramos, en que sean aplicables.
EJEMPLO
Con la regla de Simpson 3/8 integre la siguiente función:
𝑓(𝑥) = 0.2 + 25𝑥 − 200𝑥 2 + 675𝑥 3 − 900𝑥 4 + 4005
Desde a= 0 hasta b= 8
Solución
Una sola aplicación de la regla de Simpson 3/8 requiere cuatro puntos equidistantes:
𝑓(0) = 0.2 𝑓(0.2667) = 1.432724
𝑓(0.5333) = 3.487177 𝑓(0.8) = 0.232
Utilizando la ecuación:
0.2 + 3(1.432724 + 3.487177 + 0.232)
𝐼 ≅ 0.8 = 1.519170
8
5
4.1.1 REGLA TRAPEZOIDAL
La regla del trapecio es un método de integración numérica, es decir, un método para calcular
aproximadamente el valor de la integral definida.
La regla se basa en aproximar el valor de la integral de f(x) por el de la función lineal que
pasa a través de los puntos (a, f(a)) y (b, f(b)). La integral de ésta es igual al área del trapecio
bajo la gráfica de la función lineal:
𝑏
𝑓(𝑎) + 𝑓(𝑏) Donde el término error
∫ 𝑓(𝑥)𝑑𝑥 ≈ (𝑏 − 𝑎) corresponde a
𝑎 2
Siendo 𝜉 un número perteneciente al intervalo [a, b]
Geométricamente, la regla del trapecio es equivalente a aproximar el área del trapecio bajo
la línea recta que une f(a) y f(b) como se muestra en las siguientes imágenes:
La fórmula para calcular el área de un trapezoide: altura por el promedio de las bases.
Para la regla del trapecio, el concepto es el mismo, pero ahora el trapezoide está sobre su
lado; por lo tanto, la integral aproximada se representa como:
𝐼 ≅ (𝑏 − 𝑎) × altura promedio
Donde, para la regla del trapecio, la altura promedio es el promedio de los valores de la función en
los puntos extremos, o [f(a) + f(b)]/2.
6
4.2 APLICACIONES DE LA INTEGRACIÓN NÚMERICA
La integración numérica tiene diferentes aplicaciones en la rama de la ingeniería, esto permite
resolver los problemas matemáticos. En situaciones de este tipo, se aplican los métodos
numéricos usando la expresión analítica para generar una tabla de valores de la función.
En la ingeniería civil/ ambiental se utiliza las integrales numéricas para resolver problemas
donde se tiene que calcular las fuerzas del viento (f), ejercidas por pie de mástil de las velas,
varían en función de la distancia sobre la cubierta del bote (z). Para resolver problemas de
este tipo, se requiere que la fuerza distribuida f se convierta en una fuerza total equivalente
F y que se calcule su localización d sobre la cubierta. Este cálculo se complica por el hecho
de que la fuerza ejercida por pie de mástil varía con la distancia sobre la cubierta. La fuerza
total ejercida sobre el mástil se expresa como la integral de una función continua.
Por su parte en la ingeniería eléctrica se emplean cálculos del valor promedio de una
corriente eléctrica oscilante en un periodo puede ser cero, para su solución se utiliza algunas
estimaciones de la integral para varias aplicaciones de la regla del trapecio y de la regla de
Simpson 1/3.
En ingeniería muchos problemas implican el cálculo del trabajo. La fórmula general es:
Trabajo = fuerza × distancia, son algunas aplicaciones simples mediante el uso de fuerzas
que permanecían constantes durante todo el desplazamiento. Aunque ese simple cálculo es
útil para presentar el concepto, la solución de problemas reales por lo común es más
complicada. Por ejemplo, la fuerza varía durante el proceso del cálculo. En tales casos, la
ecuación para el trabajo ahora se expresa.
7
4.3 USO DE HERRAMIENTAS COMPUTACIONALES
La mayor parte de los programas computacionales que sirven para integrar una función de
una sola variable real, se basan en el método adaptativo o en fórmulas gaussianas
extremadamente exactas. Por lo general, las integrales múltiples se aproximan ampliando
buenos métodos adaptativos a dimensiones superiores.
Netlib: es una colección grande de software, documentos, bases de datos gratis que son de
interés para las comunidades científicas y de métodos numéricos. Netlib contiene software
de alta calidad que ha sido probado en forma intensiva, pero todo el software libre no tiene
garantía y poco (si existe) soporte. Para poder usar el software, primero se tiene que descargar
en su computadora y entonces compilarlo.
IMSL: La biblioteca numérica de Fortran IMSL hecha por Visual Numerics, Inc. cubre
muchas de las áreas contenidas en la biblioteca NAG. También tiene soporte para analizar y
presentar datos estadísticos en aplicaciones científicas y de negocios.
8
NUMERICAL RECIPES: son muy populares entre los ingenieros porque pueden ser
usados como libro de cocina donde se puede encontrar una “receta” para resolver algún
problema a mano. Sin embargo, el software correspondiente de Numerical Récipes no es
comparable en alcance o calidad al dado por NAG o IMSL. Es un software muy usado en
universidades, centros de investigación y por ingenieros. En los últimos años ha incluido
muchas más capacidades, como la de programar directamente procesadores digitales de
señal, crear código VHDL y otras.
GNU OCTAVE: Es un programa libre para realizar cálculos numéricos. Como indica su
nombre es parte de proyecto GNU. MATLAB es considerado su equivalente comercial; entre
varias características que comparten se puede destacar que ambos ofrecen un intérprete
permitiendo ejecutar órdenes en modo interactivo. Nótese que Octave no es un sistema de
álgebra computacional como podría ser GNU Máxima, sino que usa un lenguaje que está
orientado al análisis numérico.
9
ÍNDICE
INTRODUCCIÓN.................................................................................................................2
INTEGRACIÓN NÚMERICA
ÍNDICE................................................................................................................................10
BIBLIOFRAFÍA.................................................................................................................11
10
BIBLIOGRAFÍA
[1] Burden R. L. & Douglas J. F. (2001). Métodos Numéricos. (7ª ed.) Ceagage Learning
Latinoamérica.
[2] Chapra, C.. (2007). Métodos Numéricos para Ingenieros. (5ª ed.). McGraw Hill.
[3] Nakamura, S. (1992). Métodos numéricos aplicados con software. México. Prentice Hall
Hispanoamericana.
computo-numerico
11