Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Integracion Numerica

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

INTEGRACIÓN NUMÉRICA

Jhon Polo Vásquez

Universidad de la Costa CUC

Métodos numéricos
Motivación

Existen varias técnicas analíticas de integración, con las que


podemos resolver integrales definida, entre las más conocidas
están: sustitución e integración por partes. El problema para
usar los métodos analíticos de integración es que, es posible
que f no se pueda expresar en términos de funciones elemen-
tales, o f se conoce explicitamente pero no se puede evaluar
fácilmente. Por ejemplo,
R1 1
0 1+x5 dx
R3 1
2 ln(x) dx
R π/2
0 cos(x2 ) dx
R 5 −x2
1 e dx
R1 √3
0 1 + x3 dx

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Integración numérica

Son métodos numéricos para estimar el valor de una integral


definida de la forma
Z b
I= f (x) dx,
a

donde el intervalo [a, b] es finito y f es una función de una


variable real y valor real continua en [a, b].
Además, por el Teorema Fundamental del Cálculo, se sigue
que
Z b
f (x) dx = F (b) − F (a),
a
donde F es una antiderivada de f en [a, b].

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Fórmulas cerradas de Newton-Cotes

La regla de los trapecios.

La regla de Simpson 1/3.

La regla de Simpson 3/8.

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Discretización del dominio

Dividamos el intervalo [a, b] en N subintervalos de igual lon-


gitud,

[x0 , x1 ], [x1 , x2 ], . . . , [xk , xk+1 ], . . . , [xN −1 , xN ],

donde los N + 1 puntos x0 , x1 , . . . , xN de la partición, se obtie-


nen a partir de la fórmula

xk = a + kh, k = 0, 1, . . . , N

Notese que x0 = a, xN = b y h = xk+1 − xk .

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Regla de los trapecios

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Regla de los trapecios
Caso Compuesto (N > 1)
N −1
Z b " #
h X
f (x) dx = f (a) + f (b) + 2 f (xk )
a 2
k=1
Error local
h3
00
Ek = f (ξk ) ,

12
para algún ξ ∈ (xk , xk+1 ).
Caso simple (N = 1)
Z b
b − a 
f (x) dx = f (a) + f (b)
a 2
Error Total
h2 (b − a)
00
ET = f (ξ) ,

12
para algún ξ ∈ (a, b).
Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA
Ejemplo

Sabiendo que 0 sen(x) dx = 2. Aproximar dicha integral con
tres subintervalos y estimar el error.

Solución:
π−0
N = 3, h= 3 = π3 .

π 2π
x0 = 0, x1 = , x2 = , x3 = π
3 3
Entonces,

Z π
π/2
sen(x) dx ≈ [sen(0) + sen(π) + 2(sen(π/3) + sen(2π/3))]
0 3
≈ 1.81379

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Para estimar el error,

f 00 (x) = −sen(x)

Tomando ξ = π/2, obtenemos


(π/3)2 (π)
ET = (−sen(π/2))

12
= 0.28709

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Regla de Simpson 1/3

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Regla de Simpson 1/3
Caso Compuesto (N = 2m, con m ≥ 2)
" N −1 N −1 #
Z b 2 2
h X X
f (x) dx ≈ f (a)+f (b)+4 f (x2k+1 )+2 f (x2k )
a 3
k=0 k=1
Error local
h5
(iv)
Ek = f (ξk ) ,

90
para algún ξ ∈ (xk , xk+2 ) con k = 0, 2, . . . , N − 2.
Caso simple (N = 2)
Z b " !#
b−a a+b
f (x) dx ≈ f (a) + f (b) + 4f
a 6 2
Error Total

h4 (b − a)
ET = f (iv) (ξ) , ξ ∈ (a, b).

180
Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA
Ejemplo

Aproximar 0 sen(x) dx con N = 4 y estimar el error.

Solucion: N = 4, h = π4 .
π π 3π
x0 = 0, x1 = , x2 = , x3 = , x4 = π
4 2 4
Entonces,

Z π
π/4
sen(x) dx ≈ [sen(0) + sen(π) + 4[sen(π/4) + sen(3π/4)]
0 3
+ 2sen(π/2)]
≈ 2.00455

Tomando ξ = π/2, tenemos que el error es


(π/4)4 π
ET = 1 = 0.006641
180
Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA
Ejemplo

R π/3
Usar la regla de Simpson 1/3 simple para estimar 0 sen2 (x) dx.

Solución:

" #
Z π/3
π
sen2 (x) dx ≈ sen2 (0) + sen2 (π/3) + 4sen2 (π/6)
0 18
≈ 0.3054326

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Regla de Simpson 3/8

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Regla de Simpson 3/8
Caso Compuesto (N = 3m, con m ≥ 2)
" N −3
Z b 3
3h X
f (x) dx ≈ f (a) + f (b) + 3 f (x3k+1 )
a 8
k=0
N −3 N −3 #
X3 X3

+3 f (x3k+2 ) + 2 f (x3k )
k=0 k=1

Error local
3
Ek = h5 f (iv) (ξk ) ,

80
para algún ξ ∈ (xk , xk+3 ) con k = 0, 3, . . . , N − 3.
Caso simple (N = 3)
Z b " ! !#
b−a 2a + b a + 2b
f (x) dx ≈ f (a)+f (b)+3f +3f
a 8 3 3

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Error Total

h4 (b − a)
ET = f (iv) (ξ) , ξ ∈ (a, b).

80

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Ejemplo

Utilice la regla de Simpson 3/8 simple y compuesta con N = 6,


R π/3
para aproximar 0 sen2 (x) dx.

Solución:
Caso compuesto

Caso simple

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA


Gracias...

Jhon Polo Vásquez INTEGRACIÓN NUMÉRICA

También podría gustarte