Taller Eje 4 - Integracion y Dif Numerica 2

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 14

TALLER EJE IV

DIFERENCIACIÓN E INTEGRACIÓN
NUMÉRICA

Asignatura : ANALISIS NUMÉRICO

Docente : ANTONIO SALUSTIANO RODRIGUEZ PACHECO

Integrantes : JOHANNA MARGARITA ARAUJO ALVAREZ

JESUS ARTURO ARENAS MEJIA


MARTIN ARSENIO DOMINGUEZ PODESLEY
Introducción

En las aplicaciones prácticas, a menudo sucede que tenemos valores de la variación de una

función en algunos puntos sin tener su expresión analítica, y queremos inferir de ellos la

función derivada en algún otro punto. De la misma forma, podemos necesitar la integral

definida de esa función de la que sólo sabemos sus valores en algunos puntos; en el caso de

la integral, además, puede darse el caso de que, incluso teniendo la expresión analítica de

la función a integrar, y siendo la función integrable, no exista la función primitiva. En estos

casos la solución pasa por aproximar la derivada y la integral usando los métodos del

cálculo numérico.

Este trabajo habla de la diferenciación e integración numérica y se muestran diversos

métodos para la resolución de diversos problemas. Se explican algunas técnicas de

aproximación para la resolución de diferenciales como lo son la fórmula de diferencia

hacia adelante y hacia atrás. También se discuten métodos desarrollados específicamente

para el problema formulado como una integral definida. Hay varias razones para llevar a

cabo la integración numérica. La principal puede ser la imposibilidad de realizar la

integración de forma analítica. Los métodos que se abordan son: método del trapecio,

métodos de Simpson 1/3 y 3/8.


1. Regla de la Diferencia hacia delante

Aproximar f´(1) donde f(x) = ln (x)

a) h=0.1

𝑓(𝑥+ℎ)−𝑓(𝑥)
𝑓′ (𝑥 ) =

𝑙𝑛 (1 + 0.1) − ln(1)
𝑓′ (1) =
0.1
𝑓′ (1) = 0.95310

b) h=0.01

𝑓(𝑥+ℎ)−𝑓(𝑥)
𝑓′ (𝑥) =

𝑙𝑛(1 + 0.01) − ln(1)


𝑓′ (1) =
0.01
𝑓′ (1) = 0.99503

c) h=0.001

𝑓 (𝑥 + ℎ) − 𝑓(𝑥)
𝑓′ (𝑥 ) =

𝑙𝑛 (1 + 0.001) − 𝑙𝑛(1)
𝑓′ (1) =
0.001
𝑓′ (1) = 0.99950
2. Regla de la Diferencia centrada de tres puntos para la segunda derivada.
𝜋
Aproximar f´’( ) donde f(x) = cos (x)
4
a) h=0.1

𝑓 (𝑥 + ℎ) − 2𝑓 (𝑥 ) + 𝑓(𝑥 − ℎ)
𝑓′′ (𝑥 ) =
ℎ2
𝜋 𝜋 𝜋
𝜋 𝑐𝑜𝑠 ( + 0.1) − 2 cos ( ) + 𝑐𝑜𝑠( − 0.1)
𝑓′′ ( ) = 4 4 4
4 0.12

𝜋
𝑓 ′′ ( ) = −0.70652
4

b) h=0.01

𝑓 (𝑥 + ℎ) − 2𝑓 (𝑥 ) + 𝑓(𝑥 − ℎ)
𝑓′′ (𝑥 ) =
ℎ2
𝜋 𝜋 𝜋
𝜋 𝑐𝑜𝑠 ( + 0.01) − 2 cos ( ) + 𝑐𝑜𝑠( − 0.01)
𝑓′′ ( ) = 4 4 4
4 0.012

𝜋
𝑓 ′′ ( ) = −0.70710
4

c) h=0.001

𝑓 (𝑥 + ℎ) − 2𝑓 (𝑥 ) + 𝑓(𝑥 − ℎ)
𝑓′′ (𝑥 ) =
ℎ2
𝜋 𝜋 𝜋
𝜋 𝑐𝑜𝑠 ( + 0.001) − 2 cos ( ) + 𝑐𝑜𝑠( − 0.001)
𝑓′′ ( ) = 4 4 4
4 0.0012

𝜋
𝑓 ′′ ( ) = −0.70711
4
3. Regla del Trapecio

𝑓 (𝑥 ) = 0.2 + 25𝑥 − 200𝑥 2 + 675𝑥 3 − 900𝑥 4 + 400𝑥 5


Ocuparemos los límites 0 y 1 para realizar el cálculo:

0.8
∫ 0.2 + 25𝑥 − 200𝑥 2 + 675𝑥 3 − 900𝑥 4 + 400𝑥 5 𝑑𝑥
0

𝑓 (𝑎) = 𝑓 (0) = 0.2


𝑓 (𝑏) = 𝑓 (0.8) = 0.232

𝑏
𝑓(𝑎) + 𝑓(𝑏)
∫ 𝑓(𝑥)𝑑𝑥 = (𝑏 − 𝑎) ( )
𝑎 2

𝑏
0.2 + 0.232
∫ 𝑓(𝑥)𝑑𝑥 = (0.8 − 0) ( ) = 0,1728
𝑎 2

Por lo que:

0.8
∫ 0.2 + 25𝑥 − 200𝑥 2 + 675𝑥 3 − 900𝑥 4 + 400𝑥 5 𝑑𝑥 ≅ 0.1728
0
4. Regla del Trapecio Compuesta

Evaluar con n= 2 y n= 4
𝜋/2
∫ 6 + 3𝑐𝑜𝑠𝑥 𝑑𝑥
0

a) n= 2

x0 x1 x2
0 𝜋 𝜋
4 2
f (x0) f (x1) f (x2)
9 8.12132 6

𝑏 𝑛−1

∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑓 (𝑎) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑏))
𝑎 2
𝑖=1
𝜋
2 𝜋
∫ 6 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ (9 + 2(8.12132) + 6)
0 8
𝜋
2
∫ 6 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ 12.26896
0
b) n= 4

x0 x1 x2 x3 x4
0 𝜋 𝜋 3𝜋 𝜋
8 4 8 2
f (x0) f (x1) f (x2) f (x3) f (x4)
9 8.77164 8.12132 7.14805 6

𝑏 𝑛−1

∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑓 (𝑎) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑏))
𝑎 2
𝑖=1
𝜋
2 𝜋
∫ 6 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ (9 + 2(8.77164 + 8.12132 + 7.14805) + 6)
0 16
𝜋
2
∫ 6 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ 12.38613
0
5. Regla de Simpson 1/3 con n=4

Evalúe:

3
∫ 𝑥 2 𝑒 𝑥 𝑑𝑥
0

𝑛=1 𝑛=2
𝑏 𝑓(𝑥0 ) + 4 ∑𝑖=1,3,5 𝑓(𝑥𝑖 ) + 2 ∑𝑖=2,4 𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑛 )
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑏 − 𝑎)
𝑎 3𝑛

𝑎=0
𝑏=3
𝑛=4
𝑏−𝑎 3−0 3
ℎ= = = = 0.75
𝑛 4 4
𝑥0 = 𝑎 = 0
𝑓(𝑥0 ) = 0

𝑥1 = 𝑥0 + ℎ = 0.75
𝑓(𝑥1 ) = 𝑓 (0.75) = 1.19081

𝑥2 = 𝑥1 + ℎ = 0.75 + 0.75 = 1.5


𝑓(𝑥2 ) = 𝑓 (1.5) = 10.08380

𝑥3 = 𝑥2 + ℎ = 1.5 + 0.75 = 2.25


𝑓(𝑥3 ) = 𝑓 (2.25) = 48.03166
𝑥4 = 𝑥3 + ℎ = 2.25 + 0.75 = 3
𝑓(𝑥4 ) = 𝑓 (3) = 180.76983

Reemplazando en la fórmula de Simpson 1/3:


𝑛=1 𝑛=2
𝑏 𝑓(𝑥0 ) + 4 ∑𝑖=1,3,5 𝑓(𝑥𝑖 ) + 2 ∑𝑖=2,4 𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑛 )
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑏 − 𝑎)
𝑎 3𝑛

𝑏
0 + 4(1.19081 + 48.03166) + 2(10.08380) + 180.76983
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (3)
𝑎 3∗4

3
∫ 𝑥 2 𝑒 𝑥 𝑑𝑥 ≅ 99.4568275
0
6. Regla de Simpson 1/3 y 3/8 con n=5

Evalúe:

3
∫ 5 + 3𝑐𝑜𝑠𝑥 𝑑𝑥
0

a) Usando Simpson 1/3:


𝑛=1 𝑛=2
𝑏 𝑓(𝑥0 ) + 4 ∑𝑖=1,3,5 𝑓(𝑥𝑖 ) + 2 ∑𝑖=2,4 𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑛 )
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑏 − 𝑎)
𝑎 3𝑛

𝑎=0
𝑏=3
𝑛=4 𝑠𝑒 𝑜𝑐𝑢𝑝𝑎𝑟á 4 𝑦𝑎 𝑞𝑢𝑒 𝑛 𝑑𝑒𝑏𝑒 𝑠𝑒𝑟 𝑝𝑎𝑟
𝑏−𝑎 3−0 3
ℎ= = = = 0.75
𝑛 4 4
𝑥0 = 𝑎 = 0
𝑓(𝑥0 ) = 5 + 3 cos(0) = 8

𝑥1 = 𝑥0 + ℎ = 0.75
𝑓(𝑥1 ) = 5 + 3𝑐𝑜𝑠(0.75) = 7.19507

𝑥2 = 𝑥1 + ℎ = 0.75 + 0.75 = 1.5


𝑓(𝑥2 ) = 5 + 3𝑐𝑜𝑠(1.5) = 5.21221

𝑥3 = 𝑥2 + ℎ = 1.5 + 0.75 = 2.25


𝑓(𝑥3 ) = 5 + 3𝑐𝑜𝑠(2.25) = 3.11548
𝑥4 = 𝑥3 + ℎ = 2.25 + 0.75 = 3
𝑓(𝑥4 ) = 5 + 3𝑐𝑜𝑠(3) = 2.03002

Reemplazando en la fórmula de Simpson 1/3:


𝑛=1 𝑛=2
𝑏 𝑓(𝑥0 ) + 4 ∑𝑖=1,3,5 𝑓(𝑥𝑖 ) + 2 ∑𝑖=2,4 𝑓(𝑥𝑖 ) + 𝑓(𝑥𝑛 )
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑏 − 𝑎)
𝑎 3𝑛

𝑏
8 + 4(7.19507 + 3.11548) + 2(5.21221) + 2.03002
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (3)
𝑎 3∗4
3
∫ 5 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ 15.42416
0

b) Usando Simpson 3/8:


𝑏
𝑓(𝑥0 ) + 3𝑓(𝑥1 ) + 3𝑓(𝑥2 ) + 𝑓(𝑥3 )
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑏 − 𝑎)
𝑎 8
𝑎=0
𝑏=3
3−0
ℎ= =1
3
𝑥0 = 𝑎 = 0
𝑓(𝑥0 ) = 5 + 3 cos(0) = 8

𝑥1 = 𝑥0 + ℎ = 1
𝑓(𝑥1 ) = 5 + 3𝑐𝑜𝑠(1) = 6.62091
𝑥2 = 𝑥1 + ℎ = 1 + 1 = 2
𝑓(𝑥2 ) = 5 + 3𝑐𝑜𝑠(2) = 3.75156

𝑥3 = 𝑥2 + ℎ = 2 + 1 = 3
𝑓(𝑥3 ) = 5 + 3𝑐𝑜𝑠(3) = 2.03002

Reemplazando en la fórmula de Simpson 1/8:


𝑏 𝑓(𝑥0 ) + 3𝑓(𝑥1 ) + 3𝑓(𝑥2 ) + 𝑓(𝑥3 )
∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (𝑏 − 𝑎)
𝑎 8

𝑏 8 + 3(6.62091) + 3(3.75156) + 2.03002


∫ 𝑓 (𝑥 )𝑑𝑥 ≅ (3 − 0)
𝑎 8
3
∫ 5 + 3𝑐𝑜𝑠𝑥 𝑑𝑥 ≅ 15.43028625
0
Conclusión

La integración con métodos numéricos es una herramienta útil cuando se trata de integrar

una función muy complicada o datos tabulados.

El método de Derivación centrada es mucho más preciso que el método de diferencias hacia

atrás y el método de diferencias hacia adelante, pero por otro lado tiene más funciones de

evaluar por lo que su algoritmo es más tedioso y echo analíticamente más propenso a fallar

por el número de iteraciones.

Con el método de Simpson se puede aproximar una integral compleja a la integral de un

polinomio. Se logra convertir matemáticas superiores en aritméticas básicas. Tanto para la

regla de Simpson de 1/3 y la de 3/8 se pueden aplicar juntas sobre una misma curva para

obtener exactitudes de tercer orden sobre todo un intervalo. Además, la regla de Simpson de

1/3 y 3/8 es más exacta que otras formas numéricas de aproximar la integral como la regla

del trapecio.

Como se mostró existen diferentes fórmulas para la realización de diferenciales y pudimos

notar que algunas tienen un grado de dificultad mayor que otras, pero la utilización de ellas

depende del tipo de diferencial a resolver, es decir del tipo de problema a enfrentar, debemos

saber identificar cuál de ellas utilizar para su correcta resolución. Lo mismo sucede son los

diferentes métodos de integración numérica ya que también existen diversos

procedimientos, pero la buena resolución de un problema de integración se deberá del uso

correcto y la buena identificación del tipo de método a utilizar.


Bibliografía

https://www.ugr.es/~arobles/Telecos/Tema10.pdf

https://www.academia.edu/32875352/Diferenciaci%C3%B3n_Num%C3%A9rica

https://vdocuments.mx/metodos-de-diferenciacion-e-integracion-numerica.html

También podría gustarte